MATLAB | 1980—2022 China Meteorological Element Raster Data (Precipitation/Temperature/Wind Speed/Humidity/Pressure/Potential Evapotranspiration)

https://www.ncdc.ac.cn/portal/metadata/21691d03-bef2-4800-924e-5614e7268b87

Using monthly wind speed data as an example, this article demonstrates how to batch convert data to TIF format using MATLAB.

.rtcContent { padding: 30px; } .lineNode {font-size: 10pt; font-family: Menlo, Monaco, Consolas, "Courier New", monospace; font-style: normal; font-weight: normal; }inputFolder = 'E:\数据都在这里\巴音郭勒\ChinaMet_001deg_wind\ChinaMet_001deg_wind\';outputFolder = 'E:\数据都在这里\巴音郭勒\ChinaMet\ChinaMet_001deg_wind\';ncFiles = dir(fullfile(inputFolder, 'ChinaMet_001deg_wind_*.nc'));for k = 1:length(ncFiles)    ncPath = fullfile(inputFolder, ncFiles(k).name);    disp(['正在处理: ', ncFiles(k).name]);    lat = ncread(ncPath, 'lat');    lon = ncread(ncPath, 'lon');    wind = ncread(ncPath, 'wind');    if lat(1) > lat(end)        lat = flip(lat);        wind = flip(wind, 1);    end    wind = rot90(fliplr(wind), 1);    R = georefcells([min(lat) max(lat)], [min(lon) max(lon)], size(wind));    outputFile = fullfile(outputFolder, replace(ncFiles(k).name, '.nc', '.tif'));    geotiffwrite(outputFile, wind, R);    disp(['保存成功: ', outputFile]);enddisp('全部转换完成！');