| Re: Best routines for mapping satellite images [message #89353 is a reply to message #89352] |
Mon, 25 August 2014 16:15   |
chris_torrence@NOSPAM
Messages: 528
Registered: March 2007
|
Senior Member |
|
|
On Monday, August 25, 2014 1:13:01 PM UTC-6, Steve Super wrote:
> Does anyone have any suggestions for which routines are best for projecting satellite images? I have an un-projected image array, along with associated lat/lon arrays that I would like to map and then add some other layers/annotations.
>
>
>
> I have tried many different approaches (i.e., old IDL procedures, new IDL graphics functions, CG) with varying, but insufficient results.
>
>
>
> Thanks,
>
> Steve
I think it depends upon what you are trying to achieve. You could take the quick approach, and use IDL function graphics and see if that is what you want:
Assuming that your image coordinates are in "degrees" (since you said unprojected):
m = MAP( "Mollweide", LIMIT=[latmin, lonmin, latmax, lonmax] )
im = IMAGE(mydata, IMAGE_LOCATION=[lon0, lat0], $
IMAGE_DIMENSIONS=[widthLon, heightLat], $
GRID_UNITS="degrees", /OVERPLOT)
c = MAPCONTINENTS( )
m.Save, "map.pdf"
You will need to fill in appropriate values for latmin, etc.
Hope this helps.
-Chris
ExelisVIS
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
