利用matlab对GOCI产品进行可视化

一、综述

利用matlab对GOCI数据产品进行可视化

二、数据下载

2.1 产品数据下载

网址：http://kosc.kiost.ac.kr/index.nm
本次实验所使用的数据是2019年2月1日的叶绿素浓度产品

2.2 经纬度数据下载

GOCI数据中不自带经纬度数据，需要单独下载。
下载网址还是GOCI官网，具体下载位置如图所示：

2.3 高精度海岸线数据

之前用的是matlab自带的海岸线数据，分辨率太低，效果不是很好，如图所示：

NOAA高精度海岸线数据下载网址：https://www.ngdc.noaa.gov/mgg/shorelines/gshhs.html
下载完安装到matlab安装路径中的m_map下
结果如图所示：

可以看出来分辨率有了很大程度的提升。

三 读取显示数据

GOCI数据产品都是He5格式的，所以在读取的时候需要注意一下数据的格式问题。具体请参考：

https://blog.csdn.net/wokaowokaowokao12345/article/details/108783849?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522163487210116780274162314%2522%252C%2522scm%2522%253A%252220140713.130102334.pc%255Fall.%2522%257D&request_id=163487210116780274162314&biz_id=0&utm_medium=distribute.pc_search_result.none-task-blog-2_allfirst_rank_ecpm_v1~rank_v31_ecpm-7-108783849.first_rank_v2_pc_rank_v29&utm_term=GOCI+matlab&spm=1018.2226.3001.4187

3.1 matlab代码

%%  读取经纬度
file1='COMS_GOCI_L2P_GA_20110524031644.LON_1.he5';
LON=h5read(file1,'/HDFEOS/GRIDS/Image Data/Data Fields/Longitude Image Pixel Values');
file2='COMS_GOCI_L2P_GA_20110524031644.LAT_1.he5';
LAT=h5read(file2,'/HDFEOS/GRIDS/Image Data/Data Fields/Latitude Image Pixel Values');
%%读取chla数据
file3='‪D:\学习课件\海遥\实验一\_home_goci_nfsdb_COMS_GOCI_L2_2019_02_01_L2_COMS_GOCI_L2A_GA_20190201001640.CHL.he5\COMS_GOCI_L2A_GA_20190201001640.CHL.he5'
chla=h5read(file3,'/HDFEOS/GRIDS/Image Data/Data Fields/CHL Image Pixel Values')

%%去除异常值（此处我先统计了整体的直方图，发现>10的部分占了不到5%，所以阈值设置为了10，具体情况具体分析）
chla(chla==-999)=NaN;
chla(chla>10)==NaN;
chla1=double(chla);

%%显示其他地图信息
 m_proj('Mercator','lon',[115.5 148.6],'lat',[21.5 46.9]); %投影：墨卡托
m_pcolor(LON,LAT,chla)
m_gshhs('fc','patch',[[0.8 0.8 0.8]],'edgecolor','r');%满分辨率的海岸线
title('2019年2月1日叶绿素浓度结果图','fontsize',15)  %标题
m_grid('box','fancy')
h=colorbar('h')
 set(get(h,'title'),'string','单位：mg/m^3'); %据GOCI官网发布 叶绿素浓度数据单位是mg/m^3
 m_northarrow(-123-4.5/60,49+19.5/60,1/60,'type',4,'aspect',1.5);
%指南针 
 m_ruler([0.1 0.5],0.9,4,'color','b',...
    'linewid',5,'tickdir','out','ticklen',0.01,'fontsize',8);
set(gcf,'position',[100 100 800 600]);

注释 统计直方图如图所示：

3.2 结果