地理数据处理_矢量图处理_python_OGR

最新推荐文章于 2024-07-30 04:45:00 发布
最新推荐文章于 2024-07-30 04:45:00 发布
阅读量1.9k 收藏 29
点赞数 7
分类专栏： ArcGIS
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Mikowoo007/article/details/103880371
版权

OGR
简单要素库是地理空间数据抽象库(Geospatial Data Abstraction Library,GDAL)的一部分，它是一个非常流行的用于读写空间数据的开源库
在这里插入图片描述

环境
python 2.7

安装

conda install gdal pyproj

读取矢量数据 shapefile

在这里插入图片描述

import sys
from osgeo import ogr
fn = r'C:/Users/miko/Documents/Book_code/python地理数据处理_数据集/python地理数据处理_数据集/osgeopy-data/global/ne_50m_populated_places.shp'
ds = ogr.Open(fn,0)

if ds is None:
    sys.exit('Could not open {0}.'.format(fn))
lyr = ds.GetLayer(0)

i=0
for feat in lyr:
    pt = feat.geometry()
    x = pt.GetX()
    y = pt.GetY()
    name = feat.GetField('NAME')
    pop = feat.GetField('POP_MAX')
    print(name,pop,x,y)
    i += 1
    if i == 10:
        break
del ds

在这里插入图片描述

打开数据源

ogr.open()
通过Open函数传递文件名和一个可选更新标记，打开一个数据源。这是OGR模块中的一个独立函数，需要事先用这个模块名称定义一个函数名，以便Python可以找到它。
如果第二个变量不默认设置为0，将以只读的模式打开文件。可以通过将其设置为1或者True来以编辑的模式打开
如果shapefile文件无法打开，则Open函数返回为空

获取图层信息

ds.GetLayer()
一个数据源由一个或多个图层组成，所以打开数据源后，需要获得具体的图层
可以通过该函数来获取图层索引或者图层名，进而得到这个数据集中对应的图层
0索引
如果不给GetLayer函数提供任何参数，默认返回这个数据源中的第一个图层

获取具体图层中要素的信息

图层由一个或多个要素组成，每一个要素表示一个地理对象。地理对象的几何形状和属性值都存储在这些要素中。

for feat in lyr:
	pt = feat.geometry()
	x = pt.GetX()
    y = pt.GetY()
    name = feat.GetField('NAME')
    pop = feat.GetField('POP_MAX')

feat.geometry()
获取要素的几何形状，从而获取其坐标信息
feat.GetField()
通过属性名或索引值获取对应字段的值，该函数返回的数据类型与底层数据集中的数据类型一致。补充函数：feat.GetfieldAsString()

访问特定要素

import sys
from osgeo import ogr
fn = r'C:/Users/miko/Documents/Book_code/python地理数据处理_数据集/python地理数据处理_数据集/osgeopy-data/global/ne_50m_populated_places.shp'
ds = ogr.Open(fn,0)

if ds is None:
    sys.exit('Could not open {0}.'.format(fn))
lyr = ds.GetLayer(0)

num_features = lyr.GetFeatureCount()               #获取要素总数
last_feature = lyr.GetFeature(num_features-1)      #获取图层中的最后一个要素
print(num_features)
print(last_feature.name)

在这里插入图片描述

lyr.GetFeatureCount() 获取要素的总数
lyr.GetFeature(FIDs) 获取FID号对应的要素

ospybook 模块

环境 python2.7
安装

pip install https://github.com/cgarrard/osgeopy-code/raw/master/ospybook-latest.zip

查看属性

import ospybook as pb
fn = r'C:/Users/miko/Documents/Book_code/python地理数据处理_数据集/python地理数据处理_数据集/osgeopy-data/global/ne_50m_populated_places.shp'
pb.print_attributes(fn,3,['NAME','POP_MAX'])

在这里插入图片描述
print_attributes(lyr_or_fn,[n],[fields],[geom],[reset])
该函数将属性值信息输出到屏幕上
这个函数对于查看小数据量的属性信息运行正常，但如果使用它来输出一个大数据的所有属性，你会后悔的

