%% 读取文件信息
file='D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_443_4km.nc';
ncdisp(file);
nums = length('AQUA_MODIS.20220901_20220930.L3m.MO.RRS.Rrs_443.4km.nc');
lon = ncread(file,'lon');
lat = ncread(file,'lat');
lon_real = lon(find(lon>= 100 & lon<= 180));
lat_real = lat(find(lat >= 0 & lat <= 65));
[LON,LAT] = meshgrid(lon_real,lat_real);
%% 读取遥感反射率数据
data443 = ncread('D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_443_4km.nc','Rrs_443');
data488 = ncread('D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_488_4km.nc','Rrs_488');
data547 = ncread('D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_547_4km.nc','Rrs_547');
max1=max(data443,data488);
x = log10(max1./data547);
x_real=x(find(lon>= 100 & lon<= 180),find(lat >= 0 & lat <= 65));
data_chla = 10.^(0.2424-2.7423.*x_real+1.8017.*(x_real.^2)+0.0015.*(x_real.^3)-1.228.*(x_real.^4));
%% 筛选不好的数据
for i=1:1920
for j=1:1560
if ~isnan(data_chla(i,j)) || data_chla(i,j)>0.001
plot_chla(i,j) = log10(data_chla(i,j));
else
plot_chla(i,j) = NaN;
end
end
end
%% 遥感数据可视化
figure(1)
m_proj('miller','lon',[100 180],'lat',[0 65]);
m_pcolor(LON,LAT,plot_chla');
m_gshhs_h('color','k');
m_coast('patch',[.7 .7 .7]);
m_grid('fontsize',12,'fontname','Times New Roman');
colormap(jet(256));
ax=m_contfbar(.97,[.15 .85],[log10(0.01) log10(20)],[log10(0.01):0.001:log10(20)],'edgecolor','none','endpiece','no','ticklength',[.01 0.05]);
set(ax,'ytick',[log10(0.01) log10(0.1) log10(1) log10(10)],'yticklabel',{'0.01','0.1','1','10'}) % chla
title(ax,'Chl-a','color','k','fontweight','bold','fontname','Times New Roman')
%% 读取文件信息
file='D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_443_4km.nc';
ncdisp(file);
nums = length('AQUA_MODIS.20220901_20220930.L3m.MO.RRS.Rrs_443.4km.nc');
lon = ncread(file,'lon');
lat = ncread(file,'lat');
lon_real = lon(find(lon>= -180 & lon<= 180));
lat_real = lat(find(lat >= -90 & lat <= 90));
[LON,LAT] = meshgrid(lon_real,lat_real);
%% 读取遥感反射率数据
data443 = ncread('D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_443_4km.nc','Rrs_443');
data488 = ncread('D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_488_4km.nc','Rrs_488');
data547 = ncread('D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_547_4km.nc','Rrs_547');
max1=max(data443,data488);
x = log10(max1./data547);
x_real=x(find(lon>= -180 & lon<= 180),find(lat >= -90 & lat <= 90));
data_chla = 10.^(0.2424-2.7423.*x_real+1.8017.*(x_real.^2)+0.0015.*(x_real.^3)-1.228.*(x_real.^4));
plot_chla=log10(data_chla);
%% 遥感数据可视化
figure(1)
m_proj('miller','lon',[-180 180],'lat',[-90 90]);
m_pcolor(LON,LAT,plot_chla');
m_grid('fontsize',12,'fontname','Times New Roman');
colormap(jet(256));
ax=m_contfbar(1.05,[.15 .85],[log10(0.01) log10(20)],[log10(0.01):0.001:log10(20)],'edgecolor','none','endpiece','no','ticklength',[.01 0.05]);
set(ax,'ytick',[log10(0.01) log10(0.1) log10(1) log10(10)],'yticklabel',{'0.01','0.1','1','10'}) % chla
title(ax,'Chl-a','color','k','fontweight','bold','fontname','Times New Roman')
%% 读取文件信息
file='D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_443_4km.nc';
lon = ncread(file,'lon');
lat = ncread(file,'lat');
lon_real = lon(find(lon>= -180 & lon<= 180));
lon_real1= lon_real+180;
lat_real = lat(find(lat >= -90 & lat <= 90));
[LON,LAT] = meshgrid(lon_real1,lat_real);
%% 读取遥感反射率数据
data443 = ncread('D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_443_4km.nc','Rrs_443');
data488 = ncread('D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_488_4km.nc','Rrs_488');
data547 = ncread('D:\遥感实验\3 专题三 海洋水色遥感\A20210012021365.L3m_YR_RRS_Rrs_547_4km.nc','Rrs_547');
max1=max(data443,data488);
x = log10(max1./data547);
x_real=x(find(lon>= -180 & lon<= 180),find(lat >= -90 & lat <= 90));
data_chla = 10.^(0.2424-2.7423.*x_real+1.8017.*(x_real.^2)+0.0015.*(x_real.^3)-1.228.*(x_real.^4));
chlareal=log10(data_chla);
chlareal1 = circshift(chlareal,4320);
%% 遥感数据可视化
figure(1)
m_proj('robinson','lon',[0 360],'lat',[-90 90]);
m_pcolor(LON,LAT,chlareal1');
m_grid('fontsize',12,'fontname','Times New Roman');
colormap(jet(256));
color=m_contfbar([.2 .8],.005,[log10(0.01) log10(20)],[log10(0.01):0.001:log10(20)],'edgecolor','none','endpiece','no','ticklength',[.01 0.05]);
set(color,'ytick',[-2 -1.5 -1 -0.5 0 0.5 1 1.5],'yticklabel',{'-2','-1.5','-1','-0.5','0','0.5','1'})
title(color,'Chl-a','color','k','fontweight','bold','fontname','Times New Roman')