[GRASS-user] i.atcorr with Sentinel2

[Stefan Blumentrath]

Hi Micha,

Thanks for testing!
You would need either the release candidate for the upcoming version 7.4 or development version (7.5). GRASS 7.2 does not support Sentinel-2 in i.acorr…

Cheers
Stefan

From: Micha Silver [mailto:tsvibar at gmail.com]
Sent: tirsdag 19. desember 2017 21.58
To: Stefan Blumentrath <Stefan.Blumentrath at nina.no>; Sajid Pareeth <spareeth at gmail.com>
Cc: GRASS user list <grass-user at lists.osgeo.org>
Subject: Re: [GRASS-user] i.atcorr with Sentinel2

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:spareeth at gmail.com]
Sent: tirsdag 19. desember 2017 12.51
To: Stefan Blumentrath <Stefan.Blumentrath at nina.no><mailto:Stefan.Blumentrath at nina.no>
Cc: Moritz Lennert <mlennert at club.worldonline.be><mailto:mlennert at club.worldonline.be>; Žofie Cimburová <zoficimburova at gmail.com><mailto:zoficimburova at gmail.com>; GRASS user list <grass-user at lists.osgeo.org><mailto:grass-user at lists.osgeo.org>
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

grass-user at lists.osgeo.org<mailto:grass-user at lists.osgeo.org>

https://lists.osgeo.org/mailman/listinfo/grass-user



--

Micha Silver

Ben Gurion Univ.

Sde Boker, Remote Sensing Lab

cell: +972-523-665918
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20171219/12179807/attachment-0001.html>