编译命令:g++ main.cpp -lgdal
调用命令:./a.out 输出shp名称 操作选项
注释:操作选项(1:多边形A - 多边形B,2:B - A,3:A和B的交集部分)
#include "ogrsf_frmts.h"
#include <iostream>
using namespace std;
int main(int argc, char* argv[])
{
const char *pszDriverName = "ESRI Shapefile";
char *pShpName = argv[1];
GDALDriver *poDriver;
GDALAllRegister();
poDriver = GetGDALDriverManager()->GetDriverByName(pszDriverName );
if( poDriver == NULL )
{
printf( "%s driver not available.\n", pszDriverName );
exit( 1 );
}
GDALDataset *poDS;
char cShpName[20];
sprintf(cShpName, "%s.shp", pShpName);
poDS = poDriver->Create( cShpName, 0, 0, 0, GDT_Unknown, NULL );
if( poDS == NULL )
{
printf( "Cre