| Overplot on eastern US [message #7458] |
Tue, 19 November 1996 00:00 |
Scott Applequist
Messages: 4
Registered: November 1996
|
Junior Member |
|
|
I'm trying to plot contours of some irregularly spaced data
on a map of the eastern US. The values themselves are not a
problem, getting the grid to /overplot on the map is. All of
the following plots a map with + marks at each city.
�
;draw map
map_set, 15, -87, 0, limit=[24,-115,49,-67], /cont, /usa, /merc
;input data
zt = fltarr(3,179)
openr, 2, 'idlin.dat'
readf, 2, zt, format='(3f10.2)'
close, 2
z = transpose(zt)
x = z(*,2) ; longitude
y = z(*,1) ; latitude
c = x + y ; example data
;plot individual marks for each city
oplot, x, y, psym=1
Now as suggested in the online help, I use SPH_SCAT to
perform the spherical gridding.
r= sph_scat(x, y, c, bounds=[-115., 24., -67., 49.],gs=[0.5,0.5])
�
Now I can contour the array r, or even make a surface plot of it,
but I cannot get it to lay on top of the previously drawn map.
I have a feeling that the array coordinates must be shifted, but
I don't see how.�1 I'm using IDL. Version 4.0.1 (vms alpha).
Thanks for any help.
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
