[游戏开发] MSWEP数据nc格式转tif格式

开发: C++知识库 Java知识库 JavaScript Python PHP知识库人工智能区块链大数据移动开发嵌入式开发工具数据结构与算法开发测试游戏开发网络协议系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑笔记本显卡显示器固态硬盘硬盘耳机手机 iphone vivo oppo 小米华为单反装机图拉丁

-> 游戏开发 -> MSWEP数据nc格式转tif格式 -> 正文阅读

[游戏开发]MSWEP数据nc格式转tif格式

这里将数据从0.1°重采样到0.25°

%%
clc;
clear;

%%
fullpath = mfilename('fullpath');
[path, name] = fileparts(fullpath);

%%
ncpath = strcat(path, '\3hourly_nc\');
tifpath = strcat(path, '\tiffile_test\');
files = dir(strcat(ncpath, '*.nc'));
len = length(files);

for i = 1 : len
    filename = files(i).name;
    filename = filename(1 : size(filename, 2) - 3);
    ncfile = strcat(ncpath, filename, '.nc');
    ncinf = ncinfo(ncfile);
    
        
    if size(filename, 2) > 19
        filename = filename(1 : 19);
    end
    
    %
    stryear = filename(1 : 4);
    strdays = filename(5 : 7);
    strhour = filename(9 : 10);
    
    filedate = datetime(strcat(stryear, '-01-01')) + caldays(str2double(strdays) - 1);
    filedate = char(filedate);
    filedate = filedate([1 : 4, 6 : 7, 9 : 10]);
    filedate = strcat(filedate, strhour);
       
    OutputPath = strcat(tifpath, filedate, '.tif');
    
    if exist(OutputPath, 'file') == 2
        [num2str(i), '	', num2str(len), '	', filename, '	', filedate, '	continue']
        continue;
    end
    
    [num2str(i), '	', num2str(len), '	', filename, '	', filedate]
    
    preciSets = ncread(ncfile, 'precipitation');
    preciSets = rot90(preciSets);
    % preciSets = preciSets';
    lon = ncread(ncfile, 'lon');
    lat = ncread(ncfile, 'lat');
    
    [X, Y] = meshgrid(lon, lat);
    lat1 = 89.875 : -0.25 : -89.875;
    lon1 = -179.875: 0.25 : 179.875;
    [XI, YI] = meshgrid(lon1, lat1);
    preciSets1 = interp2(X, Y, preciSets, XI, YI, 'linear');
    
    Refference = georasterref('RasterSize', size(preciSets1), 'Lonlim', [-180 180], 'Latlim', [-90 90]);
    geotiffwrite(OutputPath, preciSets1, Refference);
end