- 积分
- 25
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2019-3-9
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
求助大家
ASCAT散射计25Km风场数据可视化,老是有莫名的条块,拉扯,想问一下是什么错误?
matlab代码:
clear
clc
ascat_path='H:\ASCAT\2018\243\';
ascat_dir=dir([ascat_path,'ascat_20180831_214500_metopa_61576_eps_o_250_2401_ovw.l2.nc']);
num_file=size(ascat_dir,1);
figure;
m_proj('stereographic','lat',90,'long',0,'radius',30); %stereographic mercator
for num=1:1 %num_file
file_name=ascat_dir(num).name;
file_name_path=fullfile(ascat_path,ascat_dir(num).name);
lon_ascat=ncread(file_name_path,'lon');
lon_ascat_t=wrapTo180(lon_ascat);
lat_ascat=ncread(file_name_path,'lat');
% lon_ascat=lon_ascat*10^-5;
% lat_ascat=lat_ascat*10^-5;
%time_ascat=ncread(file_name_path,'time');
wind_speed_ascat=ncread(file_name_path,'model_speed'); %
ascat_wvc_quality_flag=ncread(file_name_path,'wvc_quality_flag');
% wind_speed_ascat(ascat_wvc_quality_flag==64)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==128)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==256)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==512)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==1024)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==2048)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==4096)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==8192)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==16384)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==32768)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==65536)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==131072)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==262144)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==524288)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==1048576)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==2097152)=nan;
% wind_speed_ascat(ascat_wvc_quality_flag==4194304)=nan;
wind_speed_ascat(wind_speed_ascat==-32767)=nan;
% wind_speed_ascat=wind_speed_ascat*0.01;
wind_dir_ascat=ncread(file_name_path,'wind_dir');
hold on;
m_pcolor(lon_ascat_t,lat_ascat,wind_speed_ascat);
end
shading('flat');
colormap(jet)
h=colorbar('location','eastoutside');
set(h,'FontSize',13)
ylabel(h,'Wind Speed (m s^-^1)','fontsize',13);
m_coast('patch',[.7 .7 .7],'edgecolor','k');
% m_pcolor(lon_ascat,lat_ascat,wind_speed_ascat);
m_grid('linestyle','none','box','on','fontsize',13,'linest','-');
|
-
北极风场多轨数据
-
单轨数据
|