用MATLAB读取.nc文件，并绘图

学习笔记

---

前言

提示：用MATLAB简单的读取全球海表CO2分压数据，并绘制成图

这部分相对比较容易，简单易学，小白教程。

---

提示：以下是本篇文章正文内容，下面案例可供参考

一、M_map绘图工具箱是什么？

最后成图会用到一个专门用来绘制各种带有地形、海岸数据的工具包

可以到官网上进行下载，下载后添加到MATLAB安装路径文件下的toolbox文件夹中即可

M_Map: A Mapping package for Matlab (ubc.ca)---------------附下载地址，详细的安装教程在哔站都有

二、使用步骤

以我最近处理的数据为例，详细的代码如下

clc;
clear;
close;

%%读取文件
filename='C:\Users\huain\Desktop\JMA_co2map_2020.nc';
ncdisp(filename);%查看文件

lat1=ncread(filename,'lat');%提取纬度数据
lon1=ncread(filename,'lon');%提取经度数据
pCO2s=ncread(filename,'pCO2s');%提取海表CO2分压数据
time=ncread(filename,'time');%读取时间序列（1-12个月，2022年数据）

%%数据预处理
[lat,lon]=meshgrid(lat1,lon1);
pco2s=pCO2s(:,:,2);%提取第一层的数据

%%利用M_map工具画图
m_proj('robinson','lon',[120,300],'lat',[-20,20]);%将范围裁剪到你所需要的区域
m_pcolor(lon,lat,pco2s);
m_coast('patch',[.7 .7 .7],'edgecolor','none');
m_grid('tickdir','out','linew1',2);
colormap(m_colmap('jet','step',10));

h=colorbar('northoutside');
title(h,'pCO2/(uatm) at 2020 yr','fontsize',14);
set(h,'pos',get(h,'pos')+[.2 .05 -.4 0],'tickdir','out')
set(gcf,'color','w');   % Need to do this otherwise 'print' turns the lakes black
hold on

---

我截取了部分我所需要的区域，最后成图如下