[GRASS-dev] Re: [GRASS-SVN] r44252 - grass/branches/releasebranch_6_4/imagery/i.atcorr

Martin Landa landa.martin at gmail.com
Mon Nov 8 15:23:23 EST 2010


2010/11/8 Martin Landa <landa.martin at gmail.com>:
> when playing with i.atcorr, input raster maps are very small (450x450)
> in my case. Without `ialt` i.atcorr runs for few seconds, with this
> parameter it runs for ten minutes and still computes - is that normal?
> - sorry knowing very few about i.atcorr.

concretely

$ g.region rast=Brg1 -p
projection: 1 (UTM)
zone:       33
datum:      wgs84
ellipsoid:  wgs84
north:      5633223.50868776
south:      5605657.62851621
west:       497279.36486102
east:       525403.19584718
nsres:      60.05638382
ewres:      59.96552449
rows:       459
cols:       469
cells:      215271

$ i.landsat.toar -tr input_prefix=Brg output_prefix=B_rad
metfile=L4191025_02519840414_MTL.txt

$ g.region -lg
...
center_long=15.16068438
center_lat=50.72648033
...

$ r.univar --q map=demrg | grep mean:
mean: 552.711

$ cat atcorr.txt
7
04 14 9.28 15.16068045 50.72664675
0
3
50
-0.552
-1000
31 # it's MSS1

$ i.atcorr iimg=B_rad1 ialt=demrg icnd=atcorr.txt oimg=B_cor1 --o

Takes with `ialt` tens of minutes to finish computing...

Is that expected?

Martin

-- 
Martin Landa <landa.martin gmail.com> * http://geo.fsv.cvut.cz/~landa


More information about the grass-dev mailing list