将在txt中的点数据批量arcgis插值后再输出为excel格式--未完

任务：储存在txt中的点数据进行插值，然后将插值输出成带有经纬度的excel文件。

这个步骤我分解成5步

①将点数据转换成shp文件

②对shp文件进行插值后生成栅格文件

③栅格转点，还是shp文件

④对shp文件的表进行操作，增加XY坐标

⑤表转excel输出

如果是一个txt或者excel文件，

①首先在arcgis里面可以通过将“文件”——“添加数据”——“添加xy数据”的方法将点数据导入，然后在Layer图层栏，点鼠标右键“数据”——“导出”为点Shape文件；

②在“ArcToolBox”——“Spatial Analyst”——“插值工具”

Idw (in_point_features, z_field, {cell_size}, {power}, {search_radius}, {in_barrier_polyline_features})

③“转换工具箱”——“由栅格转出”——“栅格转点”

RasterToPoint_conversion (in_raster, out_point_features, {raster_field})

④“数据管理”——“要素工具集”——“添加XY坐标”

AddXY_management (in_features)

⑤“转换工具箱”——“excel”——“表转excel”

方法一：利用模型构建器ModelBuilder

第一步还是采用python转成shp，

然后卡在批量shp进行插值，因为反距离权重法要求输入是一个shp文件，不是数据集，不知道这个怎么解决，遂放弃。

对于要素选择 要素类

（下班了，此处未写完）

希望有看到的网友指导一下。

方法二：利用python

这个网址很有用http://resources.arcgis.com/zh-cn/help/main/10.1/index.html#/na/009z0000006m000000/

第一步：将txt点文件转shp文件，搜索网页发现虎贲提供了代码，

http://blog.sina.com.cn/s/blog_79c7105501015o6m.html，

帮助解决了这个问题，并打开了思路，后续步骤都可以利用python脚本来实现，再次表示感谢，谢谢无私分享！

import arcgisscripting, sys, string, os
gp = arcgisscripting.create()
try:
pathin="E:/wttdata/InputTxt"
pathout="E:/wttdata/Output/Point2Shp"
files=os.listdir(pathin)
for name in files:
in_Table = pathin+"/"+name
in_x = "nodeLon"
in_y = "nodeLat"
temp=name.strip(".txt")
out_Layer =temp.replace("-","_")
spref = r"Coordinate Systems\Geographic Coordinate Systems\World\WGS 1984.prj"
gp.MakeXYEventLayer(in_Table, in_x, in_y, out_Layer, spref)
pDSC = gp.describe(out_Layer)
print gp.getcount(out_Layer)
gp.FeatureClassToShapefile(out_Layer,pathout)
except:
print gp.GetMessages()
print "ok"

第二步：对shp文件进行插值，此处采用IDW进行插值

这里对小F福表示感谢！http://blog.csdn.net/u010603297/article/details/46535943

import arcpy
from arcpy import env
from arcpy.sa import *

arcpy.env.workspace="E:/wttdata/Output/Point2shp"
fcs=arcpy.ListFeatureClasses("*.shp")

outFeaturePath="E:/wttdata/Output/tiffiels/TIF"

for fc in fcs:
arcpy.CheckOutExtension("Spatial")

inFeature=fc
zField="prec"
cellSize=0.08
power=2
searchRadius = RadiusVariable(12,150000 )
outIdw= Idw(inFeature,zField, cellSize, power, searchRadius)
outIdw.save(outFeaturePath+fc[0:14]+".tif")

print("ok")

第三步：栅格转点

import arcpy
from arcpy import env
from arcpy.sa import *

arcpy.env.workspace="E:/wttdata/Output/Tiffiels"
rasterList=arcpy.ListRasters("*","tif")

outRasterPath="E:/wttdata/Output/Raster2Point/p"

for raster in rasterList:
arcpy.CheckOutExtension("Spatial")
inRaster=raster
outPoint=outRasterPath+raster[3:17]
arcpy.RasterToPoint_conversion(raster, outPoint, "VALUE")
print("ok")

第四步：增加XY坐标，第一步采用地理坐标系wgs_1984后，此处执行添加 XY 坐标 (数据管理)，则 POINT_X 和 POINT_Y 分别表示经度和纬度。

import arcpy
from arcpy import env
env.workspace = "E:/wttdata/Output/Raster2Point"
fcs=arcpy.ListFeatureClasses("*.shp")

for fc in fcs:
arcpy.CheckOutExtension("Spatial")
in_data= fc
in_features = fc[0:14]+"XY.shp"
arcpy.Copy_management(in_data, in_features)
arcpy.AddXY_management(in_features)
print("ok")

第五步：表转excel

import arcpy
from arcpy import env
env.workspace = "E:/wttdata/Output/Raster2Point"
fcs=arcpy.ListFeatureClasses("*XY.shp")
print fcs
for fc in fcs:

arcpy.CheckOutExtension("Spatial")
inTable= fc
outExcel = fc+".xls"

arcpy.TableToExcel_conversion(inTable,outExcel)
print("ok")