Hi, I am using Envi 5.2/IDL8.4 and try to get used to ENVIRaster functionalities instead of the old reading and writing routines but struggle with some points. 1. When I try to use ENVISubsetRaster to get a subset with its spatial references (see code below), it shows spatialref as undefined in the variables list. According to the help article, this should work. It also works when I use "spatialref = in_sub.spatialref". Is it directly implemented in ENVISubetRaster or not? Also, the function "ENVISubsetRaster" is not coloured in my compiler but it seems to do the job, except for the spatialref thing. in_raster= e.OpenRaster(in_files[i]) ; open the raster in_sub= ENVISubsetRaster(in_raster, sub_rect=[10,10,199,199], spatialref=spatialref) 2. When I read in several input files (same extent etc) and want to read out spatialref and metadata for the output, in the way I do in the code below, it looses the information in the defined variables. Why is that? for i=0, nr_in_files-1 do begin in_raster= e.OpenRaster(in_files[i]) ; open the raster in_sub= ENVISubsetRaster(in_raster, sub_rect=[10,10,199,199]) ; get the subset incluseve spatial reference of it in_data= in_sub.GetData() ; read into an array if i eq 0 then begin nb= in_raster.nbands ns= in_raster.ncolumns nl= in_raster.nrows metadata = in_raster.metadata spatialref = in_sub.spatialref stack = replicate(0.0, nr_in_files, n_elements(in_data[*,0]), n_elements(in_data[0,*])) endif in_raster.close stack[i,*,*]=in_data endfor Cheers, Kim