/nodata keyword in mapped image seems to impair the effect of LIMIT = [y1,x1,y2,x2]:

---
imj = IMAGE(prepdat, lon,lat, limit = limits , $
label_show = 0, RGB_TABLE = new_rgb, grid_units = 2, Map_Projection = "equirectangular")
--

will respect the limits. However,

---
imj = IMAGE(prepdat, lon,lat, limit = limits , $
label_show = 0, RGB_TABLE = new_rgb, grid_units = 2, Map_Projection = "equirectangular", /nodata)
---

produces a map image that is reduced in size, with blank space in the plot area. IMAGE_LOCATION does not rectify.

Don;t know if any one has noticed this but I spent most of this afternoon and at least £5 on coffee trying to figure out what was up : )

Craig Wallace