Python利用经纬度创建shpfile点图层并生成tif

最新推荐文章于 2024-07-20 03:20:08 发布

胖大狗

最新推荐文章于 2024-07-20 03:20:08 发布

阅读量3.8k

点赞数 3

分类专栏： Python 文章标签： GIS_Python

本文链接：https://blog.csdn.net/qq_39054177/article/details/100919890

版权

本文介绍如何使用Python的shapefile库，不依赖arcpy，创建shp点图层并添加投影。同时，将详细讲解如何从文本数据（如a.txt）中读取经纬度信息，生成点图层，并进一步转化为tif栅格文件。

摘要由CSDN通过智能技术生成

利用shapefile创建shp点图层并添加投影，无需使用arcpy里面的CreateFeatureclass_management或者MakeXYEventLayer_management创建点图层（另一篇博客中将讲述利用CreateFeatureclass_management工具创建点图层以及生成栅格文件）。MakeXYEventLayer_management中可以将生成的面图层，以*.lyc文件的形式存储，或将面图层放到缓存中，并不需要输出成文件，减少处理文件过程中的写入写出过程。

a.txt文件如下图所示：
TXT的存放规则（注：没有第一行头文件的文本文件，对应列内容为 x,y,value）

以下是源代码：

# -*- coding: utf-8 -*-
import shapefile
from osgeo import osr
# 创建新shp文件
shp_filename = r"C:\Users\asus\Desktop\test/shapefiles/test/point/shdemo.shp"
shp