教程 | 基于Matlab的干旱指数PDSI及CRU全球气象nc数据的处理过程-由nc文件转变为tif文件

全球干旱指数PDSI可从以下网站下载：
http://climexp.knmi.nl/selectfield_obs2.cgi?id=someone@somewher

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);

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

然后在文本中他添加经纬度信息，结果如下所示：

添加经纬度信息后，利用arcgis的ASCII码转raster功能，将该文本转变为栅格文件，并进一步输出为tif，假设文件名为样例.tif。
然后加载带有投影信息的其他栅格文件，对样例.tif文件定义投影，投影方式与加载进来的栅格文件保存一致。
经过上述步骤可得到样例数据。接下来进行批量的转换。

[aaaaa,R]=geotiffread('H:\Global\PDSI\example1.tif');%先导入纬度数据info=geotiffinfo('H:\Global\PDSI\example1.tif');
data=ncread('H:\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);
        data5=flipud(data4);
        filename=strcat('H:\Global\PDSI\年尺度pdsi\全球',int2str(year),'年PDSI.tif');
        geotiffwrite(filename,data5,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
    for mon=1:12
        data2=data1(:,:,mon);
        data4=rot90(data2);
        data5=flipud(data4);
        filename=strcat('H:\Global\PDSI\月尺度的pdsi\全球',int2str(year),'_',int2str(mon),'月PDSI.tif');
        geotiffwrite(filename,data5,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
    end
end

通过上述代码即可实现批量处理nc格式的PDSI文件了, 在新的版本中，可以不用样例文件直接进行提取，在VX公众号地学分析与算法中，结果如下：

CRU文件的处理过程

data=ncread('H:\Global\CRU\cru_ts4.01.1901.2016.tmx.dat.nc','tmx');
data3=data(:,:,1);
data4=rot90(data3);
data5(isnan(data5))=-999;
dlmwrite('样例2.txt',data5,'\t',1,1);

其余步骤同上构建一个有投影的样例tif数据，批量读取与转换

[aaaaa,R]=geotiffread('H:\Global\CRU\样例2.tif');%先导入纬度数据
info=geotiffinfo('H:\Global\CRU\样例2.tif');
data=ncread('H:\Global\CRU\cru_ts4.01.1901.2016.tmx.dat.nc','tmx');
for year=1901:2016
        data1=data(:,:,1+12*(year-1901):12*(year-1900)); %得到每年的12个月数据
        data3=sum(data1,3)/12; %对年数据求平均值，得到年平均最大气温，如果是降水，则直接去掉/12
        data4=rot90(data3);
        filename=strcat('H:\Global\CRU\tif\year\CRU',int2str(year),'_Tmx.tif');
        geotiffwrite(filename,data4,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
    for mon=1:12
        data2=data1(:,:,mon);
        data4=rot90(data2);
        filename=strcat('H:\Global\CRU\tif\month\CRU',int2str(year),'_',int2str(mon),'_Tmax.tif');
        geotiffwrite(filename,data4,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
    end
end

作者：画长空_yin
链接：https://www.jianshu.com/p/1cae03e4452e
来源：简书
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。

每日一课（免费看）

Ai尚研修最新免费教程模块：点击以下文章标题直接观看

1、每日一课 | 基于MAXENT模型的GIS空间数据剪裁、拼接与提取2、每日一课 | 基于R语言极值统计学及其相关领域中的应用3、每日一课 | 海量遥感数据处理与GEE云计算技术实践应用案例分析4、每日一课 | AQUATOX水环境与水生态模型矿化模块、植物模块详解5、每日一课 | python气象水文数据处理——Pandas中时间和日期处理6、每日一课 | 基于PyTorch深度学习无人机遥感影像——图像分割7、每日一课 | 生态环境领域中meta分析图形分析法——正态分位数图8、每日一课 | MATLAB二维曲线的基本绘图命令详解9、每日一课 | SWMM模型中如何利用GIS数据快速建立模型10、每日一课 | 地下水环评必备：地下水模拟及环评之水文地质基础与建模流程、数据要求11、每日一课 | 土壤物理模型Hydrus模型构架、典型案例分析及应用12、每日一课 | R-SEM分析（从频率学派到贝叶斯学派）13、每日一课 | 新型土地利用变化模拟模型软件研发及应用14、每日一课 | MATLAB深度学习工具箱功能详解15、每日一课 | 生态环境领域中如何用Meta分析剔除无效数据16、每日一课 | 陆面蒸散模拟净辐射（Net radiatation）计算方法17、每日一课 | MATLAB二维曲线的基本绘图命令详解18、每日一课 | SWMM模型中GIS在网络拓扑关系建立方面的应用19、每日一课 | SWMM模型中如何利用GIS数据快速建立模型20、每日一课 | 深度学习遥感图像检测案例讲解（1）——代码讲解深度学习遥感图像检测案例讲解（1）——代码讲解21、每日一课 | HEC-RAS典型应用案例：EPS读断面坐标高程点22、每日一课 | 深度学习在遥感地物分类、目标识别和图像分割实践应用中让您工作效率翻倍的实战技巧总结部分教程展示，关注Ai尚研修公众号→每日一课，海量技术资源免费学习！

上下滑动查看更多

系统性实践课程（直播课）

点击课程标题，直接进入课程详情

1、10月23日-25日 | 最新基于MATLAB编程、机器学习、深度学习在图像处理中实践技术应用高级培训班
2、10月23-24日、30日-31日 | SMOKE模型排放清单处理技术及在多模式下实践应用方法与VOCs排放量核算培训班3、10月23-24日、30日 | CASA（Carnegie-Ames-Stanford Approach）模型原理及上机操作实现4、10月23-24日、30日 | 结构方程模型（SEM）高阶应用及建模中特殊问题处理方法研修班5、10月24日 | 同源序列寻找及基因家族进化树构建技术免费直播课6、10月30日-31日 | 野外观测站生态气象数据处理分析实践应用直播课7、10月30-31日、11月6日-7日 | 基于R语言的代理模型（高斯过程、贝叶斯优化、敏感性分析、异方差性等）高级技术应用培训班8、10月30-31日、 11月6日 | 基因家族分析及SCI写作技巧实践技术培训班9、10月30-31日、11月6日-7日 | InVEST模型高级进阶及在生态系统服务供需、碳中和、论文写作等实际项目中的具体应用培训班10、10月30-31日、11月6日 | 基于HEC-RAS数值模拟与GIS水文分析技术的洪水危险性及风险评估实践培训班11、11月13日-14日、20日-21日、27日 | 高分辨率区域气候模式（WRF）实践技术学习及在生态环境、水文、能源等多领域中的应用培训班12、11月20-21日 、27日-28日 | 基于MAXENT模型的生物多样性生境模拟与保护优先区甄选、自然保护区布局优化评估及论文写作技巧高级进阶培训班13、11月20-21日、27日 | Noah-MP陆面过程模型建模方法与站点、区域模拟实践技术应用培训班14、11月20－21日、27日 | 合成孔径雷达干涉测量InSAR数据处理、地形三维重建、形变信息提取、监测等实践技术应用高级培训班15、11月20日-21日 27-28日 | 近地面无人机植被定量遥感与生理参数反演实践技术应用进阶培训班16、11月20日-21日 | 基于MATLAB植被变化趋势（变绿-变黄）分析实践应用培训班17、11月27日-28日 | 面向高校的基于算法的发明专利申请写作方法培训班

上下滑动查看更多

系统性实践课程（视频课程，可长期观看）

视频教程+课件及代码资料+云导师长期辅导

（点击数据标题，直接进入）

1、基于SMOKE模型的多模式排放清单处理技术推广课程2、大气颗粒物PMF源解析实践技术应用精品课程3、EKMA曲线及大气O3来源解析视频课程4、WRF DA资料同化系统理论、运行与与变分、混合同化新方法技术应用5、气候水文耦合WRF-Hydro模式实践技术应用精品课程6、空气质量预报模式系统（wrf-cmaq）改进与污染源排放清单建立实践技术精品课程7、区域气象-大气化学在线耦合模式（WRF/Chem）在大气环境中的应用课程8、高分辨率区域气候模式（WRF）在地学、生态环境、水文及能源领域中的实践技术应用精品课程9、基于CAMx的空气质量模拟及污染来源解析技术与案例分析精品课程10、NCL数据分析与处理实践技术应用视频课程11、R语言在气象、水文中数据处理及结果分析、绘图实践技术应用12、基于python语言气象水文数据处理及精美科研绘图实践技术视频课程进入Ai尚研修公众号，点“Ai尚课堂”，进入“全部课程”

上下滑动查看更多

数据资源（一键免费获取）

Ai尚研修海量资源（数据、课件、书籍、视频教程）一键免费领取

（点击数据标题，直接进入）

地学领域数据1、【数据分享】30米精度全球地表覆盖分类数据（2000、2010、2020）2、【数据分享】全国分省、市、县、乡镇行政区划矢量图（专业版shp、小白版ppt）免费下载3、【数据分享】中国月度1KM植被指数（NDVI）空间分布数据集4、【数据分享】中国陆地生态系统服务价值空间分布数据集5、【数据分享】中国农田熟制遥感监测数据集6、【数据分享】多个自然资源（大气、土地、水等）面板数据大放送、时间跨度长、地区分布广，快速下载7、【数据分享】中国30米分辨率坡度数据（免费获取）8、【数据分享】最新！全国人口密度数据2000-2020年9、【数据分享】中国GDP空间分布公里网格数据集10、【数据分享】建筑矢量数据shap，26个省级城市建筑轮廓11、【数据分享】中国100万地貌类型空间分布数据12、【数据分享】中国100万植被类型空间分布数据13、【数据分享】中国人口空间分布公里网格数据集14、【数据分享】自然保护区边界矢量数据下载（国家级、省级、县市级）15、【数据分享】全球100万基础地理数据（包含国家行政边界、铁路、道路、线状水系面状水系）16、【数据分享】全国电子地图SHP数据除了常有的省、市、县、乡的图层信息以外，还包括道路、水系、银行、绿地、高速出入口、医疗机构、餐饮、超市等17、【数据分享】6G超大超级详细矢量数据，来源于国内某地图导航基底数据18、【数据分享】全球冰川数据集RGI v3.219、【数据分享】2020年最新全国城市道路数据集20、【数据分享】2021年最新区县级全国行政区划数据21、【数据分享】ArcGIS矢量数据shp格式大放送：自然保护区和森林公园22、【数据分享】青海三江源草地退化空间分布数据集部分数据展示，关注Ai尚研修公众号→学术资源→数据分享，海量数据资源轻松下载！大气领域数据1、【数据分享】ERA5全球大气再分析数据和MSWEP多源融合降水数据2、【数据分享】MSWEP多源融合降水数据→免费获取3、【数据分享】CMIP6未来气候变化情景数据之一:历史时期月平均气温数据4、【数据分享】CMIP6未来气候变化情景数据之二: 历史时期降水数据集5、【数据分享】CMIP6未来气候变化情景数据之三: 历史时期风速数据集6、【数据分享】CMIP6未来气候变化情景数据之四：服务于降尺度的实测格点化气温数据集7、【数据分享】CMIP6未来气候变化情景数据之五：服务于降尺度的实测格点化降水数据集8、【数据分享】CMIP6未来气候变化情景数据之六：历史时期近地表最高气温数据集9、【数据分享】CMIP6数据分享计划之七：历史时期近地表最低气温数据集10、【数据分享】CMIP6未来气候变化情景数据之八：历史时期向下短波辐射数据集11、【数据分享】中国气象背景数据集12、【数据分享】全行业碳排放核算Excel自动计算表13、【数据分享】碳排放数据集！全国碳排放数据分省、市、县（区）14、【数据分享】1942年至今！全国气象数据：包含气温、气压、露点、风向风速、云量、降水量数据部分数据展示，关注Ai尚研修公众号→学术资源→数据分享，海量数据资源轻松下载！

土利用数据

土地利用数据1、【数据分享】2020年中国土地利用遥感监测数据丨免费下载！2、【数据分享】中国土地利用数据1980-2015 （免费获取）3、【数据分享】土地利用（覆盖）数据下载汇总4、【数据分享】全球1982～2015年逐年长时序动态土地覆盖数据集5、【数据分享】中国土壤侵蚀空间分布数据6、【数据分享】15省《中国土壤分布》土壤分布矢量数据，含使用说明7、【数据分享】世界土壤数据库（HWSD）全球土壤数据下载处理8、【数据分享】AW3D30 DSM（数字地表模型）数据9、【数据分享】2017年10米分辨率全球地表覆盖数据10、【数据分享】中国土壤质地空间分布数据11、【数据分享】中国土地市场网（土地交易数据）部分数据展示，关注Ai尚研修公众号→学术资源→数据分享，海量数据资源轻松下载！水文领域数据1、【数据分享】全国水系图SHP数据免费下载2、【数据分享】基于DEM提取的中国流域、河网数据免费下载3、【数据分享】1:25万全国三级水系流域数据集（免费下载）4、【数据分享】中国水系数据，包含长江黄河、湖泊、以及一级三级四级五级河流矢量数据5、【数据分享】中国九大流域矢量数据6、【数据分享】超详细，东亚水文流域划分数据，包括流域、子流域、实际面积等7、【数据分享】分省份水系数据8、【数据分享】全国高清水文地质图部分数据展示，关注Ai尚研修公众号→学术资源→数据分享，海量数据资源轻松下载！  年鉴类数据1、【数据分享】中国低碳年鉴2000-2016（免费下载）2、【数据分享】全Excel！中国科技统计年鉴1991-2019（免费）3、【数据分享】一个横跨35年的数据，1985-2019中国城市统计年鉴4、【数据分享】28年《中国工业统计年鉴》免费下载5、【数据分享】1985-2019年中国农村统计年鉴免费下载！6、【数据分享】中国农业统计资料1949-2019，免费下载 ！7、【数据分享】国土资源统计年鉴2005-2017&中国林业统计年鉴1992-20198、【数据分享】中国海洋统计年鉴2000-2017（免费）9、【数据分享】中国环境统计年鉴：1998-2018（免费）10、【数据分享】《中国文化和旅游统计年鉴》又是一个24年的统计年鉴集合11、【数据分享】285城空间权重矩阵&中国城乡建设统计年鉴EXCEL12、【数据分享】中国农村贫困监测报告2008-201913、【数据分享】2006-2020《中国卫生统计年鉴》免费下载14、【数据分享】最新！2000-2020年县域统计年鉴面板数据（免费）15、【数据资源】全网最全！34个省（市、自治区）2001-2020年统计年鉴（免费）16、【数据分享】中国城市统计年鉴1985-2020最新为EXCEL格式！17、【数据分享】海量数据大放送！农村研究数据合集18、【数据分享】最新EXCEL中国统计年鉴2020+经管教材16本+python学习资源部分数据展示，关注Ai尚研修公众号→学术资源→数据分享，海量数据资源轻松下载！
经济类数据1、【数据分享】1990-2020年上市公司微观数据2、【数据分享】中国及省级市场化总指数及方面指数：08-193、【数据分享】给实证研究加点料：1990-2019全球治理数据WGI和城市创新创业指数4、【数据分享】各国制度、地理、文化距离5、【数据分享】400万！人口迁徙大数据2019-20206、【数据分享】科研数据利器之二：中国土地市场网（土地交易数据）7、【数据共享】1.82G中国工业企业数据库1998-20138、【数据分享】208个地级市和31个省、市城乡泰尔指数！部分数据展示，关注Ai尚研修公众号→学术资源→数据分享，海量数据资源轻松下载！

上下滑动查看更多

扫码联系我们

扫码联系小编，视频课程最高享受7.5折优惠

END

超级会员专享

如何成为超级会员：

1.凡参加Ai尚研修收费课程即为超级会员

2.充值3500元即为Ai尚研修超级会员；

超级会员特权：

1.参加直播与现场课程专享400元会员优惠，同时成为Ai尚研修推广大使，推荐参会分享收益。
2.购买Ai尚研修平台导师视频课程（含导师随行长期指导），享受8.5折优惠。
3.免费享受Ai尚研修会议平台中会员专享下视频课程。

4.长期免费参与每月4期针对地学领域简学实用公开课（优先、无门槛）
5.免费参与任何一门课程开设的导师面对面线上答疑交流。

6.Ai尚研修资源站分享中数据免费下载；

7.会员既成为Ai尚研修合作伙伴，可以进行技术文章、视频教程投稿，分享收益红利，同时所在单位可以通过Ai尚研修新媒体平台进行合作及招聘推广。

声明： 本号旨在传播、传递、交流，对相关文章内容观点保持中立态度。涉及内容如有侵权或其他问题，请与本号联系，第一时间做出撤回。

END

Ai尚研修丨专注科研领域

技术推广，人才招聘推荐，科研活动服务

科研技术云导师，Easy  Scientific  Research