GDAL读取shp文件VS2015+GDAL2.1

因为项目需要,我们想从SHP文件中,读取给定的某个区域内的所有建筑物信息。

本来想用QGIS来实现,但是很不巧,我们的一个服务器和我的工作电脑先后驾崩,没有QGIS的环境了。

所以我想用别的办法解决,然后想到了GDAL开源库,前面我已经写了一篇关于GDAL的编译方法了,再在正好派上用场。

一、编译GDAL

编译方法我就不说了,请查看我的博客。

http://blog.csdn.net/octdream/article/details/71106474

二、GDAL环境配置

我使用的是VS2015开发环境

1.配置附加包含目录

右击项目->属性->配置属性->C/C++->常规->附加包含目录


2.配置附加库目录

右击项目->属性->配置属性->链接器->常规->附加目录


3.添加附加依赖项

右击项目->属性->配置属性->链接器->输入->附加依赖项


三、使用GDAL读取shp文件数据

在网上找了很多种版本,但是有的编译通不过,要不就是打开文件出错,后面总结出方法,不多说了直接上代码:

#include "ogrsf_frmts.h"
#include <iostream>
#include <list>
using namespace std;

int main()
{
	GDALAllRegister();
	//OGRRegisterAll();

	GDALDataset       *poDS;
	//poDS = (GDALDataset*)GDALOpenEx("E:\\Qt\\maps\\cs_map\\school_point.shp", GDAL_OF_VECTOR, NULL, NULL, NULL);
	poDS = (GDALDataset*)GDALOpenEx("E:\\Qt\\maps\\bj_map\\Building_surface.shp", GDAL_OF_VECTOR, NULL, NULL, NULL);

	if (poDS == NULL)
	{
		printf("Open failed.\n");
		exit(1);
	}
	OGRLayer  *poLayer;
	poLayer = poDS->GetLayerByName("Building_surface");

	OGRFeature *poFeature;
	poLayer->ResetReading();
	while ((poFeature = poLayer->GetNextFeature()) != NULL)
	{
		OGRFeatureDefn *poFDefn = poLayer->GetLayerDefn();
		int iField;
		int i = poFDefn->GetFieldCount(); //获得字段的数目,本实例返回5,不包括前两个字段(FID,Shape),这两个字段在arcgis里也不能被修改;
		for (iField = 0; iField < poFDefn->GetFieldCount(); iField++)
		{
			OGRFieldDefn *poFieldDefn = poFDefn->GetFieldDefn(iField);
			//根据字段值得类型,选择对应的输出
			if (poFieldDefn->GetType() == OFTInteger)
				printf("%d,", poFeature->GetFieldAsInteger(iField));
			else if (poFieldDefn->GetType() == OFTReal)
				printf("%.3f,", poFeature->GetFieldAsDouble(iField));
			else if (poFieldDefn->GetType() == OFTString)
				printf("%s,", poFeature->GetFieldAsString(iField));
			else
				printf("%s,", poFeature->GetFieldAsString(iField));
		}

		OGRGeometry *poGeometry;
		poGeometry = poFeature->GetGeometryRef();
		if (poGeometry != NULL)
		{
			OGRwkbGeometryType pGeoType = poGeometry->getGeometryType();
			if (pGeoType == wkbPoint)
			{
				OGRPoint *poPoint = (OGRPoint *)poGeometry;

				printf("%.3f,%.3f\n", poPoint->getX(), poPoint->getY());
			}
			else if (pGeoType == wkbPolygon)
			{
				OGRPolygon *pPolygon = (OGRPolygon*)poGeometry->clone();
				OGRLinearRing * poLR = pPolygon->getExteriorRing();
				//int Num = poLR->getNumPoints()

			}
			else if (pGeoType == wkbMultiPolygon)  //这里就是带空洞多边形判断
			{
				OGRMultiPolygon *pMulPolygon = (OGRMultiPolygon*)poGeometry;
				OGRPolygon *pPolygon = NULL;
				for (int i = 0; i < pMulPolygon->getNumGeometries(); i++)
				{
					pPolygon = (OGRPolygon*)pMulPolygon->getGeometryRef(i);
				}
			}
			else
			{
				printf("no point geometry\n");
			}
		
		}
		
		OGRFeature::DestroyFeature(poFeature);
	}

	//system("pause");
	getchar();
    return 0;
}

