【matlab代码】提取任意多边形内的nc数据

2023-04-18 15:50| 来源: 网络整理| 查看: 265

【matlab代码】提取任意多边形内的nc数据--以海洋温度为例子

本文来源于对象想提取一个矩形内的温度数据，从而求平均能看出时间序列变化。

由于平时我们矩形是和经纬度平行，我们可以直接使用lon，lat进行寻找。

例子如图：

这样的水平图，大家都会，直接：

find(lon>=124.5&=26.5&

大家有什么解决方案呢？可以后台私聊。

这里提供一种方法：

解决的结果如下图：

能够得到想要的数据。

至此，问题得到很好的解决。

以下内容为，matlab代码和所需要的数据函数：包括colorbar函数，命名为：colormore_66.m;

地图所需的数据为WOA2023月平均数据2月份数据：woa23_decav91C0_t02_04.nc

地图填充颜色：color.mat文件如下图：

收费不贵，安心使用，文件可联系作者，付费后，获取，；

谢谢理解，坚持分享！感谢各位厚爱！

如喜欢，多谢点赞分享

%% 本文目的写一个程序，寻找温度锋面： % 温度锋面的定义或者说计算方法：Sobel算子 % 具体计算方法见：[1]王兴宇,纪棋严,彭腾腾等.东海黑潮温度锋的三维结构特征分析[J].海洋预报,2022,39(01):67-79. % [1]王建辉. 北欧海海洋锋的时空变化研究[D].浙江海洋大学,2022.DOI:10.27747/d.cnki.gzjhy.2022.000413. % 本文作者：海洋与大气科学 % Burger_miao % 时间：2023年04月11日； % 使用数据为WOA2023的年平均、季节平均和月平均； clear;clc;close all; % 01：读取温度数据 ncdisp 'woa23_decav91C0_t02_04.nc'; file='woa23_decav91C0_t02_04.nc';%二月份温度文件 lon=double(ncread(file,'lon')); lat=double(ncread(file,'lat')); depth=double(ncread(file,'depth')); time=double(ncread(file,'time')); temp=double(ncread(file,'t_an'));%'sea_water_temperature'; % 研究范围：为了方便对比本次计算结果的准确性， % 研究范围和以往研究范围一致，120-131°E，22-33°E； % 寻找范围内的温度数据 ln=find(lon>=119&lon=21.5&lat

【本文地址】

公司简介

联系我们