WARNING: This server is unstable and will be retired in the next days. If you want to keep this forum available, please request immediately a migration on the Nabble Support forum. Forums that don't receive any migration request will be deleted forever.

 « Return to Thread: Could basemap plot in a UTM projection? - pyproj documentation

Re: Could basemap plot in a UTM projection? - pyproj documentation

by Mathew Topper :: Rate this Message:

| View in Thread

Some parts of this message have been removed. Learn more about Nabble's security policy.
On 02/04/12 23:23, Jeff Whitaker wrote:
On 4/2/12 4:02 PM, Mathew Topper wrote:
On 02/04/12 19:09, Jeff Whitaker wrote:

Mat:  I just updated git master to add the 'k_0' keyword when using projection='tmerc'.  This is the scaling factor along the central meridian, and is needed to exactly match the UTM definition.  I also added a utmtest.py example, that shows how to plot a UTM zone, and convert Basemap coordinates to UTM coordinates.

-Jeff
-- 
Jeffrey S. Whitaker         Phone  : (303)497-6313
Meteorologist               FAX    : (303)497-6449
NOAA/OAR/PSD  R/PSD1        Email  : Jeffrey.S.Whitaker@...
325 Broadway                Office : Skaggs Research Cntr 1D-113
Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg

Jeff, thanks again for all your help with this. I tried the example and it works great. My only question would be how to get the axis labels into meters rather than degrees. Any hints? I have found the customticks.py example. Is that the best route to take?
Mat:  Just use 'suppress_ticks=False' when creating the Basemap instance, and you will get axis ticks in meters. If you want to use a custom tick formatter, use the customticks.py example as a guide.

-Jeff

Jeff, I wonder whether I could run this past you? Basically, I believe the easting and northing grid lines are available just by turning on the standard matplotlib grid, so therefore I'm assuming that the easting and northing grid will be orthogonal throughout the map - is that correct?

Assuming all is well with that then I need to add the correct offset to the axes values, presumably by calculating the value, in meters, of the origin. I'm guessing I will have to use a tick formatter (like customticks.py) to do this or the map points will no longer be on the map (if, for instance, I just changed the axes ranges).

Does that make sense to you?

Thanks

Mat

--
Dr. Mathew Topper
Institute for Energy Systems
School of Engineering
The University of Edinburgh
Faraday Building
The King’s Buildings
Edinburgh EH9 3JL
Tel: +44 (0)131 650 5570
School fax: +44 (0)131 650 6554
mathew.topper@...
http://www.see.ed.ac.uk

The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.

------------------------------------------------------------------------------
Better than sec? Nothing is better than sec when it comes to
monitoring Big Data applications. Try Boundary one-second
resolution app monitoring today. Free.
http://p.sf.net/sfu/Boundary-dev2dev
_______________________________________________
Matplotlib-users mailing list
Matplotlib-users@...
https://lists.sourceforge.net/lists/listinfo/matplotlib-users

 « Return to Thread: Could basemap plot in a UTM projection? - pyproj documentation