[GRASS-user] Re: i.atcorr returns all NULL values

[Juan Benavides Duque]

Skipped content of type multipart/alternative-------------- next part --------------



##TOAR VALUES FOR LANDSAT TM5 1986 BAND1 (

Fri May 27 16:56:10 2011)                                                      
r.univar map=COL_W1_1986_atcor_1 at LANDSAT3                                       
(Fri May 27 16:56:16 2011) Command finished (6 sec)                             
(Fri May 27 16:56:30 2011)                                                      
r.univar map=COL_W1_1986_toar_1 at LANDSAT3                                        
total null and non-null cells: 53428941
total null cells: 16389576
Of the non-null cells:
----------------------
n: 37039365
minimum: -0.00298186
maximum: 0.331536
range: 0.334518
mean: 0.10724
mean of absolute values: 0.10724
standard deviation: 0.0530029
variance: 0.00280931
variation coefficient: 49.4247 %
sum: 3972086.9816170735
(Fri May 27 16:56:38 2011) Command finished (7 sec)                             
-------------- next part --------------

r.mask -r

PREFIX=COL_W1_1986_

g.region n=902700 s=695370 w=339300 e=571230 res=30


#CLOUD MASK PROCESS

#REMEMBER TO CHANGE THE FLAG FOR TM4 AQND TM5 (-t)


i.landsat.toar -s -t --overwrite input_prefix=${PREFIX} output_prefix=${PREFIX}toar_ metfile=L5009055_05519860908_MTL.txt sensor=tm5 


ICNDF=icnd_W1_1986_

FILES="${PREFIX}toar_1  ${PREFIX}toar_2  ${PREFIX}toar_3 ${PREFIX}toar_4 ${PREFIX}toar_5
       ${PREFIX}toar_7"

for k in $FILES ; do

 i=`echo $k | sed 's+'${PREFIX}toar_'++g'| sed 's+ ++g'` 

 echo "Processing ${PREFIX}toar_$i..."

i.atcorr -r -f -o --overwrite  input=$k elevation='RECLASS_SRTM3.1 at LANDSAT3' parameters=/home/dalevitt/grassdata/LANDSAT.RAW/LT50090551986251XXX04/${ICNDF}$i.txt output=${PREFIX}atcor_$i

 echo "ATMOSPHERIC  CORRECTION FILE  '${PREFIX}atcor_$i'"
done
############################################################################################################