Re: i.atcorr with Sentinel2

Micha Silver <[email protected]>
Newsgroups gmane.comp.gis.grass.user
Message-ID <[email protected]>
Dear Stefan:

I also tried your data, and failed. Any values for the geom (first
line in the 6S parameters) above 24 cause this error:

GRASS 7.2.2 (SB):~/work/tmp/Blumentrath > g.version

GRASS 7.2.2 (2017)

GRASS 7.2.2 (SB):~/work/tmp/Blumentrath > cat p6s.txt

25 # originally 25

9 5 10.6728 10.73413 59.90647

4

3

50

-0.055

-1000

165 # original 173

GRASS 7.2.2 (SB):~/work/tmp/Blumentrath > i.atcorr --v
-r --o input=S2A.B08 elev=dem output=test_atcorr param=p6s.txt
range=1,10000 rescale=1,10000

ERROR: Unsupported/unreadable format in control file (found
igeom=25)

If I put the igeom value at 24, and value for the wavelengthis above
164 (WV3 in the i.atcorr manpage) then the module segfaults.

??

Do I need a newer GRASS?

Best,

Micha

On 12/19/2017 04:24 PM, Stefan
Blumentrath wrote:

Dear Sajid,

Thanks so much for your reply and confirming
that the module works in principle!

Would you be able to provide me with the 1000 X
1000 clip of your data, so I can check if it is my
installation?

I tried also two different boxes and GRASS
versions (Ubunut 14.04 with GRASS 7.3 and 7.4 compiled from
source and Ubuntu 16.04 with GRASS 7.4 from UbuntuGIS) and
none of it worked.

Also, altering the terrain model from floating
point with negative values to integer without negative
values did not change a thing…

I still get only NULL pixels…

In case someone would be willing to test my
data (which I would be very grateful for), you can get it
directly here:

dem:
http://www.filedropper.com/dem_1

S2 band 8:
http://www.filedropper.com/s2aoperprdmsil1cpdmc20160907t044118r008v20160905t10402220160905t104245b08

62 params:
http://www.filedropper.com/p6s

I also have it on FTP and can share username
and password offlist…

Cheers

Stefan

From: Sajid Pareeth [ mailto:[email protected] ]

Sent: tirsdag 19. desember 2017 12.51

To: Stefan Blumentrath
<[email protected]>

Cc: Moritz Lennert
<[email protected]> ; Žofie Cimburová
<[email protected]> ; GRASS user list
<[email protected]>

Subject: Re: [GRASS-user] i.atcorr with Sentinel2

To add, I have used
the input as such with out applying the quantification
value.

>>

>>
----------------------------------------

>> 6s file:

>> 25- geometrical
conditions=Sentinel-2A

>> 9 5 10.6728 10.73413
59.90647- month day hh.ddd longitude
latitude

>> 4 - atmospheric
model=subarctic summer

>> 3- aerosol
model=urban

>> 50- visibility [km]
(aerosol model concentration) (this is

>> estimated, I also
tried with AOD specified)

>> -0.055- mean target
elevation above sea level [km]

>> -1000- sensor height

>> 166- Sentinel2A Blue
band B2 (440nm - 535nm)

Just noticed, the
last parameter should be 167 for Sentinel B2.

Sajid

_______________________________________________
grass-user mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/grass-user

--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

_______________________________________________
grass-user mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/grass-user
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.