<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 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 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.EmailStyle17
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
span.pl-k
{mso-style-name:pl-k;}
span.pl-s
{mso-style-name:pl-s;}
span.pl-pds
{mso-style-name:pl-pds;}
span.pl-smi
{mso-style-name:pl-smi;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 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="NO-BOK" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Hei,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">What about something like this (example converts from UTM33 to WGS84 and creates ):<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">r.stats -gn1 INPUTMAP | cs2cs -E -f "%.6f" +proj=utm +no_defs +zone=33 +a=6378137 +rf=298.257222101 +towgs84=0.000,0.000,0.000
+to_meter=1 +to +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs | tr '\t' ' ' | r.in.xyz --o --v z=6 input=- output=REPROJECTED_MAP separator=space<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Or is cs2cs to slow or does it modify also z?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">BTW, I think several modules would benefit from a -b do not build topology...<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Cheers<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Stefan<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif"> grass-dev [mailto:grass-dev-bounces@lists.osgeo.org]
<b>On Behalf Of </b>Moritz Lennert<br>
<b>Sent:</b> 30. november 2016 19:36<br>
<b>To:</b> Markus Metz <markus.metz.giswork@gmail.com><br>
<b>Cc:</b> grass-dev <grass-dev@lists.osgeo.org><br>
<b>Subject:</b> Re: [GRASS-dev] how to reproject a raster map with absolute numbers without losing data<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
<br>
<span lang="EN-US">Le 29 novembre 2016 22:33:39 GMT+01:00, Markus Metz <</span><a href="mailto:markus.metz.giswork@gmail.com"><span lang="EN-US">markus.metz.giswork@gmail.com</span></a><span lang="EN-US">> a écrit :<br>
> For<br>
>reprojection, something like </span><a href="http://r.in"><span lang="EN-US">r.in</span></a><span lang="EN-US">.xyz could work:<br>
<br>
I'm currently going through this process and it raises a series of issues:<br>
<br>
<br>
><br>
>1. convert raster cells to vector points with raster cell value as<br>
>vector z<br>
>value.<br>
<br>
</span>I use <a href="http://r.to">r.to</a>.vect for this. The -b flag allows +/- fast conversion, but it takes into account all cells. It would be great to have a flag telling it to only convert non-null cells. Or at least a hint to set a mask, which seems
to have the desired effect.<br>
<br>
<br>
>2. reproject these vector points to the target location<br>
<br>
I locally modified v.proj to add a 'Don't build topology' flag.<br>
<br>
Is there any indication against this ?<br>
<br>
>3. export the vector points with v.out.ascii format=point<br>
>4. import the exported points with <a href="http://r.in">r.in</a>.xyz method=sum<br>
<br>
This works nicely.<br>
<br>
Thanks ! If I find the time I will wrap this into a r.proj.aggregate (?) script...<br>
<br>
<br>
<br>
><br>
>The method options of <a href="http://r.in">r.in</a>.xyz could be added to r.proj for spatial<br>
>aggregation, but that would mean that most of the code of <a href="http://r.in">
r.in</a>.xyz<br>
>would<br>
>need to be copied to r.proj, then adapted. That would be an interesting<br>
>enhancement, also considering that gdalwarp offers as resampling<br>
>methods<br>
>near, bilinear, cubic, cubicspline, lanczos, average, mode, max, min,<br>
>med,<br>
>Q1, Q3.<br>
<br>
But it doesn't offer sum....<br>
<br>
Moritz<o:p></o:p></p>
</div>
</body>
</html>