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( "Creation of output file failed.\n" );exit( 1 );}OGRLayer *poLayer;poLayer = poDS->CreateLayer( pShpName, NULL, wkbPolygon, NULL );if( poLayer == NULL ){printf( "Layer creation failed.\n" );exit( 1 );}OGRFieldDefn oField( "Name", OFTString );oField.SetWidth(32);if( poLayer->CreateField( &oField ) != OGRERR_NONE )  {printf( "Creating Name field failed.\n" );exit( 1 );}char szName[33] = "hello geos";OGRFeature *poFeatureG1, *poFeatureG2, *poFeatureG3;poFeatureG3 = OGRFeature::CreateFeature(poLayer->GetLayerDefn());poFeatureG3->SetField("Name", szName);// Geometry1poFeatureG1 = OGRFeature::CreateFeature( poLayer->GetLayerDefn() );OGRPolygon* pPolygon1 = new OGRPolygon();OGRLinearRing* pLinearRing = new OGRLinearRing();pLinearRing->addPoint(0, 0);pLinearRing->addPoint(2, 0);pLinearRing->addPoint(2, 2);pLinearRing->addPoint(0, 2);pLinearRing->addPoint(0, 0);pPolygon1->addRing(pLinearRing);poFeatureG1->SetGeometry( pPolygon1 ); OGRGeometry* baseGeo = poFeatureG1->GetGeometryRef();// Geometry2poFeatureG2 = OGRFeature::CreateFeature(poLayer->GetLayerDefn());OGRPolygon* pPolygon2 = new OGRPolygon();OGRLinearRing* pLinearRing2 = new OGRLinearRing();pLinearRing2->addPoint(1, 0);pLinearRing2->addPoint(3, 0);pLinearRing2->addPoint(3, 3);pLinearRing2->addPoint(1, 3);pLinearRing2->addPoint(1, 0);pPolygon2->addRing(pLinearRing2);poFeatureG2->SetGeometry(pPolygon2);OGRGeometry* otheGeo = poFeatureG2->GetGeometryRef();int choose = atoi(argv[2]);OGRGeometry* pRet = NULL;switch(choose){case 1:pRet = baseGeo->Difference(otheGeo);   // baseGeo - otheGeobreak;case 2:pRet = otheGeo->Difference(baseGeo);   // otheGeo - baseGeobreak; case 3:pRet = baseGeo->Intersection(otheGeo); // baseGeo 和 otheGeo 相交break;default:break;}poFeatureG3->SetGeometry(pRet);if( poLayer->CreateFeature( poFeatureG3 ) != OGRERR_NONE ){printf( "Failed to create feature in shapefile.\n" );exit( 1 );}OGRFeature::DestroyFeature( poFeatureG1 );OGRFeature::DestroyFeature( poFeatureG2 );OGRFeature::DestroyFeature( poFeatureG3 );GDALClose( poDS );
}

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

像“今日头条”这样的模块滑动特效，大概咋写，求指点？---酷课堂iOS交流群问答整理（201805期）
酷课堂iOS交流群我们是一个什么样的组织: 酷课堂iOS交流群,聚集了一群热爱技术.有趣.有料,平均Q龄在10年以上的"老司机",他们遍布在全国/球各地,有知名企业iOS工程师. ...
面试题1：一个人岁数的3次方是四位数，四次方是六位数，并知道此人岁数的3次方和4次方用遍了0--9十个数字，请编程写程序求此人的岁数
代码如下: package com.javase.test; /** * @author jt *2015-10-21 下午5:14:46 * */ //一个人岁数的3次方是四位数,四次方是六 ...
【小米校招笔试】假如已知有n个人和m对好友关系（存于数字r）。如果两个人是直接或间接的好友（好友的好友的好友...），则认为他们属于同一个朋友圈，请写程序求出这n个人里一共有多少个朋友圈。
2016年小米校招笔试第三题(西安站) 3 假如已知有n个人和m对好友关系(存于数字r).如果两个人是直接或间接的好友(好友的好友的好友...),则认为他们属于同一个朋友圈,请写程序求出这n个人里一共 ...
大量的QT控件及示例发放，求泛滥
大量的QT控件及示例发放,求泛滥 2017年01月13日 14:36:22 阅读数:781 在论坛中,看到一位朋友介绍http://qt-apps.org/站上有很多好的资料和示例,现在发出来,与大 ...
创建方法求两个数的最大值max2，随后再写一个求3个数的最大值的函数max3。要求：在max3这个函数中，调用max2函数，来实现3个数的最大值计算
//创建方法求两个数的最大值max2,随后再写一个求3个数的最大值的函数max3.//要求:在max3这个函数中,调用max2函数,来实现3个数的最大值计算public static int max3 ...
写一个判断方程根的C语言程序,写出求方程ax2+bx+c=0的根的算法,画出相应的程序框图,并要求输出它的实根....
写出求方程ax2+bx+c=0的根的算法,画出相应的程序框图,并要求输出它的实根. 答案思路分析:主要考查含参问题的讨论方法及条件嵌套结构的应用.输入三个实数a,b,c,先判断a是否为零,当a≠0时 ...
python调用自己写的函数_python调用自定义函数的实例操作
在python中,想要调用自定义函数必须先声明,然后才能调用.使用函数时,只要按照函数定义的形式,向函数传递必需的参数,就可以调用函数完成相应的功能或者获得函数返回的处理结果. (1)声明函数 pyt ...
java 编辑二进制文件_java写二进制文件，不能直接编辑，该怎么操作？
java写二进制文件,不能直接编辑,该怎么操作? 关注:138 答案:3 mip版解决时间 2021-01-14 21:48 提问者醒不来的梦 2021-01-14 07:45 我现在想做写一个 ...
动态规划系列---求数组中两个元素差的最大值
题目求数组中两个元素差的最大值(后面的元素减去前面的元素):对应实际生活中的股票买卖,找出一只股票走势里面可能的最大收益: 思路类似于求数组连续和的最大值: 保存最大差值和最小值,遍历数组,如果 ...

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

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

最新文章

热门文章