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

[Yann Chemin]

Hi Martin
I am going to office now
Can you send me your input files and script
Thank you
Yann

Yann Chemin
www.csu.edu.au

On 09/11/2010, at 7:23, Martin Landa <landa.martin at gmail.com> wrote:

> 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