<P>Hello, I write some code below to compare OGRSpatialReference,</P>
<P>bool CompareWktIsUniform(const char* pszWkt1, const char* pszWkt2)<BR>{<BR> OGRSpatialReference SR1(pszWkt1);<BR> OGRSpatialReference SR2(pszWkt2);</P>
<P> int IsSame = SR1.IsSame(&SR2);</P>
<P> if (IsSame == TRUE)<BR> return true;<BR> else<BR> return false;<BR>}</P>
<P>And I have two OGRSpatialReference, one is from shapefile and the other is from Geotiff. </P>
<P>PROJCS["Albers Conical Equal Area",GEOGCS["GCS_Krasovsky_1940",DATUM["Krasovsky_1940",SPHEROID["Krasovsky_1940",6378245.0,298.3]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["longitude_of_center",105.0],PARAMETER["Standard_Parallel_1",25.0],PARAMETER["Standard_Parallel_2",47.0],PARAMETER["latitude_of_center",0.0],UNIT["Meter",1.0]]</P>
<P><BR>PROJCS["IMAGINE GeoTIFF Support<BR>Copyright 1991 - 2005 by Leica Geosystems Geospatial Imaging, LLC. All Rights Reserved<BR>@(#)$RCSfile: egtf.c $ IMAGINE 9.0 $Revision: 10.0 $ $Date: 2005/07/26 15:10:00 EST $<BR>Projection = Albers Conical Equal Area",GEOGCS["Krasovsky",DATUM["Krasovsky",SPHEROID["Krasovsky",6378245,298.3000003760163]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["standard_parallel_1",25],PARAMETER["standard_parallel_2",47],PARAMETER["latitude_of_center",0],PARAMETER["longitude_of_center",105],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]]]</P>
<P>But they are recognized the same in ArcGis and Erdas, in ArcGis they are both shown as below:</P>
<P>Projection: Albers<BR>False_Easting: 0.000000<BR>False_Northing: 0.000000<BR>Central_Meridian: 105.000000<BR>Standard_Parallel_1: 25.000000<BR>Standard_Parallel_2: 47.000000<BR>Latitude_Of_Origin: 0.000000<BR>Linear Unit: Meter (1.000000)</P>
<P>Geographic Coordinate System: GCS_Krasovsky_1940<BR>Angular Unit: Degree (0.017453292519943299)<BR>Prime Meridian: Greenwich (0.000000000000000000)<BR>Datum: D_Krasovsky_1940<BR> Spheroid: Krasovsky_1940<BR> Semimajor Axis: 6378245.000000000000000000<BR> Semiminor Axis: 6356863.018773047300000000<BR> Inverse Flattening: 298.300000000000010000</P>
<P>So, is there any way to get the same results? or what is lost?</P><br><br><span title="neteasefooter"/><hr/>
<a href="http://ym.163.com/?from=od1" target="_blank">有域名马上来，网易免费送你200户自主域名邮箱</a>
</span>