<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.hoenzb
{mso-style-name:hoenzb;}
span.E-MailFormatvorlage18
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 2.0cm 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Dear Carlos,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">I’m not an expert for projections.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">However, on Lat/Long WGS84 the actual area of cells decline from the equator towards the poles.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Thus, I would expect that cell values near the poles have “more weight” using Lat/Long WGS84 than using an equal area projection.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Near the poles I don’t understand how the values for extent and resolution should be correct (equal area), except there is a huge distortion (very likely for
a cylindrical projection)!<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Are there really 21600 cols with a nsres = 1178 m at the north and south pole. I would call this a huge distortion.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">As a test, I would calculate the statistics for a smaller area centered at the equator. I would expect that the results are very similar comparing the lat/long
and the reprojected dataset.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Regards Jörg<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b><span lang="DE" style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">Von:</span></b><span lang="DE" style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> grass-user [mailto:grass-user-bounces@lists.osgeo.org]
<b>Im Auftrag von </b>Carlos Grohmann<br>
<b>Gesendet:</b> Montag, 16. November 2015 23:32<br>
<b>Cc:</b> GRASS user list<br>
<b>Betreff:</b> Re: [GRASS-user] r.univar: different results with different projections?<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">Hello Cesar<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">That was weird, so I tested it again. The number of cells is the same for both projections, but the values differ. This must be related to reprojecting.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">To me, they shouldn't de different, since a nearest neighbor should preserve the original values. I'm not really comfortable with this, as I'm not sure I can trust the stats after projecting. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">best<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Carlos<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">GRASS 7.1.svn (latlong):~ > g.region raster=gdem_etopo1_ice -pa<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">projection: 3 (Latitude-Longitude)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">zone: 0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">datum: wgs84<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">ellipsoid: wgs84<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">north: 90N<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">south: 90S<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">west: 180W<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">east: 180E<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">nsres: 0:01<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">ewres: 0:01<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">rows: 10800<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cols: 21600<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cells: 233280000<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">GRASS 7.1.svn (latlong):~ > r.univar map=gdem_etopo1_ice -ge percentile=100<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">n=233280000<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">null_cells=0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cells=233280000<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">min=-10803<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">max=8333<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">range=19136<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">mean=-1892.40422534294<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">mean_of_abs=2644.91906490912<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">stddev=2649.98339302808<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">variance=7022411.98332463<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">coeff_var=-140.032629262802<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">sum=-441460057688<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">first_quartile=-4286<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">median=-2457<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">third_quartile=214<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">percentile_100=8333<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">GRASS 7.1.svn (eqarea):~ > g.region -p<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">projection: 99 (Equal Area Cylindrical)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">zone: 0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">datum: wgs84<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">ellipsoid: wgs84<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">north: 6363885.33192604<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">south: -6363885.33192604<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">west: -20037508.34278924<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">east: 20037508.34278924<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">nsres: 1178.49728369<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">ewres: 1855.32484655<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">rows: 10800<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cols: 21600<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cells: 233280000<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">GRASS 7.1.svn (eqarea):~ > r.univar map=gdem_etopo1_ice -ge percentile=100 <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">n=233280000<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">null_cells=0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cells=233280000<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">min=-10803<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">max=8333<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">range=19136<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">mean=-2382.28934158093<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">mean_of_abs=2845.10169015775<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">stddev=2508.93105538271<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">variance=6294735.0406638<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">coeff_var=-105.315966939504<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">sum=-555740457604<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">first_quartile=-4544<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">median=-3285<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">third_quartile=93<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">percentile_100=8333<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">On Mon, Nov 16, 2015 at 7:43 PM, César Augusto Ramírez Franco <<a href="mailto:caesarivs@gmail.com" target="_blank">caesarivs@gmail.com</a>> wrote:<o:p></o:p></p>
<div>
<div>
<p class="MsoNormal"><span style="color:#274E13">Carlos,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">2015-11-16 14:47 GMT-05:00 Carlos Grohmann <<a href="mailto:carlos.grohmann@gmail.com" target="_blank">carlos.grohmann@gmail.com</a>>:<o:p></o:p></p>
<div>
<div>
<div>
<p class="MsoNormal">GRASS 7.1.svn (base_maps):~ > g.region -p raster=gdem_etopo1_ice<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cells: 233280000<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">GRASS 7.1.svn (base_maps):~ > r.univar map=gdem_etopo1_ice -ge percentile=100<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cells=58320000<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">GRASS 7.1.svn (eqarea):~ > r.univar map=gdem_etopo1_ice -ge percentile=100<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cells=233280000<o:p></o:p></p>
</div>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><span style="color:#274E13">Notice how the number of pixels differs, that's the reason the statistics are not the same, I don't get why the region has a different number of pixels than the raster itself in the original latlong projection...
I think that's the root of the issue<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span class="hoenzb"><span style="color:#888888">-- </span><o:p></o:p></span></p>
<div>
<div>
<div>
<p class="MsoNormal"><b><span style="color:#888888">César Augusto Ramírez Franco</span></b><span style="color:#888888"><br>
Laboratorio de Sistemas Complejos Naturales<br>
Escuela de Geociencias - Facultad de Ciencias<br>
Universidad Nacional de Colombia - Sede Medellín<br>
Teléfono: (57-4) 430 9369 - 300 459 6085</span><o:p></o:p></p>
</div>
<p class="MsoNormal"><span style="color:#888888"><a href="http://labscn-unalmed.github.io/" target="_blank">http://labscn-unalmed.github.io/</a><o:p></o:p></span></p>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><br>
<br clear="all">
<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal">-- <o:p></o:p></p>
<div>
<div>
<p class="MsoNormal">Prof. Carlos Henrique Grohmann<br>
Institute of Energy and Environment - Univ. of São Paulo, Brazil<o:p></o:p></p>
<div>
<p class="MsoNormal">- Digital Terrain Analysis | GIS | Remote Sensing - <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><a href="http://carlosgrohmann.com/" target="_blank">http://carlosgrohmann.com</a><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><a href="http://orcid.org/0000-0001-5073-5572" target="_blank"><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#1155CC;background:white">http://orcid.org/0000-0001-5073-5572</span></a><o:p></o:p></p>
<div>
<p class="MsoNormal">________________<br>
Can’t stop the signal.<o:p></o:p></p>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>