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