前言
最近在帮同学写几个小程序,包括矢量数据求交(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");//矢量数据求