Hello,
When using IMAGE() with the MAP_PROJECTION and GRID_UNITS properties set, the first row and column of my image are ignored. Is this a bug, or some subtlety of map projection that I'm missing?
Cheers,
A
e.g.
;; Define an image
im=BINDGEN(3,5,4)*21b
;; And plot it normally
i0 = IMAGE(im,AXIS_STYLE=2,IMAGE_DIMENSIONS=[100,80],IMAGE_LOCATION=[-50,-40],LAYOUT=[2,2,1])
;; Now add GRID_UNITS=2 to define lat/lon and MAP_PROJECTION keywords
i1 = IMAGE(im,AXIS_STYLE=2,IMAGE_DIMENSIONS=[100,80],IMAGE_LOCATION=[-50,-40],$
LAYOUT=[2,2,2],/CURRENT,GRID_UNITS=2,MAP_PROJECTION='Geographic')
;; The problem is also here for sinusoidal (and all other grids)
i1 = IMAGE(im,AXIS_STYLE=2,IMAGE_DIMENSIONS=[100,80],IMAGE_LOCATION=[-50,-40],$
LAYOUT=[2,2,3],/CURRENT,GRID_UNITS=2,MAP_PROJECTION='Sinusoidal')
;; Fix the issue by padding the original array
im_ = BYTARR( SIZE(/DIM,im)+[0,1,1])*0b + 255b
im_[*,1:*,1:*] = im
i3 = IMAGE(im_,AXIS_STYLE=2,IMAGE_DIMENSIONS=[100,80],IMAGE_LOCATION=[-50,-40],$
LAYOUT=[2,2,4],/CURRENT,GRID_UNITS=2,MAP_PROJECTION='Sinusoidal')
|