人工增雨可播区域自动选择方法研究.pdfVIP

．236． 第十四届全国云降水物理和人工影响天气科学会议(下册) 人工增雨可播区域A动选择方法研究’ 黄毅梅1姚展予2 (1河南省人工影响天气办公室，郑州450003；2中国气象科学研究院，北京100081) 摘要：本文研究了一种方法，解决了综合利用雷达、卫星、数值模式产品和探空等各种不同类型、不 同分辨率的资料如何应用计算机软件技术自动选择可播区域。该方法使用分块法，把整个区域分解为若 干正方形小块，形成小块矩阵。通过处理不同分辨率资料，使之与一定小块之间建立相关，利用这些资料 综合判断小块可播性。然后把小块矩阵看作一个稀疏矩阵，采用正交链表存储方式，使用相邻搜索算法， 自动找出所有可播区域。本文详细描述了该方法的实现算法以及采用的相应数据结构和存储结构。该方 法实现了根据各种资料利用计算机软件技术自动选择可播区，并直接用于指挥飞机和高炮(火箭)人工增 雨，其对人工影响天气业务现代化建设具有重要的意义。 1 引言 在人工影响天气中，可播区域的选择十分重要[1’2]。过去，我们常用的方法是根据一些资料进行 综合分析，用人工的方法估计可播区域∞]，这样确定的范围往往比较粗略，并且当资料种类比较多、资 料分辨率差异较大时，人工分析的方法将无法准确判断可播区的位置，这直接影响人工增雨，特别是 高炮(火箭)人工增雨效益的充分发挥。针对这一问题，本文提出一种方法，解决如何利用雷达、卫星、 数值模式产品和探空等不同类型和不同分辨率的资料，应用计算机软件技术自动选择可播区域。 2使用单一资料选择可播区域的方法 为描述问题的方便，我们先假定使用单一的雷达回波资料来确定可播区域。假设有幅雷达回波 图，长为L，宽为D(L和D的单位均为像素)，该图每个像素点的回波强度可知，并且若强度大于 20dBZ为可播(当然可以附加其它各种条件来判断可播性，但这对本文方法研究不产生影响)，由于满 足条件的可播区可能分散地有许多块，并且每块大小和形状各异，如何应用计算机软件技术自动选择 它们? 2．1方法的设计思路 由于L，D的值较大，为了简化问题，可采用分块法进行处理。把整幅雷达回波图分解成边长 的可播条件，假设r为可播点占小块所有像素点的比率，若r大于60％，则该块可播，令PEx]Ey] 当作稀疏矩阵处理，可建立正交链表来存储小块矩阵[4]。查找可播区方法如下：遍历正交链表的 结点，先找出一个未被访问的结点，把其插入到一个单链表中，同时把正交链表中该结点标识为已 访问，然后从单链表的此结点开始，在正交链表中查找所有与该结点相邻的未被访问的结点，找出 后把其标识为已访问并插入到单链表的结尾，按上面步骤依次处理单链表中结点，直到处理完毕。 这样最后得到的单链表就是一块可播区。重复上面操作直到正交链表中所有结点均已访问，这样 · · s8一人工影响天气作业条件及指挥系统 237 就找出了所有可播区域。 2．2方法的实现 2．2．1数据结构设计 2．2．1．1稀疏矩阵数据结构的设计 本文采用正交链表存储稀疏矩阵。稀疏矩阵的每一行设置一个带表头结点的循环链表，每一列 也设置一个带表头结点的循环链表，行链表和列链表的表头可以共用一个结点。定义类Matrix— 是表头结点还是可播小块结点，down和right域用于存放稀疏矩阵列链表和行链表的最前端第一个 一个next的指针域，该域用于链结表头结点；若head一0(即为非表头结点)时，其取值为一结构体 剂量等信息，另外还包括一个next指针域，其用于链结同一可播区中其它小块结点。类MatrixNode 完整的定义见附录。 2．2．1．2线性链表数据结构的设计 个指针域，分别为link和right，其中right指向由同一可播区域所有小块组成的单链表的表首结点， 块总个数。线性链表的类的完整定义见附录。 2．2．2建立正交链表的算法 列链表和行链表的指针。在该算法中，首先建立所有的表头结点(开始时每个表头结点的next域指 向它本身)，然后，按列读取矩阵的每个可播小块P[i]EJ]，定义一个指针last，每开始读取新的一列 整last，让其指向新链入的结点，以便链结后续读入的该列的其它结点。由于表头结点H[i]的nex