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

我又来了码代码啦
这次码的是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);

这就是所有的代码啦