HI~ All ~ I am trying to a image in Polar Stereographic projection but
it is not in the right position. My code is like this: And when I
tried to change the values of lat_min, lon_min ,lat_max and lon_max,
it does not help much. And I think i have to rotate the map or the
data with some angle. How can I full that ?
Thank you very much~
lat_min=-35.3
lat_max=-26.5
lon_min=114.25
lon_max=122.45
;defind the max and min latitude and longtitude of the data area
map=MAP_PROJ_INIT('Polar Stereographic',/
GCTP,center_longitude=118.35,center_latitude=-30.9,limit=[la t_min,lon_min,lat_max,lon_max])
r= MAP_PROJ_FORWARD([lon_min,lon_max],[lat_min,lat_max],/
fill,MAP_STRUCTURE=map)
;apply the Polar Stereographic projection
result=map_proj_image(a1,
[lon_min,lat_min,lon_max,lat_max],MAP_STRUCTURE=map, missing=0)
;display the data to the map
pos = [0.05, 0.05, 0.95, 0.95]
tvimage,result, position=pos,/KEEP_ASPECT,/tv,order=1
Plot, [r[0,0],r[0,1]], [r[1,0], r[1,1]], Position=pos,/Nodata,
XStyle=5, YStyle=5, /NoErase
MAP_GRID,/label,linestyle=0, map_structure=map,
color=FSC_Color('Blue'),GLINETHICK=2,CHARSIZE=2,charthick=2
; Draw gridlines over the map
MAP_CONTINENTS, map_structure=map,
color=FSC_Color('Blue'),MLINETHICK=2
; Draw continent outlines:
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
