本文主要是介绍matlab读取中国1km分辨率逐月降水量/逐月平均气温数据集(1901-2022),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
开源数据集“中国1km分辨率逐月平均气温数据集(1901-2022)”和“中国1km分辨率逐月降水量数据集(1901-2022)”,原始数据格式为NETCDF,即.nc格式,现需要转为tif格式,采用matlab读取,代码如下:
数据的引用:
彭守璋.中国1km分辨率逐月平均气温数据集(1901-2022).国家青藏高原科学数据中心,10.11888/Meteoro.tpdc.270961[PENG Shouzhang.1-km monthly mean temperature dataset for china (1901-2022).National Tibetan Plateau / Third Pole Environment Data Center,10.11888/Meteoro.tpdc.270961]
彭守璋.中国1km分辨率逐月降水量数据集(1901-2022).时空三极环境大数据平台,10.5281/zenodo.3185722[PENG Shouzhang.1-km monthly precipitation dataset for China (1901-2022).A Big Earth Data Platform for Three Poles,10.5281/zenodo.3185722]
- 三年一期nc文件,共36个波段,
ncFilePath=['F:\pre_2000_2002.nc'];lon=ncread(ncFilePath,'lon');lat=ncread(ncFilePath,'lat');time=ncread(ncFilePath,'time');pre=ncread(ncFilePath,'pre');k=0;for y=2000:2002for j=1:12k=k+1;pre1=pre(:,:,k); data=flipud(pre1');data(data==-32768)=NaN;R = georasterref('RasterSize', size(data),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);%地理栅格数据参考对象(类)filename1=['F:\pre_',num2str(y),'_',num2str(j),'.tif'];geotiffwrite(filename1,data,R);endend
-
一年一期nc文件,共12个波段
ncFilePath=['F:\tmp_2020.nc'];lon=ncread(ncFilePath,'lon');lat=ncread(ncFilePath,'lat');time=ncread(ncFilePath,'time');tmp=ncread(ncFilePath,'tmp');for k=1:12tmp1=tmp(:,:,k); data=flipud(tmp1');data(data==-32768)=NaN;R = georasterref('RasterSize', size(data),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);%地理栅格数据参考对象(类) filename1=['F:\tmp_2020\tmp_2020','_',num2str(k),'.tif'];geotiffwrite(filename1,data,R);end
这篇关于matlab读取中国1km分辨率逐月降水量/逐月平均气温数据集(1901-2022)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!