<html><head></head><body>Ciao Nikos,<br>The pdf is available at the web link provided in the email. Thanks for giving a look.<br>Cheers,<br>Ghislain<br><br><div class="gmail_quote">Le 6 février 2019 15:19:22 GMT+01:00, Nikos Alexandris <nik@nikosalexandris.net> a écrit :<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<pre class="k9mail">* Ghislain Vieilledent | CIRAD <ghislain.vieilledent@cirad.fr> [2019-02-06 14:11:04 +0100]:<br><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 1ex 0.8ex; border-left: 1px solid #729fcf; padding-left: 1ex;">Dear GRASS GIS user,<br><br>I am using r.resamp.rst to interpolate climatic data at 30sec <br>resolution over Africa.<br><br>Grass script, data, and results are available here:<br><a href="https://nextcloud.fraisedesbois.net/index.php/s/CKzsZCSCqMTY9rq">https://nextcloud.fraisedesbois.net/index.php/s/CKzsZCSCqMTY9rq</a><br><br>#===================<br><br># Import Africa shapefile<br>v.in.ogr --overwrite input=gisdata/vectors/Africa layer=Africa output=Africa<br><br># Rasterize Africa at 30 arc-sec resolution<br>g.region -ap vector=Africa res=0:00:30<br>v.to.rast input=Africa type=area output=Africa use=val value=1<br><br># Import climate<br>r.in.gdal --overwrite input=gisdata/rasters/test1.1_044.tif output=test1_1<br><br># Resample with RST<br>g.region -ap raster=test1_1<br># Note: use maskmap=Africa to save computation time<br>r.resamp.rst --overwrite input=test1_1 elevation=test1_1_rst \\<br> ew_res=0.00833 ns_res=0.00833 maskmap=Africa<br>r.info test1_1_rst<br><br># Resample to exact resolution 0:00:30<br>g.region -ap res=0:00:30<br>r.resample --overwrite input=test1_1_rst output=test1_1_30s<br><br># Export<br>r.out.gdal --overwrite input=test1_1_30s \\<br> output=output/test1_1_30s.tif type=Float32 \\<br> createopt='compress=lzw,predictor=2'<br><br>#===================<br><br>When I use r.info on the interpolated raster, it indicates a tension <br>of 7278 in the comments which seems very high to me as the default <br>value is 40.<br><br>| Comments:<br>| tension=7278.552524<br>| dnorm=5.495598, zmult=1.000000<br>| KMAX=50, KMIN=35, errtotal=0.000535<br>| zmin_data=0.056488, zmax_data=3.756866<br>| zmin_int=-0.066779, zmax_int=3.214336<br><br>On the other side, results do not look so bad (see test.pdf figure).<br></blockquote><br>I am following this thread out of curiosity/interest on resampling. The<br>.pdf file is not attached. How big is it?<br><br>Nikos<br><br><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 1ex 0.8ex; border-left: 1px solid #729fcf; padding-left: 1ex;">Would you have any comment on this high value for the tension and do <br>you think interpolation has been performed correctly ?<br><br>Best regards,<br><br>Ghislain<br></blockquote></pre></blockquote></div><br>-- <br>Envoyé de mon appareil Android avec Courriel K-9 Mail. Veuillez excuser ma brièveté.</body></html>