基于python的grd格式批量轉(zhuǎn)為tif格式
grd格式是傳統(tǒng)的arcinfo的格式蹋砚，不適合進一步的計算，需要將grd轉(zhuǎn)換為tif格式尼荆，本文采用arcgis中自帶的python來進行批量轉(zhuǎn)換。
具體代碼如下所示：
import arcpy
arcpy.env.workspace="D:\\chazhi\\" #存放數(shù)據(jù)的文件夾
a=arcpy.ListRasters("*","grd") #得到文件夾下所有的grd名稱
for i in a:
arcpy.RasterToOtherFormat_conversion(i,"D:\\chazhi\\tif\\","TIFF")
將上述代碼復(fù)制到自帶的python上后運行即可唧垦。
基于matlab的氣象要素異常值處理
插值完后會發(fā)現(xiàn)氣象要素如降水會存在著負值捅儒，負值是由于引進了協(xié)變量DEM造成的，可以認為出現(xiàn)的負值的地方降水為0的振亮，尤其是在進入日尺度上插值時巧还，當有n多個這樣文件需要處理時，通過arcgis中的柵格計算器來一個個實現(xiàn)尤為費時坊秸，本文提供一個基于matlab的處理麸祷，同時也能夠加入投影信息。
[a,R]=geotiffread('F:\項目\dem.tif');%先導(dǎo)入投影信息
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
