Py3.8将存放经纬度的csv文件转为Shp文件（点要素和面要素）

最新推荐文章于 2024-04-11 15:17:28 发布
最新推荐文章于 2024-04-11 15:17:28 发布
阅读量1.6k 收藏 13
点赞数 2
分类专栏： python笔记 文章标签： python pandas 其他
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44814176/article/details/125445208
版权

写论文数据可视化需要用到shp文件，这里做个记录方便日后自己查看。
本文参考博主1、博主2的代码改的（故不多解释，直接上数据和代码）。并附上PythonGIS相关的知识手册。

1.数据路径和格式

在这里插入图片描述

下图中CSV文件里寸的是影像的名称和四角点以及中心点的经纬度坐标
在这里插入图片描述

2.源代码

下面给出了生成点要素和面要素的代码：

# 导入相关库
import os
from osgeo import ogr
import pandas as pd
from osgeo import osr
import glob

# 启动异常报错提示
ogr.UseExceptions()

# .shp文件保存路径
shp_path = r'F:\SHP'
# 输入的csv文件路径
csv_path = r'F:\SHP'

#点要素
def poin_shp():
    for csv_filename in glob.glob(os.path.join(csv_path, '*.csv')):

        # 读入csv文件信息，设置点几何的字段属性
        csv_df = pd.read_csv(csv_filename)

        # 利用.csv文件创建一个点shp文件

        # 获取驱动
        driver = ogr.GetDriverByName('ESRI Shapefile')

        # 创建数据源
        shp_filename = os.path.basename(csv_filename)[:-4] + '_point.shp'
        # 检查数据源是否已存在
        if os.path.exists(os.path.join(shp_path, shp_filename)):
            driver.DeleteDataSource(os.path.join(shp_path, shp_filename))
        ds = driver.CreateDataSource(os.path.join(shp_path, shp_filename))

        # 图层名
        layer_name = os.path.basename(csv_filename)[:-4]

        # 定义坐标系对象
        sr = osr.SpatialReference()
        # 使用WGS84地理坐标系
        sr.ImportFromEPSG(4326)

        # 创建点图层, 并设置坐标系
        out_lyr = ds.CreateLayer(layer_name, srs=sr, geom_type=ogr.wkbPoint)

        # 创建图层定义
        # 利用csv文件中有四个字段创建4个属性字段
        # station字段
        station_fld = ogr.FieldDefn('Image', ogr.OFTString)
        station_fld.SetWidth(6)
        out_lyr.CreateField(station_fld)
        # Latitude字段
        lat_fld = ogr.FieldDefn('latitude', ogr.OFTReal)
        lat_fld.SetWidth(9)
        lat_fld.SetPrecision(5)
        out_lyr.CreateField(lat_fld)
        # Longitude字段
        lon_fld = ogr.FieldDefn('longitude', ogr.OFTReal)
        lon_fld.SetWidth(9)
        lon_fld.SetPrecision(5)
        out_lyr.CreateField(lon_fld)

        # # pr字段
        # pr_fld = ogr.FieldDefn('pr', ogr.OFTReal)
        # pr_fld.SetWidth(5)
        # pr_fld.SetPrecision(2)
        # out_lyr.CreateField(pr_fld)

        # tas字段
        # tas_fld = ogr.FieldDefn('tas', ogr.OFTReal)
        # tas_fld.SetWidth(6)
        # tas_fld.SetPrecision(2)
        # out_lyr.CreateField(tas_fld)
        # 从layer中读取相应的feature类型，并创建feature
        featureDefn = out_lyr.GetLayerDefn()
        feature = ogr.Feature(featureDefn)

        # 设定几何形状
        point = ogr.Geometry(ogr.wkbPoint)

        # 读入csv文件信息，设置点几何的字段属性
        for i in range(len(csv_df)):
            # 设置属性值部分
            # 站点Id
            feature.SetField('Image', str(csv_df.iloc[i, 0]))
            # 纬度
            feature.SetField('latitude', float(csv_df.iloc[i, 1]))
            # 经度
            feature.SetField('longitude', float(csv_df.iloc[i, 2]))
            # # pr值
            # feature.SetField('pr', float(csv_df.iloc[i, 3]))
            # tas值
            # feature.SetField('tas', float(csv_df.iloc[i, 3]))

            # 设置几何信息部分
            # 利用经纬度创建点， X为经度， Y为纬度
            point.AddPoint(float(csv_df.iloc[i, 2]), float(csv_df.iloc[i, 1]))
            feature.SetGeometry(point)

            # 将feature写入layer
            out_lyr.CreateFeature(feature)

        # 从内存中清除 ds，将数据写入磁盘中
        ds.Destroy()

