QT 用GDAL库读写shp文件

地理信息系统离不开shp文件,自己写程序调用gdal库解析shp文件,看看shp文件里的内容。

一、GDALAllRegister()注册所有驱动

二、GDALOpenEx()打开shp文件

三、调用GDALDataset->GetLayerCount()获取图层数量

四、调用GDALDataset->GetLayer(?)获取图层

五、调用OGRLayer->GetLayerDefn()获取图层的属性表结构

六、调用OGRFeatureDefn->GetFieldCount()获取图层的字段个数

七、调用OGRFieldDefn->GetNameRef(?)字段字段名

八、调用OGRLayer->GetNextFeature()取要素

九、调用OGRFeature->GetGeometryRef()取几何对象poGeometry

十、调用OGRGeometry->getGeometryType()取几何对象类型(点线面)

1.OGRPoint:点

OGRPoint *poPoint = (OGRPoint *)poGeometry;

2.OGRLineString:线

OGRLineString* pLineGeo = (OGRLineString*)poGeometry;

3.OGRPolygon :多边形

OGRPolygon *poPolygon = (OGRPolygon *)poGeometry;

4.OGRMultiPoint :点集合

OGRMultiPoint *poMultiPoint = (OGRMultiPoint *)poGeometry;

5.OGRMultiPolygon :多边形集合

OGRMultiPoint *poMultiPoint = (OGRMultiPoint *)poGeometry

6.OGRMultiLineString :线集合

OGRMultiLineString* poMultiLine = (OGRMultiLineString *)poGeometry;

完整程序(gdal完整库和头文件):https://download.csdn.net/download/x879014419/89470363

  • 5
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是一个简单的示例代码,使用GDAL加载显示shp文件: ```python from osgeo import ogr from osgeo import gdal from osgeo import osr import matplotlib.pyplot as plt # 打开shapefile文件 driver = ogr.GetDriverByName('ESRI Shapefile') shpfile = driver.Open('path/to/your/file.shp', 0) layer = shpfile.GetLayer() # 获取图层的空间参考 spatialRef = layer.GetSpatialRef() print('空间参考:', spatialRef.ExportToWkt()) # 获取图层的范围 extent = layer.GetExtent() print('图层范围:', extent) # 获取图层的投影信息 proj4 = spatialRef.ExportToProj4() print('投影信息:', proj4) # 将图层转换为gdal格式 gdal.SetConfigOption('GDAL_FILENAME_IS_UTF8', 'NO') gdal.SetConfigOption('SHAPE_ENCODING', '') mem_drv = ogr.GetDriverByName('Memory') mem_ds = mem_drv.CreateDataSource('memData') mem_layer = mem_ds.CreateLayer('layer', srs=spatialRef, geom_type=layer.GetGeomType()) # 将原始图层的属性表结构复制到新图层中 layer_def = layer.GetLayerDefn() for i in range(layer_def.GetFieldCount()): field_def = layer_def.GetFieldDefn(i) mem_layer.CreateField(field_def) # 将原始图层中的要素复制到新图层中 for feat in layer: mem_layer.CreateFeature(feat) # 将gdal格式的图层转换为numpy数组 driver = gdal.GetDriverByName('MEM') dst_ds = driver.CreateCopy('', mem_ds, 0) arr = dst_ds.GetRasterBand(1).ReadAsArray() # 显示图层 plt.imshow(arr, cmap='gray') plt.show() ``` 说明: 1. 首先通过ogr打开shp文件,并获取图层的空间参考、范围和投影信息。 2. 然后将图层转换为gdal格式,并将其转换为numpy数组。 3. 最后使用matplotlib显示数组。 需要注意的是,该示例代码仅适用于包含单一要素类型的shp文件。如果要处理包含多个要素类型的shp文件,需要对代码进行修改。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值