| MAP_SET [message #4468] |
Wed, 21 June 1995 00:00 |
dean
Messages: 55
Registered: March 1993
|
Member |
|
|
I was just reading about the new enhanced mapping features with IDL 4.0. I
have been using IDL for several years and I am still unable to add a map to
a satellite image with one (or two) points off the globe.
Is there a way to add a map? Can you help me? Here is a sample to read in
a GOES-8 image located on anonymous FTP //castor.cira.colostate.edu/dean
called nhemi.img. I provided a mapset with the corner information and sub-
satellite point. One point is off the globe.
I also tried adding a map with just the 4 limit values. The map overlay is
to small.
If anyone has an idea to get this to work - let me know.
I am using IDL 4.0 for Window NT.
Kelly Dean
CSU/CIRA
==================================== DO_THIS.PRO============================
PRO DO_THIS
;
; Read raw image. Located at anonymous ftp//castor.cira.colostate.edu/dean.
;
xsize = 1024
ysize = 768
OPENR, unit, 'jerbe:NHEMI.IMG', /GET_LUN
image = bytarr(xsize,ysize)
scan = bytarr(xsize)
FOR J = 0,ysize-1 DO BEGIN
readu, unit, scan
image(0,(ysize-1)-J) = scan
ENDFOR
close, unit
FREE_LUN, unit
;
; Display image.
;
WINDOW, 0, XSIZE=xsize, YSIZE=ysize, Title = 'Map Test'
TV, image
;
; Add map to image.
;
Lat0 = 0.000000
Lon0 = -75.5000
; NW_lat = -9999.99
; NW_lon = -9999.99
; NE_lat = 51.2992
; NE_lon = -56.5387
; SE_lat = 13.9659
; SE_lon = -64.2103
; SW_lat = 14.3917
; SW_lon = -111.595
North = 51.2524
East = -62.5701
South = 13.9614
West =-118.715
map_set, lat0, lon0, 0, $
sat_p = [6.6107650,0,0], $
/Satellite, $
/USA, $
/GRID, $
/NOERASE, $
/NOBORDER, $
Xmargin = [0,0], $
Ymargin = [0,0], $
Limit = [South, West, North, East]
; Limit = [NW_lat, NW_lon, $
; NE_lat, NE_lon, $
; SE_lat, SE_lon, $
; SW_lat, SW_lon ]
END
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
