def test_ll2ne_bounds(self): """test of ll2ne's bounds""" from ll2ne import ll2ne #, ne2ll northing, easting, zone, h = ll2ne(-90, 0) northing, easting, zone, h = ll2ne(90, 0) northing, easting, zone, h = ll2ne(0, -180) northing, easting, zone, h = ll2ne(0, 180)
def xxxxtest_ne2ll2ne(self): """Test conversions back and forward""" from ne2ll import ne2ll from ll2ne import ll2ne eps = 1.0e-6 #Less accuracy #Northings, eastings, zone and hemisphere nes = [(6320371.85657, 716174.605757, 55, 'S'), (1505230,345745,55,'N'), (1505230,345745,55,'S')] for northing, easting, zone, hemisphere in nes: lat, lon = ne2ll(northing, easting, zone, hemisphere) northing1, easting1, zone1, hemisphere1 = ll2ne(lat, lon) print northing, northing1 print easting, easting1 assert abs(northing - northing1)/abs(northing) < eps,\ "%f %f" %(northing, northing1) assert abs(easting - easting1)/abs(easting) < eps,\ "%f %f" %(easting, easting1) assert zone == zone1 assert hemisphere == hemisphere1
def test_ll2ne(self): """Basic test of ll2ne""" from ll2ne import ll2ne #, ne2ll lat = -(37.0 + (39.0 + 10.15611/60.0)/60.0) lon = (143.0 + (55.0 + 35.383900/60.0)/60.0) northing, easting, zone, hemisphere = ll2ne(lat, lon) assert abs(northing - 5828674.33975)/abs(northing) < self.eps assert abs(easting - 758173.79727)/abs(easting) < self.eps assert zone == 54 assert hemisphere == 'S'
def test_ll2ne2ll(self): """Test conversions back and forward""" from ne2ll import ne2ll from ll2ne import ll2ne eps = 1.0e-6 #Less accuracy latlons = [(-33.234, 149.32), (53.02,2.23421), (0,0), (40, -50), (-36,0), (-50,-60), (-89.9,0), (-89,-2), (89.9,0), (0, -180), (0, 180)] for lat, lon in latlons: northing, easting, zone, hemisphere = ll2ne(lat, lon) lat1, lon1 = ne2ll(northing, easting, zone, hemisphere) #print lat, lon, zone, easting, northing, hemisphere, lat1, lon1 assert abs(lat - lat1) < eps, "%f, %f" %(lat, lat1) assert abs(lon - lon1) < eps, "%f, %f" %(lon, lon1)