新一代天气雷达产品故障识别及监测预警系统设计实现.pdf
《新一代天气雷达产品故障识别及监测预警系统设计实现.pdf》由会员分享,可在线阅读,更多相关《新一代天气雷达产品故障识别及监测预警系统设计实现.pdf(7页珍藏版)》请在咨信网上搜索。
1、Jun.2023METEORHNOLOGY3812023年6 月第51卷第3期Vol.51,No.3象技科新一代天气雷达产品故障识别及监测预警系统设计实现张东明1卢兴来1*杨豪2王晗1刘壮峰3(1浙江省大气探测技术保障中心,杭州310 0 0 0;2 浙江省宁波市气象网络与装备保障中心,宁波3150 0 0;3浙江省舟山市气象网络与装备保障中心,舟山316 0 0 0)摘要新一代天气雷达故障直接影响数据质量,异常产品的产生对短时临近预报预警带来了重大的影响,实现雷达异常产品实时监测具有一定现实意义。本系统采用WinForm3层架构十ORM模型架构,通过统计纹理分析、灰度图像特征值分析等图像处理
2、技术实现自动识别雷达异常产品。新一代天气雷达产品故障监测预警系统可以实时、高效、准确检测出异常雷达产品,满足强对流天气预警监测、天气雷达故障预警服务预判、天气雷达产品可用性服务等功能要求,从单部雷达产品数据中识别出异常回波,避免进一步影响多雷达组网拼图,对天气雷达组网产品的质量控制具有重要的意义。关键词雷达产品数据;图像处理技术;实时预警系统;异常回波检测中图分类号:P409D0l:10.19517/j.1671-6345.20220306文献标识码:A引言新一代天气雷达资料的数据质量对中小尺度灾害天气监测具有重要的影响,一旦异常气象回波产品进人模式对预报准确性影响很大且具有业务人员难以识别和
3、无法控制的特点。雷达故障造成极端数据上传,以异常数据为基础产生服务气象产品,未能及时监控并采取阻断措施,将导致严重不良后果。周红根等1,针对新一代天气雷达(CINRAD/SA),在运行过程中出现的雷达回波异常情况,从新一代天气雷达受到外界电磁波的干扰,计算机系统,接收机系统,天伺系统等方面,分析了造成雷达回波异常的故障成因。赵瑞金等2 对CINRAD/SA不同部位硬件故障的雷达基数据和回波特征进行了分析,给出不同类型雷达故障的数据特征和判定参数,指出通过对雷达扫描线位置信息逻辑检查,并结合回波图像的面积、强度、相似度检测等特征参数检测,提高硬件故障导致的错误数据识别。邵楠等3 论述基于图像处理
4、技术实现自动识别雷达异常产品,详细介绍系统实现模型及雷达数据异常自动识别算法。基于新一代天气雷达数据识别方法大多只做了个例研究,且一般针对硬件故障雷达分析或基数据研究,未开展业务化工作。本系统依托一种处理方法,该方法基于图像处理技术,进行雷达产品特征识别方法分析,包括产品数据识别、噪声特征提取、灰度图像分析4 。开发新一代天气雷达产品故障监测预警业务系统,系统自动识别异常产品,能更加高效地检测出异常雷达产品,及时采取措施阻断,提高雷达使用效率。无论是电磁环境、地物杂波干扰、超折射等外界原因造成的回波干扰,还是由设备软、硬件故障造成的产品图像异常,都会对业务正常开展造成影响5。本系统已在浙江全省
5、雷达站进行试运行,试运行期间运行平稳、可靠。1雷达产品特征分析方法当雷达出现故障时,对雷达数据主要产生如下影响:数据完整性出现异常或无数据。主要是因为雷达硬件系统出现故障,常常导致停机,对探测数据完整性会产生影响6 1。数据位置信息错误。伺http:/气象科技浙江省重点气象科技项目(项目编号:2 0 2 0 ZD12)资助作者简介:张东明,男,1990 年生,工程师主要从事气象装备保障与气象系统软件开发工作,Email:32 90 7 350 8 q q.c o m收稿日期:2 0 2 2 年7 月2 0 日;定稿日期:2 0 2 3年3月2 8 日*通信作者,Email:382第51卷象技科
6、服系统故障常常影响探测数据的方位和仰角信息。当雷达探测数据的位置信息错误时,会导致相邻体扫回波的位置、范围突变、出现丝状回波、V形等图形7。数据强度出现异常。如连续的异常高的回波强度、无数据、强度无变化或呈规律性变化等,回波上表现为相邻体扫强度、范围差别很大,出现环状、饼状回波等8 1.1产品数据完整性信息识别多普勒天气雷达数据的完整性检查包括仰角层数量、观测参量类型及数量和径向扫描数量等,数据完整性的检查基于解析后的产品数据文件信息进行判断9。首先,多普勒天气雷达有多种工作模式,各种模式的扫描层数差异较大,根据产品文件中扫描配置块的扫描层数(CutNumber)信息判断仰角层数量是否异常10
7、 ;然后,观测参量类型、数量和径向数据的检查通过产品数据文件中的径向数据块信息判断其异常情况,如径向数据头记录了径向数(RadialNumber)、方位角(Azimuth)、仰角(Elevation)以及数据类别数量(MomentNumber)等信息11。最后,根据径向数据的连续性以及各仰角层径向数据方位覆盖情况对数据完整性进行综合判定。当仰角层数为NA,每层扫描线数为NB,数据完整性异常值为:NA11或NB 2(2)1.3产品数据强度信息识别雷达故障导致数据错误会引起回波形态、范围和强度的变化,可能会出现饼状、环状等回波图14 而气象回波在演变过程中、移动速度、高度、强度等均有一定合理限度,
8、在时间和空间上发展具有一定的延续性15。因此通过相邻体扫的回波特征参数的比较可以识别数据错误。需要与相邻的正常数据进行图像特征相关才能识别出错误数据16 1.4产品数据异常信息识别基于数字图像处理对新一代天气雷达产品进行图像特征提取,形成各类判断指标,判定回波数据质量情况,推出报警或提醒等信息,部分实现雷达回波智能分析功能17 。研究上述极值出现情况,在从数据中提取相对规律,分析图像信息,设定阈值,推出报警信息,触发人工干预,避免这类异常数据连续多次传输,具体通过图像特征值提取对比分析方法实现:首先将雷达回波图每条径向上的数据解析出来,对一条径向上出现非零同值的现象进行特征值提取分析,并作出错
9、误报警。针对雷达回波图相邻点位、区块点位进行阈值对比判断,将出现严重跳变的现象进行错误预警,1.5噪声特征提取通过使用图像区域增长的方法求取所有区域,通过大量的图片实验总结出一个阈值,判定小于该阈值为噪声点,计算所有噪声点占据的像素数占总像素的比例,再经过典型图片试验对该比例找到一个优质的阈值,如果大于比例阈值则认为是噪声,否则不是噪声18 1.6灰度图像分析统计纹理分析寻找刻划纹理的数字特征,用这些特征或同时结合其他非纹理特征对图像中的区域进行分类19。对图像灰度空间分布模式提取和分析,为了能更直观地以共生矩阵描述纹理状况,从共生矩阵导出一些反映矩阵状况的参数,能更好地兼顾图像宏观性质与细微
10、结构两个方面。灰度共生矩阵中的元素(i,j)的值表示了在图像中其中一个像素的灰度值为i,另一个像素的灰度值为i。灰度共生矩阵被定义为从灰度为i的像素点出发,离开某个固定位置(相隔距离为d,方位为0/4 5/90/135)的点上灰度值为其概率,即,所有估计的值可以表示成一个矩阵的形式,以此被称为灰度共生矩阵2 0 ,如图1所示。对于纹理变化缓慢的图像,其灰度共生矩阵对角线上的数值较大;而对于纹理变化较快的图像,其灰度共生矩阵对角线上的数值较小,对角线两侧的值较大。1.6.1对比度反应了图像的清晰和纹理沟纹深浅的程度。纹理的沟纹深,其对比度越大,效果越清晰,反之,对比值小,则沟纹浅,效果模糊2 1
11、。灰度差即对比度大的像素对越多,这个值越大。灰度共生矩阵中原理383第3 期张东明等:新故障识别及监测预警系统设计实现42107厂412401220M(o)M(45)10600241LO012JLO01O1170110222L2233J图像d=1602072130704201210M(90)M(135)22223102LO2LO02图1灰度共生矩阵示意对角线的元素值越大,对比度越大。CoN=Z(i-j)P(i,j)(3)其中,对于整个画面,统计出每一种(i,j)值出现的次数,然后排列成一个方阵,在用(i,j)出现的总次数将它们归一化为出现的概率Pi,j),P(i,j)是指归一化后的灰度共生矩阵
12、。1.6.2能量能量变换反映了图像灰度分布均勾程度和纹理粗细度2 2 。若灰度共生矩阵的元素值相近,则能量较小,表示纹理细致;若其中一些值大,而其他值小,则能量值较大。能量值大表明一种较均一和规则变化的纹理模式。图像越均匀,其值越大。1.6.3熵图像包含信息量的随机性度量。当共生矩阵中所有值均相等或者像素值表现出最大的随机性时,熵最大;因此熵值表明了图像灰度分布的复杂程度,熵值越大,图像越复杂2 3ENT=-ZZp(i,j)lgp(i,j)(5)1.6.4逆方差逆方差反映了图像纹理局部变化的大小,若图像纹理的不同区域间较均匀,变化缓慢,逆方差会较大,反之较小2 4 N-1N-1p(i,j/d,
13、0)H=(6)1+(i-j)2i=0j-0其中,p(i,i ld,)表示在给定空间距离d和方向0时,灰度以i为起始点(行),出现灰度级(列)的AsM=ZZP?(i,j)(4)Z概率。1.6.5相关性用来度量图像的灰度级在行或列方向上的相似程度,当矩阵元素值均匀相等时,相关值就大;相反,如果矩阵像素值相差很大则相关值小2 5Z,2,P(i,j)-sR=(7)其中u、u y 是均值x、y是标准差。2系统设计实现2.1设计思路分析本系统总体采用基于WinForm3层架构十ORM模型的快速开发框架,使用图形、表格、声音、微信推送等方式实现系统展示、分析以及报警的功能,满足雷达站业务需求2 6 。系统主
14、要包括系统设置,产品分析、图像特征计算、报警推送等。系统整体可以分为5个子模块,分别为系统后设置、产品及时性监控、雷达图像特征值分析、产品异常监控、历史数据再分析。系统整体流程图如图2 所示。启动监控是否为网络映射路径设置映射目录设置监控目录开启监控子流程产品及时产品质量性检查检查空值产品检查收到产品绘制产品图计算回波强度极值产品检查灰度化图像份布、均值等绘制分布和趋计算图像特征根据特征值判电磁干扰检查势图值断天气状态根据检查情况根据强度分布绘制特征值趋根据判断结果报警情况报警势提示根据检查结果存储数据报警结束当前监控图2雷达产品预警系统工作流程384象第51卷科技2.2系统创新点产品质量检查
15、丰富:雷达产品质量检查流程,主要包含空值产品检查、极值产品检查、单层仰角上电磁干扰情况检查、多层仰角组合电磁干扰情况检查、绘制产品图像特征、计算产品图像特质值、雷达回波情况检查。本系统基于雷达产品数据绘制成图像,并对产品图像进行灰度二值化,计算灰度共生矩阵的能量、熵、对比度、均匀性、相关性、方差、和平均、和方差、和熵、差方差、差平均、差熵、相关信息测度以及最大相关系数,最后根据特征值判断天气状态以及雷达运行状态2.3功能模块设计采用分解的方法,将系统设计成由相对独立、功能单一的模块组成的结构,它是以系统的逻辑功能和数据流关系为基础,根据数据流程图,借助于一套标准的设计准则和图表工具,通过“自上
16、而下”和自下而上”的反复,把系统逐层划分为多个大小适当、功能明确、具有一定独立性且容易实现的模块。雷达产品分析报警系统功能结构图如图3所示。雷达产品分析报警系统系统设置产品文件监控产品文件解析产品质量检查产品图像特征分析历史数据查询报警日志记录监控相关设置报警阅值设置监控目录挂接检查及时性检查公共数据块产品头块产品数据块空值产品检查极值产品检查电磁干扰检查图像能量图像对比度图像烟数据、图像查询数据统计数据、图像导出文件记录数据库记录图3雷达产品分析报警系统功能结构(1)产品文件监控:雷达产品主要通过3种方式进行存储:直接存储在本地计算机上;雷达产品通过ftp方式存储在指定存储中,在指定计算机上
17、对映射目录进行监控;雷达产品直接存储到存储系统的共享目录中,在指定计算机中设置映射目录,对映射目录进行监控。(2)产品文件解析:按照天气雷达产品数据标准格式对新生成的产品文件进行解析,产品文件是按照二进制进行存储的。产品数据可分为公共数据块、产品头块和产品数据块3个部分,其中:公共数据块用于提供数据的站点信息、任务配置等公共信息。产品头块用于提供产品最基本的信息。数据格式中定义了数据类型和每个部分固定长度,因此通过多线程方式批量读取数据,达到快速解析一个产品数据。(3)及时性检查:定时检查与最后一次获得产品文件时间或是系统刚启动的时间对比,判断是否有超过指定时间未生成产品文件。(4)电磁干扰情
18、况检查:对单个仰角反射率产品的电磁干扰情况进行分析判断,当某一时刻多个仰角的反射率产品均获得后,对多层的反射率产品进行综合分析判断。(5)无回波或异常回波判断:对照PUP中的PPI产品,对于数据区编码值小于5代表无回波区域。对有回波的情况通过(编码值一6 6)/2 的计算方式,对超过所有数据或某方位角内的所有数据均超过55dB,则判断为异常回波。(6)绘制产品图像:根据产品文件中的数据块绘制极坐标方式的产品图像和笛卡尔坐标方式的热力图。(7)图像特征计算:对真彩图的产品图像进行灰度二值化,接着对图像的灰度共生矩阵进行统计分析2 7 ,得到14 种图像特征的统计量。图像特征值判断:根据图像特征值
19、和系统的阈值参数设置,进行分析判断。(8)检查报警模块:针对产品数据质量和图像特征值结果采用声光报警,研究标准输出控制器手机短信接口,将监控软件接人,推送告警信息,如若标准输出控制器不正常则将告警信息推送给云mas平台。雷达产品分析报警系统报警模块流程如图4 所示385第3 期张东明等:新故障识别及监测预警系统设计实现FTP方式本项目报警状态文件标准输出控制器短信猫手机服务器WRSOC图4雷达产品分析报警系统报警模块流程2.4系统功能实现系统主要模块由产品的基本信息、反射率产品图和反射率热力图、图像特征识别信息报警提示信息、多层反射率干扰综合判别报警信息、波强度分布柱状图、回波强度均值趋势柱状
20、图和图像特征值趋势图、操作控制功能区等综合模块组成。系统实时监控如图5所示。满达广品胜控系线主英时产品图最就产品的基本信息富达涯热力圈站名,NInoBo.Z9574任务,VCP21D产品代码:19仰0.48分样率,250数据开始时间:2 0 2 1/7/2 18.12:0 0能量0.934294536064相关性103384.594455数据结束时间,2 0 2 1/7/2 18 12:32对比度257.043380363C方差,3.73940272983产品生成时间,2 0 2 1/7/2 18 15:17均匀性,0.9824505483840.220609749023回波强度分布西波强度均
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 新一代 天气 雷达 产品 故障 识别 监测 预警系统 设计 实现
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。