<div dir="ltr"><div>Hi Luca,</div><br><div class="gmail_quote"><div dir="ltr">Il giorno lun 17 dic 2018 alle ore 10:39 Luca Delucchi <<a href="mailto:lucadeluge@gmail.com">lucadeluge@gmail.com</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi all,<br>
<br>
I would like to do some improvements in i.sentinel<br>
Since I'm not the main author I would like to get approval before to<br>
start to work on it.<br>
- i.sentinel.download add capability to filter out scenes with small<br>
percent (user could decide the percentage to keep) of area matching<br>
the area of interest<br>
- i.sentinel.import add capability to cut the scene on the footprint<br>
to simplify scenes patch<br>
- i.sentinel.preproc add capability to compute topographic correction<br>
using the following prcedure<br></blockquote><div><br></div><div>Of course, it's ok for me!</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
###<br>
eval `r.sunmask -sg ele=DTM year=${dat:0:4} month=$(echo $dat:4:2 |<br>
sed 's/^0*//') day=${dat:6:2} hour=${dat:9:2} minute=${dat:11:2}<br>
second=${dat:13:2} timezone=0`<br>
z=`echo "90 - $sunangleabovehorizon" | bc -l`<br>
echo "${inpmap} $z $sunazimuth"<br>
i.topo.corr -i basemap=DTM zenit=$z azimuth=$sunazimuth output=dtm_illu<br>
<br>
for m in `g.list type=raster pattern=${inpmap}*_cor`; do g.region<br>
raster=$m; r.mapcalc expr="${m}_d=double(${m})"; done<br>
<br>
i.topo.corr basemap=dtm_illu input=`g.list type=raster sep=comma<br>
pattern=${inpmap}*_cor_d` output=tcor zenith=$z method=c-factor<br>
<br>
g.remove -f type=raster pattern=${inpmap}*_cor_d<br>
###<br></blockquote><div><br></div><div>I don't have much experience with topographic correction but I'm not sure that c-factor is the best. When I tested it, I got anomalous results.</div><div><br></div><div>However, I would add the topographic correction as an option (e.g -flag "also perform  topo cor")</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>Yeah, that would be an option... What about keeping i.sentinel.import for the import part and make i.sentinel.preproc only about pre-processing, i.e., atcorr and topcorr?<br></blockquote><div><br></div><div>+1.<br></div><div> </div><div>Robi</div></div></div>