ogr_ogr python_python ogr - 腾讯云开发者社区 - 腾讯云

开发者社区

文档建议反馈控制台

最新优惠活动

文章/答案/技术大牛

发布

Shapefile属性操作之增

Show me the code. from osgeo import ogr ogr.UseExceptions() # 首先定义每个省全称到简称的映射字典 names = { '北京':...= ogr.Open('...../data/省级行政区.shp', update=True) layer: ogr.Layer = ds.GetLayer() defs: ogr.FeatureDefn = layer.GetLayerDefn...defn.GetName()} -> {defn.GetType()} -> {defn.GetWidth()}') # 添加一个省简称的字段 field: ogr.FieldDefn = ogr.FieldDefn...注意我们这里新建的字段的类型是字符串类型ogr.OFTString，当然我们还可以新建其他类型的字段，例如整形ogr.OFTInteger，实数形ogr.OFTReal，日期型ogr.OFTDate等。

6622 0

GIS空间分析之Buffer

from pathlib import Path import ogr ogr.UseExceptions() in_ds: ogr.DataSource = ogr.Open('...../data/China.shp') in_lyr: ogr.Layer = in_ds.GetLayer() fname: str = 'Buffer.shp' # 创建输出Buffer文件 driver...: ogr.Driver = ogr.GetDriverByName('ESRI Shapefile') if Path(fname).exists(): driver.DeleteDataSource...(fname) # 新建DataSource，Layer out_ds: ogr.DataSource = driver.CreateDataSource(fname) out_lyr: ogr.Layer...这里我们的空间投影直接从原始的Shapefile中进行读取，输入的空间几何体类型设置为ogr.wkbPolygon。

3.6K4 0

您找到你想要的搜索结果了吗？

是的

没有找到

Shapefile属性操作之查

from osgeo import ogr ogr.UseExceptions() ds: ogr.DataSource = ogr.Open('...../data/Provinces.shp') # 注意Layer的名称不能包含中文 layer: ogr.Layer = ds.GetLayer() # 选择出中学数量大于1万所的省份 query: str...selected: print(feature.GetField('NAME')) # 选择出中学数量最多的省份 # 我尝试使用MAX函数和嵌套的SELECT语句进行实现，但是执行报错，应该是OGR...然后使用ExecuteSQL()方法执行查询，得到的查询结果仍然是一个ogr.Layer图层类，但是图层中的要素属性只包含我们查询语句中指定的字段对于第二个案例，我们本来可以使用一个嵌套的SELECT...from osgeo import ogr ogr.UseExceptions() ds: ogr.DataSource = ogr.Open('..

9552 0

矢量数据空间查询

代码演示 import ogr ogr.UseExceptions() ds_province: ogr.DataSource = ogr.Open('...../data/Provinces.shp') l_province: ogr.Layer = ds_province.GetLayer() # 使用filter()方法找出湖北省 f_hubei: ogr.Feature...= next(filter(lambda f: '湖北' in f.GetField('NAME'), l_province)) ds_city: ogr.DataSource = ogr.Open...代码演示 import ogr ogr.UseExceptions() ds: ogr.DataSource = ogr.Open('...../data/Cities.shp') cities: ogr.Layer = ds.GetLayer() # 使用filter()方法找出武汉市 city: ogr.Feature = next(filter

1.6K3 0

GIS空间分析之Clip

代码展示 from pathlib import Path import ogr ogr.UseExceptions() # 读取被裁剪的数据 in_ds: ogr.DataSource = ogr.Open.../data/County.shp') in_lyr: ogr.Layer = in_ds.GetLayer() # 读取裁剪范围数据 method_ds: ogr.DataSource = ogr.Open...driver: ogr.Driver = ogr.GetDriverByName('ESRI Shapefile') if Path(fname).exists(): driver.DeleteDataSource...(fname) # 新建DataSource，Layer out_ds: ogr.DataSource = driver.CreateDataSource(fname) out_lyr: ogr.Layer...方法总结首先我们使用ogr.Open()函数分别读取被裁剪的数据以及裁剪范围数据；使用GetLayer()方法获取数据的图层。裁剪范围数据在GDAL中被称为Method Layer。