参数含义
lyr_or_fn可以是一个图层，也可以是一个数据源路径
如果是一个数据源，使用的是第一个图层
n是一个可选值
用于设置输出记录数目，默认输出所有数值
fields是一个可选值
用于设置输出结果中包含的属性字段列表，默认包括所有字段
geom是一个可选的布尔值
用于设置是否输出集合要素类型，默认为True
reset是一个可选的布尔值
用于设置在输出数值之前，是否重置到第一条记录，默认为True

绘制空间数据

需创建VectorPlotter类的新实例，需为构造函数传递一个布尔变量

import os
os.chdir('C:\Users\miko\Documents\Book_code\osgeopy-data\global')  # 更改工作目录，最好是英文路径
from ospybook.vectorplotter import VectorPlotter

vp = VectorPlotter(True)
vp.plot('ne_50m_admin_0_countries.shp',fill=False)  # fill参数使shp文件用一个空心多边形表示
vp.plot('ne_50m_populated_places.shp','bo')  # bo符号表示蓝色的圆圈
vp.draw()

在这里插入图片描述
ArcGIS 桌面版本打开效果：
在这里插入图片描述

plot(self,geom_or_lyr,[symbol],[name],[kwargs])

参数含义
geom_or_lyr可以是一个要素对象、图层或者指向一个数据源的路径。
如果是一个数据源，数据源中的第一个图层会被绘制显示
symbol是一个可选值
用于设置几何元素的符号样式
name是一个可选值
用于为数据设定名称、以便后期可以访问它
kwargs是一个可选值
通过关键字进行指定
kwargs经常被用作一个不确定数量的关键字参数(an indeterminate number of keyword arguments)的缩写

获取数据的元数据

获取空间范围

图层的空间范围是从上、下、左、右4个方向的最大和最小边界坐标构成的矩形
lyr.GetExtent()
获取一个图层对象的边界坐标，它会返回一个类似(min_x,max_x,min_y,max_y)的元组

import sys
from osgeo import ogr
fn = r'C:/Users/miko/Documents/Book_code/osgeopy-data/Washington/large_cities.geojson'
ds = ogr.Open(fn,0)

if ds is None:
    sys.exit('Could not open {0}.'.format(fn))
lyr = ds.GetLayer(0)

extent = lyr.GetExtent()
print(extent)

在这里插入图片描述
在这里插入图片描述

获取几何类型

lyr.GetGeomType() 返回的是一个整数值，而不是一个可理解的字符串
如果这个图层有多种几何类型，例如点和面，GetGeomType函数就会返回wkbUnknown
在ORG几何要素常量中的WKB前缀为熟知二进制(Well-KnowingBinary,WKB)的缩写，是用于不同软件程序间进行几何要素类型转换的一种二进制表示标准。因为它是二进制格式，所以人们无法直接阅读获得其表示的内容，但是熟知文本格式(Well-Known Text,WKT)可以直接阅读

import sys
from osgeo import ogr
fn = r'C:/Users/miko/Documents/Book_code/osgeopy-data/Washington/large_cities.shp'
ds = ogr.Open(fn,0)

if ds is None:
    sys.exit('Could not open {0}.'.format(fn))
lyr = ds.GetLayer(0)

print(lyr.GetGeomType())
print(lyr.GetGeomType() == ogr.wkbPoint)     #与常量进行比较，获取数据类型
print(lyr.GetGeomType() == ogr.wkbPolygon)

在这里插入图片描述

几何要素类型OGR常量
PointwkbPoint
MulitpointwkbMultiPoint
LinewkbLineString
MultilinewkbMultiLineString
PolygonwkbPolygon
MultipolygonwkbMultiPolygon
Unknown geometry typewkbUnknown
No geometrywkbNone

获取一个几何类型的可读字符串，可以从要素几何类型中进行获取

import sys
from osgeo import ogr
fn = r'C:/Users/miko/Documents/Book_code/osgeopy-data/Washington/large_cities.shp'
ds = ogr.Open(fn,0)

if ds is None:
    sys.exit('Could not open {0}.'.format(fn))
lyr = ds.GetLayer(0)

feat = lyr.GetFeature(0)
print(feat.geometry().GetGeometryName())

在这里插入图片描述

获取空间参考系统

lyr.GetSpatialRef()
得到用WKT格式描述空间参考系统的字符串

import sys
from osgeo import ogr
fn = r'C:/Users/miko/Documents/Book_code/osgeopy-data/Washington/large_cities.shp'
ds = ogr.Open(fn,0)

if ds is None:
    sys.exit('Could not open {0}.'.format(fn))
