Re: i.atcorr with Sentinel2

Thanks Daniel, I might come back to your offer!
Cheers
Stefan


From: Daniel Victoria [mailto:[email protected]]
Sent: tirsdag 23. januar 2018 15.04
To: Markus Metz <[email protected]>
Cc: Stefan Blumentrath <[email protected]>; GRASS user list <[email protected]>
Subject: Re: [GRASS-user] i.atcorr with Sentinel2

It's been a while since I wrote that script. Glad to see it's still of use :)
Anyway, I'm a bit overwhelmed so I won't be able to look at the changes asked. But I'm glad Stefan can take a look. And I can help out in any questions that might come up.
Cheers
Daniel

On Tue, Jan 23, 2018 at 11:50 AM Markus Metz <[email protected]<mailto:[email protected]>> wrote:


On Tue, Jan 23, 2018 at 2:19 PM, Stefan Blumentrath <[email protected]<mailto:[email protected]>> wrote:
>
> OK, will give it a try.
>
The first part is more important because it 1) fixes the interpolation to 2.5 nm steps, and 2) helps iwave.cpp to load the response values.
>
> For the second part I am not sure if I properly understood where min – max wavelength for each band should be printed to.
>
> Do you mean to replace
>
>     # convert limits from nanometers to micrometers
>
>     lowerlimit = wavelengths[0]/1000
>
>     upperlimit = wavelengths[-1]/1000
>
>
>
> Or just as a screen output?
Just a screen output to help update the manual

>
> You did not mean band value limits, right?

No, e.g. in the manual is

Sentinel2A Coastal blue band B1 (430nm - 457nm)

and this range could be printed out to screen based on the actual response. This range can then be copied to the manual, making life a bit easier for the person updating the manual.
>
> Not sure I am of help/useful for fixing create_iwave.py, but I`ll do my best…

Your help is welcome!

Markus M

>
> Cheers
>
> Stefan
>
>
>
>
>
> From: Markus Metz [mailto:[email protected]<mailto:[email protected]>]
> Sent: tirsdag 23. januar 2018 13.02
> To: Stefan Blumentrath <[email protected]<mailto:[email protected]>>
> Cc: Markus Neteler <[email protected]<mailto:[email protected]>>; GRASS user list <[email protected]<mailto:[email protected]>>
>
>
> Subject: Re: [GRASS-user] i.atcorr with Sentinel2
>
>
>
>
>
> On Tue, Jan 23, 2018 at 12:25 PM, Stefan Blumentrath <[email protected]<mailto:[email protected]>> wrote:
> >
> > No problem. Any chance I can contribute to maintenance of create_iwave.py?
> >
> > What needs to be done? Take your latest changes into account?
>
> nodata handling should be changed. In read_input(), missing values should be set to zero. Then for each band, leading and trailing zeros would be set to -1. This helps interpolate_band() to get the correct subrange.
>
> For each band, create_iwave.py should print out the band name and the wavelength as min - max. Care must be taken when finding reasonable min and max wavelengths: find the max response, then go back while response > threshold to find the min wavelength. Go forward while response > threshold to find the max wavelength. As threshold I suggest 0.1 to eliminate noise in the response values.
>
> Markus M
>
> >
> >
> >
> > If you think it makes sense I will try to create a patch for adding S2B to i.atcorr (would be a nice exercise and does not seem too complex)? Will open a ticket and start working on it unless you tell me different.
> >
> >
> >
> > Cheers
> >
> > Stefan
> >
> >
> >
> >
> >
> >
> >
> > From: grass-user [mailto:[email protected]<mailto:[email protected]>] On Behalf Of Markus Metz
> > Sent: tirsdag 23. januar 2018 12.03
> > To: Markus Neteler <[email protected]<mailto:[email protected]>>
> > Cc: GRASS user list <[email protected]<mailto:[email protected]>>
> > Subject: Re: [GRASS-user] i.atcorr with Sentinel2
> >
> >
> >
> >
> >
> > On Tue, Jan 23, 2018 at 11:51 AM, Markus Neteler <
> > [email protected]<mailto:[email protected]>> wrote:
> >
> > >
> > > On Tue, Jan 23, 2018 at 11:09 AM, Žofie Cimburová
> > > <
> > [email protected]<mailto:[email protected]>> wrote:
> >
> > > > Forgot to cc the list.
> > > >
> > > > The CSV and the cpp template:
> > > >
> > https://www.dropbox.com/sh/9cnhsl59kvmx4o1/AABwnEVHdxfW_PqrB5Xcpczda?dl=0
> >
> > >
> > > Great! Would you mind to also convert the S2A table? AFAIK it got
> > > updated due to recalibration compared to what we currently have.
> >
> > please wait a bit with creating the cpp template, create_iwave.py needs some maintenance
> >
> > Markus M
> >
> >
> >
> > >
> > > markusN
> > > _______________________________________________
> > > grass-user mailing list
> > > [email protected]<mailto:[email protected]>
> > > https://lists.osgeo.org/mailman/listinfo/grass-user
_______________________________________________
grass-user mailing list
[email protected]<mailto:[email protected]>
https://lists.osgeo.org/mailman/listinfo/grass-user

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