CMEMS数据（海表盐度）nc转tif代码MATLAB

最新推荐文章于 2024-03-06 19:30:00 发布

帅气土狗666

最新推荐文章于 2024-03-06 19:30:00 发布

阅读量1.5k

点赞数

分类专栏：科研笔记文章标签： matlab

本文链接：https://blog.csdn.net/qq_41651338/article/details/108230292

版权

科研笔记专栏收录该内容

6 篇文章 0 订阅

订阅专栏

我又来了码代码啦
这次码的是CMEMS盐度数据，nc格式，我又给他转了个tif，这个数据的daily数据是三维的，换算成北京时间应该是凌晨四点左右，所以还挺好用的，不用管UTC和北京时间的差距。然后空间分辨率是0.25度，转tif的时候需要转置不用翻转南北方向。

ncinfo('C:\Users\Larina\Downloads\metoffice_coupled_orca025_GL4_SAL_b20190803_dm20190801.nc');
sss=ncread('C:\Users\Larina\Downloads\metoffice_coupled_orca025_GL4_SAL_b20190803_dm20190801.nc','so');
sss=sss(:,:,1);%用的是海表盐度，所以直接调用第一层的数据就成
sss=sss';
lon=ncread('C:\Users\Larina\Downloads\metoffice_coupled_orca025_GL4_SAL_b20190803_dm20190801.nc','lon'); %读取经度
lat=ncread('C:\Users\Larina\Downloads\metoffice_coupled_orca025_GL4_SAL_b20190803_dm20190801.nc','lat'); %读取维度
R = georasterref('RasterSize', size(sss),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]); 
geotiffwrite('sss.tif',sss,R);

这就是所有的代码啦