基于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