lyr = ds.GetLayer(0)

print(lyr.GetSpatialRef())

在这里插入图片描述

获取图层自身属性字段信息

lyr.schema

import sys
from osgeo import ogr
fn = r'C:/Users/miko/Documents/Book_code/osgeopy-data/Washington/large_cities.shp'
ds = ogr.Open(fn,0)

if ds is None:
    sys.exit('Could not open {0}.'.format(fn))
lyr = ds.GetLayer(0)

for field in lyr.schema:
    print(field.name,field.GetTypeName())

在这里插入图片描述

矢量数据写入

基于已有的shp文件创建新的shapefile文件

import sys
from osgeo import ogr

ds = ogr.Open(r'C:\Users\miko\Documents\Book_code\osgeopy-data\global',1)
# 打开一个文件夹作为数据源，而不是使用shapefile文件
# 如果文件夹中大多数文件是shapefile格式文件，每个shapefile文件会被视为一个图层
# 参数 1，表示可以在文件夹中创建一个新的图层 (shapefile文件)

if ds is None:
    sys.exit('Could not open folder.')
    
in_lyr = ds.GetLayer('ne_50m_populated_places')   #获取shapefile文件
#将没有后缀.shp，传递给GetLayer函数来获取居住地shapefile文件作为一个图层

if ds.GetLayer('captial_cities'):
    ds.DeleteLayer('captial_cities')
# 因为OGR不会覆盖现有的图层，所以需要检查现有的图层是否已经存在，如果存在，就删掉它

#创建一个点图层
out_lyr = ds.CreateLayer('captial_cities',in_lyr.GetSpatialRef(),ogr.wkbPoint)
out_lyr.CreateFields(in_lyr.schema)

#创建一个空要素
out_defn = out_lyr.GetLayerDefn()
out_feat = ogr.Feature(out_defn)

#复制几何要素和属性
for in_feat in in_lyr:
    if in_feat.GetField('FEATURECLA') == 'Admin-0 capital':
        geom = in_feat.geometry()
        out_feat.SetGeometry(geom)
        for i in range(in_feat.GetFieldCount()):
            value = in_feat.GetField(i)
            out_feat.SetField(i,value)
        out_lyr.CreateFeature(out_feat)
        
del ds

新创建的为图中棕色部分点
在这里插入图片描述

创建新图层

ds.CreateLayer(name,[srs],[geom_type],[options])
创建一个新图层，用来存储输出的数据

参数含义
name表示新建的图层名称
srs表示新建的图层使用的空间参考系统
默认为空，表示没有设置任何空间参考系统
geom_type用来表示图层中要素的几何类型
默认为wkbUnknown
options是一个图层创建时的选项列表，只适合于特定的矢量数据类型
创建空要素

out_defn = out_lyr.GetLayerDefn()
创建一个要素需要获得要素定义，其包含几何类型和所有属性字段的信息，这样创建的空要素才具有相同的属性字段和几何类型
你需要获得待添加要素的图层的要素定义信息，而且要在图层做了添加、删除或更新了任何字段之后进行获取。 如果先获取要素的定义信息，然后属性字段做了任何改变，获得的信息就过时了，这意味着，基于这些过时的信息插入的要素，会与实际不符，程序崩溃

out_feat = ogr.Feature(out_defn)
为新创建的图层添加要素，需要创建一个用于存储几何要素和属性的虚拟要素，然后将其插入到图层中

要素插入图层

out_lyr.CreateFeature()
将要素插入到图层中，这个函数会在图层中保存一个要素的副本，包括添加的所有信息

图层的保存

del ds
删除ds变量强制关闭文件，并将所有的编辑成果写入到磁盘中
ds.SyncToDisk()
保持数据源打开，但编辑的成果写入磁盘中

创建新的数据源

获取驱动

方式一
是从已经打开的数据集中进行获取，这将允许创建一个新的数据源，其和已存在的数据源矢量数据格式一样

from osgeo import ogr
ds = ogr.Open(r'C:/Users/miko/Documents/Book_code/osgeopy-data/global/ne_50m_admin_0_countries.shp')
driver = ds.GetDriver()

print(driver.name)

在这里插入图片描述

from osgeo import ogr
ds = ogr.Open(r'C:/Users/miko/Documents/Book_code/osgeopy-data/Washington/large_cities.geojson')
driver = ds.GetDriver()

