<DIV>hello, everyone.</DIV>
<DIV>I need to justify whether one geotiff and one shapefile are in the same projection with IsSame function from OGR, which will compare the wkt info of OGRSpatialReference , like below.<BR><BR>bool CompareWktIsUniform(const char* pszWkt1, const char* pszWkt2) <BR>{ <BR> OGRSpatialReference SR1(pszWkt1); <BR> OGRSpatialReference SR2(pszWkt2); <BR><BR> int IsSame = SR1.IsSame(&SR2); <BR><BR> if (IsSame == TRUE) <BR> return true; <BR> else <BR> return false; <BR>} <BR><BR>And I have two OGRSpatialReference, one is from shapefile and the other is from Geotiff. <BR><BR>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]] <BR><BR><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"]]] <BR></DIV>
<DIV>In fact they should be same, as both are </DIV>
<DIV>Projection = Albers Conical Equal Area", </DIV>GEOGCS["Krasovsky
<DIV>DATUM["Krasovsky<BR></DIV>
<DIV>but in details are a little different, how to avoid that?<BR><BR> </DIV><br><br><span title="neteasefooter"/></span>