1.1K2 0

Shapefile属性操作之改

思路是我们遍历图层中的每一个Feature要素，然后通过ogr.Feature的SetField()方法更新属性值。代码展示 Talk is cheap....Show me the code. from osgeo import ogr ogr.UseExceptions() # 打开一个Shapefile ds: ogr.DataSource = ogr.Open.../data/省级行政区.shp', update=True) layer: ogr.Layer = ds.GetLayer() # 填充属性值 for feature in layer: name...feature.SetField('NAME', name) # 修改完了记得Set一下 layer.SetFeature(feature) # 关闭数据集 ds = None 方法总结首先，我们使用ogr.Open

9081 0

创建Shapefile文件并写入数据

基本思路使用GDAL创建Shapefile数据的基本步骤如下：使用osgeo.ogr.Driver的CreateDataSource()方法创建osgeo.ogr.DataSource矢量数据集...使用osgeo.ogr.DataSource的CreateLayer()方法创建一个图层使用osgeo.ogr.FieldDefn()定义Shapefile文件的属性字段创建osgeo.ogr.Feature...对象，设置每个属性字段的值，使用Feature对象的SetGeometry()定义几何属性创建Feature对象以后，使用osgeo.ogr.Layer的CreateFeature()添加Feature...) # 添加属性定义 fname = ogr.FieldDefn('Name', ogr.OFTString) fname.SetWidth(24) layer.CreateField(fname) fcx...= ogr.FieldDefn('CenterX', ogr.OFTReal) layer.CreateField(fcx) fcy = ogr.FieldDefn('CenterY', ogr.OFTReal

2.3K3 0

Python | GDAL处理矢量数据

现在OGR部分并不代表任何东西，只是历史延续。...创建和修改字段创建字段定义：ogr.FieldDefn() coord_fld = ogr.FieldDefn('X_short', ogr.OFTReal) 修改字段定义：coord_fld.SetXXXX...(ogr.wkbMultiPoint) faucet = ogr.Geometry(ogr.wkbPoint) faucet.AddPoint(67.5, 16) faucets.AddGeometry...(87.5, 24) garden1 = ogr.Geometry(ogr.wkbPolygon) garden1.AddGeometry(box1) # 第二个环 box2 = ogr.Geometry...) garden2 = ogr.Geometry(ogr.wkbPolygon) garden2.AddGeometry(box2) # 闭合 gardens = ogr.Geometry(ogr.wkbMultiPolygon

9.9K2 0

python与地理空间分析（3）shp文件操作

OGR读取数据的流程（读取中国行政区划shp文件）： #导入库 try: from osgeo import ogr except: import ogr #加载相应数据类型的驱动,相当于初始化一个对象...= ogr.FieldDefn("aqi", ogr.OFTReal) layer.CreateField(field_aqi) layer.CreateField(ogr.FieldDefn("latitude...", ogr.OFTReal)) layer.CreateField(ogr.FieldDefn("longitude", ogr.OFTReal)) #创建要素 for i in AQI.index...(0,100) outring.AddPoint(0,0) inring = ogr.Geometry(ogr.wkbLinearRing)inring = ogr.Geometry(ogr.wkbLinearRing...from osgeo import ogr except: import ogr #加载相应数据类型的驱动,相当于初始化一个对象 driver = ogr.GetDriverByName

14.8K7 1

python使用gdal对shp读取,新建和更新的实例

(oFieldID, 1) # 再创建一个叫FeatureName的字符型属性，字符长度为50 oFieldName =ogr.FieldDefn("FieldName", ogr.OFTString)...第二部：代码分析引入GDAL工具包 import osgeo.ogr as ogr import osgeo.osr as osr 注册驱动，这里是ESRI Shapefile类型，并设置shp...点 #ogr.wkbLineString 线 #ogr.wkbMultiPolygon 面这里的图层名称要与上面注册驱动的shp名称一致 layer = data_source.CreateLayer...(“ceshi”, srs, ogr.wkbLineString) 这里设置要素的属性字段，其中设置了两个字段，分别是Name、data，其中ogr.OFTString表示字符串类型，其长度都是14...field_name = ogr.FieldDefn("Name", ogr.OFTString) field_name.SetWidth(14) layer.CreateField(field_name