print(driver.name)

在这里插入图片描述
方式二
使用OGR中GetDriverByName()函数，传递给它驱动程序的简称
这些驱动程序名称都在OGR网站上有介绍

from osgeo import ogr

esri_fn = r'C:/Users/miko/Documents/Book_code/osgeopy-data/test'

esri_driver = ogr.GetDriverByName('ESRI Shapefile')
esri_ds = esri_driver.CreateDataSource(esri_fn)
if esri_ds is None:
    sys.exit('Could not create {0}.'.format(esri_fn))

新创建的数据源会打开等待写入，剩下就可以为数据源添加图层和要素信息
创建新的数据源，不能覆盖现有的数据源

from osgeo import ogr

esri_fn = r'C:/Users/miko/Documents/Book_code/osgeopy-data/test'

esri_driver = ogr.GetDriverByName('ESRI Shapefile')
esri_ds = esri_driver.CreateDataSource(esri_fn)
if esri_ds is None:
    sys.exit('Could not create {0}.'.format(esri_fn))

ds = ogr.Open(r'C:\Users\miko\Documents\Book_code\osgeopy-data\global',1)
# 打开一个文件夹作为数据源，而不是使用shapefile文件
# 如果文件夹中大多数文件是shapefile格式文件，每个shapefile文件会被视为一个图层
# 参数 1，表示可以在文件夹中创建一个新的图层 (shapefile文件)

if ds is None:
    sys.exit('Could not open folder.')
    
in_lyr = ds.GetLayer('ne_50m_populated_places')   #获取shapefile文件
#将没有后缀.shp，传递给GetLayer函数来获取居住地shapefile文件作为一个图层

if esri_ds.GetLayer('captial_cities'):
    esri_ds.DeleteLayer('captial_cities')
# 因为OGR不会覆盖现有的图层，所以需要检查现有的图层是否已经存在，如果存在，就删掉它

#创建一个点图层
out_lyr = esri_ds.CreateLayer('captial_cities',in_lyr.GetSpatialRef(),ogr.wkbPoint)
out_lyr.CreateFields(in_lyr.schema)

#创建一个空要素
out_defn = out_lyr.GetLayerDefn()
out_feat = ogr.Feature(out_defn)

#复制几何要素和属性
for in_feat in in_lyr:
    if in_feat.GetField('FEATURECLA') == 'Admin-0 capital':
        geom = in_feat.geometry()
        out_feat.SetGeometry(geom)
        for i in range(in_feat.GetFieldCount()):
            value = in_feat.GetField(i)
            out_feat.SetField(i,value)
        out_lyr.CreateFeature(out_feat)

del esri_ds

在这里插入图片描述

使用其他驱动创建数据源

from osgeo import ogr

json_fn = r'C:/Users/miko/Documents/Book_code/osgeopy-data/test'

json_driver = ogr.GetDriverByName('GeoJSON')
json_ds = esri_driver.CreateDataSource(json_fn)
if json_ds is None:
    sys.exit('Could not create {0}.'.format(json_fn))

from osgeo import ogr
driver = ogr.GetDriverByName('SQLite')
ds = driver.CreateDataSource(r'D:/osgeopy-data/global/earth.sqlite',['SPATIALITE=yes'])

新建属性字段

给新建的数据源、图层，创建属性

from osgeo import ogr

in_ds = ogr.Open(r'D:\WorkSpace\Jupyter notebook\data\osgeopy-data-washington\osgeopy-data\Washington',1)

out_esri_fn = r'D:\WorkSpace\Jupyter notebook\data\test'
out_esri_driver = ogr.GetDriverByName('ESRI Shapefile')
out_esri_ds = out_esri_driver.CreateDataSource(esri_fn)


in_lyr = in_ds.GetLayer('large_cities')
print(in_lyr.GetSpatialRef())

if out_esri_ds.GetLayer("test_create_lyr"):
    out_esri_ds.DeleteLayer("test_create_lyr")

out_lyr = out_esri_ds.CreateLayer("test_create_lyr",in_lyr.GetSpatialRef(),ogr.wkbPoint)

coord_fld = ogr.FieldDefn("name",ogr.OFTString)
coord_fld.SetWidth(6)
out_lyr.CreateField(coord_fld)

