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

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

海洋与大气科学

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

阅读量540

点赞数

文章标签：人工智能 python 算法

本文链接：https://blog.csdn.net/miaobo0/article/details/127817558

版权

【程序设计】数据海温的三维视图|三维一直是个难点，所以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');