<div class="gmail_quote"><div><div class="h5"><div class="gmail_quote">On Tue, Jan 10, 2012 at 4:10 PM, Frank Warmerdam <span dir="ltr"><<a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>On Tue, Jan 10, 2012 at 9:01 AM, Howard Butler <<a href="mailto:hobu.inc@gmail.com" target="_blank">hobu.inc@gmail.com</a>> wrote:<br>
><br>
> On Jan 10, 2012, at 10:30 AM, Howard Butler wrote:<br>
>> I'll try to dig into this some more...<br>
><br>
> GDAL isn't turning this into a proper proj.4 string.<br>
><br>
>> GTModelTypeGeoKey (Short,1): ModelTypeProjected<br>
>> GTCitationGeoKey (Ascii,14): "UTM zone 15N"<br>
>> GeogGeodeticDatumGeoKey (Short,1): Datum_North_American_Datum_1983<br>
>> GeogEllipsoidGeoKey (Short,1): Ellipse_GRS_1980<br>
>> ProjectedCSTypeGeoKey (Short,1): PCS_NAD83_UTM_zone_15N<br>
>> ProjCoordTransGeoKey (Short,1): CT_TransverseMercator<br>
>> ProjLinearUnitsGeoKey (Short,1): Linear_Meter<br>
>> VerticalCSTypeGeoKey (Short,1): VertCS_North_American_Vertical_Datum_1988<br>
>> VerticalCitationGeoKey (Ascii,19): "GEOID 2003 Applied"<br>
>> VerticalUnitsGeoKey (Short,1): Linear_Meter<br>
><br>
> is coming out of GDAL as:<br>
><br>
>> +proj=tmerc +lat_0=0 +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=NAD83 +units=m +geoidgrids=g2003conus.gtx,g2003alaska.gtx,g2003h01.gtx,g2003p01.gtx +vunits=m +no_defs<br>
<br>
</div>Folks,<br>
<br>
The problem is that ProjCoordTransGeoKey is seto to CT_TransverseMercator<br>
but there are no actual projection method parameters. When the Geokey<br>
interpreter sees ProjCoordTransGeoKey it looks for a user defined projected<br>
coordinate system and ignores the predefined one (PCS_NAD83_UTM_zone_15N).<br>
<br>
I consider this an improperly defined set of geokeys. If ProjCoordTransGeoKey<br>
is provided then the projection method parameters should also be provided.<br>
Alternatively drop the ProjCoordTransGeoKey and let the system use<br>
the ProjectedCSTypeGeoKey.<br>
<br>
Best regards,<br>
<font color="#888888">--<br>
---------------------------------------+--------------------------------------<br>
I set the clouds in motion - turn up | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a><br>
light and sound - activate the windows | <a href="http://pobox.com/~warmerdam" target="_blank">http://pobox.com/~warmerdam</a><br>
and watch the world go round - Rush | Geospatial Software Developer<br>
</font><div><div></div><div>_______________________________________________<br>
Liblas-devel mailing list<br>
<a href="mailto:Liblas-devel@lists.osgeo.org" target="_blank">Liblas-devel@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/liblas-devel" target="_blank">http://lists.osgeo.org/mailman/listinfo/liblas-devel</a><br>
<br>
</div></div></blockquote></div></div></div>Unfortunately most of the data on the internet follows this convention. Is there an easy workaround to this other than manually specifying parameters?
</div><br>