<div dir="ltr"><div>Hi Roberta,</div><div><br></div><div>thank you for your effort. I tried i.sentinel.mask just now with one Sentinel-2 granule. I have considerably better results with Level 2A Sentinel product, topographically and atmospherically corrected from Level 1C by sen2cor with default settings.</div><div><br></div><div>Regarding strange results of i.topo.corr, I also had a problem few months ago. In the first pass, if the resolution of DEM raster differs from current region settings, you can create strange illumination model. The solution was to set region from the DEM (e.g. with 25 m resolution in case of EU-DEM), then to create illumination model with "i.topo.corr -i ..."), restore region settings (e.g. with 30 m for Landsat) and finally to apply topo correction in the second pass. <br></div><div>For some unknown reason, i.topo.corr did not change region settings automatically, despite you can find "The illumination model (cos_i) with flag -i uses the actual region as limits and the resolution of the elevation map." in its manual (Notes section).  <br></div><div><br></div><div>Best regards,</div><div><br></div><div>Ivan<br></div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">2018-06-14 13:08 GMT+02:00 Stefan Blumentrath <span dir="ltr"><<a href="mailto:Stefan.Blumentrath@nina.no" target="_blank">Stefan.Blumentrath@nina.no</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div link="blue" vlink="purple" lang="NO-BOK">
<div class="m_7046544305679768337WordSection1">
<p class="MsoNormal"><span lang="EN-GB">i.topo.corr performed quite well for us in the Norwegian mountains with minnaert method…<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">What options did you try in i.topo.corr?<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">Cheers<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">Stefan<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span lang="EN-US">From:</span></b><span lang="EN-US"> Roberta Fagandini <<a href="mailto:robifagandini@gmail.com" target="_blank">robifagandini@gmail.com</a>>
<br>
<b>Sent:</b> torsdag 14. juni 2018 13:01<br>
<b>To:</b> Stefan Blumentrath <<a href="mailto:Stefan.Blumentrath@nina.no" target="_blank">Stefan.Blumentrath@nina.no</a>><br>
<b>Cc:</b> Moritz Lennert <<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>><wbr>; GRASS developers list <<a href="mailto:grass-dev@lists.osgeo.org" target="_blank">grass-dev@lists.osgeo.org</a>></span></p><div><div class="h5"><br>
<b>Subject:</b> Re: [GRASS-dev] New GRASS addon i.sentinel.mask is now available (GSoC 2018 project)<u></u><u></u></div></div><p></p><div><div class="h5">
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
<div>
<p class="MsoNormal"><span lang="EN-GB">Hi Stefan!<u></u><u></u></span></p>
<div>
<p class="MsoNormal"><span lang="EN-GB">Certainly applying the topographic correction could be very useful to avoid misclassification. <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB">I did several tests with i.topo.corr but I always got strange results especially in areas with steeper terrain. <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB">I also tested the topographic correction tool implemented in arcsi [0] but again the results were not satisfactory at all. This is why I have not further investigated this issue within my PhD research.<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB">In order to remove misclassifications, e.g. due to topographic shadows, I implemented the shadow mask cleaning procedure. It intersects the rough shadow mask with the cloud mask that is shifted according to the sun position.
 In this way, all the areas identified as shadows without a corresponding cloud are removed from the final shadow mask.<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-GB">I'm still very interested in investigating the topographic correction so if you have any suggestion please, let me know!!<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB">Ciao,<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB">Roberta<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB">[0] </span><a href="https://www.arcsi.remotesensing.info/index.html" target="_blank"><span lang="EN-GB">https://www.arcsi.<wbr>remotesensing.info/index.html</span></a><span lang="EN-GB"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
