<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
FYI,</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
If you find errors in the 3DEP data, send an email to tnm_help@usgs.gov . I have forwarded your email with the list of files to that address.</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Doug</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> gdal-dev <gdal-dev-bounces@lists.osgeo.org> on behalf of Conrad Bielski via gdal-dev <gdal-dev@lists.osgeo.org><br>
<b>Sent:</b> Thursday, March 21, 2024 8:15 PM<br>
<b>To:</b> gdal-dev@lists.osgeo.org <gdal-dev@lists.osgeo.org><br>
<b>Subject:</b> [EXTERNAL] Re: [gdal-dev] gdal-dev Digest, Vol 238, Issue 82</font>
<div> </div>
</div>
<div>
<p> </p>
<div>
<table border="2" cellspacing="4" cellpadding="0" style="border-color:black; background-color:yellow; margin-left:auto; margin-right:auto">
<tbody>
<tr>
<td>
<p style="text-align:center"><strong> This email has been received from outside of DOI - Use caution before clicking on links, opening attachments, or responding. </strong></p>
</td>
</tr>
</tbody>
</table>
</div>
<p><br>
<br>
</p>
<div>
<div class="x_ydpe723df78yahoo-style-wrap" style="font-family:Helvetica Neue,Helvetica,Arial,sans-serif; font-size:16px">
<div></div>
<div dir="ltr" data-setdir="false">Thanks <span>Even,</span></div>
<div dir="ltr" data-setdir="false"><span><br>
</span></div>
<div dir="ltr" data-setdir="false"><span>you are correct about the tag number (I sent the incorrrect one).</span></div>
<div dir="ltr" data-setdir="false"><span><br>
</span></div>
<div dir="ltr" data-setdir="false"><span>I have provided a sample tif and some other descriptions of the issue here:</span></div>
<div dir="ltr" data-setdir="false"><a href="https://drive.google.com/drive/folders/1LVJqoKrnLRf6_5FpE7a5SrUGIWBPQmqa?usp=sharing" originalsrc="https://drive.google.com/drive/folders/1LVJqoKrnLRf6_5FpE7a5SrUGIWBPQmqa?usp=sharing" shash="PQIL9ctlZb2C9NJSHXRpnsJzJP9mMqtVF9wUHBfKOyr4jmoEIxSnMWrLCY0pla5YOuSIyEvupNAYSyhgedVPgGxnUNFD6k4r2lIr2kHjKdoDmFBqjcanI+Pt8c0WiIiV3wjCf4gajP5w7bL2M7V+HsM+4tK4Hve65rOOwlEJU/c=" rel="nofollow" target="_blank" class="x_enhancr_card_1237347669">bad_geotiff
- Google Drive</a></div>
<div><br>
</div>
<div id="x_ydpf608df10enhancr_card_1237347669" class="x_ydpf608df10yahoo-link-enhancr-card x_ydpf608df10yahoo-link-enhancr-not-allow-cover x_ydpf608df10ymail-preserve-class x_ydpf608df10ymail-preserve-style" data-url="https://drive.google.com/drive/folders/1LVJqoKrnLRf6_5FpE7a5SrUGIWBPQmqa?usp=sharing" data-type="YENHANCER" data-size="MEDIUM" style="max-width:400px; font-family:YahooSans VF,YahooSans,OpenSans VF,Helvetica Neue,Segoe UI,Helvetica,Arial,sans-serif">
<a href="https://drive.google.com/drive/folders/1LVJqoKrnLRf6_5FpE7a5SrUGIWBPQmqa?usp=sharing" originalsrc="https://drive.google.com/drive/folders/1LVJqoKrnLRf6_5FpE7a5SrUGIWBPQmqa?usp=sharing" shash="Gnbeg1Tkv7CxwU6PS5L+cnPl9ObDggai7mfEy5keRA3LFU17peD/1c2LkKLDAErXzRfc5v5+nEtNM1IP7aUO5nSvu6/F2NHAy4Zo4307mnupd5NlfB4ay9cPsFjNh+ay1Bayb39F2qJNZlDi027wXpwDQoNPrWhYY4QEs11qUyQ=" class="x_ydpf608df10yahoo-enhancr-cardlink" rel="nofollow" target="_blank" style="text-decoration-line:none!important; text-decoration-style:solid!important; text-decoration-color:currentcolor!important; color:rgb(0,0,0)!important">
<table border="0" class="x_ydpf608df10card-wrapper x_ydpf608df10yahoo-ignore-table" cellpadding="0" cellspacing="0" style="max-width:400px">
<tbody>
<tr>
<td width="400">
<table border="0" class="x_ydpf608df10card x_ydpf608df10yahoo-ignore-table" cellpadding="0" cellspacing="0" width="100%" style="max-width:400px; border-width:1px; border-style:solid; border-color:rgb(224,228,233); border-radius:2px">
<tbody>
<tr>
<td>
<table border="0" class="x_ydpf608df10card-info x_ydpf608df10yahoo-ignore-table" cellpadding="0" cellspacing="0" style="background-color:rgb(255,255,255); background-repeat:repeat; background-attachment:scroll; background-image:none; background-size:auto; z-index:2; width:100%; max-width:400px; border-radius:0px 0px 2px 2px; border-top:1px solid rgb(224,228,233)">
<tbody>
<tr>
<td style="background-color:#ffffff; padding:16px 0 16px 12px; vertical-align:top; border-radius:0 0 0 2px">
<img class="x_ydpf608df10card-object-1 x_ydpf608df10yahoo-ignore-inline-image x_ydpf608df10ymail-preserve-class" height="36" style="min-width:36px; margin-top:3px" src="https://s.yimg.com/lb/brands/80x80_google.png"></td>
<td style="vertical-align:middle; padding:12px 24px 16px 12px; width:99%; font-family:YahooSans VF,YahooSans,OpenSans VF,Helvetica Neue,Segoe UI,Helvetica,Arial,sans-serif; border-radius:0 0 2px 0">
<h2 class="x_ydpf608df10card-title" style="font-size:14px; line-height:19px; margin:0px 0px 6px; font-family:YahooSans VF,YahooSans,OpenSans VF,Helvetica Neue,Segoe UI,Helvetica,Arial,sans-serif; color:rgb(29,34,40); max-width:314px">
bad_geotiff - Google Drive</h2>
<p class="x_ydpf608df10card-description" style="font-size:12px; line-height:16px; margin:0px; color:rgb(151,158,168)">
</p>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
</a></div>
</div>
<div class="x_ydpe723df78yahoo-style-wrap" style="font-family:Helvetica Neue,Helvetica,Arial,sans-serif; font-size:16px">
<div><br>
</div>
<div dir="ltr" data-setdir="false">The three files there are the following:</div>
<div dir="ltr" data-setdir="false">
<div>
<div>1. A bad tiff (this is slightly smaller; they are all pretty big).</div>
<div>2. A GDAL info dump (survey feet is almost at the bottom)</div>
<div>3. A PDF with an explanation of the issue</div>
</div>
</div>
<div dir="ltr" data-setdir="false"><br>
</div>
<div dir="ltr" data-setdir="false">thanks for your help.</div>
<div dir="ltr" data-setdir="false">Conrad<br>
</div>
<div><br>
</div>
<div><br>
</div>
</div>
<div id="x_ydp394a6135yahoo_quoted_1765526266" class="x_ydp394a6135yahoo_quoted">
<div style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif; font-size:13px; color:#26282a">
<div>On Thursday, March 21, 2024 at 04:50:19 PM EDT, gdal-dev-request@lists.osgeo.org <gdal-dev-request@lists.osgeo.org> wrote:
</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div dir="ltr">Send gdal-dev mailing list submissions to<br>
</div>
<div dir="ltr"> <a href="mailto:gdal-dev@lists.osgeo.org" rel="nofollow" target="_blank">
gdal-dev@lists.osgeo.org</a><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">To subscribe or unsubscribe via the World Wide Web, visit<br>
</div>
<div dir="ltr"> <a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" originalsrc="https://lists.osgeo.org/mailman/listinfo/gdal-dev" shash="baefZ/6p7KeysNbMfDcyU3gU4PJ0B+cZx9kdLC6zBjK6+EAzlK5IvXU55HbHrgHgsB5+t1Q2DO1Ad6/xYRnDWbPHZBvowleKe6rs9tF9Vgo65Yk8eVF5ssRxCkKPthh1URY0B3JkTXSZjjWDdNZgC6sgw1cZCwfGGRFqtjk/yFk=" rel="nofollow" target="_blank">
https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</div>
<div dir="ltr">or, via email, send a message with subject or body 'help' to<br>
</div>
<div dir="ltr"> <a href="mailto:gdal-dev-request@lists.osgeo.org" rel="nofollow" target="_blank">
gdal-dev-request@lists.osgeo.org</a><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">You can reach the person managing the list at<br>
</div>
<div dir="ltr"> <a href="mailto:gdal-dev-owner@lists.osgeo.org" rel="nofollow" target="_blank">
gdal-dev-owner@lists.osgeo.org</a><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">When replying, please edit your Subject line so it is more specific<br>
</div>
<div dir="ltr">than "Re: Contents of gdal-dev digest..."<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Today's Topics:<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr"> 1. USGS 3DEP (3D Elevation Program) - feet should be metres, how<br>
</div>
<div dir="ltr"> can I fix this? (Conrad Bielski)<br>
</div>
<div dir="ltr"> 2. Re: USGS 3DEP (3D Elevation Program) - feet should be metres,<br>
</div>
<div dir="ltr"> how can I fix this? (Even Rouault)<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">----------------------------------------------------------------------<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Message: 1<br>
</div>
<div dir="ltr">Date: Thu, 21 Mar 2024 20:45:27 +0000 (UTC)<br>
</div>
<div dir="ltr">From: Conrad Bielski <<a href="mailto:conrad_bielski@yahoo.com" rel="nofollow" target="_blank">conrad_bielski@yahoo.com</a>><br>
</div>
<div dir="ltr">To: "gdal-dev@lists.osgeo.org" <<a href="mailto:gdal-dev@lists.osgeo.org" rel="nofollow" target="_blank">gdal-dev@lists.osgeo.org</a>><br>
</div>
<div dir="ltr">Subject: [gdal-dev] USGS 3DEP (3D Elevation Program) - feet should be<br>
</div>
<div dir="ltr"> metres, how can I fix this?<br>
</div>
<div dir="ltr">Message-ID: <<a href="mailto:1045195440.383478.1711053927068@mail.yahoo.com" rel="nofollow" target="_blank">1045195440.383478.1711053927068@mail.yahoo.com</a>><br>
</div>
<div dir="ltr">Content-Type: text/plain; charset="utf-8"<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Hello GDALers,<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">I have a question about reading USGS 3DEP (3D Elevation Program) data. Inside of this data, a GEOTIFF tag 42114 is provided which is causing problems with datum shifts.<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">So when I use GDAL to compute the datum shifts, the tag is read and interprets that the DEM is showing elevation in 'feet' while the DEM is actually in metres. The DEMs are in fact in meter elevations and meter UTM horizontal coordinates. This
is obviously erroneously integrated into the tag.<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">So my question is whether anyone has come across this issue and found a solution? Is there a way to edit the offending tag so that it is correctly interpreted as metres instead of feet? Other potential solutions so that GDAL interprets the elevation
correctly?<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Thanks in advance for your help. I provide below the output of gdalinfo which I hope could help with the offending info in bold.<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Have a great day,<br>
</div>
<div dir="ltr">Conrad<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Driver: GTiff/GeoTIFF<br>
</div>
<div dir="ltr">Files: I:\v3_us_3dep\asheville\ashevilleUSGS_1M_17_x34y392_NC_Phase5_2018_A18.tif<br>
</div>
<div dir="ltr">? ? ? ?I:\v3_us_3dep\asheville\ashevilleUSGS_1M_17_x34y392_NC_Phase5_2018_A18.tif.aux.xml<br>
</div>
<div dir="ltr">Size is 10012, 10012<br>
</div>
<div dir="ltr">Coordinate System is:<br>
</div>
<div dir="ltr">PROJCRS["NAD83 / UTM zone 17N",<br>
</div>
<div dir="ltr">? ? BASEGEOGCRS["NAD83",<br>
</div>
<div dir="ltr">? ? ? ? DATUM["North American Datum 1983",<br>
</div>
<div dir="ltr">? ? ? ? ? ? ELLIPSOID["GRS 1980",6378137,298.257222101,<br>
</div>
<div dir="ltr">? ? ? ? ? ? ? ? LENGTHUNIT["metre",1]]],<br>
</div>
<div dir="ltr">? ? ? ? PRIMEM["Greenwich",0,<br>
</div>
<div dir="ltr">? ? ? ? ? ? ANGLEUNIT["degree",0.0174532925199433]],<br>
</div>
<div dir="ltr">? ? ? ? ID["EPSG",4269]],<br>
</div>
<div dir="ltr">? ? CONVERSION["UTM zone 17N",<br>
</div>
<div dir="ltr">? ? ? ? METHOD["Transverse Mercator",<br>
</div>
<div dir="ltr">? ? ? ? ? ? ID["EPSG",9807]],<br>
</div>
<div dir="ltr">? ? ? ? PARAMETER["Latitude of natural origin",0,<br>
</div>
<div dir="ltr">? ? ? ? ? ? ANGLEUNIT["degree",0.0174532925199433],<br>
</div>
<div dir="ltr">? ? ? ? ? ? ID["EPSG",8801]],<br>
</div>
<div dir="ltr">? ? ? ? PARAMETER["Longitude of natural origin",-81,<br>
</div>
<div dir="ltr">? ? ? ? ? ? ANGLEUNIT["degree",0.0174532925199433],<br>
</div>
<div dir="ltr">? ? ? ? ? ? ID["EPSG",8802]],<br>
</div>
<div dir="ltr">? ? ? ? PARAMETER["Scale factor at natural origin",0.9996,<br>
</div>
<div dir="ltr">? ? ? ? ? ? SCALEUNIT["unity",1],<br>
</div>
<div dir="ltr">? ? ? ? ? ? ID["EPSG",8805]],<br>
</div>
<div dir="ltr">? ? ? ? PARAMETER["False easting",500000,<br>
</div>
<div dir="ltr">? ? ? ? ? ? LENGTHUNIT["metre",1],<br>
</div>
<div dir="ltr">? ? ? ? ? ? ID["EPSG",8806]],<br>
</div>
<div dir="ltr">? ? ? ? PARAMETER["False northing",0,<br>
</div>
<div dir="ltr">? ? ? ? ? ? LENGTHUNIT["metre",1],<br>
</div>
<div dir="ltr">? ? ? ? ? ? ID["EPSG",8807]]],<br>
</div>
<div dir="ltr">? ? CS[Cartesian,2],<br>
</div>
<div dir="ltr">? ? ? ? AXIS["(E)",east,<br>
</div>
<div dir="ltr">? ? ? ? ? ? ORDER[1],<br>
</div>
<div dir="ltr">? ? ? ? ? ? LENGTHUNIT["metre",1]],<br>
</div>
<div dir="ltr">? ? ? ? AXIS["(N)",north,<br>
</div>
<div dir="ltr">? ? ? ? ? ? ORDER[2],<br>
</div>
<div dir="ltr">? ? ? ? ? ? LENGTHUNIT["metre",1]],<br>
</div>
<div dir="ltr">? ? USAGE[<br>
</div>
<div dir="ltr">? ? ? ? SCOPE["Engineering survey, topographic mapping."],<br>
</div>
<div dir="ltr">? ? ? ? AREA["North America - between 84?W and 78?W - onshore and offshore. Canada - Nunavut; Ontario; Quebec. United States (USA) - Florida; Georgia; Kentucky; Maryland; Michigan; New York; North Carolina; Ohio; Pennsylvania; South Carolina;
Tennessee; Virginia; West Virginia."],<br>
</div>
<div dir="ltr">? ? ? ? BBOX[23.81,-84,84,-78]],<br>
</div>
<div dir="ltr">? ? ID["EPSG",26917]]<br>
</div>
<div dir="ltr">Data axis to CRS axis mapping: 1,2<br>
</div>
<div dir="ltr">Origin = (339993.999958705157042,3920005.999981591943651)<br>
</div>
<div dir="ltr">Pixel Size = (1.000000000000000,-1.000000000000000)<br>
</div>
<div dir="ltr">Metadata:<br>
</div>
<div dir="ltr">? AREA_OR_POINT=Area<br>
</div>
<div dir="ltr">Image Structure Metadata:<br>
</div>
<div dir="ltr">? COMPRESSION=LZW<br>
</div>
<div dir="ltr">? INTERLEAVE=BAND<br>
</div>
<div dir="ltr">? LAYOUT=COG<br>
</div>
<div dir="ltr">? PREDICTOR=3<br>
</div>
<div dir="ltr">Corner Coordinates:<br>
</div>
<div dir="ltr">Upper Left ?( ?339994.000, 3920006.000) ( 82d45'43.98"W, 35d24'38.21"N)<br>
</div>
<div dir="ltr">Lower Left ?( ?339994.000, 3909994.000) ( 82d45'36.92"W, 35d19'13.36"N)<br>
</div>
<div dir="ltr">Upper Right ( ?350006.000, 3920006.000) ( 82d39' 7.17"W, 35d24'43.82"N)<br>
</div>
<div dir="ltr">Lower Right ( ?350006.000, 3909994.000) ( 82d39' 0.55"W, 35d19'18.95"N)<br>
</div>
<div dir="ltr">Center ? ? ?( ?345000.000, 3915000.000) ( 82d42'22.15"W, 35d21'58.63"N)<br>
</div>
<div dir="ltr">Band 1 Block=512x512 Type=Float32, ColorInterp=Gray<br>
</div>
<div dir="ltr">? Min=669.509 Max=1604.949 <br>
</div>
<div dir="ltr">? Minimum=669.509, Maximum=1604.949, Mean=988.295, StdDev=159.804<br>
</div>
<div dir="ltr">? NoData Value=-999999<br>
</div>
<div dir="ltr">? Overviews: 5006x5006, 2503x2503, 1251x1251, 625x625, 312x312<br>
</div>
<div dir="ltr">? Unit Type: US survey foot<br>
</div>
<div dir="ltr">? Metadata:<br>
</div>
<div dir="ltr">? ? STATISTICS_APPROXIMATE=YES<br>
</div>
<div dir="ltr">? ? STATISTICS_MAXIMUM=1604.9490966797<br>
</div>
<div dir="ltr">? ? STATISTICS_MEAN=988.29463693706<br>
</div>
<div dir="ltr">? ? STATISTICS_MINIMUM=669.50915527344<br>
</div>
<div dir="ltr">? ? STATISTICS_STDDEV=159.80361228109<br>
</div>
<div dir="ltr">? ? STATISTICS_VALID_PERCENT=100<br>
</div>
<div dir="ltr">-------------- next part --------------<br>
</div>
<div dir="ltr">An HTML attachment was scrubbed...<br>
</div>
<div dir="ltr">URL: <<a href="http://lists.osgeo.org/pipermail/gdal-dev/attachments/20240321/68938236/attachment-0001.htm" originalsrc="http://lists.osgeo.org/pipermail/gdal-dev/attachments/20240321/68938236/attachment-0001.htm" shash="U6fey0TfC72MwYDFpjrJZym/P1c4sRXCfUjInXEUi7qLMFkf5D9Dy6Lb3hpaYvBUSjqlNsyaRuJC/oXrQhGdxB0ZtqUqNywu1480ufOwPjXYKUW6dcaf+K7k0K/PMnJ7n20q5yDUDQC2UmMWzR+ljGV5hnwQzCOEGa0MeHdpIZg=" rel="nofollow" target="_blank">http://lists.osgeo.org/pipermail/gdal-dev/attachments/20240321/68938236/attachment-0001.htm</a>><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">------------------------------<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Message: 2<br>
</div>
<div dir="ltr">Date: Thu, 21 Mar 2024 21:49:55 +0100<br>
</div>
<div dir="ltr">From: Even Rouault <<a href="mailto:even.rouault@spatialys.com" rel="nofollow" target="_blank">even.rouault@spatialys.com</a>><br>
</div>
<div dir="ltr">To: Conrad Bielski <<a href="mailto:conrad_bielski@yahoo.com" rel="nofollow" target="_blank">conrad_bielski@yahoo.com</a>>,<br>
</div>
<div dir="ltr"> "gdal-dev@lists.osgeo.org" <<a href="mailto:gdal-dev@lists.osgeo.org" rel="nofollow" target="_blank">gdal-dev@lists.osgeo.org</a>><br>
</div>
<div dir="ltr">Subject: Re: [gdal-dev] USGS 3DEP (3D Elevation Program) - feet should<br>
</div>
<div dir="ltr"> be metres, how can I fix this?<br>
</div>
<div dir="ltr">Message-ID: <<a href="mailto:7832c1d9-61ce-41b0-b524-2f00859c817a@spatialys.com" rel="nofollow" target="_blank">7832c1d9-61ce-41b0-b524-2f00859c817a@spatialys.com</a>><br>
</div>
<div dir="ltr">Content-Type: text/plain; charset="utf-8"; Format="flowed"<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Le 21/03/2024 ? 21:45, Conrad Bielski via gdal-dev a ?crit?:<br>
</div>
<div dir="ltr">> Hello GDALers,<br>
</div>
<div dir="ltr">><br>
</div>
<div dir="ltr">> I have a question about reading USGS 3DEP (3D Elevation Program) data.
<br>
</div>
<div dir="ltr">> Inside of this data, a GEOTIFF tag 42114 is provided which is causing
<br>
</div>
<div dir="ltr">> problems with datum shifts.<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">There's no such thing as a GEOTIFF tag 42114.<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">The closest tags are:<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">#define TIFFTAG_GDAL_METADATA 42112<br>
</div>
<div dir="ltr">#define TIFFTAG_GDAL_NODATA 42113<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">It would be helpuful if you could provide the file, or the result of a
<br>
</div>
<div dir="ltr">"tiffdump -m 1000 the.tif" and "listgeo the.tif" on it.<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">And was the .tif.aux.xml provided with the TIFF, or computed by GDAL.
<br>
</div>
<div dir="ltr">Because the information might come from it too.<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">><br>
</div>
<div dir="ltr">> So when I use GDAL to compute the datum shifts, the tag is read and
<br>
</div>
<div dir="ltr">> interprets that the DEM is showing elevation in 'feet' while the DEM
<br>
</div>
<div dir="ltr">> is actually in metres. The DEMs are in fact in meter elevations and
<br>
</div>
<div dir="ltr">> meter UTM horizontal coordinates. This is obviously erroneously <br>
</div>
<div dir="ltr">> integrated into the tag.<br>
</div>
<div dir="ltr">><br>
</div>
<div dir="ltr">> So my question is whether anyone has come across this issue and found
<br>
</div>
<div dir="ltr">> a solution? Is there a way to edit the offending tag so that it is
<br>
</div>
<div dir="ltr">> correctly interpreted as metres instead of feet? Other potential
<br>
</div>
<div dir="ltr">> solutions so that GDAL interprets the elevation correctly?<br>
</div>
<div dir="ltr">><br>
</div>
<div dir="ltr">> Thanks in advance for your help. I provide below the output of <br>
</div>
<div dir="ltr">> gdalinfo which I hope could help with the offending info in bold.<br>
</div>
<div dir="ltr">><br>
</div>
<div dir="ltr">> Have a great day,<br>
</div>
<div dir="ltr">> Conrad<br>
</div>
<div dir="ltr">><br>
</div>
<div dir="ltr">><br>
</div>
<div dir="ltr">><br>
</div>
<div dir="ltr">><br>
</div>
<div dir="ltr">> Driver: GTiff/GeoTIFF<br>
</div>
<div dir="ltr">> Files: <br>
</div>
<div dir="ltr">> I:\v3_us_3dep\asheville\ashevilleUSGS_1M_17_x34y392_NC_Phase5_2018_A18.tif<br>
</div>
<div dir="ltr">> ?I:\v3_us_3dep\asheville\ashevilleUSGS_1M_17_x34y392_NC_Phase5_2018_A18.tif.aux.xml<br>
</div>
<div dir="ltr">> Size is 10012, 10012<br>
</div>
<div dir="ltr">> Coordinate System is:<br>
</div>
<div dir="ltr">> PROJCRS["NAD83 / UTM zone 17N",<br>
</div>
<div dir="ltr">> ? ? BASEGEOGCRS["NAD83",<br>
</div>
<div dir="ltr">> ? ? ? ? DATUM["North American Datum 1983",<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ELLIPSOID["GRS 1980",6378137,298.257222101,<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ? ? LENGTHUNIT["metre",1]]],<br>
</div>
<div dir="ltr">> ? ? ? ? PRIMEM["Greenwich",0,<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ANGLEUNIT["degree",0.0174532925199433]],<br>
</div>
<div dir="ltr">> ? ? ? ? ID["EPSG",4269]],<br>
</div>
<div dir="ltr">> ? ? CONVERSION["UTM zone 17N",<br>
</div>
<div dir="ltr">> ? ? ? ? METHOD["Transverse Mercator",<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ID["EPSG",9807]],<br>
</div>
<div dir="ltr">> ? ? ? ? PARAMETER["Latitude of natural origin",0,<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ANGLEUNIT["degree",0.0174532925199433],<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ID["EPSG",8801]],<br>
</div>
<div dir="ltr">> ? ? ? ? PARAMETER["Longitude of natural origin",-81,<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ANGLEUNIT["degree",0.0174532925199433],<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ID["EPSG",8802]],<br>
</div>
<div dir="ltr">> ? ? ? ? PARAMETER["Scale factor at natural origin",0.9996,<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? SCALEUNIT["unity",1],<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ID["EPSG",8805]],<br>
</div>
<div dir="ltr">> ? ? ? ? PARAMETER["False easting",500000,<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? LENGTHUNIT["metre",1],<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ID["EPSG",8806]],<br>
</div>
<div dir="ltr">> ? ? ? ? PARAMETER["False northing",0,<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? LENGTHUNIT["metre",1],<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ID["EPSG",8807]]],<br>
</div>
<div dir="ltr">> ? ? CS[Cartesian,2],<br>
</div>
<div dir="ltr">> ? ? ? ? AXIS["(E)",east,<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ORDER[1],<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? LENGTHUNIT["metre",1]],<br>
</div>
<div dir="ltr">> ? ? ? ? AXIS["(N)",north,<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? ORDER[2],<br>
</div>
<div dir="ltr">> ? ? ? ? ? ? LENGTHUNIT["metre",1]],<br>
</div>
<div dir="ltr">> ? ? USAGE[<br>
</div>
<div dir="ltr">> ? ? ? ? SCOPE["Engineering survey, topographic mapping."],<br>
</div>
<div dir="ltr">> ? ? ? ? AREA["North America - between 84?W and 78?W - onshore and
<br>
</div>
<div dir="ltr">> offshore. Canada - Nunavut; Ontario; Quebec. United States (USA) -
<br>
</div>
<div dir="ltr">> Florida; Georgia; Kentucky; Maryland; Michigan; New York; North <br>
</div>
<div dir="ltr">> Carolina; Ohio; Pennsylvania; South Carolina; Tennessee; Virginia;
<br>
</div>
<div dir="ltr">> West Virginia."],<br>
</div>
<div dir="ltr">> ? ? ? ? BBOX[23.81,-84,84,-78]],<br>
</div>
<div dir="ltr">> ? ? ID["EPSG",26917]]<br>
</div>
<div dir="ltr">> Data axis to CRS axis mapping: 1,2<br>
</div>
<div dir="ltr">> Origin = (339993.999958705157042,3920005.999981591943651)<br>
</div>
<div dir="ltr">> Pixel Size = (1.000000000000000,-1.000000000000000)<br>
</div>
<div dir="ltr">> Metadata:<br>
</div>
<div dir="ltr">> ? AREA_OR_POINT=Area<br>
</div>
<div dir="ltr">> Image Structure Metadata:<br>
</div>
<div dir="ltr">> ? COMPRESSION=LZW<br>
</div>
<div dir="ltr">> ? INTERLEAVE=BAND<br>
</div>
<div dir="ltr">> ? LAYOUT=COG<br>
</div>
<div dir="ltr">> ? PREDICTOR=3<br>
</div>
<div dir="ltr">> Corner Coordinates:<br>
</div>
<div dir="ltr">> Upper Left ?( ?339994.000, 3920006.000) ( 82d45'43.98"W, 35d24'38.21"N)<br>
</div>
<div dir="ltr">> Lower Left ?( ?339994.000, 3909994.000) ( 82d45'36.92"W, 35d19'13.36"N)<br>
</div>
<div dir="ltr">> Upper Right ( ?350006.000, 3920006.000) ( 82d39' 7.17"W, 35d24'43.82"N)<br>
</div>
<div dir="ltr">> Lower Right ( ?350006.000, 3909994.000) ( 82d39' 0.55"W, 35d19'18.95"N)<br>
</div>
<div dir="ltr">> Center ? ? ?( ?345000.000, 3915000.000) ( 82d42'22.15"W, 35d21'58.63"N)<br>
</div>
<div dir="ltr">> Band 1 Block=512x512 Type=Float32, ColorInterp=Gray<br>
</div>
<div dir="ltr">> ? Min=669.509 Max=1604.949<br>
</div>
<div dir="ltr">> ? Minimum=669.509, Maximum=1604.949, Mean=988.295, StdDev=159.804<br>
</div>
<div dir="ltr">> ? NoData Value=-999999<br>
</div>
<div dir="ltr">> ? Overviews: 5006x5006, 2503x2503, 1251x1251, 625x625, 312x312<br>
</div>
<div dir="ltr">> *Unit Type: US survey foot*<br>
</div>
<div dir="ltr">> ? Metadata:<br>
</div>
<div dir="ltr">> ? ? STATISTICS_APPROXIMATE=YES<br>
</div>
<div dir="ltr">> ? ? STATISTICS_MAXIMUM=1604.9490966797<br>
</div>
<div dir="ltr">> ? ? STATISTICS_MEAN=988.29463693706<br>
</div>
<div dir="ltr">> ? ? STATISTICS_MINIMUM=669.50915527344<br>
</div>
<div dir="ltr">> ? ? STATISTICS_STDDEV=159.80361228109<br>
</div>
<div dir="ltr">> ? ? STATISTICS_VALID_PERCENT=100<br>
</div>
<div dir="ltr">><br>
</div>
<div dir="ltr">><br>
</div>
<div dir="ltr">> _______________________________________________<br>
</div>
<div dir="ltr">> gdal-dev mailing list<br>
</div>
<div dir="ltr">> <a href="mailto:gdal-dev@lists.osgeo.org" rel="nofollow" target="_blank">
gdal-dev@lists.osgeo.org</a><br>
</div>
<div dir="ltr">> <a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" originalsrc="https://lists.osgeo.org/mailman/listinfo/gdal-dev" shash="hOKjXY959JjRDnyrcylAV/V14X3BJSiKog9UyM4/o+G8lMh/dS3mUqjiUZBedustpwUKuwhipMY8aFgBQLWHFgSKKBMumJUJQ0wLD8Z6t71X3kS6FRkb7ojiBJgY3Uhici6IjFqjz0Fm/S4UgyO771oRhHeBcX2LJ2VN+CGh2So=" rel="nofollow" target="_blank">
https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">-- <br>
</div>
<div dir="ltr"><a href="http://www.spatialys.com/" originalsrc="http://www.spatialys.com/" shash="KpcnsYs6BT5Cn9z2yww+LRc1KPAjNTC7V7OVPcvPGdL2PI5ZBN9TrFhuKQxECo/bRot4P9ngfnkyh3tDIhZdbB6ESFwhhIJTwvPKRSdK+ifoD/cGKlAam6hIhbaT515mr/tMsohkf+d8bdYSpheqJZry+zEAkY42DyYHpqjzeh8=" rel="nofollow" target="_blank">http://www.spatialys.com</a><br>
</div>
<div dir="ltr">My software is free, but my time generally not.<br>
</div>
<div dir="ltr">-------------- next part --------------<br>
</div>
<div dir="ltr">An HTML attachment was scrubbed...<br>
</div>
<div dir="ltr">URL: <<a href="http://lists.osgeo.org/pipermail/gdal-dev/attachments/20240321/1e68a1d4/attachment.htm" originalsrc="http://lists.osgeo.org/pipermail/gdal-dev/attachments/20240321/1e68a1d4/attachment.htm" shash="ABtLXLRHguvia4tooavuz+7eusAIojx8eDQjruOA+zgYsjsR6g4nN4EtcelGfICHxbi2q05Qw3ubgl9805R4fNiSRk1bFu+8wZ/8/jXaNdkh1S78XHwsv9OK5876ZH9vX1/mEG/tGnbrd0ljSevE1Xg80EEvJVYdScOYrbHXooI=" rel="nofollow" target="_blank">http://lists.osgeo.org/pipermail/gdal-dev/attachments/20240321/1e68a1d4/attachment.htm</a>><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">------------------------------<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Subject: Digest Footer<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">_______________________________________________<br>
</div>
<div dir="ltr">gdal-dev mailing list<br>
</div>
<div dir="ltr"><a href="mailto:gdal-dev@lists.osgeo.org" rel="nofollow" target="_blank">gdal-dev@lists.osgeo.org</a><br>
</div>
<div dir="ltr"><a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" originalsrc="https://lists.osgeo.org/mailman/listinfo/gdal-dev" shash="o10vhFrM5Jw28ExuCr+QGQ9VP2oATQ/LzjUpZfU4ag43M8kF9qtJSLF1qwhJUhptl0QMju+duSk7IrOR0c+bg0s3OruQalNdO++KHLFZyzfGZXzIZFkmkncj0phXwCq57bIzSzyRo45fXuGGY49yatTbobnH9F4/BS2L1chD0xM=" rel="nofollow" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">------------------------------<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">End of gdal-dev Digest, Vol 238, Issue 82<br>
</div>
<div dir="ltr">*****************************************<br>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>