<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Hello,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Hello devlist,<o:p></o:p></p>
<p class="MsoNormal">I have small patch of alaska elevation data in utm4/navd88 and I would like to convert this into wgs84 ellipsoidal. I know this data uses navd88,but gdalinfo does not report this vertical datum in the WKT. To reproject this data into
wgs84 ellipsoidal elevation, I’m assuming it should be possible with the following gdalwarp command:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">gdalwarp -s_srs "+proj=utm +zone=4 +datum=WGS84 +units=m +no_defs +geoidgrids=g2009alaska.gtx" -t_srs "+proj=longlat +datum=WGS84 +units=m +no_def" alaska_sample_utm4_navd88.tif alaska_sample_wgs84.tif<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I’m seeing the following errors reported from this command:<o:p></o:p></p>
<p class="MsoNormal">ERROR 1: point not within available datum shift grids<o:p></o:p></p>
<p class="MsoNormal">ERROR 1: GDALWarperOperation::ComputeSourceWindow() failed because the pfnTransformer failed.<o:p></o:p></p>
<p class="MsoNormal">Warning 1: Unable to compute source region for output window 0,0,351,106, skipping.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I’m using the geoidgrids from <a href="http://trac.osgeo.org/proj/wiki/VerticalDatums">
http://trac.osgeo.org/proj/wiki/VerticalDatums</a>, moved into the gdal/projlib directory<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The source file is in utm4:<o:p></o:p></p>
<p class="MsoNormal">PROJ.4 string is:<o:p></o:p></p>
<p class="MsoNormal">'+proj=utm +zone=4 +datum=WGS84 +units=m +no_defs '<o:p></o:p></p>
<p class="MsoNormal">Origin = (412011.848695503660000,7085064.815277445100000)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">For reference, this is about 100msurrounding the following lat/long:<o:p></o:p></p>
<p class="MsoNormal">Center (-160.7909106, 63.8815138) (160d47'27.28"W, 63d52'53.45"N)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">A simple gdalinfo on g2009alaska.gtx reports the following:<o:p></o:p></p>
<p class="MsoNormal">PROJ.4 string is:<o:p></o:p></p>
<p class="MsoNormal">'+proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs '<o:p></o:p></p>
<p class="MsoNormal">Origin = (171.991666666666670,72.008333333333340)<o:p></o:p></p>
<p class="MsoNormal">Pixel Size = (0.016666666666667,-0.016666666666667)<o:p></o:p></p>
<p class="MsoNormal">Corner Coordinates:<o:p></o:p></p>
<p class="MsoNormal">Upper Left ( 171.9916667, 72.0083333) (171d59'30.00"E, 72d 0'30.00"N)<o:p></o:p></p>
<p class="MsoNormal">Lower Left ( 171.9916667, 48.9916667) (171d59'30.00"E, 48d59'30.00"N)<o:p></o:p></p>
<p class="MsoNormal">Upper Right ( 234.008, 72.008) (234d 0'30.00"E, 72d 0'30.00"N)<o:p></o:p></p>
<p class="MsoNormal">Lower Right ( 234.008, 48.992) (234d 0'30.00"E, 48d59'30.00"N)<o:p></o:p></p>
<p class="MsoNormal">Center ( 203.000, 60.500) (203d 0' 0.00"E, 60d30' 0.00"N)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Does anyone see anything wrong with my gdalwarp/proj4 usage? Are there better ways to handle these sorts of transformations?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks,<o:p></o:p></p>
<p class="MsoNormal">Alex<o:p></o:p></p>
</div>
</body>
</html>