Overlap: Red and Green drawings from HDFLook, white drawing from IDL
Click on the image to enlarge.
| |
; IDL SHELL READ A GEOTIFF FILE ( LAMBERT AZIMUTHAL PROJECTION ) AND ADD A MAP
;
image = READ_TIFF('~/IDL/DATA_HDFLOOK/GEOTIFF_LAMBERT.tif', GEOTIFF = geotag, /VERBOSE)
sz=reverse(SIZE(image, /DIM))
limit_meters=DBLARR(8)
limit_meters[0:1]=[geotag.modeltiepointtag[3], geotag.modeltiepointtag[4]]
scalex=geotag.modelpixelscaletag[0]
limit_meters[2]=limit_meters[0]+sz[0]*scalex
limit_meters[3]=limit_meters[1]
limit_meters[4]=limit_meters[0]+sz[0]*scalex
limit_meters[5]=limit_meters[1]-sz[1]*scaley
limit_meters[6]=limit_meters[0]
limit_meters[7]=limit_meters[1]-sz[1]*scaley
window, xsize=500, ysize=500
tv, reverse(image, 3), true=1
print, 'Center Longitude', geotag.PROJCENTERLONGGEOKEY
print, 'Center Latitude', geotag.PROJCENTERLATGEOKEY
res = map_proj_init('lambert azimuthal', /gctp, center_latitude = geotag.PROJCENTERLATGEOKEY, center_longitude = geotag.PROJCENTERLONGGEOKEY, false_northing = 0, false_easting = 0, sphere_radius = geotag.GEOGSEMIMAJORAXISGEOKEY)
limit = map_proj_inverse(limit_meters, map_structure=res)
limit=reform(limit,2,4)
limit=REVERSE(limit)
map_set,geotag.PROJCENTERLATGEOKEY, geotag.PROJCENTERLONGGEOKEY, pos=[0,0,1,1], /lambert, limit=limit, /hires, /noborder, /noerase
map_grid, londel=1, latdel=1
map_continents, /hires, /coast