Python读写矢量数据（2）矢量数据写入（属性数据）——Python地理数据处理学习分享

已于 2022-10-19 11:05:12 修改
阅读量1.8k 收藏 13
点赞数 2
分类专栏： Python地理数据处理学习分享 文章标签： python
于 2022-10-19 00:06:38 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/remote_giser/article/details/127377213
版权

        这一节主要介绍矢量数据的写入（只有属性数据，无几何），如果有读者没有读取的基础建议先看一下上一篇文章，需要对矢量数据读取有一定的了解才能继续学习本节。在这里我们用到的数据仍为goble文件夹下的数据，数据获取可参考下方链接中的附录获取。Python读写矢量数据（1）针对读取矢量数据——Python地理数据处理学习分享_Cishazhu_ooook的博客-CSDN博客https://blog.csdn.net/remote_giser/article/details/127351424?spm=1001.2014.3001.5502

1.矢量数据的复制

        当我们想新建矢量文件（shapefile）比复制以前的矢量数据时，可以采用以下方法。

# 导入ogr模块
from osgeo import ogr
# 设置路径为global文件夹，该目录下存在很多shapefile文件，并打开文件夹
path="D:\\Python_study\\osgeopy-data\\global"
doc=ogr.Open(path,1)
# 打开图层
in_layer=doc.GetLayer("ne_50m_populated_places")
# 查询图层，存在删除图层
if doc.GetLayer('capital_cities'):
	doc.DeleteLayer('capital_cities')
# 创建一个新图层，空间参考信息与in_layer一致，并继承其属性字段
our_layer=doc.CreateLayer('capital_cities',in_layer.GetSpatialRef(),ogr.wkbPoint)
our_layer.CreateFields(in_layer.schema)
# 创建一个用于存储属性和几何的虚拟存储要素，插入图层中
our_defn=our_layer.GetLayerDefn()
our_feat=ogr.Feature(our_defn)
# 给要素赋予属性
for in_feat in in_layer:
	if in_feat.GetField('FEATURECLA')=='Admin-0 capital':
		geom=in_feat.geometry()
		our_feat.SetGeometry(geom)
		for i in range(in_feat.GetFieldCount()):
			value=in_feat.GetField(i)
			our_feat.SetField(i,value)
		our_layer.CreateFeature(our_feat)
del doc

        这个程序打开方式就与众不同，首先我们打开的不是某个shapefile文件，而是将路径定义在某个文件夹下（这个文件夹包含大量shapefile文件），并用编辑的方式打开（Open(x,1/0)设为1为编辑模式，0为只读）。打开这个文件夹后，这个文件夹中的所有shapefile文件都可被看作是一个单独的图层。因此，在我们可以创建一个新图层，就会为我们默认创建一个shapefile文件。接着我们使用GetLayer（）函数获取了"ne_50m_populated_places"文件当作图层in_layer。

        在创建新图层capital_cities（shapefile文件）之前，我们要确认文件夹中有无这个图层，如果有将其删除，重新创建，这一步的目的是后面我对程序多次运行对图层进行修改（没有这一步，多次运行将会持续报错有同名文件）。用if语句来判断文件夹时候能获取到（GetLayer（））这个图层，如果存在将用DeleteLayer()函数将其删除。

        随后我们创建了我们的图层our_layer。在这里我们使用的是CreateLayer()函数在文件夹下创建了新的图层our_layer。这里需要注意CreateLayer()函数有四个参数。CreateLayer(name,[srs],[geo_type],[options])如下：

  1. name，代表新建的图层名称。在这里我们使用的是'capital_cities'。
  2. srs,代表图层使用的空间参考，默认为空（如不输入将表示没有任何空间参考）。这里继承了in_layer的空间参考，采用的是GetSpatialRef()函数（上节介绍过）。
  3. geom_type代表一个几何常量，意思是什么类似的几何要素。这里我们使用了点要素ogr.wkbPoint，上节也介绍过了（wkb你可以理解为给电脑看的格式，wkt是给人看的直白的话就是点point\线line\面polygon，在wkb中点对应wkbPoint、线对应wkbLineString、面对应wkbPolygon）想要更多了解可看上节的对应表。
  4. options，一个图层创建的的选项列表，只适用于特定的矢量数据类型。这里我未作设置，只有前三个参数（个人没搞懂这个参数，大概不进行输入的以options=value形式存储字符串）。

        创建图层完成后，我们接下来要给图层创建属性，利用CreateFields（）的方法继承in_layer的属性（schema属性会获得图层FieldDefn对象列表）。到这里我们创建图层属性的任务基本完成。如果这时候运行程序会发现会获得继承文件的属性，可在GIS软件中查看。

