from pylab import *
from geospatial_data import *
from pmesh.create_mesh import convert_points_from_latlon_to_utm
from coordinate_transforms.redfearn import degminsec2decimal_degrees
from os import sep
from utilities.polygon import inside_polygon
import project


datadir = project.datadir
ion()
hold(False)
    
file1 = datadir+'onslow_onshore_30m_dted.pts'
file2 = datadir+'onslow_offshore_points.xya'

G1 = Geospatial_data(file_name = file1)
G2 = Geospatial_data(file_name = file2)

pts1 = G1.get_data_points(absolute = True)
pts2 = G2.get_data_points(absolute = True)
z1 = G1.get_attributes()
z2 = G2.get_attributes()

minz1 = min(z1)
maxz1 = max(z1)
minz2 = min(z2)
maxz2 = max(z2)

print ''
print 'Offshore elevation in range [%.2f %.2f]   ' %(minz1, maxz1)
print 'Onshore elevation in range [%.2f %.2f]   ' %(minz2, maxz2)
print ''

minelev = min(minz1,minz2)
maxelev = max(maxz1,maxz2)

# bounding polygon for Onslow scenario
polygon_bound = project.polyAll

x_bound = []
y_bound = []
n = len(polygon_bound)
for i in range(n+1):
    if i == n:
        thispt = polygon_bound[0]
    else:
        thispt = polygon_bound[i]
    x_bound.append(thispt[0])
    y_bound.append(thispt[1])

# interior regions
figure(1)
plot(pts1[:,0],pts1[:,1],'g.',
     pts2[:,0],pts2[:,1],'b.',
     x_bound,y_bound,'r-')
#title('Offhore (g) and onshore (b) data: Onslow scenario. \n Elevation in range [%.2f %.2f]' %(minelev, maxelev))
xlabel('x')
ylabel('y')
savefig('onslow_data_extent')