Contour plots of lat/lon data

*To*: mathgroup at smc.vnet.net*Subject*: [mg132072] Contour plots of lat/lon data*From*: amannucci <Anthony.J.Mannucci at jpl.nasa.gov>*Date*: Sun, 1 Dec 2013 08:23:04 -0500 (EST)*Delivered-to*: l-mathgroup@mail-archive0.wolfram.com*Delivered-to*: l-mathgroup@wolfram.com*Delivered-to*: mathgroup-outx@smc.vnet.net*Delivered-to*: mathgroup-newsendx@smc.vnet.net

I am interested in using Mathematica to plot functions of latitude/longitude (e.g. on Earth). I am thinking of ContourPlot type of functionality, but with latitude/longitude coordinates. I am familiar how this is done in matplotlib/python. It works like this: One initializes the map projection thus: m = Basemap(projection='npaeqd',boundinglat=10,lon_0=270,resolution ='l') This projection is "North Polar Azimuthal Equidistant Projection". Then, given an array of data in lat, lon coordinates z(lat,lon), one converts to projection coordinates thus: x,y = m(lon,lat) And then executes the contour method of the map projection object: m.contourf(x,y,z) Is there a polar projection in Mathematica? I cannot find one. If not, I can probably get reasonably close with the existing set of projections. More broadly, I don't see how to create a geographically registered contour map in Mathematica. Is there a function that plots contour data taking into account a map projection that is already defined? I cannot find any examples of contour plotting on a user-specified geographic projection. Does one manually convert all the coordinates and then use ContourPlot? Thanks for your help! -Tony