eef4d921e30c4959854a59381e310566.png

        从代码的14行开始，我们将进行继承每一个要素对应的属性（几何和文本）添加到文件中。首先我们要创建一个用于存储属性和几何的虚拟要素（变量our_defn），将图层对象属性信息传递给这个虚拟要素，这里采用GetLayerDefn()的方法来传递，让out_defn接受in_layer的图层对象属性。 

        随后我们需要使用ogr.Feature把这个虚拟要素包含的属性真正给一个变量，在这里我们将虚拟要素our_defn的特征给了our_feat（这里称变量吧）。

        接下来我们将进行要素对应的属性赋值，这也是最繁琐的一步。在被继承的图层in_layer中遍历每个要素（for循环中的feat）。如果某个要素的'FEATURECLA'等于'Admin-0 capital'，建立一个继承要素的遍量geom，使用geometry()方法继承in_feat要素（循环中指每个要素）的几何。再使用SetGeometry（）的方法让我们变量our_feat继承原来图层要素的几何。

        几何继承完后我们就可以继承属性值，当然两者继承没有先后（只是我觉着有了几个位置再向内部插值比较直观）。就一次循环过程中而言（最外层的for），使用for循环遍历要素in_feat的属性数量，让变量value获取in_field的属性(有for循环，采用的函数为GetField())，最后用变量our_feat使用SetField（）逐一获取单一value，直到读完整个要素的value。这样我们的特征都被存入再our_feat中了，接下来使用CreateFeatyre()让我们图层全部继承其几何和属性。

2.创建新的数据源

2.1创建正确的驱动

        在创建新的数据源之前我们要了解到，新的数据源创建需要正确驱动。例如，GeoJSON驱动程序并不能创建一个shapefile文件，即使创建后扩展名为.shp，其核心还是一个GeoJSON文件。

        下面有两者方法来创建所需的驱动。

        1.继承从打开的数据中进行获取驱动（该方法创建的数据源将会和打开的数据源一致）。

# 导入模块，打开文件
from osgeo import ogr
doc=ogr.Open('D:\\Python_study\\osgeopy-data\\global\\ne_50m_populated_places.shp')
# 获取驱动
driver=doc.GetDriver()

        在这里使用.GetDriver（）的方式继承了文件doc的驱动。

        2.使用OGR的GetDriverByName（）函数，传递给他驱动名字。

# 导入模块
from osgeo import ogr
# 设置驱动
new_driver=ogr.GetDriverByName("ESRI Shapefile")

        这里我们使用是驱动名为ESRI Shapefile,其他驱动的访问可以参考下代码（该代码与创建驱动无关，但我们可能需要特定驱动才能正确创建数据源）：

# 导入模块
from osgeo import ogr
# 查看驱动数量
num=ogr.GetDriverCount()
print(num)
# 查看驱动
for i in range(int(num)):
	print(ogr.GetDriver(i).name)

cdb0599ceed64b9f8246fb34ed61155e.png

        从图可以看到这里可以获取我们常用的格式 。如果不清楚对应文件可取OGR网站查询每种格式对应的文件。

2.2创建新的数据源

        只有在正确的驱动下我们才可以进行创建正确的数据源。创建数据源代码如下：

# 导入模块
from osgeo import ogr
# 设置驱动
new_driver=ogr.GetDriverByName("ESRI Shapefile")
nwe_doc=new_driver.CreateDataSource('D:\\Python_study\\osgeopy-data\\global\\new_shp')

 40a1362c3f414b49a22a977010127b0d.png

         从图片可以看到，我们在文件夹目录下创建了一个new_shp文件（数据源是个文件夹）。注意，创建新的数据源时，不能覆盖现有的数据源，那么我们就需要在创建新数据源之前删除旧的数据源，目前采用的方法是Python的os.path.exists()函数处理。

# 导入模块
from osgeo import ogr
import os
# 设置驱动
new_driver=ogr.GetDriverByName("ESRI Shapefile")
# 检查是否已经存在
if os.path.exists("D:\\Python_study\\osgeopy-data\\global\\new_shp"):
	new_driver.DeleteDataSource("D:\\Python_study\\osgeopy-data\\global\\new_shp")
