Hi All,
I have a pretty basic question.
I have a list of MODIS images in a folder, all beginning with the
prefix "MOD09GA" followed by a string of other characters, i.e.
'.A2009138.h09v05.005.2009140163311_Grid_2D_reprojected.img'
I have written a very basic little program to reorder the bands in the
MODIS image, but I've only gotten it to work one image at a time.
Ideally, I'd have a loop that takes each image from the folder,
reorders the bands, then writes each file with a unique name.
_________________________________
pro MOD09GA_bandreorder
ns=4578
nl=2367
nb=7
image = fltarr(ns,nl,nb)
morder=fltarr(ns,nl,nb)
;outfile=fltarr(ns,nl,nb)
for i=1,11 do begin ; is this opening multiple files, or the same
over and over?? figure it out AB!
openr, i, '/users/u0609216/documents/modis/2009/Geo_2009/MOD09GA.*'
readu, i, image
close, i
b1 = image(*,*,0)
b2 = image(*,*,1)
b3 = image(*,*,2)
b4 = image(*,*,3)
b5 = image(*,*,4)
b6 = image(*,*,5)
b7 = image(*,*,6)
; band order should be (3,4,1,2,5,6,7)
morder[*,*,0]=b3
morder[*,*,1]=b4
morder[*,*,2]=b1
morder[*,*,3]=b2
morder[*,*,4]=b5
morder[*,*,5]=b6
morder[*,*,6]=b7
print, "Now processing:" ,'file', i
get_lun,unit
openw,unit,'/users/u0609216/documents/modis/2009/reorder_200 9/
reorder.bip'
writeu,unit,morder
close,unit
endfor
print, 'Done!'
end
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
