2301_80740374的博客

原创 arcpy批量重分类，分区统计

output_raster = os.path.join(output_path, f'风险等级_{year}.tif')pop_raster = os.path.join(pop_path, f'学术严选_china_{year}.tif')print(f"\n处理第 {index}/{total_files} 个文件 ({progress:.1f}%)")print(f"\n处理第 {index}/{total_files} 个年份: {year}")

原创 gis中批量归一化

模糊隶属度中的线性函数归一化，批量处理，模型构建器。

原创 栅格修改行列数保持一致可以叠加运算

前提先把gis大环境的处理范围设置为标准栅格文件，捕捉栅格也是同样设置，然后批量修改。

原创 模型构建器批量分区统计

print(f'所有文件已合并并保存到 {output_file}')批量汇总，表格里进行筛选、提取各县级单位的指数序列，制图。# 将合并后的数据保存到一个新的 Excel 文件中。# 创建一个空的 DataFrame 用来合并数据。# 获取所有 xls 文件的路径列表。# 遍历所有 xls 文件并合并数据。# 按照文件名顺序排序。

原创 gee下载tvdi之后有缺失值nodata

（Nibble 仅适用于掩膜栅格中的 NoData 值。请注意，掩膜栅格中所有非NoData 像元均会接收输入栅格中的值。这些像元的值及其位置均将被用于向掩膜栅格各 NoData 位置处分配值。NoData 位置将接收输入栅格中被识别为掩膜栅格中最近非 NoData 像元的像元值。gee下载ndvi、lst并用python计算tvdi，出现缺失值nodata，在gis中先复制原始tif，将nodata设置为-9999，tvdico.tif。栅格计算器，给nodata赋值。

原创 镶嵌至新栅格

原因是在像素类型没有修改，如下图改为32位有符号即可。（注：可以参考原始tif的属性来确定这里的输入）使用镶嵌至新栅格工具镶嵌相邻两省的人口时，结果值错误。

原创 SPEI指数处理

col_6_data = df.iloc[:, 5] # 注意：列索引从0开始，所以第6列的索引是5。模型构建器进行分区统计，并转出为Excel表格。分区统计导入研究区矢量，并采用MEAN统计。注意：分区统计结果保存在文件夹，则要加dbf后缀；解析路径中要选择NAME。# 获取文件名的第8至第13个字符作为列标题，，9-14位。用python代码提取数据并汇总为总的excel。）数据量大，设置FTP下载。批量裁剪（模型构建器）。

原创 土壤nc文件处理

IDW插值，用全球矢量掩膜（arcgis自带的continents矢量文件，desktop10.8-templatedata）创建netcdf要素图层，nc文件全英文路径！

原创 arcgis Pro 环境设置