[GRASS-dev] i.sentinel improvements

Luca Delucchi lucadeluge at gmail.com
Mon Dec 17 01:38:55 PST 2018


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

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

let me know.

-- 
ciao
Luca

www.lucadelu.org


More information about the grass-dev mailing list