用泰森多边形法求子流域气象数据 (matlab)
需要先用Arcgis画泰森多边形并求出各面积占比(权重),然后用该matlab脚本计算。
-
面积占比所在Excel文件格式:三列,分别为站点编号、子流域编号、权重
-
把各站点数据放在一个目录里,各文件名需要包含站点编号,把目录作为输入。
-
按注释确定输入,运行。
%%
clc;clear;close all;
%此脚本用于通过泰森多边形方法求子流域气象数据,输入数据和格式包括:
%% 输入
ThiesenResult = xlsread('D:\气象处理\泰森多边形法结果.xlsx','Sheet1'); %三列,分别为站点编号、子流域编号、权重
stationDir = 'D:\evp\'; %站点数据所在的目录(最后有斜杠)
subbasinNum = 224; % 子流域数目
dataColNum = 4; % 站点文件中数据所在的列号(在第几列)
outputFile = 'Evp_daily.txt';
%% 读取所有气象文件信息&#x