gdal--矢量求交

前言

最近在帮同学写几个小程序，包括矢量数据求交（Intersection）。网上相关文章很少，所以我参考GDAL官网写了一个完整的程序（裁剪、合并、擦除、更新等功能只要修改代码的求交函数部分，下面会有说明）。

---

代码

#include "gdal_priv.h"  
#include "ogrsf_frmts.h" //for ogr 

bool VectorIntersection(const char *pszSrcShp, const char *pszMethodShp, const char *pszDstShp, const char* pszFormat);

int main()
{
    const char *pszSrcFile = "C:\\Users\\liuwei\\Desktop\\CreateIm\\shpManage\\SrcLayer.shp";//原始文件
    const char *pszMethodFile = "C:\\Users\\liuwei\\Desktop\\CreateIm\\shpManage\\MethodLayer.shp";//用来求交的文件
    const char *pszOutFile = "C:\\Users\\liuwei\\Desktop\\shp\\Intersection.shp";//求交后的结果文件
    VectorIntersection(pszSrcFile, pszMethodFile, pszOutFile, "ESRI Shapefile");//矢量数据求