charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=3DContent-Type content=3D"text/html; charset=3Diso-8859-1"=
>
<META content=3D"MSHTML 6.00.2712.300" name=3DGENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=3D#ffffff>
<DIV><FONT face=3DArial size=3D2>Frank,</FONT></DIV>
<DIV><FONT face=3DArial size=3D2></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2>I have added the following few lines to=20
gt_wkt_srs.cpp. It seemed a shame not to take advantage of the EPSG data if=
it=20
was present. It would be nice if you could include them. The idea could be=
=20
extended to look for a GEOGCS string in the same file if the PROJCS one=20
failed.</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=3DArial size=3D2></FONT> </DIV>
<DIV><FONT face=3DArial=20
size=3D2>/* &nbs=
p; =
=20
GTIFSetFromOGISDefn() =
&nb=
sp; =20
*/<BR>/* &=
nbsp; &nbs=
p; &=
nbsp; &nbs=
p; &=
nbsp; =20
*/<BR>/* Write GeoTIFF projection tags from a=
n OGC=20
WKT definition. =20
*/<BR>/********************************************************************=
****/</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=3DArial size=3D2>int GTIFSetFromOGISDefn( GTIF * psGTIF, co=
nst char=20
*pszOGCWKT )</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=3DArial size=3D2>{<BR>............</FONT></DIV>
<DIV><FONT face=3DArial size=3D2>/*=20
--------------------------------------------------------------------=20
*/<BR>/* Handle the projection=20
transformation. =
&nb=
sp; =20
*/<BR>/* ------------------------------------------------------------------=
--=20
*/<BR> const char *pszProjection =3D poSRS->GetAttrVal=
ue(=20
"PROJECTION" );<BR>// RFJ<BR> // Find the EPSG projected cs code from =
the=20
wkt proj cs name<BR> if (nPCS =3D atoi(CSVGetField(=20
CSVFilename("horiz_cs.csv"),<BR> &=
nbsp; &nbs=
p; &=
nbsp;=20
"HORIZCS_EPSG_NAME", poSRS->GetAttrValue("PROJCS"),=20
CC_ApproxString,<BR> &=
nbsp; &nbs=
p; =20
"HORIZCS_CODE")))<BR> {<BR> GTIFKeySet(psGTIF,=20
ProjectedCSTypeGeoKey, TYPE_SHORT, 1, nPCS );<BR> GTIFKeySet(psG=
TIF,=20
GTModelTypeGeoKey, TYPE_SHORT, 1,=20
ModelTypeProjected);<BR> }<BR> else</FONT></DIV>
<DIV><FONT face=3DArial size=3D2>// RFJ<BR> if( pszProjec=
tion =3D=3D=20
NULL )</FONT></DIV>
<DIV><FONT face=3DArial size=3D2></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2>Roger</DIV>
<DIV><BR></DIV></FONT></BODY></HTML>