爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 10749|回复: 26

[程序设计] 反射率基数据去噪程序

[复制链接]

新浪微博达人勋

发表于 2013-4-17 13:52:57 | 显示全部楼层 |阅读模式

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

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

x
论文中用到雷达反射率资料,需要进行去噪预处理,不知道怎么编程序,求高手解答

评分

参与人数 1金钱 +2 收起 理由
又是那隻貓 + 2 很给力!

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-5-18 11:50:54 | 显示全部楼层

function ReadDataHead(fid)

global state;%径向数据状态
global azimuth%方位角
global elevation%仰角
global nstep%体扫内的仰角数
global firstdistanceR%反射率数据的第一个距离库的实际距离
global firstdistanceV%多普勒数据的第一个距离库的实际距离
global gatewidR%反射率数据的距离库长
global gatewidV%多普勒的距离库长
global gatenumR%反射率的距离库数
global gatenumV%多普勒的距离库数
global Rpoint%反射率数据指针(偏离雷达数据信息头的字节数),表示第一个反射率数据的位置
global Vpoint%速度数据指针(偏离雷达数据信息头的字节数),表示第一个速度数据的位置
global Wpoint%谱宽数据指针(偏离雷达数据信息头的字节数),表示第一个谱宽数据的位置
global mode%体扫(VCP)模式
global resolution;%多普勒速度分辨率
global Nyquistvelo;%Nyquist速度(最大不模糊速度)

head=fread(fid,128,'uchar');%读文件头的重要信息
head=head';
state=head(41)+head(42)*256;
azimuth=(head(37)+head(38)*256)/8.*(180./4096);
elevation=(head(43)+head(44)*256)/8.*(180./4096);
nstep=head(45)+head(46)*256;
firstdistanceR=head(47)+head(48)*256;
firstdistanceV=head(49)+head(50)*256;
gatewidR=head(51)+head(52)*256;
gatewidV=head(53)+head(54)*256;
gatenumR=head(55)+head(56)*256;
gatenumV=head(57)+head(58)*256;
Rpoint=head(65)+head(66)*256;
Vpoint=head(67)+head(68)*256;
Wpoint=head(69)+head(70)*256;
mode=head(73)+head(74)*256;
resolution=head(71)+head(72)*256;%2表示表示0.5米/秒,4表示表示1.0米/秒
Nyquistvelo=(head(89)+head(90)*256)/100.;%最大不模糊速度
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2013-4-18 14:38:34 | 显示全部楼层
嗯,也是我在处理基数据的时候问同学的,去噪要看是指孤立回波还是晴空回波,如果只是孤立回波,用5*5的模板判断一下,就能去除噪声了。晴空回波不太好去,要做回波识别或者简单点的设个阈值,比如把小于多少dBz的都识别为晴空回波。地物回波就必须做回波识别了。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-4-19 06:53:04 | 显示全部楼层

回帖也是一种美德
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-4-19 13:04:08 | 显示全部楼层
爱哭贪吃睡不着 发表于 2013-4-18 14:38
嗯,也是我在处理基数据的时候问同学的,去噪要看是指孤立回波还是晴空回波,如果只是孤立回波,用5*5的模板 ...

谢啦  我那个已经编好了,由于第一次用matlab编程,程序中间出了点问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-4-19 13:04:28 | 显示全部楼层
wagxx1205 发表于 2013-4-19 06:53
回帖也是一种美德

八错。。。。。。。。。。。。。。。。。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-4-19 22:04:39 | 显示全部楼层
学习学习。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-5-13 23:34:20 | 显示全部楼层

回帖奖励 +2 金钱

请问能不能将matlab读取雷达基数据的程序共享一下呢 ,以及一些相关的程序,非常感谢,目前正在学习中
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-5-18 11:59:53 | 显示全部楼层
searching· 发表于 2013-5-18 11:50
function ReadDataHead(fid)

global state;%径向数据状态

非常感谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-6-2 15:14:02 | 显示全部楼层
回帖也是一种美德 请问有没有雷达估测降水方面的同行 介绍一下具体获取雷达反射率的步骤 以及Z-I关系中Z获取的方法
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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