<div dir="ltr"><div style="font-size:12.8px">Ik probeer om de 3-uur neerslagaccumulaties van het KNMI (<a 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"><br clear="all"><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"><br></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"><br></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"><br></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><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"></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></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">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"><br></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"><br></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"><br></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"><br><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 href="http://www.topx-geo-ict.nl" target="_blank"><img src="https://docs.google.com/uc?id=0B8UeIGgvyx2VUEhWSzJiZWZUXzhrbWlPVkNiY0tVNTZRMFpF&export=download" width="78" height="96" alt="TopX Geo-ICT"></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"><br></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"><br></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 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 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 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"><br></p></span>
</div></div></div></div></div></div></div></div></div></div></div>
</div>