OGR示例:写shp，求面与面的交和差操作

编译命令：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