3.9K3 0

Shapefile属性操作之删

Show me the code. from osgeo import ogr ogr.UseExceptions() # 从给定图层中读取字段的定义，根据给定字段名称找到该字段的索引编号 def...get_field_index_by_name(layer: ogr.Layer, name: str): defs: ogr.FeatureDefn = layer.GetLayerDefn(...).GetName(): return i raise ValueError(f'{name} not found') # 打开一个Shapefile文件 ds: ogr.DataSource...= ogr.Open('...../data/省级行政区.shp', update=True) layer: ogr.Layer = ds.GetLayer() # 删除Abbr字段 index = get_field_index_by_name

7471 0

MongoDB导入Shapefile数据

两种解决方案：一、将整个shapefile转为GeoJSON然后直接导入mongoDB数据库中首先，将shapefile数据转为WGS84地理坐标，然后使用GDAL的命令行工具ogr2ogr...进行格式的转换，转换命令如下： ogr2ogr -f geoJSON continents.json continents.shp 删除生成JSON文件的前两行{ "type": "FeatureCollection

1.9K2 0

矢量数据投影转换

方法介绍跟栅格数据投影转换一样，使用GDAL库，我们有两种方法进行矢量数据的重投影：使用命令工具及其对应的命令行API接口进行转换（简单，准确，实践中一定要用这种方法） GDAL提供了ogr2ogr...命令行工具进行矢量数据投影转换，命令如下：ogr2ogr -t_srs "+proj=aea +lat_1=25 +lat_2=47 +lat_0=30 +lon_0=105 +x_0=0 +y_0=0...units=m +no_defs """ gdal.VectorTranslate(dst_file, src_file, dstSRS=srs_def, reproject=True) src_ds = ogr.Open...src_layer = src_ds.GetLayer(0) src_srs = src_layer.GetSpatialRef() # 输入数据投影调用基本API函数实现 from osgeo import ogr...src_feature: geometry = src_feature.GetGeometryRef() geometry.Transform(ctx) dst_feature = ogr.Feature

1.8K3 0

生成面要素的中心点-GDAL

PointOnSurface(Geometry self) -> Geometry OGRGeometryH OGR_G_PointOnSurface(OGRGeometryH hGeom) Returns...If OGR is built without the GEOS library, this method will always fail, issuing a CPLE_NotSupported error...OGR 1.10 Centroid(Geometry self) -> Geometry int OGR_G_Centroid(OGRGeometryH hGeom, OGRGeometryH hCentroidPoint...If OGR is built without the GEOS library, this function will always fail, issuing a CPLE_NotSupported...import ogr p = ogr.CreateGeometryFromWkt('POLYGON((1 1,5 1,5 5,1 5,1 1),(2 2,2 4,4 4,4 2,2 2)) ') #

2.5K2 0

在Python中用GDAL实现矢量对栅格的切割实例

8 -*- """ @author lzugis @date 2017-06-02 @brief 利用shp裁剪影像 """ from osgeo import gdal, gdalnumeric, ogr...world file) info srcImage = gdal.Open(raster_path) geoTrans = srcImage.GetGeoTransform() # Create an OGR...layer from a boundary shapefile shapef = ogr.Open(shapefile_path) lyr = shapef.GetLayer( os.path.split...raster_path = 'world.tif' main( shapefile_path, raster_path ) 补充知识：Python+GDAL | 读取矢量并写出txt 这篇文章主要描述了如何使用GDAL/OGR...代码 import ogr,sys,os import numpy as np os.chdir(r'E:\') #设置driver,并打开矢量文件 driver = ogr.GetDriverByName

1.1K1 0

GDAL的java环境配置以及将shp转换为json

实现的关键代码如下： package com.lzugis.gdal; import org.gdal.ogr.*; import org.gdal.ogr.Driver; import org.gdal.gdal....*; public class GdalShpTest { public static void main(String[] args) { // 注册所有的驱动 ogr.RegisterAll...String strVectorFile = "D:\\data\\china\\capital.shp"; //打开数据 DataSource ds = ogr.Open...); Driver dv = ogr.GetDriverByName("GeoJSON"); if (dv == null) {

2.4K3 0

GIS项目中的开源架构分析

4 GDAL/OGR功能特点· 多格式支持：GDAL/OGR支持多种栅格和矢量地理空间数据格式，包括ESRI Shapefiles、S-57、SDTS、PostGIS、Oracle Spatial等。...· 多语言支持：GDAL/OGR提供了多种编程语言的绑定，如C、C++、Perl、Python、Java、C#等。· 高性能：GDAL/OGR使用面向对象的C++语言编写，具有很高的执行效率。...应用场景· 遥感数据处理：GDAL/OGR适用于遥感数据的处理和分析，支持多种卫星数据和地理信息系统数据格式。...· 地理空间数据转换：GDAL/OGR提供了强大的数据转换功能，适用于需要跨平台数据处理的场景。...QGIS、GeoServer、PostgreSQL with PostGIS、GDAL/OGR、MapServer和CesiumJS等开源项目在各自的应用领域中表现出色，广泛应用于政府、企业、科研等多个领域

1141 0

基于8211lib库对s57电子海图的解析和存储

ISO 8211lib是C++编写的开源库，专门用于读取符合ISO 8211封装标准的数据文件，如S-57电子海图数据格式和空间数据传输标准数据集（SDTS）；GDAL是一个开源栅格空间数据转换库，其中OGR...OGR中关于S-57读取的官方文档在http://www.gdal.org/drv_s57.html可以看到，其中需要注意OGR的S-57驱动模块将处理S-57文件里所有的feature要素，其中S-57...在OGR中，从S-57文件读取的所有要素根据物标类型（OBJL）归属于相应的图层，可以类比于GIS里的shape格式，对应的图层下有要素（feature）,要素下面有字段（field）和几何对象（geometry...directory with the s57objectclasses.csv and related files,具体可以参见 http://osgeo-org.1560.x6.nabble.com/OGR-to-access-s57

3.7K4 0

如何在Ubuntu 14.04上安装和配置PostGIS

当我们安装GDAL时，我们还将获得OGR（OpenGIS简单功能参考实现）和命令ogr2ogr。这是一个矢量数据转换库，我们将把它用于将Shapefile转换为PostGIS可以使用的数据。...使用包管理器安装GDAL： sudo apt-get install gdal-bin 现在再次切换到用户postgres： sudo -i -u postgres 现在将使用ogr2ogr把Natural...Earth的Shapefile转换为PostGIS表，如下所示： ogr2ogr -f PostgreSQL PG:dbname=test1 -progress -nlt PROMOTE_TO_MULTI...ogr2ogr命令将基于文件中的前几个特征对几何类型进行假设。我们导入的数据包含多边形类型和多部分多边形或MultiPolygons的混合。...最后，我们指定输入文件的路径： /home/sammy/nedata/ne_110m_admin_0_countries.shp 访问ogr2ogr网站查看完整的选项。

1.6K0 0

python与地理空间分析（一）

可供选择的矢量格式数目很饿的，开源矢量库OGR支持超过86种矢量格式，FME支持超过188种矢量格式。...库，python中解析Shapefile文件的模块大都是基于OGR库的。...库 OGR库是python处理矢量数据的通用库，也是许多处理矢量数据包的依赖项。...from osgeo import ogr shp=ogr.Open("point.shp") layer=shp.GetLayer() for feature in layer: geometry...哭的数据访问功能提供了一套简洁的python api，大大提高了OGR库在python中的易用性。

8K5 2

点击加载更多

扫码

添加站长进交流群

领取专属 10元无门槛券

手把手带您无忧上云

扫码加入开发者社群

相关资讯

热门标签

活动推荐

运营活动

活动名称

广告关闭