今天有些不一样，发这篇文章并不是项目需要。单纯的想到有这个功能没使用Python实现，所以就去研究了一下，第一时间就和大家分享。如何使用Python的osgeo库实现面矢量数据与线矢量数据的互相转换。
一、导入所需库
import os
from osgeo import ogr二、面转线
代码中注释都给了，就不讲解了。入参就是输入面矢量的路径，输出线矢量的路径。有一个小问题就是这段代码只适用于多个单面，即一个面为一个要素的情况。如果多个面合并成一个面的也能转换，不过要素会缺失！
# -*- coding: utf-8 -*-
"""
@Time ： 2023/11/21 11:02
@Auth ： RS迷途小书童
@File ：Vector Face To Line.py
@IDE ：PyCharm
@Purpose：面矢量、线矢量相互转换
"""def face_to_line(int_path, out_path):""":param int_path: 输入面矢量路径:param out_path: 输出线矢量路径:return: None"""driver = ogr.GetDriverByName('ESRI Shapefile')ds = ogr.Open(int_path, 0)layer = ds.GetLayer()# 打开面矢量数据，并获取其图层src_proj = layer.GetSpatialRef()# 获取其源坐标信息if os.path.exists(out_path):driver.DeleteDataSource(out_path)# 如果目标文件已存在，则删除它ds_result = driver.CreateDataSource(out_path)layer_result = ds_result.CreateLayer(out_path, srs=src_proj, geom_type=ogr.wkbLineString)defn_result = layer_result.GetLayerDefn()# 创建目标资源、目标图层、目标要素for feature in layer:# 遍历面矢量中的所有要素geom = feature.GetGeometryRef()# 获取该要素的地理空间范围line_geom = geom.GetGeometryRef(0)# 获取线格式地理空间范围feature_result = ogr.Feature(defn_result)# 创建一个新的要素。要素是Shapefile中的数据实体，它们有几何形状和属性。feature_result.SetGeometry(line_geom)# 将图形赋值到要素上layer_result.CreateFeature(feature_result)# 创建该要素，写入layer_result = Noneif __name__ == '__main__':os.chdir(r'G:\彭俊喜')face_to_line('1.shp', '样本_line.shp')# 参数：输入面矢量，输出线矢量三、线转面
这里的逻辑就和面转线不一样，不能够直接读取要素范围然后写入，只能通过线的范围创建wkt格式的面数据再写入。同样入参为线矢量路径和面矢量路径。
# -*- coding: utf-8 -*-
"""
@Time ： 2023/11/21 11:02
@Auth ： RS迷途小书童
@File ：Vector Face To Line.py
@IDE ：PyCharm
@Purpose：面矢量、线矢量相互转换
"""def line_to_face(int_path, out_path):""":param int_path: 输入线矢量路径:param out_path: 输出面矢量路径:return: None"""driver = ogr.GetDriverByName('ESRI Shapefile')ds = ogr.Open(int_path, 0)layer = ds.GetLayer()# 打开面矢量数据，并获取其图层src_proj = layer.GetSpatialRef()# 获取其源坐标信息if os.path.exists(out_path):driver.DeleteDataSource(out_path)# 如果目标文件已存在，则删除它ds_result = driver.CreateDataSource(out_path)layer_result = ds_result.CreateLayer(out_path, srs=src_proj, geom_type=ogr.wkbPolygon)# 创建一个数据资源，格式为面矢量，坐标系为src_projdefn_result = layer_result.GetLayerDefn()# 创建目标资源、目标图层、目标要素for feature in layer:# 遍历面矢量中的所有要素geom = feature.GetGeometryRef()# 获取该要素的地理空间范围feature_result = ogr.Feature(defn_result)# 创建一个新的要素。要素是Shapefile中的数据实体，它们有几何形状和属性。polygon = ogr.CreateGeometryFromWkt("Polygon(%s)" % str(geom)[str(geom).find("("):])# print("Polygon"+str(geom)[11:])feature_result.SetGeometry(polygon)# 将图形赋值到要素上layer_result.CreateFeature(feature_result)# 创建该要素，写入layer_result = Noneif __name__ == '__main__':os.chdir(r'G:\彭俊喜')line_to_face('样本_line.shp', '123312.shp')# 参数：输入线矢量，输出面矢量
四、总结
网上貌似还没有线转面的教程，我应该算是首发了。当然我说的是使用osgeo库，arcpy除外，毕竟这玩意相当于ARCGIS一样太变态了。
本文章主要是分享个人在学习Python过程中写过的一些代码。有些部分借鉴了前人以及官网的教程，如有侵权请联系作者删除，大家有问题可以随时留言交流，博主会及时回复。
)