只上了大概的代码,后面的功能大学自己拓展吧!


版权声明:本文为octdream原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/octdream/article/details/74327887

智能推荐

Ubuntu 14.04 下,安装 Java8

下载Java http://www.oracle.com/technetwork/java/javase/downloads/jdk8-downloads-2133151.html 打开上述链接, , 下载 特定的 Java8 版本。笔者选择的是 jdk-8u172-linux-x64.tar.gz。 - 配置Java环境变量 将 jdk-8u172-linux-x64.tar.gz 解压至指定目...

OpenCV学习之路(五)图像的几何变换

在这一章将要学习图像的移动、旋转,仿射变换等 扩展缩放 我们如果想要改变图像的大小,我们就需要对图像进行扩展缩放,opencv提供给我们控制扩展缩放的函数: 参数解释: src:进行扩展缩放的原图片 dst:可以在此处设置缩放因子,也可手动设置尺寸 interpolation:在缩放时我们推荐使用cv2.INTER_AREA, 在扩展时我们推荐使用cv2.INTER_CUBIC(慢) 和 cv2....

2018.8.27

2018.8.27...

HTML 表单元素的基本样式

HTML 表单元素的基本样式 原创 ixygj197875 发布于2018-02-22 17:48:53 阅读数 2296 收藏 更新于2018-05-20 15:35:58 分类专栏: 揭秘 CSS 揭秘 CSS 收起 表单元素主要包括 label、input、textarea、select、datalist、******、progress、meter、output等,以及对表单元素进行分组的 ...

php输出语句

php输出语句 常见的输出语句 echo(): 可以一次输出多个值,多个值之间用逗号分隔。echo是语言结构(language construct),而并不是真正的函数,因此不能作为表达式的一部分使用。 print(): 函数print()打印一个值(它的参数),如果字符串成功显示则返回true,否则返回false。 print_r(): 可以把字符串和数字简单地打印出来,而数组则以括起来的键和值...

猜你喜欢

工厂模式

简介 常见的实例化对象模式。 用工厂方法替代new操作的一种模式。 当我们使用new操作实例化对象时,调用构造函数完成初始化。若初始化仅是进行赋值等简单的操作,写入构造函数即可。但如果初始化时需要执行一长串复杂的代码,将多个工作装入一个方法,是不妥的。 创建实例与使用实例分离。将创建实例所需的大量初始化工作从基类的构造函数中分离出去。 简单工厂模式、工厂方法模式针对的是一个产品等级结构;而抽象工厂...

B1105 Spiral Matrix (画图)

B1105 Spiral Matrix (25分) //第一次只拿了21分 矩阵的长和宽,求最大因子,从sqrt(num)开始枚举. 每次循环一次,s++,t--,d--,r++ 测试点四运行超时,是因为输入一个数字的时候,需要直接输出这个数字。//1分 测试点二运行超时,最后一个数字不必再while循环一次,直接输出即可。//3分 最后一个测试点卡了好久/(ㄒoㄒ)/~~ 螺旋矩阵...

Java基础=>String,StringBuffer与StringBuilder的区别

字符串常量池 什么是字符串常量池? JVM为了减少字符串对象的重复创建,其维护了一块特殊的内存,这段内存被称为字符串常量池(存储在方法区中)。 具体实现 当代码中出现字符串时,JVM首先会对其进行检查。 如果字符串常量池中存在相同内容的字符串对象,如果有,则不再创建,直接返回这个对象的地址返回。 如果字符串常量池中不存在相同内容的字符串对象,则创建一个新的字符串对象并放入常量池,并返回新创建的字符...

java调用其他java项目的Https接口

项目中是这样的: 用户拿出二维码展示,让机器识别二维码, 机器调用开门的后台系统接口, 然后开门的后台系统接口需要调用管理系统的接口, 管理系统需要判断能不能开门.这两个系统是互相独立的.当时使用http调用是没有问题的.当时后来要求必须用https.废话不说,直接代码: 我的项目中调用的是 HttpsUtils.Get(utlStr) 这个接口 开门系统接口如下图:   管理系统的接口...

Hadoop1.2.1全分布式模式配置

一 集群规划 主机名            IP                               安装的软件 &nbs...