一种基于视觉边缘检测的实验室空间波浪实时测量方法与流程

本发明属于海洋工程领域，尤其是涉及一种基于视觉边缘检测的实验室空间波浪实时测量方法。

背景技术：

防波堤的破坏通常用护面块体的移动、断裂等条件来衡量的，因此，对实验室块体运动检测装置的研究是十分有必要的。常用的块体运动检测方法有目视检查、特写摄影以及从固定位置摄影测量等。目视检查对于检查特定的损坏是有用的，每段防波堤的破碎块体的数量和移动块体的数量可以进行目视检查，但这种方法耗时，不适合检测整个防波堤；特写摄影只是记录视觉检测结果，对于检查局部损伤的详细状况是有用的；从固定位置进行摄影测量，以制作覆盖防波堤整个水面状况的重叠照片，是防波堤检测中最有用和最具成本效益的方法；通过试验前后采集到的图像来分析防波堤的状态是一个可行的方法，然而，对图像中目标的检测和提取对计算要求很高，因此不能用于计算能力有限的系统。但上述方法均未对块体的运动和位移只能进行定性地分析，均未实现定量的测量和实时测量。

技术实现要素：

有鉴于此，本发明提出一种基于视觉边缘检测的实验室空间波浪实时测量方法以解决现有检测方法无法对块体的运动和位移只能进行定性地分析，实现定量的测量和实时测量的问题。

为达到上述目的，本发明的技术方案是这样实现的：

一种基于视觉边缘检测的实验室空间波浪实时测量方法，包括以下步骤：

s1：对摄像机进行标定；

s2：通过摄像机拍摄防波堤表面块体的数字图像，将数字图像发送到基于现场可编程逻辑门阵列的数字信号处理系统；

s3：数字信号处理系统对图像中的块体进行特征点检测；

s4：特征点检测完成后进行坐标转换；

s5：比较试验前后块体的特征点位置的变化，对试验前后的两张图像的坐标作差得到特征点变化值，从而得到块体的位移量；

s6：将位移计算结果显示在数据处理器上。

进一步的，步骤s3所述的特征点检测包括以下步骤：

s31：使用hessian矩阵生成数字图像的边缘点，图像中每个边缘点都设置一个hessian矩阵

s32：利用数字图像构建高斯金字塔；

s33：对hessian矩阵处理过的每个像素点与其三维邻域内点的大小进行比较，如果该像素点是邻域内像素点的最大值或最小值，则保留下来，作为初步的特征点；

s34：统计特征点邻域内的harr小波特征；

s35：根据harr小波特征生成特征点描述子；

s36：通过计算两个特征点间的距离判断两个特征点的匹配度，两个特征点的距离越短匹配度越高；

s37：对与每个块体对应的特征点进行筛选，保留匹配度最高的特征点来表示该块体，从而完成特征点的检测。

进一步的，步骤s31利用的使用hessian矩阵生成数字图像的边缘点方法如下：

其中，f(x,y)为图像的像素值；

hessian矩阵的判别式为：

当hessian矩阵的判别式取得局部极大值时，可判定当前点是比周围邻域内其它点更亮或更暗的点，则该点为特征点的位置。

进一步的，步骤s32利用的构建高斯金字塔，图像的大小不变，只改变高斯模糊模板的尺寸和尺度大小。

进一步的，步骤s34利用的统计像素点邻域内的harr小波特征的步骤如下：

s341：以特征点为中心，计算邻域内所有点在水平和垂直方向的haar小波响应总和；

s342：给haar小波响应值赋高斯权重系数，使得靠近特征点的响应贡献比远离特征点的响应贡献大；

s343：给邻域内haar小波响应相加以形成新的矢量；

s344：遍历整个区域，选择最长矢量的方向为该特征点的主方向。

进一步的，步骤s35利用的生成特征点描述子是在特征点周围沿着特征点的主方向取4*4的矩形区域块，统计每个子区域像素的水平方向和垂直方向的harr小波特征，该haar小波特征包括水平方向值之和、水平方向绝对值之和、垂直方向之和、垂直方向绝对值之和。

进一步的，步骤s4利用的坐标转换方法如下：

建立o0uv像素坐标系，o0为像素坐标系的原点，(u0,v0)为图像平面中心的像素坐标，建立o1xy为物理坐标系，o1为物理坐标系的原点；

dx为每个像素在u轴方向上的物理尺寸，dy为每个像素在v轴方向上的物理尺寸。

进一步的，步骤s1利用的对摄像机进行标定采用张正友标定法。

相对于现有技术，本发明具有以下有益效果：

本发明提出的一种基于特征点识别的块体运动实时检测装置，采用远程控制的数码摄像机发送灰度像素和基于现场可编程逻辑门阵列的数字信号处理系统中实现的块体运动测量算法对灰度像素进行分析来测量实验室块体的运动。块体运动测量算法构建了高斯金字塔，不同组间图像的尺寸都是一致的，不同的是不同组间使用的盒式滤波器的模板尺寸逐渐增大，同一组不同层图像使用相同尺寸的滤波器，但是滤波器的尺度空间因子逐渐增大，以此来实现尺度不变性；所述块体运动测量算法在生成特征描述子前将图片旋转到主方向上，保证了一个特征点生成描述子用的是同一块图像的信息，实现了旋转不变性，对光照、阴影及焦点丢失时发生的尺度变化有更好的效果。最后，在现场可编程逻辑门阵列中对块体运动测量算法进行了硬件实现，大大降低了系统对计算能力的要求，能够实时计算试验前后图像中块体的位移，响应速度快。此外，本发明所述的实验室块体运动测量装置还具有造价低、安装方便、不易损坏、维修成本低的优点。

附图说明

构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1为本发明实施例所述的一种基于视觉边缘检测的实验室空间波浪实时测量方法流程示意图；

图2为本发明实施例所述的防波堤示意图。

附图标记说明：

1-数据处理器；2-数字信号处理系统；3-数码摄像机；4-图像采集控制器；5-水面；6-防波堤；7-块体。

具体实施方式

需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。

在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”等仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”等的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。

在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以通过具体情况理解上述术语在本发明中的具体含义。

下面将参考附图并结合实施例来详细说明本发明。

如图1所示，一种基于特征点识别的块体7运动实时检测方法，包括以下步骤：

s1：对摄像机进行标定；

s2：通过摄像机拍摄防波堤6表面块体7的数字图像，将数字图像发送到基于现场可编程逻辑门阵列的数字信号处理系统；

s3：数字信号处理系统对图像中的块体7进行特征点检测；

s4：特征点检测完成后进行坐标转换；

s5：比较试验前后块体7的特征点位置的变化，对试验前后的两张图像的坐标作差得到特征点变化值，从而得到块体7的位移量；

s6：将位移计算结果显示在数据处理器上。

数码摄像机设置于水槽上方且垂直于模型护面，所述图像采集控制器直接和数码摄像机相连来控制连续拍摄图像，所述数字信号处理系统将采集的图像进行处理分析，并将分析结果传递给数据处理器。

数字信号处理系统用于对图像进行实时处理，数码摄像机被配置为可以发送灰度像素，该摄像机使用但不限于usb总线与数字信号处理通信，图像采集控制器在不接触摄像机即可远程控制摄像机拍摄连续的图像。

所述现场可编程逻辑门阵列模块包括：视频采集模块、图像存储模块、数据处理模块和图像显示模块。

如图1所示，步骤s3所述的特征点检测包括以下步骤：

s31：使用hessian矩阵生成数字图像的边缘点，图像中每个边缘点都设置一个hessian矩阵

s32：利用数字图像构建高斯金字塔；

s33：对hessian矩阵处理过的每个像素点与其三维邻域内点的大小进行比较，如果该像素点是邻域内像素点的最大值或最小值，则保留下来，作为初步的特征点；

s34：统计特征点邻域内的harr小波特征；

s35：根据harr小波特征生成特征点描述子；

s36：通过计算两个特征点间的距离判断两个特征点的匹配度，两个特征点的距离越短匹配度越高；

s37：对与每个块体7对应的特征点进行筛选，保留匹配度最高的特征点来表示该块体7，从而完成特征点的检测。

步骤s31利用的使用hessian矩阵生成数字图像的边缘点方法如下：

其中，f(x,y)为图像的像素值；

hessian矩阵的判别式为：

当hessian矩阵的判别式取得局部极大值时，可判定当前点是比周围邻域内其它点更亮或更暗的点，则该点为特征点的位置。

步骤s32利用的构建高斯金字塔，图像的大小不变，只改变高斯模糊模板的尺寸和尺度大小。

步骤s34利用的统计像素点邻域内的harr小波特征的步骤如下：

s341：以特征点为中心，计算邻域内所有点在水平和垂直方向的haar小波响应总和；

s342：给haar小波响应值赋高斯权重系数，使得靠近特征点的响应贡献比远离特征点的响应贡献大；

s343：给邻域内haar小波响应相加以形成新的矢量；

s344：遍历整个区域，选择最长矢量的方向为该特征点的主方向。

步骤s35利用的生成特征点描述子是在特征点周围沿着特征点的主方向取4*4的矩形区域块，统计每个子区域像素的水平方向和垂直方向的harr小波特征，该haar小波特征包括水平方向值之和、水平方向绝对值之和、垂直方向之和、垂直方向绝对值之和。

步骤s4利用的坐标转换方法如下：

建立o0uv像素坐标系，o0为像素坐标系的原点，(u0,v0)为图像平面中心的像素坐标，建立o1xy为物理坐标系，o1为物理坐标系的原点；

dx为每个像素在u轴方向上的物理尺寸，dy为每个像素在v轴方向上的物理尺寸。

进一步的，步骤s1利用的对摄像机进行标定采用张正友标定法。

如图2所示，当水面5完全平静之后使用数码摄像机3在图像采集控制器4的远程控制下分别对试验前后的防波堤6护面进行数字图像的拍摄，然后将灰度像素传输到基于现场可编程逻辑门阵列的数字信号处理系统2，在数字信号处理系统中对图像中的块体7进行特征点检测、坐标转换以及块体7位移的计算，最后将结果显示在数据处理器1上。

本发明所述的一种基于特征点识别的实验室块体运动实时检测装置，采用远程控制的数码摄像机发送灰度像素和基于现场可编程逻辑门阵列的数字信号处理系统中实现的块体运动测量算法对灰度像素进行分析来测量实验室块体的运动。所述的块体运动测量算法构建了高斯金字塔，不同组间图像的尺寸都是一致的，不同的是不同组间使用的盒式滤波器的模板尺寸逐渐增大，同一组不同层图像使用相同尺寸的滤波器，但是滤波器的尺度空间因子逐渐增大，以此来实现尺度不变性；所述块体运动测量算法在生成特征描述子前将图片旋转到主方向上，保证了一个特征点生成描述子用的是同一块图像的信息，实现了旋转不变性，对光照、阴影及焦点丢失时发生的尺度变化有更好的效果。最后，在现场可编程逻辑门阵列中对块体运动测量算法进行了硬件实现，大大降低了系统对计算能力的要求，能够实时计算试验前后图像中块体的位移，响应速度快。此外，本发明所述的实验室块体运动测量装置还具有造价低、安装方便、不易损坏、维修成本低的优点。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。