Hi Luca, Il giorno lun 17 dic 2018 alle ore 10:39 Luca Delucchi <lucadel...@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
_______________________________________________ grass-dev mailing list grass-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-dev