【程序设计】数据海温的三维视图|三维一直是个难点,所以AI很香,奈何我不会.

【程序设计】数据海温的三维视图|三维一直是个难点,所以AI很香,奈何我不会.

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

%% 本脚本给你们说清楚三维涡旋的事情:
% 上次发布之后,因为数据是我自己处理好的你们没懂;
%这次给你们讲一讲;
%% 微信公众号:海洋与大气科学
%% 时间:2022年11月12日
%% 作者:mapm
ncdisp 'SODA_2.2.4_201012.cdf'
clear;close all;clc;
file='SODA_2.2.4_201012.cdf';
lon=double(ncread(file,'lon'));   %%% 读取lon变量,读取经度;
lat=double(ncread(file,'lat'));   %%% 读取lat变量,读取纬度;
u=double(ncread(file,'u')); %%% 读取u变量,读取速度;
v=double(ncread(file,'v'));%%% 读取v变量,读取速度;
temp=double(ncread(file,'temp'));%%% 读取v变量,读取速度;
salt=double(ncread(file,'salt'));%%% 读取v变量,读取速度;
ssh=double(ncread(file,'ssh'));   %%% 读取lon变量,读取经度;
depth=double(ncread(file,'depth'));
%%
area1=[30 110;-60 -20;0 1500];%32.25
%%%% 根据实际经纬度范围;确定起始和结束位置;
ln1=find(lon>=area1(1,1)&lon<=area1(1,2));
la1=find(lat>=area1(2,1)&lat<=area1(2,2));
ld=find(depth>=area1(3,1)&depth<=area1(3,2));
lonx1=lon(ln1);laty1=lat(la1);depth1=-depth(ld);
temp1=temp(ln1,la1,ld);
temp1 = permute(temp1,[2 1 3]);
[X1,Y1,Z1]=meshgrid(lonx1,laty1,depth1);% depth是深度;
%%
figure
hslicer=slice(X1,Y1,Z1,temp1,[30],[-32.25],[5.01000000000000])%(1:5)画5层的盐度图;
shading flat%(去网格) 
set(hslicer(1),'MarkerEdgeColor','k','LineWidth',2)
set(hslicer(2),'MarkerEdgeColor','k','LineWidth',2)
set(hslicer(3),'MarkerEdgeColor','k','LineWidth',2)
get(hslicer(1))
%% 温度三维视图1
figure
slice(X1,Y1,Z1,temp1,[],[-32.25 -32.25] ,[0:1500])%(1:5)画5层的盐度图;
shading flat%(去网格) 
colormap jet
caxis([0 30])
hold on
axis([30 110 -60 -20 -1500 0])
view(50,45)
xlabel('lon')
ylabel('lat')
zlabel('depth\m')
title('温度三维视图')
set(gca,'xtick',[30:20:110],'xticklabel',{'30°E','50°E','70°E','90°E','110°E'},...
    'ytick',[-60:20:-20],'yticklabel',{'-60°S','-40°S','-20°S'},'box','on',...
    'fontsize',12,'fontweight','bold','fontname','time news roman','linewidth',2);
export_fig 温度三维视图1.png -r300
%% 温度三维视图2
figure
slice(X1,Y1,Z1,temp1,[90],[-32.25 -32.25] ,[0:1500])%(1:5)画5层的盐度图;
shading flat%(去网格) 
colormap jet
caxis([0 30])
hold on
axis([30 110 -60 -20 -1500 0])
view(50,45)
xlabel('lon')
ylabel('lat')
zlabel('depth\m')
title('温度三维视图')
set(gca,'xtick',[30:20:110],'xticklabel',{'30°E','50°E','70°E','90°E','110°E'},...
    'ytick',[-60:20:-20],'yticklabel',{'-60°S','-40°S','-20°S'},'box','on',...
    'fontsize',12,'fontweight','bold','fontname','time news roman','linewidth',2);
export_fig 温度三维视图2.png -r300
%% 温度三维视图3
figure
slice(X1,Y1,Z1,temp1,[109.75],[-59.75] ,[-5.01])%(1:5)画5层的盐度图;
shading flat%(去网格) 
colormap jet
caxis([0 30])
hold on
axis([30 110 -60 -20 -1500 0])
view(50,45)
xlabel('lon')
ylabel('lat')
zlabel('depth\m')
title('温度三维视图')
set(gca,'xtick',[30:20:110],'xticklabel',{'30°E','50°E','70°E','90°E','110°E'},...
    'ytick',[-60:20:-20],'yticklabel',{'-60°S','-40°S','-20°S'},'box','on',...
    'fontsize',12,'fontweight','bold','fontname','time news roman','linewidth',2);
export_fig 温度三维视图3.png -r300
%% 温度三维视图4
figure
slice(X1,Y1,Z1,temp1,[100.75],[-50] ,[-500])%(1:5)画5层的盐度图;
shading flat%(去网格) 
colormap jet
caxis([0 30])
hold on
axis([30 110 -60 -20 -1500 0])
view(50,45)
xlabel('lon')
ylabel('lat')
zlabel('depth\m')
title('温度三维视图')
set(gca,'xtick',[30:20:110],'xticklabel',{'30°E','50°E','70°E','90°E','110°E'},...
    'ytick',[-60:20:-20],'yticklabel',{'-60°S','-40°S','-20°S'},'box','on',...
    'fontsize',12,'fontweight','bold','fontname','time news roman','linewidth',2);
export_fig 温度三维视图4.png -r300
%% 温度三维视图5
figure
slice(X1,Y1,Z1,temp1,[40:100],[-60:-20] ,[-1500:0])%(1:5)画5层的盐度图;
shading flat%(去网格) 
colormap jet
colorbar
caxis([0 30])
hold on
axis([30 110 -60 -20 -1500 0])
view(30,35)
xlabel('lon')
ylabel('lat')
zlabel('depth\m')
title('温度三维视图')
set(gca,'xtick',[30:20:110],'xticklabel',{'30°E','50°E','70°E','90°E','110°E'},...
    'ytick',[-60:20:-20],'yticklabel',{'-60°S','-40°S','-20°S'},'box','on',...
    'fontsize',12,'fontweight','bold','fontname','time news roman','linewidth',2);
set(gcf,'color','w','position',[50 50 1000 750])
export_fig 温度三维视图5.png -r300
%% 以下colormore使用
cmap=(colormore_9(100));
num=9;
figure
slice(X1,Y1,Z1,temp1,[109.75],[-59.75] ,[-5.01])%(1:5)画5层的盐度图;
shading flat%(去网格) 
colormap(cmap)
colorbar
caxis([0 30])
hold on
axis([30 110 -60 -20 -1500 0])
view(50,45)
xlabel('lon')
ylabel('lat')
zlabel('depth\m')
title('温度三维视图')
set(gca,'xtick',[30:20:110],'xticklabel',{'30°E','50°E','70°E','90°E','110°E'},...
    'ytick',[-60:10:-20],'yticklabel',{'-60°S','-50°S','-40°S','-30°S','-20°S'},'box','on',...
    'fontsize',12,'fontweight','bold','fontname','time news roman','linewidth',2);
set(gcf,'color','w','position',[50 50 1000 750])
export_fig(['温度三维视图colormore_',num2str(num),'.png'],'-r300');






  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 3
    评论
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

海洋与大气科学

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值