[GRASS-user] i.atcorr with Sentinel2

[Stefan Blumentrath]

Thanks again Sajid,

Tested now with your data and parameters and the command you send earlier:

r.in.gdal input=~/sen2data/elev_clip.tif output=elev –o
r.in.gdal input=~/sen2data/S2A_OPER_MSI_L1C_TL_MTI_20151012T101024457Z_T32TPR_B02_clip.tif output=S2A_OPER_MSI_L1C_TL_MTI_20151012T101024457Z_T32TPR_B02 -o
g.region -p raster=S2A_OPER_MSI_L1C_TL_MTI_20151012T101024457Z_T32TPR_B02
i.atcorr -r input=S2A_OPER_MSI_L1C_TL_MTI_20151012T101024457Z_T32TPR_B02 elevation=elev range=0,10000 output=test_atcorr rescale=0,10000 parameters=6Sparameters.txt –v --o

but I get only NULL values also for those data and settings and when running i.atcorr in GRASS 7.4.0RC1.
I stumbled over the following line in i.atcorr output (see full verbose output below):
*    view zenith angle:     0.00 deg  view azimuthal angle:         0.00 deg   *
Do you get that as well?

Given that the i.atcorr code did not change the last ½ year and given that you succeeded with current trunk as well (also with my data), I am afraid I have to conclude this is an issue of my local installation.

Does anyone know if an underlying library can be causing trouble or what parts of my setup I would have to look for?

Sajid, can you send me the output of your
grass75 --config
so I can build trunk just like you did…

Again, thanks for all your support on trouble shooting…

Cheers
Stefan

P.S.: Output from i.atcorr

* ****************************** 6s version 4.2b ****************************** *
*                        geometrical conditions identity                        *
*                        -------------------------------                        *
*                        sentinel2a observation                                 *
*                                                                               *
*    month: 10 day: 29                                                         *
*    solar zenith angle:   61.03 deg  solar azimuthal angle:      166.87 deg   *
*    view zenith angle:     0.00 deg  view azimuthal angle:         0.00 deg   *
*    scattering angle:    118.97 deg  azimuthal angle difference: 166.87 deg   *
*                                                                               *
*                        atmospheric model description                          *
*                        -----------------------------                          *
*            atmospheric model identity :                                       *
*                midlatitude winter  (uh2o=.853g/cm2,uo3=.395cm-atm)            *
*                                                                               *
*            aerosols type identity :                                           *
*                    Continental aerosols model                                 *
*                                                                               *
*            optical condition identity :                                       *
*                 visibility :   15.00 km  opt. thick. 550nm :   0.3158       *
*                                                                               *
*                        spectral condition                                     *
*                        ------------------                                     *
*            sentinel2a blue b2                                                 *
*                value of filter function :                                     *
*                 wl inf=    0.300 mic   wl sup=    2.600 mic                  *
*                                                                               *
*                        target type                                            *
*                        -----------                                            *
*            homogeneous ground                                                 *
*              constant reflectance over the spectra     0.000                 *
*                                                                               *
*                        target elevation description                           *
*                        ----------------------------                           *
*            ground pressure  [mb]        961.80                               *
*            ground altitude  [km]        0.450                                *
*                 gaseous content at target level:                              *
*                 uh2o=    0.711 g/cm2        uo3=    0.396 cm-atm             *
*                                                                               *
*                         atmospheric correction activated                      *
*                         --------------------------------                      *
Atmospheric correction...




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

Hi




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?

here it is, including the dem and parameters text file: (epsg:32632)

http://scientificpubs.com/data_share/sen2data.tar.gz


Sajid
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20171219/87b12d6a/attachment-0001.html>