<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<p>Paul,</p>
<p>Ik gebruik voor conversie van KNMI naar RD de volgende argumenten
in GDALWARP (niet gdal_translate):</p>
<p>D:\gdal64\bin\gdalwarp -of netCDF -dstnodata -999 -te 9000 29900
296000 651000 -tr 1000 1000 -s_srs "+proj=stere +x_0=0 +y_0=0
+lat_0=90 +lon_0=0 +lat_ts=60 +a=6378.388 +b=6356.906" -t_srs
"EPSG:28992" <br>
</p>
Ik zie wat kleine verschillen met jouw projectiestring.<br>
<br>
mvg,<br>
Siebe Bosch<br>
<br>
<br>
<div class="moz-cite-prefix">On 12-12-16 12:17, Paul Meems wrote:<br>
</div>
<blockquote
cite="mid:CAHNf2YQs5HvRQAyU5QkAUcrKjZk-CJQxbHNFPmFDFmz2KViK+Q@mail.gmail.com"
type="cite">
<div dir="ltr">
<div style="font-size:12.8px">Ik probeer om de 3-uur
neerslagaccumulaties van het KNMI (<a moz-do-not-send="true"
href="https://data.knmi.nl/datasets/radar_tar_corr_accum_03h/2.0"
target="_blank">https://data.knmi.nl/<wbr>datasets/radar_tar_corr_accum_<wbr>03h/2.0</a>)
te gebruiken.</div>
<div style="font-size:12.8px"><br>
</div>
<div style="font-size:12.8px">Al onze geo-data is in WGS84, dus
mijn eerste stap is het converteren van het HDF5-bestand van
KNMI naar een andere projectie en dan gelijk ook naar GeoTiff
wat toch wat makkelijker te verwerken is.</div>
<div style="font-size:12.8px">Hier gebruiken we de GDAL tools
voor, v2.0.</div>
<div style="font-size:12.8px">Als ik rechtstreeks GDAL_Warp
gebruik krijg ik een foutmelding over ontbrekende
transformatie matrix. Via Google heb ik soortgelijke problemen
gevonden en daar was de oplossing om eerst GDAL_Translate te
doen en dan de <i>georeferenced bounds</i> toe te voegen.</div>
<div style="font-size:12.8px"><br>
</div>
<div style="font-size:12.8px">Dat doe ik moet dit commando:</div>
<div style="font-size:12.8px">
<pre class="gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-java" style="white-space:pre-wrap;margin-top:0px;margin-bottom:0px;padding:0px;max-height:30em;overflow:auto;word-wrap:normal;color:rgb(51,51,51);font-size:12px;background-color:rgb(245,245,245)">gdal_translate -of GTiff -a_nodata 65535 -a_srs <span class="gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote" style="color:rgb(0,145,0)">"+proj=stere +lat_0=90 +lon_0=0 +lat_ts=60 +a=6378.14 +b=6356.75 +x_0=0 y_0=0"</span> -a_ullr 0 55.973602 9.0092793 48.895298 <span class="gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote" style="color:rgb(0,145,0)">"HDF5:\"</span>RAD_NL25_RAC_03H_20161204<wbr>1200.h5\<span class="gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote" style="color:rgb(0,145,0)">":<span class="gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-comment" style="color:rgb(128,128,128)">//image1/image_data"</span> <span class="gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote">"test-translate.tif"</span></span></pre></div><div style="font-size:12.8px">
<div><div class="gmail-m_6090844043269395631gmail-m_4197702424755525384gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">-a_ullr moet dan de bounds zijn. Die haal ik uit de meta-data van het bestand.</span></p><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">De projectiestring komt ook uit de meta-data.</span></p><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">
</span></p><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">Vervolgens kan ik dan een reprojectie doen met:</span></p><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap"></span></p><pre class="gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-java" style="white-space:pre-wrap;margin-top:0px;margin-bottom:0px;padding:0px;max-height:30em;overflow:auto;word-wrap:normal;color:rgb(51,51,51);font-size:12px;background-color:rgb(245,245,245)">gdalwarp -overwrite -s_srs <span class="gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote" style="color:rgb(0,145,0)">"+proj=stere +lat_0=90 +lon_0=0 +lat_ts=60 +a=6378.14 +b=6356.75 +x_0=0 y_0=0"</span> -t_srs EPSG:3857 -dstnodata 65535 -of GTiff <span class="gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote" style="color:rgb(0,145,0)">"test-translate.tif"</span> <span class="gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote" style="color:rgb(0,145,0)">"test-warp.tif"</span></pre><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">
</span></p><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">Het resultaat had dan in NL moeten liggen met coördinaten van 52 07' 25'', 5 15' 46'' (Amersfoort)</span></p><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">Maar ik krijg 89 29' 23'', 173 44' 12''</span></p><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">
</span></p><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">Volgens mij zijn er drie mogelijke oorzaken:</span></p><ol><li style="margin-left:15px"><font face="arial"><span style="font-size:12.6667px;white-space:pre-wrap">Ik doe de conversie en reprojectie niet goed</span></font></li><li style="margin-left:15px"><font face="arial"><span style="font-size:12.6667px;white-space:pre-wrap">Ik geef niet de juiste bounds op in -a_ullr</span></font></li><li style="margin-left:15px"><font face="arial"><span style="font-size:12.6667px;white-space:pre-wrap">Ik gebruik niet de juiste projectiestring</span></font></li></ol><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">Ik ben vast niet de eerste die dit wil doen, weet iemand de juiste workflow?</span></p><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">
</span></p><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;vertical-align:baseline;white-space:pre-wrap">Alvast bedankt.</span></p></div></div></div></div></div></div></div></div></div><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><span><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;color:rgb(34,34,34);background-color:rgb(255,255,255);font-weight:400;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;text-decoration:none;vertical-align:baseline;white-space:pre-wrap">
</span></p><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;color:rgb(34,34,34);background-color:rgb(255,255,255);font-weight:400;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;text-decoration:none;vertical-align:baseline;white-space:pre-wrap">Met vriendelijke groet,</span></p><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;color:rgb(34,34,34);font-weight:400;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;text-decoration:none;vertical-align:baseline;white-space:pre-wrap;background-color:rgb(255,255,255)"><font face="arial, helvetica, sans-serif">Paul </font></span></p><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><b style="font-weight:normal">
</b></p><div dir="ltr" style="margin-left:0pt"><table style="border:none;border-collapse:collapse"><colgroup><col width="467"></colgroup><tbody><tr style="height:0px"><td style="border-width:1px;border-style:solid;border-color:rgb(204,204,204);vertical-align:top;padding:7px">
<div dir="ltr" style="margin-left:0pt"><table style="border:none;border-collapse:collapse"><colgroup><col width="145"><col width="439"></colgroup><tbody><tr style="height:0px"><td style="border-width:0px;border-style:solid;border-color:rgb(0,0,0);vertical-align:top;padding:7px"><p dir="ltr" style="line-height:1.2;margin-top:0pt;margin-bottom:0pt"><span style="font-size:12.6667px;font-family:arial;color:rgb(34,34,34);background-color:rgb(255,255,255);font-weight:400;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;text-decoration:none;vertical-align:baseline;white-space:pre-wrap"><a moz-do-not-send="true" href="http://www.topx-geo-ict.nl" target="_blank"><img moz-do-not-send="true" src="https://docs.google.com/uc?id=0B8UeIGgvyx2VUEhWSzJiZWZUXzhrbWlPVkNiY0tVNTZRMFpF&export=download" alt="TopX Geo-ICT" height="96" width="78"></a></span></p></td><td style="border-width:0px;border-style:solid;border-color:rgb(0,0,0);vertical-align:top;padding:7px"><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:6pt"><span style="font-size:13.3333px;font-family:verdana;color:rgb(0,0,0);background-color:rgb(248,248,248);font-weight:700;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;text-decoration:none;vertical-align:baseline;white-space:pre-wrap">Paul Meems</span><span style="font-size:13.3333px;font-family:verdana;color:rgb(0,0,0);background-color:rgb(248,248,248);font-weight:700;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;text-decoration:none;vertical-align:baseline;white-space:pre-wrap">
</span><span style="font-size:13.3333px;font-family:verdana;color:rgb(128,128,128);background-color:rgb(248,248,248);font-weight:400;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;text-decoration:none;vertical-align:baseline;white-space:pre-wrap">Senior GIS consultant</span><span style="font-size:13.3333px;font-family:verdana;color:rgb(128,128,128);background-color:rgb(248,248,248);font-weight:400;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;text-decoration:none;vertical-align:baseline;white-space:pre-wrap">
</span><span style="font-family:verdana;color:rgb(34,34,34);font-weight:400;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;text-decoration:none;vertical-align:baseline;white-space:pre-wrap;background-color:rgb(248,248,248)"><font size="1">06-53989481</font></span></p><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:6pt"><span style="color:rgb(0,0,0);font-family:verdana;font-size:13.3333px;white-space:pre-wrap;line-height:1.38;background-color:rgb(248,248,248)">Wij bieden ondersteuning
</span><span style="font-size:13.3333px;font-family:verdana;color:rgb(0,0,0);vertical-align:baseline;white-space:pre-wrap;background-color:rgb(248,248,248)">voor </span><a moz-do-not-send="true" href="http://www.mapwindow.org/" style="line-height:1.38;font-size:12.8px;text-decoration:none" target="_blank"><span style="font-size:13.3333px;font-family:verdana;color:rgb(17,85,204);text-decoration:underline;vertical-align:baseline;white-space:pre-wrap;background-color:rgb(248,248,248)">MapWindow GIS</span></a></p></td></tr></tbody></table></div><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:6pt"><a moz-do-not-send="true" href="https://github.com/MapWindow/MapWinGIS/releases" style="text-decoration:none" target="_blank"><span style="font-size:13.3333px;font-family:verdana;color:rgb(17,85,204);background-color:rgb(248,248,248);font-weight:400;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;text-decoration:underline;vertical-align:baseline;white-space:pre-wrap">Download the latest MapWinGIS mapping engine.</span></a></p><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:6pt"><a moz-do-not-send="true" href="https://github.com/MapWindow/MapWindow5/releases" style="text-decoration:none" target="_blank"><span style="font-size:13.3333px;font-family:verdana;color:rgb(17,85,204);background-color:rgb(248,248,248);font-weight:400;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;text-decoration:underline;vertical-align:baseline;white-space:pre-wrap">Download the latest MapWindow 5 open source desktop application.</span></a></p></td></tr></tbody></table></div><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:0pt">
</p></span>
</div></div></div></div></div></div></div></div></div></div></div>
</div>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre wrap="">_______________________________________________
Dutch mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Dutch@lists.osgeo.org">Dutch@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="http://lists.osgeo.org/mailman/listinfo/dutch">http://lists.osgeo.org/mailman/listinfo/dutch</a></pre>
</blockquote>
<div class="moz-signature">--
<img src="cid:part7.DE14C3F4.12E5D247@hydroconsult.nl" border="0"></div></body></html>