Matlab利用shp文件裁剪netCDF（nc）数据及将数据写为nc文件

Matlab利用shp文件裁剪区域数据

---

前言

Matlab利用shp文件裁剪区域数据

一、第一种方法

%读取shp文件
shape = shaperead('shapefile.shp');
% 读取 netCDF 文件
nc = ncinfo('data.nc');
data = ncread('data.nc', 'data');
% 确定裁剪区域
lon = ncread('data.nc', 'lon');
lat = ncread('data.nc', 'lat');
lonlim = [min(shape.X) max(shape.X)];
latlim = [min(shape.Y) max(shape.Y)];
ilat = find(lat >= latlim(1) & lat <= latlim(2));
ilon = find(lon >= lonlim(1) & lon <= lonlim(2));
% 裁剪数据
data_new = data(ilon, ilat);

二、第二种方法（通过调用函数）

data_new = maskregion(lon,lat,var,shape);

function varout = maskregion(lon,lat,var,shp_file)
% lon 为原数据的lon
% lat 为原数据的lat
% var 为原始数据（变量）
% shp_file为要裁剪的shp file
[xx,yy]=meshgrid(lon,lat);
index=zeros(size(var,2),size(var,1));
for kk=1:length(shp_file)
    out=[W(kk).X ;W(kk).Y]';
    aa = inpolygon(xx,yy,out(:,1),out(:,2)) +0;
    try
        index=index+aa;
    catch
        index=index+aa';
    end
end
index(index==0)=nan;
for kk=1:size(var,3)
    varout(:,:,kk)=var(:,:,kk).*index';
end
end

将数据写为nc

% 创建一个新的 netCDF 文件并写入裁剪后的数据
ncwriteatt('data_cropped.nc', '/', 'Conventions', 'CF-1.6');
ncwriteatt('data_cropped.nc', '/', 'title', 'Cropped Data');
ncwriteatt('data_cropped.nc', '/', 'history', ['Created on ' datestr(now)]);
ncwrite('data_cropped.nc', 'lon', lon(ilon));
ncwrite('data_cropped.nc', 'lat', lat(ilat));
ncwrite('data_cropped.nc', 'data_new', data_new);