# 创建新数据源
nwe_doc=new_driver.CreateDataSource("D:\\Python_study\\osgeopy-data\\global\\new_shp")

2.3新建属性字段

        OGR定义字段与数据库类型时，首先需查看有哪些字段类型。这里统计了一些常见的字段数据类型在OGR中的对应。

枚举值

类型

位数

编码值

OFTInteger

整型

32位整型

0

OFTIntegerList

整型列表

32位整型

1

OFTReal

实数型

双精度浮点型

2

OFTRealList

实数型列表

双精度浮点型

3

OFTString

字符串

ASCII字符

4

OFTStringList

字符串列表

ASCII字符串数组

5

OFTDate

日期型

9

OFTTime

时间型

10

OFTDateTime

日期时间型

11

        要将一个字段添加到图层中，首先需要一个包含字段名称、数据类型、字段宽度、精度等重要信息的对象（使用ogr.FieldDefn()函数）。然后基于这个对象将其添加到图层。

# 导入模块
from osgeo import ogr
import os
# 设置驱动
new_driver=ogr.GetDriverByName("ESRI Shapefile")
# 检查是否已经存在
if os.path.exists("D:\\Python_study\\osgeopy-data\\global\\new_shp"):
	new_driver.DeleteDataSource("D:\\Python_study\\osgeopy-data\\global\\new_shp")
# 创建新数据源
new_doc=new_driver.CreateDataSource("D:\\Python_study\\osgeopy-data\\global\\new_shp")
# 创建新图层
layer=new_doc.CreateLayer("new_layer")
# 定义字段
add_1=ogr.FieldDefn("add_the_1",ogr.OFTInteger)
add_1.SetWidth(8)
add_1.SetPrecision(3)
# 将定义字段添加到图层
layer.CreateField(add_1)

        别忘了，我们在这里创建完数据源后记得添加新的图层（即shape file文件，在这里我忘了添加坐标系和类型。），然后用变量add_1指向新建的字段并输入字段名称和类型（在对应表中寻找）。随后我们设置了新字段的宽度和精度（分别用到了SetWidth()、SetPrecision()）函数。最后使用CreatField（）添加到了图层中，并在GIS软件中可以到如下结果。

76e291329a65481b8d733aba2603ab83.png

         当然有时我们可能想要创建仅仅是字段名不同的两个字段，这里会提供SetName()函数来改变字段名而其他不发生改变，但需要重新将第二个字段添加到属性表中（实际修改了变量然后增进行了第二次增加）。

# 导入模块
from osgeo import ogr
import os
# 设置驱动
new_driver=ogr.GetDriverByName("ESRI Shapefile")
# 检查是否已经存在
if os.path.exists("D:\\Python_study\\osgeopy-data\\global\\new_shp"):
	new_driver.DeleteDataSource("D:\\Python_study\\osgeopy-data\\global\\new_shp")
# 创建新数据源
new_doc=new_driver.CreateDataSource("D:\\Python_study\\osgeopy-data\\global\\new_shp")
# 创建新图层
layer=new_doc.CreateLayer("new_layer")
# 定义字段
add_1=ogr.FieldDefn("add_the_1",ogr.OFTInteger)
add_1.SetWidth(8)
add_1.SetPrecision(3)
# 将定义字段添加到图层
layer.CreateField(add_1)
# 仅名称不一样的字段
add_1.SetName("add_the_2")
layer.CreateField(add_1)

        查看结果如下：

375da677458441c9b966baf17c55e0bf.png

2.4修改字段信息 

        对于实例中的add_the_1名称不满意，在这里我们可以对字段就行修改（包括精度、类型等）。采用AlterFieldDefn(iField,field_def,nFlags)进行修改。

  1. iField，索引的字段名，本例去索引的为add_the_1
  2. field_def，新属性字段定义
  3. nFlags，改变的属性，下表代表属性对应的代码

         dcc5bbb7945c46deb63364bbfc05d73a.png

# 导入模块
from osgeo import ogr
import os
# 设置驱动
new_driver=ogr.GetDriverByName("ESRI Shapefile")
# 检查是否已经存在
if os.path.exists("D:\\Python_study\\osgeopy-data\\global\\new_shp"):
	new_driver.DeleteDataSource("D:\\Python_study\\osgeopy-data\\global\\new_shp")
