matlab滴m_map学习心得3

乾坤毅志

已于 2024-03-29 17:03:20 修改

阅读量135

点赞数 1

分类专栏： m_map实战文章标签：算法

于 2023-09-26 17:08:38 首次发布

本文链接：https://blog.csdn.net/qq_58518586/article/details/133316632

版权

m_map实战专栏收录该内容

8 篇文章 1 订阅

订阅专栏

%% 读取文件信息
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')