<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.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_Mercator"],</font></div><div><font face="monospace, monospace"> PARAMETER["latitude_of_origin",0],</font></div><div><font face="monospace, monospace"> PARAMETER["central_meridian",9],</font></div><div><font face="monospace, monospace"> PARAMETER["scale_factor",0.9996],</font></div><div><font face="monospace, monospace"> PARAMETER["false_easting",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,5860662.620930000200000)</font></div><div><font face="monospace, monospace">Pixel Size = (0.125010000000000,-0.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.644180297852</font></div><div><font face="monospace, monospace"> STATISTICS_MEAN=18.604999564386</font></div><div><font face="monospace, monospace"> STATISTICS_MINIMUM=-0.17713606357574</font></div><div><font face="monospace, monospace"> STATISTICS_STDDEV=7.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="gmail_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>