| Plotting data on orthographic globe projection [message #79422] |
Tue, 28 February 2012 13:32 |
Sir Loin Steak
Messages: 42
Registered: January 2012
|
Member |
|
|
Hi all
I have some data in a (72,36) array, corresponding to 72 longitudes
and 36 latitudes. I am wanting to plot them onto a globe, and so have
used IDLs map_set and contour routines. However, depending upon which
angle I choose to view the image, there are problems with the filling
of the contours.
As an example, see the four images at http://www.physics.open.ac.uk/~lsteele/globe.png.
I have plotted the exact same data in each image, the only difference
being I changed the value of 'loncentre' in the plotting routine
below:
;---------------------------------------
lat = findgen(36)*180/35-90
lon = findgen(72)*360/71-180
latcentre=30
loncentre=-70
angle=0
map_set,latcentre,loncentre,angle,/orthographic,/isotropic
loadct, 13
contour,reverse(vapour(*,*,250),2),lon,lat,/overplot,
levels=findgen(30), /fill
map_grid, latdel=5, londel=5, glinestyle=0, glinethick=0.5
;---------------------------------------
Has anyone used these routines and experienced the same problem? Or
can anyone suggest any other routines that produce something similar?
I'm stumped!
Many thanks,
Liam
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
