<div dir="ltr">To answer my own question:<div><br></div><div>GdalWarp is smart enough to internally reproject if needed.</div><div>So I can use my tiff in UTM and my border in 28992 and I get a perfectly clipped output.</div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><br>
<br>
Paul<br>
<br>
<div style="font-size:13.3px;font-family:Verdana,Arial,Helvetica,sans-serif"><div style="border:1px solid #e5e5e5;background-color:#f8f8f8;padding:0.8em 0.8em 0.3em 0.8em"><div style="margin:0 0 8px 0"><p style="margin:0"><span><b><img style="float:left;margin-left:5px;margin-right:5px" src="http://www.bontepaarden.nl/bontepaarden/images/newButton.png" alt="" height="120" width="160">Paul Meems </b><br><span>Release manager, configuration manager</span><br><span>and forum moderator of MapWindow GIS.</span><br><a href="http://www.mapwindow.org/" target="_blank"><span>www.mapwindow.org</span></a><br><br><span>Owner of MapWindow.nl - Support for</span><br><span>Dutch speaking users.</span><br><a href="http://www.mapwindow.nl/" target="_blank"><span>www.mapwindow.nl</span></a><br></span></p>
<p style="margin:0px"><b><br></b></p><p style="margin:0px"><b>The MapWindow GIS project has moved to <a href="https://github.com/MapWindow" target="_blank">GitHub</a>!</b></p><p style="font-size:13.3px;margin:0px"><br><a href="https://github.com/MapWindow/MapWinGIS/releases" target="_blank">Download the latest MapWinGIS mapping engine.</a><br></p><p style="font-size:13.3px;margin:0px"><span style="font-size:13.3px"><a href="https://github.com/MapWindow/MapWindow5/releases" target="_blank">Download the latest MapWindow 5 open source desktop application.</a></span></p></div></div><img src="http://p1.wisestamp.com/pixel.png?p=mozilla&v=3.9.12&t=1322134056890&u=482737fd0c4ca423" height="1" width="1"></div>
</div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">2017-06-12 16:14 GMT+02:00 Paul Meems <span dir="ltr"><<a href="mailto:bontepaarden@gmail.com" target="_blank">bontepaarden@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hi,</div><div><br></div>I'm trying to clip a tiff-file with a shapefile.<div>I'm using GDAL v2.1.3 with the SWIG-bindings in my C# application.</div><div><br><div>The tiff-file has this info:</div><div><div><font face="monospace, monospace">Driver: GTiff/GeoTIFF</font></div><div><font face="monospace, monospace">Files: </font></div><div><font face="monospace, monospace">Size is 5140, 6065</font></div><div><font face="monospace, monospace">Coordinate System is:</font></div><div><font face="monospace, monospace">PROJCS["WGS 84 / UTM zone 32N",</font></div><div><font face="monospace, monospace"> GEOGCS["WGS 84",</font></div><div><font face="monospace, monospace"> DATUM["WGS_1984",</font></div><div><font face="monospace, monospace"> SPHEROID["WGS 84",6378137,298.257223563,</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","7030"]],</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","6326"]],</font></div><div><font face="monospace, monospace"> PRIMEM["Greenwich",0,</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","8901"]],</font></div><div><font face="monospace, monospace"> UNIT["degree",0.<wbr>0174532925199433,</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","9122"]],</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","4326"]],</font></div><div><font face="monospace, monospace"> PROJECTION["Transverse_<wbr>Mercator"],</font></div><div><font face="monospace, monospace"> PARAMETER["latitude_of_origin"<wbr>,0],</font></div><div><font face="monospace, monospace"> PARAMETER["central_meridian",<wbr>9],</font></div><div><font face="monospace, monospace"> PARAMETER["scale_factor",0.<wbr>9996],</font></div><div><font face="monospace, monospace"> PARAMETER["false_easting",<wbr>500000],</font></div><div><font face="monospace, monospace"> PARAMETER["false_northing",0],</font></div><div><font face="monospace, monospace"> UNIT["metre",1,</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","9001"]],</font></div><div><font face="monospace, monospace"> AXIS["Easting",EAST],</font></div><div><font face="monospace, monospace"> AXIS["Northing",NORTH],</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","32632"]]</font></div><div><font face="monospace, monospace">Origin = (360810.191610000040000,<wbr>5860662.620930000200000)</font></div><div><font face="monospace, monospace">Pixel Size = (0.125010000000000,-0.<wbr>125010000000000)</font></div><div><font face="monospace, monospace">Metadata:</font></div><div><font face="monospace, monospace"> AREA_OR_POINT=Area</font></div><div><font face="monospace, monospace"> TIFFTAG_SOFTWARE=pix4dmapper</font></div><div><font face="monospace, monospace">Image Structure Metadata:</font></div><div><font face="monospace, monospace"> COMPRESSION=LZW</font></div><div><font face="monospace, monospace"> INTERLEAVE=BAND</font></div><div><font face="monospace, monospace">Corner Coordinates:</font></div><div><font face="monospace, monospace">Upper Left ( 360810.192, 5860662.621) ( 6d55'54.04"E, 52d52'39.50"N)</font></div><div><font face="monospace, monospace">Lower Left ( 360810.192, 5859904.435) ( 6d55'55.21"E, 52d52'14.98"N)</font></div><div><font face="monospace, monospace">Upper Right ( 361452.743, 5860662.621) ( 6d56'28.39"E, 52d52'40.10"N)</font></div><div><font face="monospace, monospace">Lower Right ( 361452.743, 5859904.435) ( 6d56'29.55"E, 52d52'15.58"N)</font></div><div><font face="monospace, monospace">Center ( 361131.467, 5860283.528) ( 6d56'11.80"E, 52d52'27.54"N)</font></div><div><font face="monospace, monospace">Band 1 Block=256x256 Type=Float32, ColorInterp=Gray</font></div><div><font face="monospace, monospace"> Min=-0.177 Max=33.644 </font></div><div><font face="monospace, monospace"> Minimum=-0.177, Maximum=33.644, Mean=18.605, StdDev=7.894</font></div><div><font face="monospace, monospace"> NoData Value=-10000</font></div><div><font face="monospace, monospace"> Metadata:</font></div><div><font face="monospace, monospace"> STATISTICS_MAXIMUM=33.<wbr>644180297852</font></div><div><font face="monospace, monospace"> STATISTICS_MEAN=18.<wbr>604999564386</font></div><div><font face="monospace, monospace"> STATISTICS_MINIMUM=-0.<wbr>17713606357574</font></div><div><font face="monospace, monospace"> STATISTICS_STDDEV=7.<wbr>8939827089249</font></div><div><br></div><div>The border can be in a different projection. Because I create this file earlier I know its EPSG-code.</div><div>Before clipping I want to check if both files are in the same projection.</div><div>How can I do this? I know I can use <font face="monospace, monospace">SpatialReference</font>. But how to convert the projection string in the metadata of my tiff to EPSG code?</div><div>If they have different projections I need to reproject my border shapefile first using VectorTranslate.</div><div><br></div><div>I was think of using these options:</div><div><div><font face="monospace, monospace">var options = new[]</font></div><div><font face="monospace, monospace">{</font></div><div><font face="monospace, monospace"> "-s_srs", "EPSG:" + sridInput,</font></div><div><font face="monospace, monospace"> "-t_srs", "EPSG:" + sridOutput,</font></div><div><font face="monospace, monospace"> "-f", "ESRI Shapefile"</font></div><div><font face="monospace, monospace">};</font></div><div><font face="monospace, monospace">VectorTranslate(input, output, options);</font></div></div><div><br></div><div>But then I need to get the EPSG code of the tiff file.</div><div><br></div><div>Please advice.</div><div><br></div><div>Thanks,</div><div><div class="m_1744896948388436183gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">
<br>
Paul<br>
</div></div></div></div></div></div>
</div></div></div>
</blockquote></div><br></div>