基于ANUSPLIN的批量气象插值-从数据处理到最终结果(6)

基于python的grd格式批量转为tif格式

grd格式是传统的arcinfo的格式,不适合进一步的计算,需要将grd转换为tif格式,本文采用arcgis中自带的python来进行批量转换。
具体代码如下所示:

import arcpy
arcpy.env.workspace="D:\\chazhi\\"  #存放数据的文件夹
a=arcpy.ListRasters("*","grd")  #得到文件夹下所有的grd名称
for i in a:
    arcpy.RasterToOtherFormat_conversion(i,"D:\\chazhi\\tif\\","TIFF")

将上述代码复制到自带的python上后运行即可。

基于matlab的气象要素异常值处理

插值完后会发现气象要素如降水会存在着负值,负值是由于引进了协变量DEM造成的,可以认为出现的负值的地方降水为0的,尤其是在进入日尺度上插值时,当有n多个这样文件需要处理时,通过arcgis中的栅格计算器来一个个实现尤为费时,本文提供一个基于matlab的处理,同时也能够加入投影信息。

[a,R]=geotiffread('F:\项目\dem.tif');%先导入投影信息
info=geotiffinfo('F:\项目\dem.tif');
e=dir('*.tif'); 
for i=1:size(e,1)
      data=importdata(e(i).name);
      data(data<0)=0;
      geotiffwrite(e(i).name,data,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
end

推荐阅读更多精彩内容