<div dir="ltr">I'm trying to use the precipitation accumulations from our national weather agency: <a href="https://data.knmi.nl/datasets/radar_tar_corr_accum_03h/2.0" target="_blank" style="font-size:12.8px">https://data.knmi.nl/datasets<wbr>/radar_tar_corr_accum_03h/2.0</a><div>The data is in HDF5 with two subsets.</div><div>I only need one subset and I need the data in WGS84.<br></div><div><br></div><div>I'm using GDAL v2.2 from GisInternals.com</div><div>I tried using GDALWarp but then I get an error about missing the transformation matrix.</div><div>Using Google I figured out I need to do GDALTranslate first to set the georeferenced bounds.</div><div><br></div><div>This is the result of GDALInfo:</div><div><div><font face="monospace, monospace">Driver: HDF5/Hierarchical Data Format Release 5</font></div><div><font face="monospace, monospace">Files: RAD_NL25_RAC_03H_201612041200.h5</font></div><div><font face="monospace, monospace">Size is 512, 512</font></div><div><font face="monospace, monospace">Coordinate System is `'</font></div><div><font face="monospace, monospace">Metadata:</font></div><div><font face="monospace, monospace"> geographic_geo_column_offset=0 </font></div><div><font face="monospace, monospace"> geographic_geo_dim_pixel=KM,KM</font></div><div><font face="monospace, monospace"> geographic_geo_number_columns=700 </font></div><div><font face="monospace, monospace"> geographic_geo_number_rows=765 </font></div><div><font face="monospace, monospace"> geographic_geo_par_pixel=X,Y</font></div><div><font face="monospace, monospace"> geographic_geo_pixel_def=LU</font></div><div><font face="monospace, monospace"> geographic_geo_pixel_size_x=1.0000013 </font></div><div><font face="monospace, monospace"> geographic_geo_pixel_size_y=-1.0000055 </font></div><div><font face="monospace, monospace"> geographic_geo_product_corners=0 49.362053 0 55.973602 10.856429 55.388977 9.0092793 48.895298 </font></div><div><font face="monospace, monospace"> geographic_geo_row_offset=3649.9792 </font></div><div><font face="monospace, monospace"> geographic_map_projection_projection_indication=Y</font></div><div><font face="monospace, monospace"> geographic_map_projection_projection_name=STEREOGRAPHIC</font></div><div><font face="monospace, monospace"> geographic_map_projection_projection_proj4_params=+proj=stere +lat_0=90 +lon_0=0 +lat_ts=60 +a=6378.14 +b=6356.75 +x_0=0 y_0=0</font></div><div><font face="monospace, monospace"> image1_calibration_calibration_flag=Y</font></div><div><font face="monospace, monospace"> image1_calibration_calibration_formulas=GEO = 0.010000 * PV + 0.000000</font></div><div><font face="monospace, monospace"> image1_calibration_calibration_missing_data=0 </font></div><div><font face="monospace, monospace"> image1_calibration_calibration_out_of_image=65535 </font></div><div><font face="monospace, monospace"> image1_image_bytes_per_pixel=2 </font></div><div><font face="monospace, monospace"> image1_image_data_CLASS=IMAGE</font></div><div><font face="monospace, monospace"> image1_image_data_VERSION=1.2</font></div><div><font face="monospace, monospace"> image1_image_geo_parameter=PRECIP_[MM]</font></div><div><font face="monospace, monospace"> image1_image_number_accumulated_products=37 </font></div><div><font face="monospace, monospace"> image1_image_product_name=RAD_NL25_RAC_H1.5_03H</font></div><div><font face="monospace, monospace"> image1_image_size=535500 </font></div><div><font face="monospace, monospace"> image1_statistics_stat_max_value=9.96 </font></div><div><font face="monospace, monospace"> image1_statistics_stat_min_value=0 </font></div><div><font face="monospace, monospace"> image2_calibration_calibration_flag=Y</font></div><div><font face="monospace, monospace"> image2_calibration_calibration_formulas=GEO=0.001*PV+-32.768</font></div><div><font face="monospace, monospace"> image2_calibration_calibration_missing_data=0 </font></div><div><font face="monospace, monospace"> image2_calibration_calibration_out_of_image=65535 </font></div><div><font face="monospace, monospace"> image2_image_bytes_per_pixel=2 </font></div><div><font face="monospace, monospace"> image2_image_data_CLASS=IMAGE</font></div><div><font face="monospace, monospace"> image2_image_data_VERSION=1.2</font></div><div><font face="monospace, monospace"> image2_image_geo_parameter=ADJUSTMENT_FACTOR_[DB]</font></div><div><font face="monospace, monospace"> image2_image_number_accumulated_products=37 </font></div><div><font face="monospace, monospace"> image2_image_product_name=RAD_NL25_RAC_H1.5_03H</font></div><div><font face="monospace, monospace"> image2_image_size=535500 </font></div><div><font face="monospace, monospace"> image2_statistics_stat_max_value=0 </font></div><div><font face="monospace, monospace"> image2_statistics_stat_min_value=0 </font></div><div><font face="monospace, monospace"> overview_hdftag_version_number=3.6</font></div><div><font face="monospace, monospace"> overview_number_image_groups=2 </font></div><div><font face="monospace, monospace"> overview_number_radar_groups=1 </font></div><div><font face="monospace, monospace"> overview_number_satellite_groups=0 </font></div><div><font face="monospace, monospace"> overview_number_station_groups=0 </font></div><div><font face="monospace, monospace"> overview_product_datetime_end=04-DEC-2016;12:00:00.000</font></div><div><font face="monospace, monospace"> overview_product_datetime_start=04-DEC-2016;09:00:00.000</font></div><div><font face="monospace, monospace"> overview_product_group_name=RAD_NL25_RAC_03H</font></div><div><font face="monospace, monospace"> overview_products_missing=NA</font></div><div><font face="monospace, monospace"> radar1_radar_adjustment=F=-0.00dB</font></div><div><font face="monospace, monospace"> radar1_radar_location=5.17834 52.101681 </font></div><div><font face="monospace, monospace"> radar1_radar_name=DeBilt</font></div><div><font face="monospace, monospace">Subdatasets:</font></div><div><font face="monospace, monospace"> SUBDATASET_1_NAME=HDF5:"RAD_NL25_RAC_03H_201612041200.h5"://image1/image_data</font></div><div><font face="monospace, monospace"> SUBDATASET_1_DESC=[765x700] //image1/image_data (16-bit unsigned integer)</font></div><div><font face="monospace, monospace"> SUBDATASET_2_NAME=HDF5:"RAD_NL25_RAC_03H_201612041200.h5"://image2/image_data</font></div><div><font face="monospace, monospace"> SUBDATASET_2_DESC=[765x700] //image2/image_data (16-bit unsigned integer)</font></div><div><font face="monospace, monospace">Corner Coordinates:</font></div><div><font face="monospace, monospace">Upper Left ( 0.0, 0.0)</font></div><div><font face="monospace, monospace">Lower Left ( 0.0, 512.0)</font></div><div><font face="monospace, monospace">Upper Right ( 512.0, 0.0)</font></div><div><font face="monospace, monospace">Lower Right ( 512.0, 512.0)</font></div><div><font face="monospace, monospace">Center ( 256.0, 256.0)</font></div><div> </div></div><div>This is my GDALTranslate command, the projectionstring and a_ullr come from the meta data:</div><div><pre class="gmail-m_-2999790335521197642gmail-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_-2999790335521197642gmail-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_-2999790335521197642gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote" style="color:rgb(0,145,0)">"HDF5:\"</span>RAD_NL25_RAC_03H_<wbr>201612041200.h5\<span class="gmail-m_-2999790335521197642gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote" style="color:rgb(0,145,0)">":<span class="gmail-m_-2999790335521197642gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-comment" style="color:rgb(128,128,128)">//image1/<wbr>image_data"</span> <span class="gmail-m_-2999790335521197642gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote">"test-translate.tif"</span></span></pre></div><div><br></div><div>This is my GDALWarp command:</div><div><pre class="gmail-m_-2999790335521197642gmail-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_-2999790335521197642gmail-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_-2999790335521197642gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote" style="color:rgb(0,145,0)">"test-translate.tif"</span> <span class="gmail-m_-2999790335521197642gmail-m_6090844043269395631gmail-m_4197702424755525384gmail-code-quote" style="color:rgb(0,145,0)">"test-warp.tif"</span></pre></div><div><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>The warped tif is not where I expect it to be.</div><div>I expect a origin around 52 en 5 degrees (The Netherlands)</div><div>I'm getting: Origin = (169.559973821725920,89.530840368222059)</div><div><br></div><div>I've checked with the data provider and the projectionstring is correct.</div><div>What am I doing wrong?</div><div dir="ltr"><br>
Thanks,</div><div dir="ltr"><br>
Paul<br>
</div></div></div></div></div></div>
</div></div>