OGRGeometry和geos::Geometry之间的相互转换

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的时候是不是会出现，后来使用熟悉了就没有问题了。