#面要素
def Polygon_shp():
    for csv_filename in glob.glob(os.path.join(csv_path, '*.csv')):

        # 读入csv文件信息，设置点几何的字段属性
        csv_df = pd.read_csv(csv_filename)

        # 利用.csv文件创建一个点shp文件

        # 获取驱动
        driver = ogr.GetDriverByName('ESRI Shapefile')

        # 创建数据源
        shp_filename = os.path.basename(csv_filename)[:-4] + '_Polygon.shp'
        # 检查数据源是否已存在
        if os.path.exists(os.path.join(shp_path, shp_filename)):
            driver.DeleteDataSource(os.path.join(shp_path, shp_filename))
        ds = driver.CreateDataSource(os.path.join(shp_path, shp_filename))

        # 图层名
        layer_name = os.path.basename(csv_filename)[:-4]

        # 定义坐标系对象
        sr = osr.SpatialReference()
        # 使用WGS84地理坐标系
        sr.ImportFromEPSG(4326)

        # 创建面图层, 并设置坐标系
        out_lyr = ds.CreateLayer(layer_name, srs=sr, geom_type=ogr.wkbPolygon)

        # 创建图层定义
        # 利用csv文件中有四个字段创建4个属性字段
        # Image字段
        station_fld = ogr.FieldDefn('Image', ogr.OFTString)
        station_fld.SetWidth(6)
        out_lyr.CreateField(station_fld)
        # Latitude字段
        lat_fld = ogr.FieldDefn('latitude', ogr.OFTReal)
        lat_fld.SetWidth(9)
        lat_fld.SetPrecision(5)
        out_lyr.CreateField(lat_fld)
        # Longitude字段
        lon_fld = ogr.FieldDefn('longitude', ogr.OFTReal)
        lon_fld.SetWidth(9)
        lon_fld.SetPrecision(5)
        out_lyr.CreateField(lon_fld)

        # # pr字段
        # pr_fld = ogr.FieldDefn('pr', ogr.OFTReal)
        # pr_fld.SetWidth(5)
        # pr_fld.SetPrecision(2)
        # out_lyr.CreateField(pr_fld)

        # tas字段
        # tas_fld = ogr.FieldDefn('tas', ogr.OFTReal)
        # tas_fld.SetWidth(6)
        # tas_fld.SetPrecision(2)
        # out_lyr.CreateField(tas_fld)
        # 从layer中读取相应的feature类型，并创建feature
        featureDefn = out_lyr.GetLayerDefn()
        feature = ogr.Feature(featureDefn)

        # 读入csv文件信息，设置点几何的字段属性
        for i in range(len(csv_df)):
            # 设置属性值部分
            # 站点Id
            feature.SetField('Image', str(csv_df.iloc[i, 0]))
            # 纬度
            # feature.SetField('tl_lat', float(csv_df.iloc[i, 3]))
            # feature.SetField('tr_lat', float(csv_df.iloc[i, 5]))
            # feature.SetField('bl_lat', float(csv_df.iloc[i, 7]))
            # feature.SetField('br_lat', float(csv_df.iloc[i, 9]))
            # # 经度
            # feature.SetField('tl_lon', float(csv_df.iloc[i, 4]))
            # feature.SetField('tr_lon', float(csv_df.iloc[i, 6]))
            # feature.SetField('bl_lon', float(csv_df.iloc[i, 8]))
            # feature.SetField('br_lon', float(csv_df.iloc[i, 10]))

            # 设置几何信息部分
            # 利用经纬度创建点， X为经度， Y为纬度
            ring = ogr.Geometry(ogr.wkbLinearRing)
            ring.AddPoint(float(csv_df.iloc[i, 4]), float(csv_df.iloc[i, 3]))  #左上
            ring.AddPoint(float(csv_df.iloc[i, 6]), float(csv_df.iloc[i, 5]))   #右上
            ring.AddPoint(float(csv_df.iloc[i, 10]), float(csv_df.iloc[i, 9]))  # 右下
            ring.AddPoint(float(csv_df.iloc[i, 8]), float(csv_df.iloc[i, 7]))   #左下


            # 设定几何形状
            Polygon = ogr.Geometry(ogr.wkbPolygon)
            Polygon.AddGeometry(ring)

            # ds = ogr.Open(shp_path, True)
            # layer = ds.GetLayer(0)
            # feature = ogr.Feature(layer.GetLayerDefn())
            feature.SetGeometry(Polygon)


            # 将feature写入layer
            out_lyr.CreateFeature(feature)

        # 从内存中清除 ds，将数据写入磁盘中
        ds.Destroy()