coord_fld = ogr.FieldDefn("X_default",ogr.OFTReal)
out_lyr.CreateField(coord_fld)

coord_fld = ogr.FieldDefn("Y_default",ogr.OFTReal)
out_lyr.CreateField(coord_fld)

coord_fld = ogr.FieldDefn("X_precision",ogr.OFTReal)
coord_fld.SetWidth(8)
coord_fld.SetPrecision(3)
out_lyr.CreateField(coord_fld)

coord_fld = ogr.FieldDefn("Y_precision",ogr.OFTReal)
coord_fld.SetWidth(8)
coord_fld.SetPrecision(3)
out_lyr.CreateField(coord_fld)

del out_esri_ds

在这里插入图片描述

字段类型
字段数据类型OGR常量
IntegerOFTInteger
List of integersOFTIntegerList
Floating point numberOFTReal
List of floating point numbersOFTRealList
StringOFTString
List of StringOFTStringList
DateOFTDate
Time of dayOFTTime
Date and timeOFTDateTime

定义属性字段的简单写法

coord_fld = ogr.FieldDefn('X_default', ogr.OFTReal)
out_lyr.CreateField(coord_fld)
coord_fld.SetName('Y_default')
out_lyr.CreateField(coord_fld)

给属性表赋值

from osgeo import ogr

in_ds = ogr.Open(r'D:\WorkSpace\Jupyter notebook\data\osgeopy-data-washington\osgeopy-data\Washington',1)

out_esri_fn = r'D:\WorkSpace\Jupyter notebook\data\test'
out_esri_driver = ogr.GetDriverByName('ESRI Shapefile')
out_esri_ds = out_esri_driver.CreateDataSource(esri_fn)


in_lyr = in_ds.GetLayer('large_cities')
print(in_lyr.GetSpatialRef())

if out_esri_ds.GetLayer("test_create_lyr"):
    out_esri_ds.DeleteLayer("test_create_lyr")

out_lyr = out_esri_ds.CreateLayer("test_create_lyr",in_lyr.GetSpatialRef(),ogr.wkbPoint)

coord_fld = ogr.FieldDefn("name",ogr.OFTString)
coord_fld.SetWidth(6)
out_lyr.CreateField(coord_fld)

coord_fld = ogr.FieldDefn("X_default",ogr.OFTReal)
out_lyr.CreateField(coord_fld)

coord_fld = ogr.FieldDefn("Y_default",ogr.OFTReal)
out_lyr.CreateField(coord_fld)

coord_fld = ogr.FieldDefn("X_Short",ogr.OFTReal)
coord_fld.SetWidth(8)
coord_fld.SetPrecision(3)
out_lyr.CreateField(coord_fld)

coord_fld = ogr.FieldDefn("Y_Short",ogr.OFTReal)
coord_fld.SetWidth(8)
coord_fld.SetPrecision(3)
out_lyr.CreateField(coord_fld)

out_feat = ogr.Feature(out_lyr.GetLayerDefn())

for in_feat in in_lyr:
    pt = in_feat.geometry()
    name = in_feat.GetField("name")
    out_feat.SetGeometry(in_feat.geometry())
    out_feat.SetField("name",name)
    out_feat.SetField("X_default",pt.GetX())
    out_feat.SetField("Y_default",pt.GetY())
    out_feat.SetField("X_Short",pt.GetX())
    out_feat.SetField("Y_Short",pt.GetY())
    out_lyr.CreateFeature(out_feat)
del out_esri_ds

在这里插入图片描述

更改图层属性字段定义
需要更改的字段属性OGR常量
Field name onlyALTER_NAME_FLAG
Field type onlyALTER_TYPE_FLAG
Field width and/or precision onlyALTER_WIDTH_PRECISION_FLAG
All of the aboveALTER_ALL_FLAG

更新字段的单个属性信息

from osgeo import ogr

ds = ogr.Open(r'D:\WorkSpace\Jupyter notebook\data\test',1)

lyr = ds.GetLayer("test_create_lyr")

i = lyr.GetLayerDefn().GetFieldIndex("name")
print(i)

fld_defn = ogr.FieldDefn('City_Name',ogr.OFTString)
lyr.AlterFieldDefn(i,fld_defn,ogr.ALTER_NAME_FLAG)

del ds

在这里插入图片描述
更新字段的多个属性信息

from osgeo import ogr

