Hi,

I need to find out if two WKTs have the same horizontal coordinate systems.  I 
found OSRIsSameVertCS() but I need some like it for horizontal CS.

I'm thinking something like this should do the trick:

   ref1 = OSRNewSpatialReference(wkt1);
   if (OSRIsCompound(ref1))
      ref1 = <get horizontal cs>(ref1);

   ref2 = OSRNewSpatialReference(wkt2);
   if (OSRIsCompound(ref2))
      ref2 = <get horizontal cs>(ref2);

   hasSameHorizontalCS = OSRIsSame(ref1, ref2);

However, I don’t see a way to extract the horizontal cs out of a compound cs.

Questions:
1) Is there a way to extract the horizontal coordinate system?  I’m using GDAL 
1.9.
2) Is there a better way to solve my problems?

-john
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to