<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="NO-BOK" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">OK, will give it a try.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">For the second part I am not sure if I properly understood where min – max wavelength for each band should be printed
to. <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Do you mean to replace<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"> # convert limits from nanometers to micrometers<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"> lowerlimit = wavelengths[0]/1000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"> upperlimit = wavelengths[-1]/1000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Or just as a screen output?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">You did not mean band value limits, right?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Not sure I am of help/useful for fixing create_iwave.py, but I`ll do my best…<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Cheers<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Stefan<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif"> Markus Metz [mailto:markus.metz.giswork@gmail.com]
<br>
<b>Sent:</b> tirsdag 23. januar 2018 13.02<br>
<b>To:</b> Stefan Blumentrath <Stefan.Blumentrath@nina.no><br>
<b>Cc:</b> Markus Neteler <neteler@osgeo.org>; GRASS user list <grass-user@lists.osgeo.org><br>
<b>Subject:</b> Re: [GRASS-user] i.atcorr with Sentinel2<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span lang="EN-GB"><br>
<br>
On Tue, Jan 23, 2018 at 12:25 PM, Stefan Blumentrath <</span><a href="mailto:Stefan.Blumentrath@nina.no"><span lang="EN-GB">Stefan.Blumentrath@nina.no</span></a><span lang="EN-GB">> 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?<o:p></o:p></span></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span lang="EN-GB">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.<o:p></o:p></span></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span lang="EN-GB">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.<o:p></o:p></span></p>
</div>
<p class="MsoNormal"><span lang="EN-GB">Markus M<o:p></o:p></span></p>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span lang="EN-GB">><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:</span><a href="mailto:grass-user-bounces@lists.osgeo.org"><span lang="EN-GB">grass-user-bounces@lists.osgeo.org</span></a><span lang="EN-GB">] On Behalf Of Markus Metz<br>
> Sent: tirsdag 23. januar 2018 12.03<br>
> To: Markus Neteler <</span><a href="mailto:neteler@osgeo.org"><span lang="EN-GB">neteler@osgeo.org</span></a><span lang="EN-GB">><br>
> Cc: GRASS user list <</span><a href="mailto:grass-user@lists.osgeo.org"><span lang="EN-GB">grass-user@lists.osgeo.org</span></a><span lang="EN-GB">><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>
> </span><a href="mailto:neteler@osgeo.org"><span lang="EN-GB">neteler@osgeo.org</span></a><span lang="EN-GB">> wrote:<br>
><br>
> ><br>
> > On Tue, Jan 23, 2018 at 11:09 AM, Žofie Cimburová<br>
> > <<br>
> </span><a href="mailto:zoficimburova@gmail.com"><span lang="EN-GB">zoficimburova@gmail.com</span></a><span lang="EN-GB">> wrote:<br>
><br>
> > > Forgot to cc the list.<br>
> > ><br>
> > > The CSV and the cpp template:<br>
> > ><br>
> </span><a href="https://www.dropbox.com/sh/9cnhsl59kvmx4o1/AABwnEVHdxfW_PqrB5Xcpczda?dl=0"><span lang="EN-GB">https://www.dropbox.com/sh/9cnhsl59kvmx4o1/AABwnEVHdxfW_PqrB5Xcpczda?dl=0</span></a><span lang="EN-GB"><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>
> > </span><a href="mailto:grass-user@lists.osgeo.org"><span lang="EN-GB">grass-user@lists.osgeo.org</span></a><span lang="EN-GB"><br>
> > </span><a href="https://lists.osgeo.org/mailman/listinfo/grass-user"><span lang="EN-GB">https://lists.osgeo.org/mailman/listinfo/grass-user</span></a><span lang="EN-GB"><o:p></o:p></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>