# 创建新数据源
new_doc=new_driver.CreateDataSource("D:\\Python_study\\osgeopy-data\\global\\new_shp")
# 创建新图层
layer=new_doc.CreateLayer("new_layer")
# 定义字段
add_1=ogr.FieldDefn("add_the_1",ogr.OFTInteger)
add_1.SetWidth(8)
add_1.SetPrecision(3)
# 将定义字段添加到图层
layer.CreateField(add_1)
# 修改字段信息
i=layer.GetLayerDefn().GetFieldIndex('add_the_1')
fld_defn=ogr.FieldDefn("change_1",ogr.OFTString)
fld_defn.SetWidth(10)
fld_defn.SetPrecision(6)
layer.AlterFieldDefn(i,fld_defn,ogr.ALTER_ALL_FLAG)

        其实这串代码也十分好解释，首先让 i 获取到要修改的字段（寻找字段名），重新定义一个字段信息并设置精度和宽度。最后使用上述提过的函数来完成字段的改变，当然我们也可以针对单一进行改变。

10e0ede815474bf185c6aff2fba668a1.png

 要素添加、更新和删除

        经过上方的学习，我们已经能修改数据源、图层、属性表的字段了，但我们现在还需要对属性表的要素进行更改。（下方代码用到的数据是50m人口）

# 导入模块
from osgeo import ogr
# 读取文件路径
fn="D:\\Python_study\\osgeopy-data\\global\\ne_50m_populated_places.shp"
# 打开文件
ds=ogr.Open(fn,1)
# 打开图层
lyr=ds.GetLayer(0)
# 添加新字段
add_1=ogr.FieldDefn("ID_add",ogr.OFTInteger)
add_1.SetWidth(8)
add_1.SetPrecision(3)
lyr.CreateField(add_1)
# 要素添加
n=1
for feat in lyr:
	feat.SetField("ID_add",n)
	lyr.SetFeature(feat)
	n+=1
	if n==10:
		break
# 关闭文件
del ds

