1. 数据的准备
data.txt
点数据：（纬度，经度，值）
data.png
2. shp模板准备（可以直接通过arcgis生成）
我把自己做好的放到了百度盘中，有需要的可以下载链接：https://pan.baidu.com/s/1PhVJBGinK1GcCKQSRvSynQ
提取码：as6x
3. 将数据写入到shp中
demo.py
import sys
import arcpy
shp=r'E:\\module\\DZS.shp'
txt=r'E:\\data.txt'
cur=arcpy.da.InsertCursor(shp,["shape@","z"])
f = open(txt,"r")
inter=1 # 需要跳几个点，不跳点就是1，跳几个点依次加1
count=0
while True:
line = f.readline()
if not line:
break
lineSplits = line.split(' ')
length=len(lineSplits)
if length != 3:
continue
count=count+1
if count%inter!=0:
continue
flon = float(lineSplits[1])
flat = float(lineSplits[0])
z = str(lineSplits[2])
row = ((flon, flat, z), z)
cur.insertRow(row)
f.close()
print('success')
在arcgis中运行demo.py
按F5运行
image.png
4. 生成shape文件
在ide中运行完成后生成的结果：
shape.png
5. 生成栅格数据
打开arcgis的工具箱arctoolbox
选择Conversion tools -> to Raster -> point to Raster
arctoolbox.png
设置.png
6. 生成的结果
image.png
