往期推荐
图片
【python海洋专题一】查看数据nc文件的属性并输出属性到txt文件
【python海洋专题二】读取水深nc文件并水深地形图
【python海洋专题三】图像修饰之画布和坐标轴
【Python海洋专题四】之水深地图图像修饰
【Python海洋专题五】之水深地形图海岸填充
【Python海洋专题六】之Cartopy画地形水深图
【python海洋专题】测试数据
【Python海洋专题七】Cartopy画地形水深图的陆地填充
【python海洋专题八】Cartopy画地形水深图的contourf填充间隔数调整
【python海洋专题九】Cartopy画地形等深线图
【python海洋专题十】Cartopy画特定区域的地形等深线图
【python海洋专题十一】colormap调色
【python海洋专题十二】年平均的南海海表面温度图
【python海洋专题十三】读取多个nc文件画温度季节变化图
【python海洋专题十四】读取多个盐度nc数据画盐度季节变化图
【python海洋专题十五】给colorbar加单位
【python海洋专题十六】对大陆周边的数据进行临近插值
【python海洋专题十七】读取几十年的OHC数据,画四季图
【python海洋专题十八】读取Soda数据,画subplot的海表面高度四季变化图
【python海洋专题十九】找范围的语句进阶版本
【python海洋专题二十】subplots_adjust布局调整
【python海洋专题二十一】subplots共用一个colorbar
【python海洋专题二十二】在海图上text
【python海洋专题二十三】共用坐标轴
【python海洋专题二十四】南海年平均海流图
【python海洋专题二十五】给南海年平均海流+scale
【python海洋专题二十六】南海海流流速图
【python海洋专题二十七】南海四季海流图
【python海洋专题二十八】南海四季海流流速图
【python海洋专题二十九】读取CTD文件数据并画温度点剖面图
【python海洋专题三十】画南海115°E的温度剖面图
【python海洋专题三十一】画南海115°E的地形温度剖面图
【python海洋专题三十二】画南海115°E的地形温度流速剖面图
【python海洋专题三十三】画海洋表面的风场分布
【python海洋专题三十四】调用自己的colormore
【python海洋专题三十五】加密数据–二维插值
【python海洋专题三十六】两个一维数组的相关系数–为海洋指数作准备
【python海洋专题三十七】海洋指数画法–折线图样式一
【python海洋专题三十八】海洋指数画法–折线图样式二
【python海洋专题三十九】海洋指数画法–折线图样式三–不同颜色的线条
【python海洋专题四十】海洋指数画法–单色填充图
【python海洋专题四十一】海洋指数画法–渐变填色图
【python海洋专题四十二】海洋指数画法–双色柱状图
【python海洋专题四十三】海洋指数画法–单色渐变柱状图
【python海洋专题四十四】海洋指数画法–多色渐变柱状图
【python海洋专题四十五】海洋研究区域示意图
【python海洋专题四十六】研究区域示意放大图
【python海洋专题海洋指数画法】大气与海洋指数画法汇总
【MATLAB海洋专题】历史汇总
【matlab程序】(1-5)五坐标轴的精细修饰
【matlab程序】图片平面制作||文末点赞分享||海报制作等
大佬推荐一下物理海洋教材吧?
【matlab海洋专题】高级玫瑰图–风速风向频率玫瑰图–此图细节较多
【上千种颜色包|全平台可用】收集自Matlab、python、R、NCL等颜色包
R语言_RColorBrewer包–全平台可用
海洋专用cmocean颜色包_共22种–全平台可用
【matlab教程】matlab不规则区域的外围填充
代码分享:
% 本次程序处理txt文件,并读取台风数据,并画图。
% 时间:20231123
% 海洋与大气科学
clear;clc;close all;
%% 数据处理
%01 第一个台风文件
data_typhoon1=importdata(‘typhoon1.txt’); % 将生成的typhoon1.txt文件导入工作空间,变量名为data_typhoon1,实际上它不显示出来
lon1=data_typhoon1(:,4);% 台风经度
lat1=data_typhoon1(:,3);% 台风经度
data_typhoon2=importdata(‘typhoon2.txt’); % 将生成的typhoon2.txt文件导入工作空间,变量名为data_typhoon2,实际上它不显示出来
lon2=data_typhoon2(:,4);% 台风经度
lat2=data_typhoon2(:,3);% 台风经度
data_typhoon3=importdata(‘typhoon3.txt’); % 将生成的typhoon3.txt文件导入工作空间,变量名为data_typhoon3,实际上它不显示出来
lon3=data_typhoon3(:,4);% 台风经度
lat3=data_typhoon3(:,3);% 台风经度
data_typhoon4=importdata(‘typhoon4.txt’); % 将生成的typhoon4.txt文件导入工作空间,变量名为data_typhoon4,实际上它不显示出来
lon4=data_typhoon4(:,4);% 台风经度
lat4=data_typhoon4(:,3);% 台风经度
data_typhoon5=importdata(‘typhoon5.txt’); % 将生成的typhoon5.txt文件导入工作空间,变量名为data_typhoon5,实际上它不显示出来
lon5=data_typhoon5(:,4);% 台风经度
lat5=data_typhoon5(:,3);% 台风经度
%% 画图
close all
figure % 建立画布,改为白色背景和画布大小。
set(gcf,‘color’,‘w’,‘position’,[100 50 800 650])
m_proj(‘miller’,‘lon’,[100 120],‘lat’,[10 24])
m_gshhs_i(‘patch’,[.65 .65 .65],‘edgecolor’,‘k’);
% 添加国界
World = shaperead(‘china_country.shp’);
wx = [World(.X];wy = [World(.Y];%读shp文件
hold on
m_plot(wx,wy,‘k’,‘linewidth’,1,‘color’,[0 0 0]);%绘制
% 添加省份
World = shaperead(‘diquJie_polyline.shp’);
wx = [World(.X];wy = [World(.Y];%读shp文件
hold on% 使用经纬度画台风
m_plot(wx,wy,‘k’,‘linewidth’,1,‘color’,[0 0 0]);%绘制
hold on
h1=m_plot(lon1-10,lat1,‘linewidth’,2,‘color’,[0,0,0]./255);
hold on
h2=m_plot(lon2-8,lat2,‘linewidth’,2,‘color’,[0,0,0]./255);
hold on
h3=m_plot(lon3-6,lat3,‘linewidth’,2,‘color’,[139,0,139]./255);
hold on
h4=m_plot(lon4-4,lat4,‘linewidth’,2,‘color’,[255,0,255]./255);
hold on
h5=m_plot(lon5-1,lat5,‘linewidth’,2,‘color’,[220,20,60]./255);
hold on% 加上观测点
m_scatter(109.566,21.583,‘Marker’,‘o’,‘MarkerFaceColor’,‘r’,‘MarkerEdgeColor’,‘r’,‘sizedata’,20)
hold on
m_text(109.666,21.583,‘验证点’,‘fontweight’,‘bold’,‘fontsize’,15)
hold on% 加上图例
m_text(102,18,‘图例’,‘fontweight’,‘bold’,‘fontsize’,15)
hold on
m_text(102,15,‘验证站点’,‘fontweight’,‘bold’,‘fontsize’,15)
hold on
m_scatter(103,15.5,‘Marker’,‘o’,‘MarkerFaceColor’,‘r’,‘MarkerEdgeColor’,‘r’,‘sizedata’,18)
m_grid(‘ytick’,0:2:24,‘xtick’,100:2:120,‘linewidth’,1,‘linestyl’,‘none’,‘box’,‘fancy’,‘fontsize’,15);