33707567f88844c7b17d2ee924fa931f.png

         首先我们要找到我们需要添加的字段名称，然后使用for循环来改变要素。SetField("字段名“，更改内容)是用来确定某一要素对应的字段和要改变的内容。SetFeature()是对这个要素对应字段进行改变。

改变单一要素对应值：

# 导入模块
import sys
from osgeo import ogr
# 读取文件路径
fn="D:\\Python_study\\osgeopy-data\\global\\ne_50m_populated_places.shp"
# 打开文件
ds=ogr.Open(fn,1)
if ds is None:
	sys.exit('Could not open{0}.'.format(fn))
# 打开图层
lyr=ds.GetLayer(0)
# 改变单一要素字段
for feat in lyr:
	if feat.GetField('ID_add')==4:
		feat.SetField("ID_add",99)
		lyr.SetFeature(feat)
# 关闭文件
del ds

7df066cbec6248468397f53242b106c7.png

        从上面我们可以知道，要获取某要素时我们往往通过for来进行遍历，然后让其等于某字段的内容。（常采用FID或者名称）这是有一种常用的矢量写入的方法。

        对于要素的删除，则更简单。我们只需要知道要素某个字段的值，便可删除。 这里使用DeleteFeature()的方法将循环参数对应FID的要素删除。（FID的获取使用GetFID()）

# 导入模块
import sys
from osgeo import ogr
# 读取文件路径
fn="D:\\Python_study\\osgeopy-data\\global\\ne_50m_populated_places.shp"
# 打开文件
ds=ogr.Open(fn,1)
if ds is None:
	sys.exit('Could not open{0}.'.format(fn))
# 打开图层
lyr=ds.GetLayer(0)
# 要素删除
for feat in lyr:
	if feat.GetField('ID_add')==5:
		lyr.DeleteFeature(feat.GetFID())
# 关闭文件
del ds

a53cb2935d434c13a879459a5f8f0f20.png

Cishazhu_ooook
关注
  • 2
    点赞
  • 踩
  • 13
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
专栏目录
读取指定_Python GDAL 地学分析（三） 读取矢量数据
weixin_39862871的博客
01-01 489
矢量数据读取OGR提供对矢量数据格式的读写支持，它所支持的文件格式包括：ESRI Shapefiles，PostGIS，Oracle Spatial， Mapinfo mid/mif ， Mapinfo TAB等。OGR包括如下几部分：Geometry：类Geometry (包括OGRGeometry等类)封装了OpenGIS的矢量数据模型，并提供了一些几何操作，WKB(Well Knows Bi...
Caltech 256数据集
不务正业的猿
10-08 3119
Caltech 256数据集是加利福尼亚理工学院收集整理的数据集，该数据集选自google image数据集，并手工去除了不符合其类别的图片。在该数据集中，图片被分为256类，每个类别的图片超过80张。 Caltech 256数据集是 caltech-101 数据集的改进版图片数据集，该数据集收集了 256 个类的 20607 张图片。有几处改进： a）类别数量增加一倍以上 b）任何类别中图像的最小数量从 31 增加到 80 c）避免因图像旋转造成的伪影 d）引入了一个新的更大的杂波类别来测试背景
矢量数据文件解析（shapefile、geojson、kml）
qq_19800029的博客
04-11 492
1、org.geotools.data.geojson.GeoJSONReader解析给geojson,但是它不会读取坐标系，所以不管什么坐标系的geojson文件读出来都是默认的坐标系WGS84.KML(Keyhole Markup Language,Keyhole 标记语言)是一种基于XML 的标记语言，利用XML 语法格式描述地理空间数据(如点、线、面、多边形和模型等)。2、org.geotools.geojson.feature.FeatureJSON解析geojson可以读到文件对应的坐标系。
python 矢量数据的读取
m0_56180742的博客
08-19 887
import os import sys from osgeo import ogr from osgeo import gdal driver=ogr.GetDriverByName('ESRI Shapefile') filename='C:\\Users\\admin\\Desktop\\空间信息处理与系统开发资料\\第六周\\Lesson6_osgeopy-data\\global\\ne_50m_populated_places.shp' dataSource=driver.Open(filen.
Python读写矢量数据（1）针对读取矢量数据——Python地理数据处理学习分享
remote_giser的博客
10-16 2307
分享Python地理数据处理学习内容，仅有矢量读取部分，输入部分在下篇。大家加油！！！
【Python&GIS】Python处理矢量数据的基本操作（查询、修改、删除、新建）
Remote Sensing迷途小书童
08-10 1334
本文章分享了如何使用Python对矢量数据进行查询、修改、删除、新建等操作
Python地学分析 — 矢量数据集介绍
一叶知秋的博客
10-31 2336
欢迎关注博主的微信公众号：“智能遥感”。 该公众号将为您奉上Python地学分析、爬虫、数据分析、Web开发、机器学习、深度学习等热门源代码。 什么叫矢量数据，在这里，小编先给大家普及一下！ 地图上的要素都绘制为几何对象——点，线(包含多边形)、面。这些类型的数据（地理要素都是不同的对象）称为矢量数据集。 在最基本的矢量数据中，地理要素以离散几何形式表示，特别是点，线、面。具有不同边界...
GDAL 矢量属性数据修改方式(python)
12-20
Case：需要给一个现有的shp数据创建一个字段，并将属性表中原有的一个文本类型的属性转换为整型后填入新创建的字段。 Problem：新字段创建成功，但是赋值操作无效，即无法成功给字段写入值。 solution：对字段进行...
基于Python的地理信息数据处理.rar
最新发布
05-07
2. **GDAL/OGR库**：GDAL（Geospatial Data Abstraction Library）和OGR（Open Geospatial Library）是Python中处理地理信息数据的核心库。GDAL主要用于栅格数据，OGR用于矢量数据。它们提供了读写多种地理空间文件...
使用Python自动化管理和分析空间数据
02-19
2. 属性数据分析：利用pandas的统计函数对地理数据的属性进行计算，如平均值、中位数等。 3. 空间关系判断：如点在多边形内、线交叉等，这些可通过Geopandas和Shapely联合完成。 四、空间数据可视化 1. Matplotlib...
Fiona：Fiona读写地理数据文件
01-30
Fiona 是一个基于 Python 的库，用于处理地理空间数据，它提供了一种简单而高效的方式来读取和写入多种矢量数据格式。Fiona 借助于 GDAL/OGR 库，这是一个强大的地理空间数据处理框架，广泛应用于地理信息系统（GIS...
Shp矢量数据批量处理
10-16
批量处理矢量数据属性表，新增字段，计算几何面积；掩膜提取；重分类等，欢迎一起交流学习，yang811901@163.com
python 读取shp文件
03-14
Reading and writing vector data with OGR
Python 遥感数据 地理空间数据 分析常用库 适合Python3.9
09-25
Python在遥感和地理空间数据分析领域扮演着重要角色，提供了许多强大的库来处理、分析以及可视化这类数据。这里我们将深入探讨标题和描述中提到的一些关键库，并了解它们如何助力于Python 3.9版本的遥感和地理空间...
python读取矢量图形
weixin_42577243的博客
02-13 706
要读取矢量图形，你可以使用一些 Python 库，例如： Matplotlib：它是一个非常强大的可视化工具，可以读取并绘制多种格式的图形，包括 SVG。 Pycairo：它是 Cairo 图形库的 Python 封装，可以读取并绘制 SVG 和其他格式的图形。 svg.path：它是一个纯 Python 库，可以读取和解析 SVG 路径字符串。 如果你需要读取特定的图形格式，可以查找与...
Python的标量与矢量运算分析
昊大侠的博客
05-02 2758
文章目录技术NumPyPytorchTensorFlowNumbaCython标量运算矢量运算总结 技术 Python 的官方解释器是用 C 实现的 CPython 。其他解释器有用 Java 实现的 Jython、C# 实现的 IronPython 和 Python 自身实现的 PyPy 。用 C 语言编写的 CPython 利于使用 C 语言接口的外部库。一定程度可以解决 Python 执行慢的问题，但是用 C 编写必要的代码与 Python 通信相当麻烦 NumPy NumPy 是 Python 中科
python实现最短路径、线要素矢量地图自动生成图
Fan-Yan的博客
11-16 628
只有线要素矢量地图，想求任意两点之间的最短路径，arcgis输出的路径只有长度信息，没有详细的线路信息，于是利用python地理处理包实现最短路径。核心步骤就是两步，首先生成图，再基于图实现最短路径。利用momepy和osmnx生成带编号的MultiGraph。生成的最短路径是以图的结点为起终点，要精确到任意输入点的位置，需要重构路网图。
[Python] GeoPandas操作矢量数据（shp、GeoJSON）
热门推荐
geodoer
12-30 2万+
字段类型常量 字段数据类型 OGR常量 Integer OFTInteger List of integers OFTIntegerList Floating point number OFTReal List of floating point numbers OFTRealList String OFTString List of strings OFTSt...
chatgpt赋能python：PythonList:如何将数存入列表
a1sssssb的博客
06-06 1505
本文由chatgpt生成，文章没有在chatgpt生成的基础上进行任何的修改。以上只是chatgpt能力的冰山一角。作为通用的Aigc大模型，只是展现它原本的实力。对于颠覆工作方式的ChatGPT，应该选择拥抱而不是抗拒，未来属于“会用”AI的人。🧡AI职场汇报智能办公文案写作效率提升教程 🧡专注于AI+职场+办公方向。下图是课程的整体大纲下图是AI职场汇报智能办公文案写作效率提升教程中用到的ai工具。
Python地理数据处理 四：矢量数据读写
07-28
你好！对于Python地理数据处理中的矢量数据读写，可以使用一些常见的库来实现，比如`geopandas`、`fiona`和`shapely`等。这些库提供了丰富的功能，可以读取和写入各种矢量数据格式，如Shapefile、GeoJSON等。 首先，你需要安装相应的库。可以使用pip命令来安装它们： ``` pip install geopandas fiona shapely ``` 安装完成后，你可以使用`geopandas`库来读取和写入矢量数据。下面是一个简单的示例： ```python import geopandas as gpd # 读取Shapefile文件 data = gpd.read_file('path/to/shapefile.shp') # 查看数据的前几行 print(data.head()) # 写入GeoJSON文件 data.to_file('path/to/output.geojson', driver='GeoJSON') ``` 在上面的示例中，我们使用`gpd.read_file()`函数来读取Shapefile文件，并使用`print()`函数来查看数据的前几行。然后，我们使用`data.to_file()`函数将数据写入GeoJSON文件，指定`driver='GeoJSON'`来指定输出格式为GeoJSON。 除了Shapefile和GeoJSON格式，`geopandas`还支持其他常见的矢量数据格式，如CSV、KML等。你可以根据具体需求选择合适的格式进行读写操作。 希望这个简单示例对你有所帮助！如果你有更多的问题，请随时提问。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Cishazhu_ooook

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值