- 积分
 - 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','-'); 
 
 |   
- 
北极风场多轨数据 
 
 
 
- 
单轨数据 
 
 
 
 
 
 
 
 |