The output of r.info is the following
r.info map=immagine@mapset
+----------------------------------------------------------------------------+
| Map: immagine@mapset Date: Tue Nov 30 12:51:12 2021 |
| Mapset: mapset Login of Creator: [...] |
| Location: WGS84_UTM32 |
| DataBase: C:\Users\[...]\grassdata |
| Title: immagine |
| Timestamp: none |
|----------------------------------------------------------------------------|
| |
| Type of Map: raster Number of Categories: 0 |
| Data Type: CELL |
| Rows: 10980 |
| Columns: 10980 |
| Total Cells: 120560400 |
| Projection: UTM (zone 32) |
| N: 4900020 S: 4790220 Res: 10 |
| E: 609780 W: 499980 Res: 10 |
| Range of data: min = 0 max = 18876 |
| |
| Data Description: |
| generato da r.in.gdal |
| |
| Comments: |
| r.in.gdal -e input="C:\Users\[...]\T32TNP_2020062\ |
| 2T100559_B04.jp2" output="immagine" memory=300 offset=0 num_digits=0 |
| |
+----------------------------------------------------------------------------+
Do you mean the sixth line of the file? From the i.atcorr help it is the "mean
target elevation above sea level [km]", and it is correct that it is negative.
I'm now wondering if the problem can be having put "0" for visibility (aerosol
model concentration, fifth line). What do you think about it?
-----Messaggio originale-----
Da: Moritz Lennert <[email protected]>
Inviato: giovedì 27 gennaio 2022 14:08
A: [email protected]; 'Maris Nartiss' <[email protected]>
Cc: 'GRASS user list' <[email protected]>
Oggetto: Re: R: [GRASS-user] NULL values as output of i.atcorr
On 27/01/22 14:01, [email protected] wrote:
> Dear Moritz,
> thank you for your reply.
> I tried with and without the -r flag and the results is the same.
I reformulate this part of the question then: what exactly is the content of
immagine@mapset ? Please provide the output of
r.info map=immagine@mapset
> Please find attached the parameters of the file 22giugno2020_b4.txt:
>
> 26
> 6 22 17.05 9.68728523 43.75514951
> 2
> 1
> 0
> -0.293
Are you sure this should be a negative number ? This should be visibility in km
IIRC.
> -1000
> 182
Moritz
>
> -----Messaggio originale-----
> Da: Moritz Lennert <[email protected]>
> Inviato: giovedì 27 gennaio 2022 13:22
> A: Ilaria Ferrando <[email protected]>; Maris Nartiss
> <[email protected]>
> Cc: GRASS user list <[email protected]>
> Oggetto: Re: [GRASS-user] NULL values as output of i.atcorr
>
> On 26/01/22 19:50, Ilaria Ferrando wrote:
>> Dear Māris,
>> Thank you for your reply.
>> Yes, the image is inside the computational region.
>> Any other ideas?
>
>
> 1) Are you sure your input data is reflectance (-r flag used) ?
> 2) Could you provide us with the contents of the 22giugno2020_b4.txt file ?
>
> Moritz
>
>>
>> Il mer 26 gen 2022, 19:45 Maris Nartiss <[email protected]
>> <mailto:[email protected]>> ha scritto:
>>
>> One thing to always check – does your computational region contain
>> data you are analysing? Try to run `g.region rast=immagine` before
>> running i.atcorr.
>>
>> Māris.
>>
>> trešd., 2022. g. 26. janv., plkst. 18:13 — lietotājs
>> <[email protected] <mailto:[email protected]>>
>> rakstīja:
>> >
>> > Dear all,
>> >
>> > I am having troubles in processing a Sentinel-2 image with i.atcorr.
>> >
>> > The command I use is the following:
>> >
>> > i.atcorr -r -input=immagine@mapset range=0,10000
>> parameters=C:\...\22giugno2020_b4.txt output=immagine_atcorr rescale=0,1
>> >
>> > and I obtain this output:
>> >
>> > wavelength less than 0.25 micron:
>> >
>> > let's take s(l)=s(0.25)
>> >
>> > Atmospheric correction...
>> >
>> > Atmospheric correction complete.
>> >
>> >
>> >
>> > The map I obtain is constituted by pixels with 6.669033e-34 values.
>> >
>> > Can anyone explain the meaning of the output and why I obtain a
>> NULL value corrected map?
>> >
>> > Thank you for your support.
>> >
>> >
>> >
>> > _______________________________________________
>> > grass-user mailing list
>> > [email protected] <mailto:[email protected]>
>> > https://lists.osgeo.org/mailman/listinfo/grass-user
>> <https://lists.osgeo.org/mailman/listinfo/grass-user>
>>
>>
>> _______________________________________________
>> grass-user mailing list
>> [email protected]
>> https://lists.osgeo.org/mailman/listinfo/grass-user
>
>
_______________________________________________
grass-user mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/grass-user
