pdsi matlab,基于Matlab的干旱指数PDSI及CRU全球气象nc数据的处理

image.png

从该网站上可直接下载全球1901-2016年的逐月PDSI数据，文件为nc格式，为进一步和其他栅格数据进行计算，需要将nc文件转变为tif文件，因此本文提供一种能够批量转换的处理方式。

首先准备个样例数据

ncdisp('H:\Global\PDSI\scPDSI.cru.3.25.bams2017.GLOBAL.1901.2016.nc');

data1=ncread('H:\Global\PDSI\scPDSI.cru.3.25.bams2017.GLOBAL.1901.2016.nc','scpdsi');

data3=data1(:,:,1);

data4=rot90(data3);

data5=flipud(data4);

data5(isnan(data5))=-999;

dlmwrite('样例1.txt',data5,'\t',1,1);

经过上述转换后可得到文本格式的样例数据，结果如下：

image.png

然后在文本中他添加经纬度