<div dir="ltr"><div><div><div>It's been a while since I wrote that script. Glad to see it's still of use :)<br><br></div>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.<br><br></div>Cheers<br></div>Daniel <br></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Jan 23, 2018 at 11:50 AM Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com">markus.metz.giswork@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><br><br>On Tue, Jan 23, 2018 at 2:19 PM, Stefan Blumentrath <<a href="mailto:Stefan.Blumentrath@nina.no" target="_blank">Stefan.Blumentrath@nina.no</a>> wrote:<br>><br>> OK, will give it a try.<br>><br></div></div><div dir="ltr"><div>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.<br></div></div><div dir="ltr"><div>><br>> For the second part I am not sure if I properly understood where min – max wavelength for each band should be printed to.<br>><br>> Do you mean to replace<br>><br>>     # convert limits from nanometers to micrometers<br>><br>>     lowerlimit = wavelengths[0]/1000<br>><br>>     upperlimit = wavelengths[-1]/1000<br>><br>>  <br>><br>> Or just as a screen output?<br><br></div></div><div dir="ltr">Just a screen output to help update the manual</div><div dir="ltr"><br><div>><br>> You did not mean band value limits, right?<br></div><div><br></div></div><div dir="ltr"><div>No, e.g. in the manual is</div><div><br></div><div>Sentinel2A Coastal blue band B1 (430nm - 457nm)<br></div><div><br></div><div>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.<br></div></div><div dir="ltr"><div>><br>> Not sure I am of help/useful for fixing create_iwave.py, but I`ll do my best…<br></div><div><br></div></div><div dir="ltr"><div>Your help is welcome!</div><div><br></div><div>Markus M<br></div></div><div dir="ltr"><div><br></div><div>><br>> Cheers<br>><br>> Stefan<br>><br>>  <br>><br>>  <br>><br>> From: Markus Metz [mailto:<a href="mailto:markus.metz.giswork@gmail.com" target="_blank">markus.metz.giswork@gmail.com</a>]<br>> Sent: tirsdag 23. januar 2018 13.02<br>> To: Stefan Blumentrath <<a href="mailto:Stefan.Blumentrath@nina.no" target="_blank">Stefan.Blumentrath@nina.no</a>><br>> Cc: Markus Neteler <<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>>; GRASS user list <<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a>><br>><br>><br>> Subject: Re: [GRASS-user] i.atcorr with Sentinel2<br>><br>>  <br>><br>><br>><br>> On Tue, Jan 23, 2018 at 12:25 PM, Stefan Blumentrath <<a href="mailto:Stefan.Blumentrath@nina.no" target="_blank">Stefan.Blumentrath@nina.no</a>> wrote:<br>> ><br>> > No problem. Any chance I can contribute to maintenance of create_iwave.py?<br>> ><br>> > What needs to be done? Take your latest changes into account?<br>><br>> 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.<br>><br>> 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.<br>><br>> Markus M<br>><br>> ><br>> >  <br>> ><br>> > 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.<br>> ><br>> >  <br>> ><br>> > Cheers<br>> ><br>> > Stefan<br>> ><br>> >  <br>> ><br>> >  <br>> ><br>> >  <br>> ><br>> > From: grass-user [mailto:<a href="mailto:grass-user-bounces@lists.osgeo.org" target="_blank">grass-user-bounces@lists.osgeo.org</a>] On Behalf Of Markus Metz<br>> > Sent: tirsdag 23. januar 2018 12.03<br>> > To: Markus Neteler <<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>><br>> > Cc: GRASS user list <<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a>><br>> > Subject: Re: [GRASS-user] i.atcorr with Sentinel2<br>> ><br>> >  <br>> ><br>> ><br>> ><br>> > On Tue, Jan 23, 2018 at 11:51 AM, Markus Neteler <<br>> > <a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>> wrote:<br>> ><br>> > ><br>> > > On Tue, Jan 23, 2018 at 11:09 AM, Žofie Cimburová<br>> > > <<br>> > <a href="mailto:zoficimburova@gmail.com" target="_blank">zoficimburova@gmail.com</a>> wrote:<br>> ><br>> > > > Forgot to cc the list.<br>> > > ><br>> > > > The CSV and the cpp template:<br>> > > ><br>> > <a href="https://www.dropbox.com/sh/9cnhsl59kvmx4o1/AABwnEVHdxfW_PqrB5Xcpczda?dl=0" target="_blank">https://www.dropbox.com/sh/9cnhsl59kvmx4o1/AABwnEVHdxfW_PqrB5Xcpczda?dl=0</a><br>> ><br>> > ><br>> > > Great! Would you mind to also convert the S2A table? AFAIK it got<br>> > > updated due to recalibration compared to what we currently have.<br>> ><br>> > please wait a bit with creating the cpp template, create_iwave.py needs some maintenance<br>> ><br>> > Markus M<br>> ><br>> >  <br>> ><br>> > ><br>> > > markusN<br>> > > _______________________________________________<br>> > > grass-user mailing list<br>> > > <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>> > > <a href="https://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br><br></div></div>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a></blockquote></div>