matlab滴m_map学习心得3

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

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值