def xyp2latlon(x, y, panel, rotate_lat=RLAT, rotate_lon=RLON, R=1): xyz = xyp2xyz(x, y, panel) xr, yr, zr = xyz_rotate_reverse(xyz, rotate_lat, rotate_lon) lat, lon = xyz2latlon(xr, yr, zr, R) """ print 'x=%.18f, y=%.18f'%(x,y) print 'x=%.18f, y=%.18f, z=%.18f'%xyz print 'xr=%.18f, yr=%.18f, zr=%.18f'%(xr,yr,zr) print 'lat=%.18f, lon=%.18f'%(lat,lon) """ return (lat, lon)
def test_xyp_rotate_reverse(): ''' xyp_rotate() -> xyz_reverse_rotate() : check consistency, repeat 1000 times ''' from cart_rotate import xyz_rotate, xyz_rotate_reverse from cart_ll import latlon2xyz N = 1000 for i in xrange(N): lat = pi*rand() - pi/2 lon = 2*pi*rand() x, y, z = latlon2xyz(lat, lon) rlat = pi*rand() - pi/2 rlon = 2*pi*rand() xr, yr, zr = xyz_rotate((x, y, z), rlat, rlon) x2, y2, z2 = xyz_rotate_reverse((xr, yr, zr), rlat, rlon) aa_equal((x,y,z), (x2,y2,z2), 15)