|
I think what you are going to need here is to write out to a text file and then format the outputs accordingly. The formatting can be a bit tedious so the best approach is to print to the output console first to see how it looks and then print to the file. Try it with one image first to make each test quicker. Consider the following snippets from some code, please note that the xmap/ymap may or may not exactly apply for your situation:
;Set up an output *.txt file to print the results to...
out_file =dialog_pickfile(/write, filter='*.txt', /default_extension)
if (out_file eq '') then return
openw, lun, out_file, /get_lun
;Now go to each point (lon, lat)
for i=0, n_elements(xf)-1 do begin
print,''
printf,lun,''
print,'Location',[i+1],'Lon/Lat:',ixmap[i], iymap[i]
printf,lun,'Location',[i+1],' Lon/Lat:',ixmap[i], iymap[i]
;Get the pixel values at each point
for k=0, nb-1 do begin
dims = ([0 ,xf[i] , xf[i] , yf[i] , yf[i]])
p_value= envi_get_data(fid=fid, pos=pos[k], dims=dims)
print,k+1,p_value, format='("Band number:",i5," ","Pixel value: ",3f0)'
printf,lun,k+1,p_value, format='("Band number:",i5," ","Pixel value: ",3f0)'
endfor ; for pixel values
endfor ;for xf file coordinates
|