用泰森多边形法求子流域气象数据

原创

已于 2022-02-22 15:01:15 修改 · 1.1k 阅读

7 ·

CC 4.0 BY-SA版权

文章标签：

#算法 #几何学 #动态规划

于 2022-02-22 14:59:38 首次发布

该脚本详细说明了如何利用泰森多边形法结合ArcGIS处理后的数据，在Matlab环境中计算各子流域的气象数据。首先在ArcGIS中获取站点与子流域的权重信息，保存到Excel文件。然后读取各站点的气象数据，通过站点编号匹配权重，逐时计算子流域的气象指标，最终将结果输出到txt文件。

用泰森多边形法求子流域气象数据 (matlab)

需要先用Arcgis画泰森多边形并求出各面积占比（权重），然后用该matlab脚本计算。

面积占比所在Excel文件格式：三列，分别为站点编号、子流域编号、权重
把各站点数据放在一个目录里，各文件名需要包含站点编号，把目录作为输入。
按注释确定输入，运行。

%%
clc;clear;close all;

%此脚本用于通过泰森多边形方法求子流域气象数据，输入数据和格式包括：


%% 输入
ThiesenResult = xlsread('D:\气象处理\泰森多边形法结果.xlsx','Sheet1');  %三列，分别为站点编号、子流域编号、权重
stationDir = 'D:\evp\';     %站点数据所在的目录（最后有斜杠）
subbasinNum = 224;  % 子流域数目
dataColNum = 4;     % 站点文件中数据所在的列号（在第几列）
outputFile = 'Evp_daily.txt';

%% 读取所有气象文件信息，并建立对应键值对
h=waitbar(