Hi,
I am trying to plot thevertices of a polygon over an existing raster in an IMAGE() window. The snippetis as follows:
;-----------------------------------------
shp_and_raster = IMAGE(REFORM(ortho2[0,*,*]),GEOTIFF=gtiff_VI_l_shrinked,RGB_TABLE=8,/ORDER, $
SEMIMAJOR_AXIS=6378137.0d, SEMIMINOR_AXIS=6356752.314245d, ZONE=17,ELLIPSOID='WGS 84', GRID_UNITS="m", $
IMAGE_LOCATION=[0,754],$
IMAGE_DIMENSIONS=[new_ns*gtiff_VI_l_shrinked.MODELPIXELSCALETAG[0],new_nl*gtiff_VI_l_shrinked.MODELPIXELSCALETAG[1]],$
MAP_PROJECTION='UTM')
POLYFILL,(*RESULT.VERTICES)[0,*],(*result.VERTICES)[1,*],/data
;--------------------------------------------
The POLYFILL procedure opens a new IDL WINDOW instead of plotting it in the existing IMAGE() window. Unfortunately there’s no /overlay and /current keywords in thePOLYFILL procedure. Expecting a solution to this. Thanks,
Regards,
Puneeth
|