爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 9764|回复: 4

[讨论] ASCAT散射计数据风场资料可视化问题

[复制链接]

新浪微博达人勋

发表于 2020-7-16 22:26:09 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

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','-');

北极风场多轨数据

北极风场多轨数据

单轨数据

单轨数据
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-9-28 16:05:41 | 显示全部楼层
画图方法的问题;你应该把这些数据看成点数据,先把plot上点,然后对这些点数据的值赋上颜色。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-9-28 16:21:47 | 显示全部楼层
还有一种方式是把数据重采样成标准的网格格式,使用m_pcolor就不会出现这种问题了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-9-28 16:35:08 | 显示全部楼层
或者可以试一试 m_grid.
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-12-22 23:39:36 | 显示全部楼层
请问这个数据是怎么下载的呀  我只能找到图片
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表