Matlab读取葵花卫星NC数据后显示图像

记

初学者，记录下来。虽然很简单的东西，但是真是开头难。希望能够稍微帮助大家。

读取文件

因为我需要 3.9um 和 11.2um

clear all;clc;

file='NC_H08_20190401_0200_R21_FLDK.02401_02401.nc_'
ncdisp(file)
longitude=ncread(file,'longitude');
latitude=ncread(file,'latitude');
tbb_07=ncread(file,"tbb_07");
tbb_14=ncread(file,"tbb_14");

对于不同通道图显示

imagesc(tbb_14);
colormap(gray);
colorbar;

采用mmap进行显示。

m_proj('Mercator','longtitude',[117,123],'latitude',[34,41]);
m_contourf(longitude,latitude,tbb_03);
hold on
m_coast('color',[0 0 0],'linewidth',2);  %绘制海岸线，填充陆地
m_grid('box','fancy') 
hold off

整体代码：

1. 第一种

%%读取nc文件

clear all;clc;

file='NC_H08_20190401_0200_R21_FLDK.02401_02401.nc_'
ncdisp(file)
longitude=ncread(file,'longitude');
latitude=ncread(file,'latitude');
tbb_07=ncread(file,"tbb_07");
tbb_14=ncread(file,"tbb_14");

%%数据保存形式
longitude=longitude(741:881);
latitude=latitude(381:521);

%%选择显示区域的范围
tbb_07=tbb_07(520:-1:380,1880:2020); %%数据保存形式
tbb_14=tbb_14(520:-1:380,1880:2020);

imagesc(tbb_14);
colormap(gray);
colorbar;

图像的数据可以通过 窗口调节

2. 第二种采用mmap

%%读取nc文件

clear all;clc;

file='NC_H08_20190401_0200_R21_FLDK.02401_02401.nc_'
ncdisp(file)
longitude=ncread(file,'longitude');
latitude=ncread(file,'latitude');
tbb_07=ncread(file,"tbb_07");
tbb_14=ncread(file,"tbb_14");

%%数据保存形式
longitude=longitude(741:881);
latitude=latitude(381:521);

%%选择显示区域的范围
tbb_07=tbb_07(520:-1:380,1880:2020); %%数据保存形式
tbb_14=tbb_14(520:-1:380,1880:2020);

m_proj('Mercator','longtitude',[117,123],'latitude',[34,41]);
m_contourf(longitude,latitude,tbb_07);
hold on
m_coast('color',[0 0 0],'linewidth',2);  %绘制海岸线，填充陆地
m_grid('box','fancy') 
hold off