ds = ogr.Open(r'D:\WorkSpace\Jupyter notebook\data\test',1)

lyr = ds.GetLayer("test_create_lyr")

i = lyr.GetLayerDefn().GetFieldIndex("X_default")
print(i)

width = lyr.GetLayerDefn().GetFieldDefn(i).GetWidth()
fld_defn = ogr.FieldDefn('X_coord',ogr.OFTReal)
fld_defn.SetWidth(width)
fld_defn.SetPrecision(4)
flag = ogr.ALTER_NAME_FLAG + ogr.ALTER_WIDTH_PRECISION_FLAG
lyr.AlterFieldDefn(i,fld_defn,flag)

del ds

在这里插入图片描述

复制数据源与文件

import os
import ospybook as pb

original_fn = r'D:\WorkSpace\Jupyter notebook\data\osgeopy-data-washington\osgeopy-data\Washington\large_cities.shp'
new_fn = r'D:\WorkSpace\Jupyter notebook\data\test\large_cities2.shp'

pb.copy_datasource(original_fn, new_fn)

在这里插入图片描述

查看图层属性

import os
import ospybook as pb

original_fn = r'D:\WorkSpace\Jupyter notebook\data\osgeopy-data-washington\osgeopy-data\Washington\large_cities.shp'
new_fn = r'D:\WorkSpace\Jupyter notebook\data\test\large_cities2.shp'

pb.copy_datasource(original_fn, new_fn)

# Open the copied shapefile for writing.
ds = ogr.Open(new_fn, 1)
if ds is None:
    sys.exit('Could not open {0}.'.format(new_fn))
lyr = ds.GetLayer(0)

# Take a look at the attributes before you change anything.
print('Original attributes')
pb.print_attributes(lyr, geom=True)

在这里插入图片描述
在这里插入图片描述