if __name__ == "__main__":
     # poin_shp()  #点要素
     Polygon_shp()   #面要素

3.结果展示

转成功后便生成如下文档。
在这里插入图片描述
想要看到图形，用ArcMap打开即可。步骤为：ArcMap工具栏选择“文件”——"添加数据“——选择相应的”.SHP“文件即可。
在这里插入图片描述
这里就只展示图的一部分（点要素）。

在这里插入图片描述

是谁的兴啊
关注
  • 2
    点赞
  • 踩
  • 13
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
python把经纬度生成shp,Python GeoPandas 文本经纬度转换为点要素、线要素
weixin_32165219的博客
03-26 2342
需求目的：将以CSV格式存储的经纬度点坐标文件，转换为shapefile的point类型和line类型的矢量数据，主要依赖geopandas和shapely包。核心思想：采用shapely里面的Point空间数据模型来构建Point对象，并将其作为参数传入构建GeoPandas的主要类型之一GeoDataFrame的geometry属性中。在线要素的创建中，需要采用shapely的LineStri...
CSV文件转Shp文件
追求灵魂代码
03-20 5259
准备工作 环境：Python3.7，需要引入shapefile、csv两个库 开发工具：Pycharm 示例文件：链接：https://pan.baidu.com/s/16xbjqLMaTiVy2hB8FGUh4w 提取码：ch2d 原始教程网站(Pyshp文档)：https://pypi.org/project/pyshp/ 代码实现 import shapefile import csv ...
CSV文件转SHP文件
02-08
不需要安装Arcgis等Gis软件，就能将csv文件转shp文件，比较适合管网管线数据入库等方面。目前借助shape-tools中DOS命令创建shp文件，如有二次开发需要，可将直接采用二进制读写方式，去掉第三方软件，软件进一步精简，通用性更强。压缩包内有测试数据，可按照测试数据格式自行编辑。
如何使用Arcgis将csv文件转化成shp文件
m0_48901857的博客
11-05 1万+
（1）添加数据（csv文件格式）——右键——显示XY数据(选择经纬度坐标)——确定——右键——数据——导出（dbf格式） （2）点击dbf文件（勾选，可以在图上展示出散点）——右键——数据——导出数据（shp文件） 以上为csv文件转化成点要素文件的方法，面要素的还在探索中。。。 ...
【Arcpy】将包含坐标的CSV文件批量转为矢量点
最新发布
withowl的博客
04-11 181
在这里我们用代码的方式实现批量转换。（1）文件-添加数据-添加XY数据。（2）右键导出数据成矢量。
Python使用gdal/ogr包以及pandas包批量将带经纬度信息的站点csv文件转为点shape文件
lidahuilidahui的博客
11-24 3294
近日，接到一个博士师兄的要求，希望我写一个Python程序,将其2006-2099年各年份的带经纬度的站点模拟数据.csv文件转为点shape文件（.shp）。网上有利用arcpy转的，也有利用shapely包转的，但是，作为开源GIS/RS爱好者，gdal/ogr包作为开源GIS官方（osgeo）社区维护的包，很多较为底层的操作都能实现，当然，是最好的选择。不过，由于操作较为底层也相对较为繁琐一...
中国地质公园空间分布矢量数据点文件 .shp
02-14
中国地质公园空间分布矢量数据 .shp
python OSM(.map)数据转shp
qq_20373723的博客
02-18 1945
1.OSM数据下载 下载就不多说了，去官网地址OpenStreetMaphttps://www.openstreetmap.org/需要登陆账号，没有的自己注册下，登录后如下，手动选择想要下载的区域，不能太大，否则无法下载 下载好后的文件是map.osm数据 2.OSM数据转化 这里用gdal来转化，参考链接：Vector Layers — Python GDAL/OGR Cookbook 1.0 documentationhttps://pcjericks.github.io/py-gda..
QGIS将CSV格式数据导出为shp图层的方法
qq_43342037的博客
06-14 1767
QGIS将CSV格式数据导出为shp图层的方法。
利用ArcGIS软件将csv文件转换为shp格式
weixin_44848172的博客
11-21 6029
利用ArcGIS软件将csv文件转换为shp格式
中国世界遗产地空间分布 点文件 .shp
02-14
中国世界遗产地空间分布 点文件 .shp
34个省市直辖市省会驻地点要素 shp文件（最新版）.rar
11-11
客官您好！ 文件中包含了我国最新的34省、直辖市、特别行政区首府驻地点要素文件，也就是矢量数据shp文件，可用GIS软件打开。精心制作，良心数据，安心购买！
森林公园空间分布点文件 .shp
02-14
森林公园空间分布点文件 .shp
如何利用ArcGis把经纬度转成shp数据
热门推荐
老马的程序人生
09-20 1万+
这段时间在处理一批地理信息数据，由于部分数据是经纬度坐标，如下图所示： 这样，面对的第一个问题，就是把这批数据转换成shp格式。下面做一个记录，与大家分享。 工具： ArcGIS 软件 Step1：用Excel把经纬度的坐标换算为小数的形式。 “分”“秒”换算成小数的规则为：“分”位除以60，“秒”位除以3600，比如 37°15′54′′ = 37 + 15/60 + 54/3600 ...
ArcGIS 已知某点经纬度，如何制以该点为中心的10km的矢量shp文件
zhangxiaohua18的博客
03-16 1463
文章目录前言一、将某点经纬度输出到csv文件中二、将csv文件添加到ArcGIS中三、将经纬度位置到处为shp四、对该点建立10km区域缓冲区总结 前言 ArcGIS 已知某点经纬度，制以该点为中心的10km的矢量shp文件 一、将某点经纬度输出到csv文件中 首先将度分秒的格式全部转换为度的格式。 然后将转换后的经纬度输出到excel文件中，然后保存为csv格式的文件。如下图所示： 二、将csv文件添加到ArcGIS中 文件——添加数据——添加XY数据 选择要上述的csv文件。X字段选择经度；Y字
【Python Geopandas】csv 转 shp 文件
Ricce的博客
11-03 1355
csv 转 shp
python3操作shp文件
weixin_43931979的博客
06-12 1021
python3操作shp文件
Qgis软件编写python脚本实现批量将txt经纬度文件转换成shp文件
06-07
好的，这是一个关于QGIS和Python编程的问题。我会尽力回答。 首先，需要安装QGIS软件，并确保已经安装了Python插件。然后，可以使用以下Python脚本实现批量将txt经纬度文件转换成shp文件： ```python import os from qgis.core import * import qgis.utils # 设置QGIS应用程序路径 QgsApplication.setPrefixPath("C:/Program Files/QGIS 3.10", True) qgs = QgsApplication([], False) qgs.initQgis() # 设置经纬度txt文件路径 input_folder = "C:/input_folder" # 设置输出shp文件路径 output_folder = "C:/output_folder" # 遍历输入文件夹中的所有txt文件 for file_name in os.listdir(input_folder): if file_name.endswith(".txt"): # 加载txt文件 uri = "file:///" + os.path.join(input_folder, file_name) + "?delimiter=," layer = QgsVectorLayer(uri, file_name[:-4], "delimitedtext") # 设置坐标系 crs = QgsCoordinateReferenceSystem(4326, QgsCoordinateReferenceSystem.EpsgCrsId) layer.setCrs(crs) # 保存为shp文件 output_file = os.path.join(output_folder, file_name[:-4] + ".shp") QgsVectorFileWriter.writeAsVectorFormat(layer, output_file, "UTF-8", crs, "ESRI Shapefile") print("Converted {} to {}".format(file_name, output_file)) # 退出QGIS应用程序 qgs.exitQgis() ``` 在这个脚本中，首先设置了QGIS应用程序的路径，并初始化了QGIS应用程序。然后，设置了输入文件夹和输出文件夹的路径，并遍历了输入文件夹中的所有txt文件。对于每个txt文件，首先使用QGIS的delimitedtext数据提供程序加载文件，并设置其坐标系为WGS 84。然后，使用QGIS的ESRI Shapefile格式将图层保存为shp文件。最后，退出QGIS应用程序。 希望这个脚本能够帮助到你。如果你有任何问题，请随时问我。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值