[GRASS-dev] i.sentinel improvements

[Roberta Fagandini]

Hi Luca,

Il giorno lun 17 dic 2018 alle ore 10:39 Luca Delucchi <lucadeluge at gmail.com>
ha scritto:

> Hi all,
>
> I would like to do some improvements in i.sentinel
> Since I'm not the main author I would like to get approval before to
> start to work on it.
> - i.sentinel.download add capability to filter out scenes with small
> percent (user could decide the percentage to keep) of area matching
> the area of interest
> - i.sentinel.import add capability to cut the scene on the footprint
> to simplify scenes patch
> - i.sentinel.preproc add capability to compute topographic correction
> using the following prcedure
>

Of course, it's ok for me!


>
> ###
> eval `r.sunmask -sg ele=DTM year=${dat:0:4} month=$(echo $dat:4:2 |
> sed 's/^0*//') day=${dat:6:2} hour=${dat:9:2} minute=${dat:11:2}
> second=${dat:13:2} timezone=0`
> z=`echo "90 - $sunangleabovehorizon" | bc -l`
> echo "${inpmap} $z $sunazimuth"
> i.topo.corr -i basemap=DTM zenit=$z azimuth=$sunazimuth output=dtm_illu
>
> for m in `g.list type=raster pattern=${inpmap}*_cor`; do g.region
> raster=$m; r.mapcalc expr="${m}_d=double(${m})"; done
>
> i.topo.corr basemap=dtm_illu input=`g.list type=raster sep=comma
> pattern=${inpmap}*_cor_d` output=tcor zenith=$z method=c-factor
>
> g.remove -f type=raster pattern=${inpmap}*_cor_d
> ###
>

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.

However, I would add the topographic correction as an option (e.g -flag
"also perform  topo cor")


>
> 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?
>

+1.

Robi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20181219/9649214f/attachment.html>