基于ANUSPLIN的批量气象插值 |
您所在的位置：网站首页 › 最终的数据 › 基于ANUSPLIN的批量气象插值 |
|
引言
气象数据是进行气候变化研究必不可少的基础数据，对其空间插值是实现面源数据的重要步骤，空间插值的方式是多种的，如常见的Kring，IDW等，上述两种通过arcgis既能够快速的完成，但当我们的研究时间尺度较短，需要的数据量特别大时，即需要我们进行批量的气象数据插值，这里假设我们需要同时插值全国2000-2015年日气压，日相对湿度，日风速，日照时数，空间分辨率为8km，时间分辨率为1km，将会产生366416=23424个文件，AUSPLIN是专业的气象插值软件，主要是通过薄样条函数进行插值，协变量的引入（通常是高程，气象要素的分布一般与高程相关）能够使得插值精度更高，插值结果更为平滑。 基于matlab的文本气象数据初处理 基于matlab的文本批量拼接首先从中国气象数据网上下载中国地面气候资料日值数据集（V3.0）。如下图所示： 我这里下载的是1951-2016年的，包括了日气压，日相对湿度，日风速，日照时数等 image.png以下是下载下来的气压的文本文件，可以看到每个文本文件是以月为单位的，我们需要将月的全部拼接到年上去。 image.png 代码如下所示： % filename='H:\Day1951-2012yuanshi\';%%根目录，自行修改 e=dir(); %得到文件夹的名称 for i=5:10 file=strcat(filename,e(i).name,'\'); %得到子文件夹 if i==5 id='-13011-'; %对应不同的数据 end if i==6 id='-10004-'; end if i==7 id='-13003-'; end if i==8 id='-14032-'; end if i==9 id='-12001-'; end if i==10 id='-11002-'; end for m=2000:2015 value=[];value1=[];value2=[]; for n=1:12 if n |
| CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |