[发明专利]一种基于地基云图的云团移动预测方法

摘要

本发明公开了一种基于地基云图的云团移动预测方法，对经过云团区域识别的图像依次进行云团匹配、移动预测和面积校正。本发明的云团移动预测方法，基于云团计算位移矢量和云团面积伸缩系数，对云团的运动趋势提供了一种更加准确的预测方法，为光伏电站发电功率超短期准确预测提供了前提条件。

主权项

一种基于地基云图的云团移动预测方法，其特征在于，包括下列步骤：步骤S1，云团匹配，包括以下步骤：步骤S11，首先对图像进行预处理，将当前帧和前一帧经过云团区域识别的图像二值化，使云团区域用1表示，背景区域用0表示；然后对二值化图像用边长为3个像素的方块进行开运算，以将较小的云块和细小的连接消除，根据8连通规则，对开运算后的图像进行连通分量的检测，每检测完一个连通分量就相应作一标记，每个连通分量里的像素被分配给一个唯一的整数，该整数的范围为从1到连通分量的总数,每个连通分量被认为是一个云团，同时可以获得该云团的几何特征，所述该云团的几何特征包括该云团的面积以及包含该云团的最小矩形的位置和大小；步骤S12，对于当前帧中的每一个云团在前一帧某一给定搜索区域内根据累积灰度差异值最小原则搜寻最佳匹配区域，所述累积灰度差异值的计算公式如下：C(Δx,Δy)=ΣΔx=-NNΣΔy=-NN|ft(xi,yj)-ft-1(xi+Δx,yj+Δy)|---(1)]]>其中，C(Δx,Δy)为累积灰度差异值的矩阵，Δx,Δy表示搜索框相对于待匹配云团中心的偏移；ft(xi,yj)为当前帧中，中心在(xi,yj)，搜索框大小为包含待匹配云团的最小矩形的子图区域，ft‑1(xi+Δx,yj+Δy)为前一帧中，中心在(xi+Δx,yj+Δy)，大小与ft(xi,yj)相等的子图区域；N表示横纵方向上的搜索范围；所述累积灰度差异值越小，表明两幅子图像之间的差异越小，可以近似认为这两幅子图像中的云团为同一云团在不同时刻的形态；对某一云团，令其相应的搜索框的大小为包含该云团的最小矩形，标记该云团的位置为该矩形的中心，并设置横纵方向最大的搜索半径为20个像素；首先按照搜索框的大小和位置对当前帧待识别云图像取相应的子图像即为包含该云团的最小矩形区域，然后在前一帧图像中，提取所有中心在由该云团中心和搜索半径限定的范围的搜索框区域，并计算它们与当前帧提取出的子图像的累积灰度值差异；找到满足累积灰度值差异为最小累积灰度值差异1.1倍的搜索框，计算它们的中心相对于该云团位置偏移的平均距离，以及该云团的移动矢量；步骤S13，对所述S11中标注的每一个云团重复所述S12过程，得到所有云团的移动矢量；步骤S2，面积校正：云团在运动的过程中，除了位置的偏移还有面积的膨胀和缩小，对云团进行面积校正是很有必要的，当前帧与前一帧中相匹配云团的面积比定义为面积伸缩系数Kt，Kt的表示如下：Kt=StSt-1---(2)]]>其中，St表示当前帧中待匹配云团的面积，St‑1为前一帧中匹配的云团的面积，若Kt>1时，表示云团面积膨胀，若Kt<1时，表示面积压缩，对标注的每一个云团，令当前帧和前一帧相匹配的云团面积相比，得到所有云团的面积校正系数，预测时，假定云团面积的膨胀或者压缩是均匀变化的，即云团面积的变化服从：St0+t＝Kt0*St0   (3)其中，St0为当前时刻云团的面积，St0+t为预测t时刻后云团的面积，Kt0为当前时刻求得该云团的面积伸缩系数，这样通过面积校正得到对云团面积变化更加准确的描述；步骤S3，移动预测：对上述面积校正之后的云团，根据获得的移动矢量，对它的运动进行预测，假定云团的运动为匀速直线运动，即云团质心的运动服从：(xt0+t,yt0+t)＝(xt0+dx*t,yt0+dy*t)   (4)其中，(xt0,yt0)为当前时刻云团质心，(xt0+t,yt0+t)为预测t时刻后云团质心，dx、dy一一对应地为该云团在横纵方向上单位时间内的位移。