Hi:
> I've hacked together some scripts that let me plot data on
> top of world maps with different map projections in
> matplotlib. It's quite slow, but it works! There's an
> example at
> http://www.cdc.noaa.gov/people/jeffrey.s.whitaker/plotmap.py
> (the plot should look like
> http://www.cdc.noaa.gov/people/jeffrey.s.whitaker/plotmap.png).
> To run the example you'll need the data files
> http://www.cdc.noaa.gov/people/jeffrey.s.whitaker/wcl.txt
> (world coastlines) and
> http://www.cdc.noaa.gov/people/jeffrey.s.whitaker/topodata.pickle
> (data to plot), as well as a module I wrote to compute
> cartographic transformations using proj4
> http://www.cdc.noaa.gov/people/jeffrey.s.whitaker/python/proj.py.
> If anyone knows of a faster way to do this, I'd love to hear
> about it.
I am having trouble running this - proj.py is making the system call
'stdin,stdout=os.popen2(cmd,mode='b')' where cmd is
proj +R=63712000 +lon_0=-107 +proj=lcc +lat_1=50 +lat_2=50 -b
but I don't have proj on my system. Am I missing something? Is this
a tool I can easily install?
Also, while I have your ear, what version of matplotlib are you
working with and what backend do you typically use?
Nice map, by the way!
JDH