<div>
<p class="MsoNormal"><span lang="EN-GB">2018-06-14 10:55 GMT+02:00 Stefan Blumentrath <</span><a href="mailto:Stefan.Blumentrath@nina.no" target="_blank"><span lang="EN-GB">Stefan.Blumentrath@nina.no</span></a><span lang="EN-GB">>:<u></u><u></u></span></p>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<div>
<div>
<p class="MsoNormal"><span lang="EN-GB">Hi Roberta,
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"> <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">What about topographic correction (i.topo.corr [1])?<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">Should that be run before as well?<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">I would assume it has some impact on shadow detection in steeper terrain…<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"> <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">Cheers<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">Stefan<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"> <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">1:
</span><a href="https://grass.osgeo.org/grass74/manuals/i.topo.corr.html" target="_blank"><span lang="EN-GB">https://grass.osgeo.org/<wbr>grass74/manuals/i.topo.corr.<wbr>html</span></a><span lang="EN-GB"><u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"> <u></u><u></u></span></p>
<p class="MsoNormal"><b><span lang="EN-US">From:</span></b><span lang="EN-US"> grass-dev <</span><a href="mailto:grass-dev-bounces@lists.osgeo.org" target="_blank"><span lang="EN-US">grass-dev-bounces@lists.<wbr>osgeo.org</span></a><span lang="EN-US">>
<b>On Behalf Of </b>Roberta Fagandini<br>
<b>Sent:</b> torsdag 14. juni 2018 10:29<br>
<b>To:</b> Moritz Lennert <</span><a href="mailto:mlennert@club.worldonline.be" target="_blank"><span lang="EN-US">mlennert@club.worldonline.be</span></a><span lang="EN-US">><br>
<b>Cc:</b> GRASS developers list <</span><a href="mailto:grass-dev@lists.osgeo.org" target="_blank"><span lang="EN-US">grass-dev@lists.osgeo.org</span></a><span lang="EN-US">><br>
<b>Subject:</b> Re: [GRASS-dev] New GRASS addon i.sentinel.mask is now available (GSoC 2018 project)</span><span lang="EN-GB"><u></u><u></u></span></p>
<div>
<div>
<p class="MsoNormal"><span lang="EN-GB"> <u></u><u></u></span></p>
<div>
<p class="MsoNormal"><span lang="EN-GB">Hi Moritz!<u></u><u></u></span></p>
<div>
<p class="MsoNormal"><span lang="EN-GB"> <u></u><u></u></span></p>
<div>
<p class="MsoNormal"><span lang="EN-GB">2018-06-14 10:03 GMT+02:00 Moritz Lennert <</span><a href="mailto:mlennert@club.worldonline.be" target="_blank"><span lang="EN-GB">mlennert@club.worldonline.be</span></a><span lang="EN-GB">><wbr>:<u></u><u></u></span></p>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<p class="MsoNormal"><span lang="EN-GB">On 13/06/18 17:35, Roberta Fagandini wrote:<u></u><u></u></span></p>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<p class="MsoNormal"><span lang="EN-GB">Hi all!<br>
The first version of i.sentinel.mask is now available and can be installed from the official svn repository.<br>
i.sentinel.mask is a module for cloud and shadow masks computation and it is the first step of my GSoC project.<u></u><u></u></span></p>
</blockquote>
<p class="MsoNormal"><span lang="EN-GB"><br>
Thank you, Roberta. Great job !<u></u><u></u></span></p>
</blockquote>
<div>
<p class="MsoNormal"><span lang="EN-GB"> <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB">Thank you! <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB"> <u></u><u></u></span></p>
</div>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<p class="MsoNormal"><span lang="EN-GB"><br>
One question: ISTR that you mentioned that the input to this module has to be atmospherically corrected images, i.e. land surface reflectance. Is this true ? If yes, then this needs to be mentioned in the man page.<u></u><u></u></span></p>
</blockquote>
<div>
<p class="MsoNormal"><span lang="EN-GB"> <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB">Yes, it's true! I'm going to add it. <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB"> <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal">Roberta<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="color:#888888"><br>
<br>
<span class="m_7046544305679768337gmail-m-7064033895571345447hoenzb">Moritz</span></span><u></u><u></u></p>
</blockquote>
</div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div></div></div>
</div>

<br>______________________________<wbr>_________________<br>
grass-dev mailing list<br>
<a href="mailto:grass-dev@lists.osgeo.org">grass-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/grass-dev</a><br></blockquote></div><br></div></div>