全球干旱指数PDSI可从以下网站下载:
http://climexp.knmi.nl/selectfield_obs2.cgi?id=someone@somewhere
下载下来后是nc格式,将nc格式转化为tif格式后,可以非常容易的在arcgis中用矢量文件裁剪出来自己的研究区域。
本文提供matlab代码来将nc转为tif从而直接提取全部的PDSI数据集,代码如下:
GeoRef = georasterref('Rastersize',[360,720],'Latlim',[-90,90],'Lonlim',[-180,180]);data=ncread('F:\Global\PDSI\scPDSI.cru.3.25.bams2017.GLOBAL.1901.2016.nc','scpdsi');for year=1901:2016 data1=data(:,:,1+12*(year-1901):12*(year-1900)); %得到每年的12个月数据 data3=sum(data1,3)/12; data4=rot90(data3); filename=strcat('F:\Global\PDSI\年尺度pdsi\1全球',int2str(year),'年平均PDSI.tif'); %注意放在自己的路径中 geotiffwrite(filename,data4,GeoRef) for mon=1:12 data2=data1(:,:,mon); data4=rot90(data2); filename=strcat('F:\Global\PDSI\月尺度的pdsi\全球',int2str(year),'_',int2str(mon),'月PDSI.tif');%注意放在自己的路径中 geotiffwrite(filename,data4,GeoRef) endend
将上述代码直接复制到matlab编辑器中的新建脚本中,修改pdsi文件的存放路径,就可以直接运行得到年均值和月值的PDSI数据集。
空气质量指数
Sen+MK趋势分析
VOSviewer文献综述
MATLAB一元线性回归趋势
Hurst指数---判断未来趋势
社区发现算法--FN算法Python实现
蚁群算法求解旅行商(TSP)问题
Manner-Kendall(M-K)---突变检验
ArcPy批量定义投影和批量投影转换
GIMMS NDVI数据处理-nc格式转tif
机器人局部规划算法--DWA算法原理
旅行商(TSP)问题的动态规划解法
ArcGIS时间滑块实现车辆轨迹动态展示
世界各国GDP排名变化--Matlab动图实现
![9b1f4e58647ac78f37e65969fc55d731.png](https://i-blog.csdnimg.cn/blog_migrate/b3de0174485d7c9e4ba150024a67fc87.png)
![5a581971a5b2124922cf9a741f5b324a.png](https://i-blog.csdnimg.cn/blog_migrate/ad6290b9d856108859d07246d2ea93a2.jpeg)
![573b9c4f25a250fea1a34364bef16339.png](https://i-blog.csdnimg.cn/blog_migrate/4643d6509ae5732389bb6f296498577d.png)