| sec : U Re: Fitting a map to an image [message #31169 is a reply to message #31023] |
Thu, 06 June 2002 16:22  |
Andrew Cool
Messages: 219
Registered: January 1996
|
Senior Member |
|
|
Sean Raffuse wrote:
>
> Hello.
>
> I have an image of known size and flat geocoordinates. I want to overlay
> the US states vector in a display window on top of the image. My problem is
> that I cannot get the two to match up properly. There seems to be a border
> on the vector, but the image draws to the edge of the window. I tried using
> the noborder keyword on my call to map_set and also tried setting the !x and
> !y margins to zero. No luck. Can you help? Here is the relevant code. . .
I may as well throw in the example...
PRO USA
Device, decomposed=0
window
LoadCT, 40
!P.position=[0.1,0.1,0.9,0.9]
lowerlat = 25
upperleftlat = 50
upperleftlon = 360-140
rightlon = 360-60
map_set,40,-110,0,limit=[lowerlat,upperleftlon,upperleftlat, rightlon],/merc
image = DIST(800,250)
geo_warped_image = map_image(image,startx,starty,$
LATMIN = lowerlat,$
LATMAX = upperleftlat,$
LONMIN = upperleftlon,$
LONMAX = rightlon,$
COMPRESS = 1)
TV, BYTSCL(geo_warped_image), startx,starty,/device
map_continents, /usa, color = 200, mlinethick=2
MAP_GRID,/BOX,color = 200
END
------------------------------------------------------------ -----------------
Andrew D. Cool .->-.
Electromagnetics & Propagation Group '-<-'
Intelligence, Surveillance & Reconnaissance Division Transmitted on
Defence Science & Technology Organisation 100% recycled
PO Box 1500, Edinburgh electrons
South Australia 5111
Phone : 061 8 8259 5740 Fax : 061 8 8259 6673
Email : andrew.cool@crapdsto.defence.gov.au
Please remove the crap from my email address to reply ;-)
------------------------------------------------------------ -----------------
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
