OGRGeometry属于OGR库，OGR库通常是和GDAL一起封装供用户使用的。GDAL&OGR库支持了大多数地理数据格式的读写，通常用户使用OGR库可以完成大部分矢量相关的任务。但是OGR并没有对特征之间的空间关系计算提供原生支持。尤其是空间矢量之间的拓扑关系计算，其内部是采用的第三方库geos。所以有时候我们干脆直接使用geos进行矢量运算。这时候就有可能遇到OGRGeometry和geos::Geometry之间相互转换的问题。
本文只列出了OGRGeometry转geos::Geometry。反过来原理是一样的。
两者之间可以通过wkt和wkb相互转换。
第二个方案采用wkb，这方法稍显麻烦，因为geos的WKBReader类是这样定义的：
class GEOS_DLL WKBReader { public: geom::Geometry* read(std::istream &is);输入参数是std::istream，我整了一下午不知道如何构造std::istream。当然可以落盘，形成磁盘文件，但那不是我想要的。于是采用第二个方案：采用wkb，geos的WKTReader定义read函数如下：
class GEOS_DLL WKTReader { public: geom::Geometry* read(const std::string &wellKnownText);我们发现，传入的是string。
而OGRGeometry恰好有这么个函数：
class CPL_DLL OGRGeometry { public: virtual OGRErr exportToWkt( char ** ppszDstText, OGRwkbVariant=wkbVariantOldOgc ) const = 0;这就给我们了一个思路，顺着思路实现：
geos::io::WKTReader wkt_read; OGRFeature* fea = poLayer->GetFeature(i); char * buff = nullptr; fea->GetGeometryRef()->exportToWkt(&buff); wkt_read.read(buff); geoms = GeometryPtr(wkt_read.read(buff));经过测试，该方法是可行的。顺着这个思路就可以实现geos::Geometry转到OGRGeometry。
多说一句，有时候使用geos库会出现莫名其妙的崩溃问题。这些问题无非是内存重复销毁引起的。新手多看看geos的代码注释，非常有营养。尤其是在涉及到指针传递的时候，看看注释里有没有写使用者获得指针所有权等等。
这种情况在我刚接触geos的时候是不是会出现，后来使用熟悉了就没有问题了。
