1. 数据的准备
data.txt
点数据:(纬度,经度,值)
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运行
4. 生成shape文件
在ide中运行完成后生成的结果:
5. 生成栅格数据
打开arcgis的工具箱arctoolbox
选择Conversion tools -> to Raster -> point to Raster