python学习——csv文件转换为shp(Python实现)_csv点数据怎么转成点shp-程序员宅基地

python——csv文件转换为shp

转载自https://blog.csdn.net/qq_23926575/article/details/82316965

学地信的经常会有这样的需求，即将csv格式的文件转换为shapefile格式加载到arcMap中进行显示与分析，arcMap本身提供了这样的功能，即

文件->添加数据->添加XY数据：选择文件路径，指定X字段(longitude)、Y字段(latitude)即可，默认坐标系为WGS84坐标系。

arcMap的这个功能非常方便，但因为我的数据量较大，有接近100万条数据，故通过代码进行格式转换效率更高点。
下面的代码参考如下链接：Using pyshp to convert .csv file to .shp?
需要用到pyshp模块、csv模块、codecs模块，pip install 安装即可。
下面贴上代码：
目前只能转换point类型数据

#-*-coding:utf-8-*-
import shapefile as shp
import csv
import codecs
import os

def trans_point(folder, fn, idlng, idlat, delimiter=','):
    # create a point shapefile
    output_shp = shp.Writer(shp.POINT)
    # for every record there must be a corresponding geometry.
    output_shp.autoBalance = 1
    # create the field names and data type for each.you can omit fields here
    # output_shp.field('id','N') # number    
    output_shp.field('longitude', 'F', 10, 8) # float
    output_shp.field('latitude', 'F', 10, 8) # float
    output_shp.field('locname','C',100) # string, max-length
    # access the CSV file
    with codecs.open(folder + fn, 'rb', 'utf-8') as csvfile:
        reader = csv.reader(csvfile, delimiter=delimiter)
        # skip the header
        next(reader, None)
        #loop through each of the rows and assign the attributes to variables
        for row in reader:
            # idx = row[0]
            locname= row[1]
            lng= float(row[idlng])
            lat = float(row[idlat])
            print lng, lat
            # create the point geometry
            output_shp.point(lng, lat)
            # add attribute data
            output_shp.record(lng, lat, locname)
    output_shp.save(folder + "%s.shp"%fn.split('.')[0]) # save the Shapefile

if __name__ == '__main__':
    folder = 'C:\Users\MaMQ\Desktop' + os.sep
    fn = 'test.csv'
    trans_point(folder, fn, 2, 3)

示例数据

id	locname	x	y
1	Haryana, Haryana, India	76.9806	29.6161
2	Mumbai, Maharashtra, India	72.8258	18.975

结果
这里写图片描述

参考资料
[1] Using pyshp to convert .csv file to .shp?

本文链接：https://blog.csdn.net/knkn123/article/details/85601607

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

如何使用SonarQube Scanner扫描器对代码进行扫描-程序员宅基地

文章浏览阅读3.2k次。2019独角兽企业重金招聘Python工程师标准>>> ..._代码扫描工具对转测前代码进行扫描

SQL常用语法( WITH 语句)-程序员宅基地

文章浏览阅读1.3w次，点赞6次，收藏46次。WITH语句

解决MySQL登录出现错误码10061_2008服务器登陆错误10061-程序员宅基地

文章浏览阅读2.5k次，点赞3次，收藏5次。原因：没有启动MySQL服务解决方案：以管理员身份打开cmd命令行窗口，使用cd命令将目录切换为你MySQL的bin目录下，输入 net start mysql 启动MySQL服务即可解决此问题！！！电脑中搜索服务找到MySQL服务，启动它！管理员身份命令行窗口打开Windows快捷键： Win + X + A如果找不到MySQL服务，需要安装MySQL服务。解决的方法是用管理员的方式打开mysql安装目录bin。执行mysqld.exe -install，安装后就可以启动服务了！..._2008服务器登陆错误10061

STM32移植工程教程包括解决Keil 一些常见的错误等等_s32ds 移植 keil startup error-程序员宅基地

文章浏览阅读2.9k次，点赞6次，收藏55次。然后打开Keil 在左边的能看到即说明添加成功出现如图则说明添加成功编译一下，把多余的头文件删去有6个错误我们一般从第一个错误开始看这里的错误是TIM定时器的结构体没定义，但其实这些在工程里面的库函数这些都已经包含了这时候，找到stm32f10x_conf.h打开，可以看到stm32f10x_tim.h被注释掉了，这就是导致编译错误的原因，只需要删掉注释重新编译一下即可重新编译一下错误就没有了有时候会遇到引脚冲突，解决方法一修改文件里面的宏定义更换为别的引脚解决方法二使用重定义功能_s32ds 移植 keil startup error

networkx学习与使用——（6）图划分与介数计算_边介数计算方法-程序员宅基地

文章浏览阅读6.8k次，点赞7次，收藏33次。networkx学习与使用——（5）图划分与介数计算摘要图划分例子生成介数定义及计算定义networkx计算边介数通过networkx的最短路算法实现使用networkx的内置函数计算结果分析参考摘要图划分按照一定规则将一个连通图划分成几个连通分量，看上去有点像聚类的感觉。从网络的角度，会根据一些重要的节点或边来进行划分，这里介绍划分图的指标——边介数。图划分图划分一般有两种方法，“删边法"和"聚集法”。删边法通过删除某条"重要"的边进行划分。聚集法通过将最"接近"的节点聚集起来构成不同的区域。这里_边介数计算方法

Python3读取Excel表格数据时报错:ImportError: Missing optional dependency ‘xlrd‘. Install xlrd ＞= 1.0.0 for Exc_install xlrd >= 1.0.0 for excel support-程序员宅基地

文章浏览阅读388次。ImportError: Missing optional dependency 'xlrd'. Install xlrd >= 1.0.0 for Excel support Use_install xlrd >= 1.0.0 for excel support

随便推点

系统丢失ucrtbased.dll文件一些软件无法运行问题-程序员宅基地

文章浏览阅读202次。其实很多用户玩单机游戏或者安装软件的时候就出现过这种问题，如果是新手第一时间会认为是软件或游戏出错了，其实并不是这样，其主要原因就是你电脑系统的该dll文件丢失了或没有安装一些系统软件平台所需要的动态链接库，这时你可以下载这个ucrtbased.dll文件(挑选合适的版本文件)把它放入到程序或系统目录中，当我们执行某一个.exe程序时，相应的DLL文件就会被调用，因此将缺失的文件放回到原目录之后就能打开你的软件或游戏了.那么出现ucrtbased.dll丢失要怎么解决？_ucrtbased.dll