Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2013

[Date Index] [Thread Index] [Author Index]

Search the Archive

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



  • Next by Date: Creating an ordered list of pairs
  • Next by thread: Creating an ordered list of pairs