Contents:
Planetary mapping routines.
phi = 0 faces toward the observer phi = pi thus faces away from the observer theta=pi/2 is the z-axis or ‘north pole’ theta=-pi/2 is the ‘south pole’
Give the chi-squared to an input timeseries to get the wedge coefficients and inclination angle.
param – [inclination, phi0, wedge coefficients] phi,theta – from MAKEGRID rangle – rotation angles at which flux was measured meas – measured flux err – error on measurement (one-sigma)
Make grids of phi and theta values with the specified number of points in each direction. Phi ranges from 0 to 2pi, and theta ranges from -pi/2 to pi/2.
Returns (phi, theta)
Make a spot-map.
EXAMPLE: import maps inc, rot = pi/2., 0. phi, theta = maps.makecoord(120,60) phi2, theta2 = maps.rotcoord(phi,theta,inc,rot) vmap = maps.visiblemap(phi,theta,inc,rot)
maps.makespot(phi,theta,inc=pi/2,rot=0.,lat=pi/4,long=1.*pi) maps.makespot(phi2,theta2,vmap=vmap,lat=pi/4,long=1.*pi)
Return a 2D map of the projected area. Note that you need to determine for yourself which parts of the areal map are actually visible to the observer.
rotate coordinate system from local (planetary) coordinates into observer-oriented coordinages. rangle and iangle are the rotation and inclination angles of the planet in radians.
returns (phi2,theta2)
phi2 will be in the range (0,2pi) and theta2 will be in the range (-pi/2,pi/2)
Return a 2D boolean map that’s True for the planetary latitude, longitude values visible from an observer. For rangle (rotation angle) zero and iangle (inclination angle) equal to pi/2, phi=pi is toward the observer.
Return a set of basis functions for the flux from each wedge.
phi,theta are the observer-centered grids from MAPS.MAKEGRID iangle is the inclination of the system (0 = pole-on) rangle (seq.) is a sequence of (0,2pi) rotation values at which
the wedge-based flux is evaluated.
nwedge is the number of wedges phi0 defines the center of wedge zero.
fwedge – an optional sequence to set individual flux values for each wedge.