Arcgis pro中Arcpy对Excel文件进行克里金插值

具体流程：excel→table（投影）→shp→克里金插值（掩膜）

import arcpy
from arcpy import env
from arcpy.sa import *
import pandasarcpy.env.overwriteOutput = True  # 允许覆盖之前产生的数据集，便于测试'''获取excel中sheet名称'''
file_name = "G:/Arcgis/ChangshaXian/Rainfall data.xlsx"
excel_total = pandas.ExcelFile(file_name)
excel_sheet = excel_total.sheet_names'''excel→表→shp→投影'''
arcpy.env.workspace = "G:/Arcgis/ChangshaXian"  # 设置工作空间
out_table = "G:/Arcgis/ChangshaXian/CSX_table"  # 表输出路径
out_shp = "G:/Arcgis/ChangshaXian/CSX_shp"  # shp输出路径
i = 0
while i < len(excel_sheet):arcpy.ExcelToTable_conversion("Rainfall data.xlsx", out_table + "/" + excel_sheet[i], excel_sheet[i])with arcpy.EnvManager(outputCoordinateSystem='PROJCS["WGS_1984_UTM_Zone_50N",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]'):arcpy.management.XYTableToPoint(out_table + "/" + excel_sheet[i] + ".dbf", out_shp + "/" + excel_sheet[i],"经度", "纬度", "rainfall")i += 1'''克里金插值'''
z_field = "rainfall"
semiVariogram_props = KrigingModelOrdinary("SPHERICAL")
cell_size = 20
search_radius = "VARIABLE 12"  # 搜索半径（可选）
out_variance_prediction_raster = None  # 预测栅格的输出方差（可选）
mask_region = "G:/Arcgis/ChangshaXian/CSX_basemap/研究区域.shp"  # 掩膜区域限制插值范围out_kriging = "G:/Arcgis/ChangshaXian/CSX_kriging"  # 克里金插值输出路径
j = 0
while j < len(excel_sheet):in_point_features = out_shp + "/" + excel_sheet[j] + ".shp"with arcpy.EnvManager(  # 设置输出坐标系、掩膜范围outputCoordinateSystem='PROJCS["WGS_1984_UTM_Zone_50N",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]',extent=mask_region):outKriging = Kriging(in_point_features, z_field, semiVariogram_props, cell_size,search_radius, out_variance_prediction_raster)outKriging.save(out_kriging + "/" + excel_sheet[j])j += 1print("已完成%d个插值" % j)

Arcgis pro中Arcpy对Excel文件进行克里金插值相关推荐

克里金插值中重要参数变量
ArcGIS中,克里金插值是地统计向导中地统计插值创建表面的一个重要模块,其中包括普通克里金.简单克里金.通用克里金.指示器克里金.概率克里金.析取克里金.经验贝叶斯克里金和面插值. 涉及到克里金插值 ...
在ArcGIS中利用降雨量数据进行克里金插值
今天来分享下如何在ArcGIS中插值,使用的源数据为降水量数据,数据集名称为中国国家级地面气象站基本气象要素日值数据集(V3.0),来源于国家科学气象信息中心(http://data.cma.cn/) ...
ArcGIS中使用协同克里金插值（co-kriging interplotation ）对气象数据插值
ArcGIS中如何使用协同克里金插值(co-kriging interplotation )对气象数据插值 ANUSPLIN气象站点数据插值局限性百度搜索ArcGIS 克里金插值搭建梯子搜索Arc ...
ArcGIS之克里金插值教学
本文来自:GIS科研实验室基本概念 1.什么是克里金插值? 克里金插值又称空间局部插值法,是以半变异函数理论和结构分析为基础,在有限区域内对区域化变量进行无偏最优估计的一种方法,是地统计学的主要内容 ...
matlab 克里金插值,克里金插值(arcgis克里金插值步骤)
1. 克里格方法概述克里格方法(Kriging)又称空间局部插值法,是以变异函数理论和结构分析为基础, 在有限区域内对区域化变量进行无偏最优估计的一种方法,是地. 克里金差值最后的出来的克里金误差有 ...
C语言给出点坐标进行克里金插值,Arcgis笔记之克里金插值——求助surfer8.0
arcgis的插值方法比较丰富,但是有一点,那就是克里金插值.arcgis中克里金插值的结果不平滑,看起来非常难看.如下图,是对图中点数据某个属性值的插值结果. arcgis克里金插值结果使用su ...
利用ArcGIS的模型构建器实现批量投影栅格和批量克里金插值
博主要批量投影几百个栅格,但是ArcGIS中的投影栅格的批处理工具需要手动一条条设置参数,还是麻烦,因此学习了使用模型构建器(Model Builder)来批量投影栅格,可以导入需投影栅格的文件夹,然 ...
ArcGIS克里金插值ERROR999999解决办法汇总
1.插值点shp文件坐标系与插值范围坐标系不一致解决办法:两种shp数据其中修改一个坐标系,直接投影转换就行. 2.用csv文件插值的时候通常是提前生成一个Event事件,但是Event事件有时候不 ...
克里金插值（Kriging）在MATLAB中的实现（克里金工具箱）
一,直接献上克里金插值MATLAB工具箱链接:https://pan.baidu.com/s/1p0lt2G8KQ-els7Dyw5jPFg 提取码:wcss 下载后将该程序添加到MATLAB安装文 ...

Arcgis pro中Arcpy对Excel文件进行克里金插值

Arcgis pro中Arcpy对Excel文件进行克里金插值相关推荐

最新文章

热门文章