Mikowoo007
关注
  • 7
    点赞
  • 踩
  • 29
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
基于Python的地理信息数据处理.rar
05-07
本压缩包"基于Python的地理信息数据处理.rar"提供了关于如何利用Python进行地理信息数据处理与可视化的实例。以下是其中涉及的关键知识点： 1. **Python编程基础**：Python是这个项目的基础，了解Python的基本语法...
Python地理空间分析指南（第2版）_地理_地理空间分析_python_
10-02
GDAL/OGR提供了一个强大的接口来处理矢量和栅格数据，而Fiona则简化了与GDAL/OGR的交互。Shapely库允许用户创建、操作和分析几何对象，如点、线和多边形。Geopandas是基于Pandas扩展的，为地理数据提供了更高效的...
SuperMap GIS 三维方向矢量数据处理详解
小客服的博客
11-21 2785
SuperMap GIS 三维方向矢量数据处理详解 webgl
GDAL python教程基础篇（2）——用OGR写入矢量数据
qq_45373844的博客
03-08 1751
GDAL python教程基础篇（2）——用OGR写入矢量数据
osgeo库中的ogr模块说明
学习历程记录和分享
05-31 1009
OSGeo库中的ogr模块是GDAL库的一部分，用于处理矢量数据。ogr是OpenGIS Simple Features Reference Implementation的简称，它提供了简单的接口来读取和写入多种矢量地理数据格式。读取和写入矢量数据创建和编辑矢量数据的几何形状查询矢量数据属性执行空间和属性的查询支持的矢量数据格式包括：Shapefile, GeoJSON, KML, ESRI File Geodatabase等。DataSource: 表示整个矢量数据文件或数据集。
Python地理数据处理 七：使用OGR处理几何要素
amyniez的博客，欢迎交流讨论
02-04 2805
目录1. 点集1.1 创建编辑一个点1.2 创建编辑多个点 1. 点集 1.1 创建编辑一个点   OGR中带有额外z坐标的几何要素被认为是2.5D，而不是3D。再OGR进行空间操作时，不考虑z值。   OGR常量表示的不同几何类型： firepit = ogr.Geometry(ogr.wkbPoint) firepit.AddPoint(10, 21) x, y = firepit.GetX(), firepit.GetY() print('{}, {}'.format(x, y)) print(
Python OGR 矢量操作
qq_40292148的博客
08-22 400
矢量数据的读取顺序：文件(datasource) → 图层(layer) → 特征(feature)、字段(field)
Python地理数据处理 六：使用OGR过滤数据
amyniez的博客，欢迎交流讨论
02-02 1896
目录写在前面1. 属性过滤条件 写在前面   过滤条件可以将不想要的要素抛弃，通过过滤条件可以选出符合特定条件的要素，也可以通过空间范围限定要素，这样就可以简单地处理感兴趣的数据。 1. 属性过滤条件   过滤条件需要一个条件语句，类似于SQL语句中的Where子句。如： ‘Population < 50000’ ‘Population = 50000’ ‘Name = “Tokyo”’   注：比较字符串时，需要用引号将字符串包住，并且保证它们与包住整个查询字符串的引号不同。否则，将会导致语法错误
世界道路_道路矢量数据.rar
05-15
为了处理这些数据，你需要GIS软件或编程库，比如ArcGIS、QGIS、GDAL/OGR或者Python的geopandas库。在编程环境中，你可以读取、操作和可视化这些道路数据，甚至与其他数据源结合，进行复杂的分析任务。 总的来说，...
yao'gan.zip_GaN_geotiff python_python tiff_python envi_图像处理软件
09-20
OGR是GDAL的一部分，主要处理矢量数据。 4. **PIL/Rasterio**：PIL是Python中的图像处理库，可以读取和操作各种图像格式，但对GeoTIFF的支持有限。而Rasterio是专门针对地理空间数据的，它基于GDAL，提供了一种更...
GIS数据下载处理系统_Python_PLSQL_下载.zip
04-26
例如，`GDAL/OGR`库用于处理栅格和矢量数据，` Fiona `用于读写GIS矢量数据，`Shapely`用于几何对象操作，`geopandas`提供了Pandas-like接口来处理地理数据，`Cartopy`则用于创建地图投影和地理元素。这些库使得...
GDAL python教程基础篇（7）OGR空间计算
qq_45373844的博客
03-16 1081
GDAL python教程基础篇（7）OGR空间计算
python(ogr)处理geojson为本地shp文件
weixin_73810008的博客
03-09 694
快别用GIS桌面软件处理geojson数据了！！
GDAL python教程基础篇（1）——用OGR读矢量数据
qq_45373844的博客
02-16 1406
利用ogr读取矢量数据
【Python的GDAL/OGR库】
最新发布
qq_36253366的博客
07-30 451
定义：GDAL是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式，并提供了一系列命令行工具来进行数据转换和处理。功能：GDAL主要提供对栅格数据的支持，包括读取、写入、转换和处理等操作。它支持多种栅格数据格式，如GeoTIFF、JPEG2000、HDF5、NetCDF等。类与结构：GDAL使用一系列类来管理和处理栅格数据，如GDALDataset类用于表示一个栅格数据集，GDALDriver类用于管理文件格式驱动等。应用。
wxpython pyqt_分别使用wxPython与PyQt显示Mapnik的渲染结果出现问题
weixin_39637049的博客
11-25 247
Python与开源GIS教程的内容，开发了单独的内容，请打开https://www.osgeo.cn/pygis/查看。本页面的内容不再更新。Mapnik 是一个开源的 Python/C++ 地图渲染引擎，他的功能是把数据形式的地图，包含一些地理对象，如地图、层、数据源、特征和地理几何等，通过一个样式表的定义，渲染成位图格式，用来提供 WMS 服务。其核心是一个 C++的共享库提供空间数据访问和可...
Python+GDAL/OGR矢量数据读写
03-02 9099
       常见的矢量数据格式有Shapefile、GeoJSON、CSV，及文件数据库gdb和空间数据库PostGIS，不论是何种格式的数据或如何存储，一旦打开数据源、获取矢量图层后（详情参考OGR操作矢量数据的类结构图），对数据的操作都一样。下面对矢量数据的读写进行详细的介绍。    一、打开不同的矢量文件    1、定义打开数据源的函数，并遍历所有的图层，输出他们的名字和图层''' ...
Python——使用OGR操作矢量数据
qq_45751615的博客
10-07 930
OGR库
GDAL与OGR：地理空间数据处理的核心库
学习GDAL不仅有助于理解地理空间数据处理的基本原理，还能掌握如何在不同编程语言中高效地操作和转换各种数据格式。对于任何从事GIS或遥感分析的人来说，掌握GDAL和OGR是必不可少的技能，它们是地理信息系统的基石，...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值