裂纹分析装置和方法与流程

本发明大体涉及一种装置和方法。特别地，但非排他地，本发明涉及用于分析车辆的玻璃窗板中的裂纹的装置和方法。进一步特别地，但非排他地，本发明涉及用于分析玻璃——特别是车辆的玻璃窗板——中的裂缝的方法和设备。

背景技术：

在驾驶过程中，如果道路上有碎片和其他物质，会导致这些物质转移到行车道中，当这些物质碰撞到挡风玻璃上时，会导致挡风玻璃出现裂缝，裂纹或其他损坏，这就需要修复或更换车辆的玻璃窗板。

出于安全原因和经济原因，修复或更换工作必须尽快进行，因为这些裂缝可能由于寒冷天气的影响而进一步扩大到挡风玻璃中，例如，寒冷天气会导致裂缝增大。这可能导致裂缝从需要小规模修复变为需要完全更换挡风玻璃。

需要评估是否可以通过修复来补救车辆的玻璃窗板的损坏。如果评估表明修复不可行，则需要更换玻璃窗板。

基于对前述内容的考虑，构思了各方面和实施例。

技术实现要素：

从第一方面来看，把发明提供了一种用于分析车辆玻璃窗板中的裂纹的裂纹分析方法，该方法包括：捕获车辆玻璃窗板中的裂纹的图像；处理裂纹的图像。从第二方面来看，本发明提供了一种用于分析车辆玻璃窗板中的裂纹的裂纹分析装置，该装置包括：相机，用于捕获车辆玻璃窗板中的裂纹的图像；处理模块，用于处理裂纹的图像。

可选地，该装置可包括移动计算设备，该移动计算设备包括相机。移动计算设备为用于捕获图像的电子设备，可以包括移动电话(例如智能电话)，笔记本计算机，平板电脑，平板手机或照相机。移动计算设备包括用于捕获图像的相机。

移动计算设备还可以包括处理模块。

根据第一和第二方面的方法或设备使得能够使用表面中的裂纹的图像来分析裂纹。这就使技术人员不必亲自到场对裂纹进行任何分析。

可以基于图像的处理确定是否需要更换玻璃窗板。因此，该方法可包括基于图像的处理来确定是否需要更换玻璃窗板，和/或玻璃窗板是否适合修复的步骤。

可以以相对于车辆玻璃窗板倾斜的角度来捕获裂纹的图像。

该图像可由与玻璃窗板的表面保持接触的移动计算设备捕获，其中移动计算设备包括相机。移动计算设备可以例如是移动电话，例如设置有相机的智能电话。本发明可以通过软件构件实现，该软件构件用于处理来自相机的图像数据，以便确定是否可以修复裂纹或者是否需要更换玻璃窗板。可以将软件构件下载到设备，例如作为独立应用程序(app-application)等。软件元件可包括用于进行确定的算法，并且还优选地显示指示如何实现确定方法的用户指令。

因此，根据另一方面，本发明包括一种用于确定车辆玻璃窗板中存在的裂纹的技术信息的计算机实现的方法，其中该方法包括将软件构件下载到具有相机的移动计算设备的步骤，该软件构建用于处理通过相机捕获的图像信息，以确定优选是修复裂纹还是更换。

移动计算设备可以是倾斜的，以在相机的视场中的预定位置查看裂纹。可以通过确保计算设备向玻璃窗板倾斜并且边缘接触玻璃窗板来实现基准位置。

可以通过在相机上或在包括相机的移动计算设备上(例如在屏幕上)显示的标记来指示预定位置，在该预定位置移动计算设备可倾斜来观看裂纹。

附加地或替代地，标记可以指示视场的中心(例如，在屏幕上)。

可选地，相机或移动计算设备可最初平放在玻璃窗板的表面上。然后，移动计算设备可枢转离开玻璃窗板或倾斜于玻璃窗板，其中至少部分移动计算设备保持与玻璃窗板接触。

可选地，在移动计算设备枢转之前，相机定位在相对于裂纹的预定位置。例如，该方法可以包括将移动计算设备的特征对准在相对于裂纹的预定位置。

在一些实施例中，图像捕获模块的边缘——例如顶部边缘——可与裂纹的最低点紧邻(或紧挨着其下方)处对齐。这意味着可以使用移动计算设备的几何结构来确定相机和裂纹的最低点之间的初始距离。

然后，图像捕获模块或移动计算设备可以枢转离开玻璃窗板，其中移动计算设备的底部边缘保持与玻璃窗板接触。然后捕获裂纹的图像。

该方法可以包括枢转或倾斜移动计算设备，直到在相机的视场中的预定位置处可以看到裂纹。

可选地，该方法包括枢转移动计算设备，直到裂纹位于相机的视场的中心。

该方法可以包括使用移动计算设备的几何参数和相机的镜头参数来估计裂纹的参数。裂纹的参数包括指示裂纹尺寸的一个或多个空间维度。

例如，使用上述方法，移动计算设备的几何参数和相机的镜头参数可用于确定裂纹的一条或多条支路的长度和/或裂纹中心的尺寸(例如，宽度/直径)。

移动计算设备的几何参数可用于确定移动计算设备和玻璃窗板之间的旋转角度(或枢转角度)。

如果估计的裂纹尺寸(例如裂纹的一个支路的长度)超过给定阈值，则该方法可确定需要更换玻璃窗板。如果估计的裂纹的尺寸(例如裂纹的一个支路的长度)小于给定阈值，则该方法可以确定玻璃窗板适合修复。

图像的处理可基于移动计算设备的几何参数来捕获裂纹的图像。

图像的处理还可基于相机和/或移动计算设备的芯片参数。

图像的处理可以生成一组数据点，这些数据点用于生成裂纹的直径，该直径可以用于确定更换车辆玻璃窗板的要求。

该方法可包括输出指示玻璃窗板是否需要修复或更换的信号或指示。

玻璃窗板中的裂纹可包括中心及一条或多条支路。当小石头或其他物体撞击玻璃窗板时，这种裂纹形式很常见。所述一条或多条支路(或裂缝(cracks))通常从裂纹(break)中心向周围辐射。

裂纹中心的形状可基本上为圆形。

确定更换或修复玻璃窗板的需要可以包括生成指示裂纹中心及裂纹的一条或多条支路的数据。

该方法可以包括识别裂纹中心及裂纹的一条或多条支路。

该方法可以包括生成比例因子，该比例因子指示与裂纹中心的尺寸相比，裂纹的一条或多条支路的相对长度。裂纹中心的尺寸可以是裂纹中心的直径、宽度和/或长度。

该方法可以包括使用比例因子估计裂纹的一条或多条支路的长度(即，以cm或mm等为单位的绝对长度)。例如，可以通过将比例因子乘以预定值来估算一条或多条支路的长度。

预定值可为对裂纹中心的实际(即绝对)尺寸的估计，优点是不需要校准对象，这对于用户来说至少是更方便。

因此，本发明的方法可包括确定图像的比例值以估计裂纹的一条或多条支路的长度。

如果一条或多条支路的估计长度超过给定阈值，则该方法可确定需要更换玻璃窗板。如果一条或多条支路的估计长度小于给定阈值，则该方法可以确定玻璃窗板中的裂纹适合修复。

该方法可以包括：如果一条支路或多条支路的估计长度小于给定阈值，则输出指示需要修复玻璃窗板的信号。

该方法可以包括：如果一条支路或多条支路的估计长度超过给定阈值，则输出指示需要更换玻璃窗板的信号。

已经发现，当比较璃窗板中的不同裂纹时，裂纹中心的尺寸通常是不同的且小于裂纹的一条或多条支路的长度。这样，预定值可以是玻璃窗板中的裂纹中心的测量尺寸的平均值或众数。

可选地，裂纹中心的实际宽度(或直径)的估计值(即预定值)可以在1mm和3mm之间。裂纹中心的特别优选的预定宽度(或直径)可以是2mm。这些范围/值已经由申请人进行的裂纹研究确定。

这样的效果是，裂纹中心的估计尺寸，例如裂纹中心的估计直径，可用于估计裂纹支路的长度，因为预定值可用于在相对于裂纹中心的尺寸一条或多条支路的相对长度，和裂纹的一条或多条支路的估计的实际长度之间进行变化。

例如，如果我们知道裂纹中心总是宽约(或直径为)2mm，并且生成的图像数据表明支路长度是裂纹中心直径的两倍，那么该方法可以包括将2mm乘以比例因子2。这就估计出支路的长度为4mm，有助于在数据中构建裂纹维度的图片。

一条或多条支路的生成的(或估计的)长度可用于指示裂纹的估计的尺寸。可以将裂纹的尺寸与阈值参数进行比较，以确定需要更换玻璃窗板或者修复玻璃窗板。

如果估计的裂纹尺寸超过给定阈值，则该方法可确定需要更换玻璃窗板。如果估计的裂纹尺寸小于给定阈值，则该方法可以确定玻璃窗板适合修复。

可以对裂纹阈值参数和裂纹的最大跨度距离进行比较。

可选地，裂纹中心的尺寸的预定的估计可取决于一个或多个参数。参数可以由用户输入和/或预先设置到设备或处理模块中。例如，参数可包括：玻璃窗板的一个或多个属性(例如类型，尺寸等)，和/或当发生破裂时车辆行进的速度。

图像的处理可包括过滤图像以去除背景部分以识别裂纹。

可以对图像应用形态细化以从图像中去除任何杂质，并改善图像数据的质量，该图像数据用作确定是否需要更换玻璃窗板的基础。

该方法可包括在捕获裂纹的图像之前清洁玻璃窗板。这可以帮助去除任何可能影响图像处理的污垢。例如，污垢有可能被图像处理软件构建为裂纹的风险。

该方法可包括在捕获裂纹的图像之前，禁用图像捕获模块或设备的闪光功能。如果使用闪光摄影，则光可能不利地影响图像处理软件的准确性。例如，闪光可以在玻璃窗板中反射，这可能影响对裂纹的识别或分析。

该方法可使用计算机实现的指令来实现，当安装到存储器中时，所述指令指示处理器实施如上定义的方法，优选为可下载的软件构建(例如app)。

应当理解，可使用本发明的装置执行该方法的任何特征。

参考本文所述的实施例，本发明的这些方面和其他方面将变得清楚明了。

附图说明

现在将仅通过示例并参考附图来描述本发明的第一和第二实施例。其中：

图1示出了一个带有裂缝的挡风玻璃；

图2示出了使用相机捕获裂缝的图像；

图3示出了可用于分析图1的挡风玻璃中的裂缝的处理模块；

图4示出了详细描述使用图3的系统评估挡风玻璃中的裂缝所涉及的步骤的流程图；

图5示出了裂缝的傅立叶变换的图像；以及

图6示出了裂缝的轮廓图像；

图7示出了可用于模拟裂缝102的布置；

图8示出了生成裂纹参数所涉及的步骤；

图9a示意性地示出了相机对挡风玻璃的俯视视场；

图9b示意性地示出了相机对挡风玻璃的倾斜视场；

图10示出了如何对相机进行建模以便使用相机的参数来分析裂缝；以及

图11示出了从系统输出的，用以确定挡风玻璃中的裂缝的尺寸的图像。

具体实施方式

在第一实施例中，图1示出了挡风玻璃100，其具有由弹击到挡风玻璃上的石头引起的裂缝102。然后，挡风玻璃100封装的汽车的驾驶员使用移动电话104捕获裂缝102的图像，移动电话104包含用于捕获裂缝102的图像的相机106。这种设置在图2中以侧面视图示出。相机106的焦距固定在100毫米以下，以使相机的焦点近距离紧密聚焦。

然后响应用户输入，捕获裂缝102的图像，移动电话104用于向用户提供提示，以请求用户指明他们希望将图像从相机106传输到图像处理模块108。现在我们将参考图3描述图像处理模块108。该步骤使用户能够自己评价图像的质量，从而，如果他们认为图像不清晰——例如在恶劣天气条件下，可能会有水珠凝结在相机106的镜头上——他们可以选择捕获另一张图像。

使用任何合适的方法，相机106将捕获的图像转换为图像数据阵列。相机106可将图像数据保存为可交换图像文件(exif)，其中相机的镜头参数也被存储。

也就是说，相机106是图像捕获模块的示例，其可运行以捕获裂缝102的图像，并且以捕获的图像数据的形式将捕获的图像传输到图像处理模块108，在图像处理模块108中该图像可进一步被处理，以提取裂缝的细节。

图像处理模块108可形成为移动电话104的一部分，或者它可地理上相对远离移动电话104。图像数据通过任何合适的手段——例如，数据总线或互联网——传输到图像处理模块108。

在优选实施例中，图像处理模块108作为下载到移动电话的软件构建得以实现。这可以通过将软件构建下载为app来实现。软件构建能够处理来自相机的图像数据，以确定裂纹是否易于修复，或者是否需要将更换玻璃窗板作为优选解决方案。

裂纹分析模块112可为下载到移动电话的软件构建，优选地，作为与图像处理模块108结合的单个下载。优选地，单个下载的软件构建i用于处理来自电话相机的图像数据，并使用软件中实现的一个或多个算法分析裂纹。

在一个实施例中，捕获的图像数据由图像处理模块108在数据输入接口110处接收。然后将图像数据传输到裂纹分析模块112，裂纹分析模块112用于访问例程库114，其中可存储例程，以在对捕获的图像数据的分析过程中，对捕获的数据执行操作。

裂纹分析模块112还用于访问设备参数数据库116，其中存储与移动电话104相关的参数。

与移动电话104相关的参数包括定义相机106的图像捕获性能的芯片参数，例如，焦距，镜头的传感器尺寸，以及移动电话104的维度参数，例如，移动电话104的长度，以及移动电话104的顶部边缘与相机106的图像中心之间的距离。

裂纹分析模块112还可用于与显示模块118连接，显示模块118可运行以在显示器上显示从裂纹分析模块112传输的图像数据，并在显示器上显示从裂纹分析模块112传输的参数数据。

现在，我们将参考图4描述使用裂纹分析模块112对裂缝102进行的分析。

在步骤s400中，裂纹分析模块112接收图像数据。然后，在步骤s402中，裂纹分析模块112从例程库114调用傅里叶变换例程，并使用傅里叶变换例程将二维离散傅立叶变换(discretetwo-dimensionalfouriertransform)应用于图像数据以产生如图5a所示的已变换的图像。

在图5a中，我们可以看到已变换的图像。图5a根据相应的傅立叶分量的量级绘制了的空间频率。可以看出，低空间频率占据已变换的图像中心，当我们远离已变换的图像中心时，可以看到更高的空间频率。

对图像进行傅里叶变换，使得裂纹分析模块112能够对图像的分量空间频率(componentspatialfrequencies)和相位方面进行图像分析。如现在将描述的，这去除了我们不感兴趣的空间频率并且使我们能够通过保留感兴趣的空间频率来重建我们感兴趣的图像。

在步骤s404中，然后，通过裂纹分析模块112将巴特沃斯带通滤波器(butterworthbandpassfilter)应用于已变换的图像。由巴特沃斯带通滤波器实现的掩模(mask)在图5b中示出。巴特沃斯带通滤波器在图5a所示的已变换的图像上实现了一个掩模，并去除了低空间频率(图5b中心的黑点所示)和非常高的空间频率(图5b中图像中的暗边界)，非常高的空间频率代表图像上的灰尘和污垢斑点。

在步骤s406中，然后，通过从例程库114调用傅里叶变换反演例程来反转图像数据的傅里叶变换，以对已变换的图像数据执行逆离散二维傅里叶变换(inversediscretetwo-dimensionalfouriertransform)。

对已变换的图像数据的逆傅立叶变换过程，将已变换的图像数据从傅里叶域变换到实域，以生成实域图像数据。得到的实域图像数据在图6a和6b中示出。

使用傅立叶变换产生图6b中所示图像，具有将裂缝从背景隔离的效果。

在步骤s408中，将实域图像数据与阈值强度4进行比较，以更清楚地描绘感兴趣区域。

图6a示出了不使用巴特沃斯带通滤波器的实域图像数据。图6b示出了在将巴特沃斯带通滤波器应用于图5a中所示的已变换的数据之后的实域图像数据，并且应用阈值处理以使用阈值强度4生成二值图像。该示例中的巴特沃斯带通滤波器具有衰减(rolloff)值3。

可对巴特沃斯带通滤波器的上截止频率和下截止频率进行建模，因为它们与图像最长边的像素个数(表示为m)线性地相关，可分别表示为：

可以使用标准试验和数值实验来改变这种关系。

图6b中所示的图像是不只包括裂缝102的图像。该图像还可包括已经通过步骤s400至s408的图像数据，但是由挡风玻璃上的污垢斑点和通过裂纹分析模块112所执行的处理的其他伪影(artefacts)所造成的。

使用阈值强度4来生成图6b中所示的二值图像有助于更清楚地显示感兴趣的区域。图6b中所示的实域图像数据所示的图像突出了裂缝——包括中心裂缝区域——这是一个低空间频率的区域。

可以看出，傅里叶方法做了一个非常巧妙的工作，即将裂缝区域从杂乱的背景中隔离，从而假设它处于焦点而背景没有处于焦点。

然后，裂纹分析模块112可以从例程库114调用形态例程，以在步骤s410从图6b中所示的图像中去除任何杂质。

形态学例程对图6b中所示的图像执行若干操作。该图像是二值图像。黑色区域为零值，白色区域值为非零值。像素的值存储在裂纹分析模块112的存储器中，并且是步骤s400至s408中处理的结果。

这些操作中的第一个是填充操作，该操作使用形态重建，用白色来填充由白色区域包围的像素大小的黑色区域，通过根据[1]所列出的过程将零值替换为非零值来进行填充。

这些操作中的第二个是清理操作，其丢弃非常小的非零值区域。非常小的非零值区域被定义为占据小于(图像的最大维度/500)平方的区域的非零值区域。图像的最大维度可通过裂纹分析模块简单地通过比较图像的宽与图像的长来确定。

然后重复第一形态学操作以填充由第二形态学操作产生的被白色区域围绕的任何像素大小的黑色区域。这是第三种形态学操作。

然后实施第四形态学操作以连接裂缝102的图像中的任何支路，所述支路具有间隙。这是使用[2]中描述的形态学闭合操作实现的。通过利用半径为(图像的最大维度\5312)乘以20的盘形结构元件进行侵蚀，然后进行扩张。值20是凭经验确定的，并且可以改变。可以确定该值，而对不同的图像分辨率不造成任何不适当的负担。

然后重复第一形态学操作以填充由第四形态学操作产生的被白色区域围绕的任何像素大小的黑色区域。这是第五次形态学操作。

然后执行第六形态学操作以丢弃任何非零的小区域。小区域定义为面积等于(图像的最大维度/100)平方的区域。

然后执行第七形态学操作以移除图像中的任何断开连接的对象。感兴趣的断开连接的对象是比最接近图像中心的最大对象的半径的3/4要远的对象。这意味着包括仍不连贯的裂缝的支路，但包括多余的伪影。通过为图像中的每个剩余区域找到质心(即图像的质心)和该区域的长轴的长度来实现第七形态学操作。基于质心与图像中心的接近程度，为每个区域分配额外的加权。

加权w＝1/d²，其中d是质心和图像中心之间的欧几里德距离(euclideandistance)。选择最接近图像中心的最大区域，并且其长轴长度用于设置半径(或者从其质心算起的长轴长度的3/4)，丢弃在该半径之外的任何区域。也就是说，形态例程和质心到边界计算用于将所有“斑点”保持在一定距离内，该距离为距图像中的最大对象的半径的裂缝中心的距离加上该半径的一半，以确保不会丢失任何裂缝102中的不连续部分。

在应用形态学以细化图像数据之后，图像数据可以用于确定裂缝102的维度。

裂纹分析模块112应用进一步的边缘检测、形态学、模糊处理和阈值处理来确定裂缝102中心。

通过实验已经观察到裂缝中心直径通常约为2mm。裂纹分析模块112用于使用细化的图像数据和由裂缝102中心的确定而得到的数据，估计裂缝102的支路的长度，并确定比例值，该比例值表征与裂缝102中心的直径相比的支路的长度，即与裂缝102中心相比的支路的比例因子。使用观察到的裂缝中心通常为2mm，则比例因子可以用来确定支路的长度。这就提供了裂缝102的尺寸的未校准分析。

然后所确定的支路的长度可以用来接近裂缝102的尺寸，并且使裂纹分析模块112能够输出是否需要更换挡风玻璃，或者是否修复挡风玻璃就足够了，因为裂缝的尺寸是决定做出哪种决定的重要因素。通过将裂缝102的尺寸与修复/更换阈值进行比较，裂纹分析模块112可以自动化该决定。裂纹分析模块112将该决定输出到显示模块118。

在步骤s412中，然后，使用显示模块118显示输出，即是否需要更换挡风玻璃。

使用观察到的裂缝中心的估计来估计裂缝支路的尺寸，这取决于对裂缝径向的尖锐程度的假定，这就意味着可以拍摄裂缝的图像，并用于在现场不借助任何标度的情况下，分析裂缝以提供裂缝102的尺寸。

该方法能够在各种条件下进行裂缝分析，而无需技术人员的参与。

在第二实施例中，我们现在描述如何使用移动电话104的参数和相机106的镜头来获得裂缝102的参数。这可以帮助校正角度对图像的任何影响。

在不偏离本发明的情况下，第二实施例可以与第一实施例组合。

图2中所示的布置使得使用相机106的芯片参数和移动电话104的几何参数来估计裂缝的维度成为可能。

为了计算移动电话104相对于挡风玻璃的旋转角度(或枢转角度或倾斜角度)，我们可以使用移动电话104的几何参数。

在将裂缝102定位在相机106透镜的视场的中心时，能够绘制直角三角形。这是参考图7进行描述的。

在发现裂缝102之后，移动电话104平放在挡风玻璃上，顶部边缘位于裂缝102的底部。这意味着，移动电话104的底部边缘与裂缝的底部之间的距离等于移动电话104的长度。移动电话104从移动电话104的底部边缘倾斜，直到裂缝102位于相机106的视场的中心。可以在移动电话104的显示区上作标记，以指示视场的中心。

可以从设备参数数据库116中获得移动电话104的底部边缘和相机106的镜头之间的距离。因此，形成由移动电话104的底部边缘和挡风玻璃100之间的旋转角度，相机镜头的z轴，以及底部边缘与裂缝底部之间形成的距离限定的直角三角形。

我们现在描述如何使用移动电话104的几何参数和镜头参数来估计裂缝的参数。

捕获裂缝的图像与上述过程一致，其中旋转移动电话104，直到裂缝102位于相机106的视场的中心。

这就使得能够通过相机镜头的z轴、所述底部边缘和裂缝底部之间形成的距离、以及所述底部边缘和相机镜头之间的长度形成直角三角形。

我们参考图8来描述如何使用移动电话104的几何形状和镜头参数来估计裂纹的参数。

在步骤s800中，裂纹分析模块112从设备参数数据库116获得移动电话104的底部边缘与裂缝底部之间的距离(即移动电话104的长度)，以及移动电话104的底部边缘与相机镜头之间的长度。然后，可以在步骤s802中使用底部边缘和裂缝底部之间的距离与底部边缘和相机镜头之间的长度之间的余弦关系来计算移动电话104的旋转角度。

然后，我们需要使用相机参数来推导出相机的像素与图像的真实世界空间维度之间的平面-到-平面单应映射(plane-planehomographicmapping)。然后在步骤s804中从例程库114调用平面-到-平面单应映射例程，以导出图像的真实世界空间维度。

单应映射的推导以提供图像的真实世界空间维度是单应映射基于“针孔相机模型”，其中相机认为相机的视图区域是相对于相机106的镜头展开的底部为矩形的锥体，如图9a和9b中所示。

图9a仅用于说明，其示出了移动电话104位于挡风玻璃100正上方的情况。也就是说，相机106提供了挡风玻璃100的俯视图。在这种情况下，视场区域a1是矩形，且每个像素占用相同数量的真实世界空间(以毫米为单位)。

在这种情况下，如图9b所示，移动电话104相对于挡风玻璃100成一角度。在步骤s802中计算该角度。然后，视场区域a2变为梯形，这意味着靠近相机的像素代表的距离比距离更远的像素代表的更少的毫米。

我们描述了如何推导出平面-到-平面单应映射的理论基础，但是应当理解，这将使用例程以数值表示而实现，使用例程库114，该例程将对裂纹分析模块112可用。

考虑形成相机106的一部分的矩形图像传感器和从平面倾斜一个旋转角θ的传感器，传感器观察到的区域映射到等腰梯形。该梯形的基部宽度直接取决于θ。使用平面-到-平面单应映射例程，我们可以使用该原理使用相机106上的像素知识对裂缝102的参数进行数值估计。

我们定义围绕x轴的3d旋转矩阵，作为θ的函数，如下所示：

应当理解，θ是移动电话104相对于挡风玻璃的角度。我们可以在(0，0，0)处定义笛卡尔(cartesian)x，y和z维度的原点，即世界原点。这是与x轴对齐的移动电话104的底边的中间的点。然后，该坐标系的y轴从底部垂直指向电话顶部。为了简单而不失普遍性，如果我们假设相机位于y轴上，距离手机底部一定距离dc，则相机中心定义为：

然后，在步骤s806中，可以从设备参数数据库116中获得相机106的镜头的焦距及垂直和水平传感器的尺寸。这些参数可以称为芯片参数。这使我们能够从相机计算视野区域。视野区域由两个量定义，称为水平和垂直视角(分别表示为αh和αv)，由以下等式定义：

其中，sx和sy为水平和垂直传感器的尺寸，f是焦距。

步骤s808，在计算了水平和垂直视角之后，裂纹分析模块112使用平面-到-平面单应映射例程来计算视场金字塔的边缘，以在挡风玻璃100上为我们提供视场。这就提供了图9b中所示的梯形，即我们需要校正的梯形，以补偿距离透镜更远的像素相对于更靠近透镜的像素所占据的不同空间量。也就是说，我们需要缩放梯形以确保正在进行的计算为每个像素分配相等的真实世界空间量。

这是在由裂纹分析模块112使用的平面-到-平面单应映射例程中由线——即射线——建模的，该线从透镜沿着透镜和裂缝102之间的视线延伸。该线将与由挡风玻璃代表的平面相交——即平面-到-平面的单应映射例程建模为平面。

在步骤s810中，平面-到-平面的单应映射例程从例程库114调用数值求解器例程以求解定义挡风玻璃的平面的联立方程以及从透镜沿着透镜和裂缝102之间的视线延伸的线。假定限定挡风玻璃100的平面是平的并且相机106相对于该平面倾斜的情况下，编程平面-到-平面单应映射例程。这就提供了从镜头沿着视线延伸的线与挡风玻璃100的平面之间的交点。

理论上，上述内容可以表示为对射线的计算：所述射线从相机的中心点放射，经由传感器/图像平面的各个角到挡风玻璃上，形成上述梯形。

我们首先在给定如上所述水平和垂直视角αh和αv的情况下，获得射线与平面的交点，该平面以单位距离与图像平面平行。

有四条射线，每条用于矩形传感器的每个角。最小和最大x值可以定义为：

同理，我们可以将最小和最大y值定义为：

然后我们可以将矩形传感器的角定义为：

按照这些坐标的量级对其进行标准化，为我们提供了4条射线的方向。我们将每个坐标的射线方向定义为：

如果我们假设手机在x轴上旋转θ的角度，我们可以计算出相机中心的位置为：

这使我们能够将射线的方向定义为：

这就给我们提供了笛卡尔坐标系中的射线，具有一个与平行于图像平面的平面相交的已知的点，我们知道这种相交只发生一次。这就提供了指示现实世界中的视场的梯形。

我们将梯形的角定义为：

vi，i∈(tl，tr，br，bl)

我们使用[3]中描述的线平面交点公式计算梯形的顶点。

我们知道挡风玻璃平面的法线是向量n＝(0，0，-1)，位于世界原点上，这意味着交点公式简化为：

其中，所述点vi，i∈(tl，tr，br，bl）为梯形的顶点，我们需要使用梯形顶点和图像坐标之间的四点对应技术，来定义从图像平面到现实世界中的平面的单应映射h，图像坐标：

其中w是图像的宽度，h是图像的高度。在[4]中讨论涉及如何获得该单应映射的算法。

挡风玻璃上方的相机的高度可以由裂纹分析模块112使用毕达哥拉斯(pythagoras)定理计算，因为在步骤s800中，移动电话104的底部边缘与裂缝的基部之间形成的距离(即移动电话104的长度)和移动电话104的底部边缘与相机镜头之间的长度已经从设备参数数据库116获取并且仍然在裂纹分析模块112的存储器中。

步骤s810的输出是现实世界中的梯形视图(x1，x2，x3，x4)。参数(x1，x2，x3，x4)与挡风玻璃上捕获的图像的各个角之间的比较(在步骤s812中由裂纹分析模块112执行)提供了将相机106的像素的位置映射到挡风玻璃100上的视场上以毫米为单位的位置所需的标定(scaling)。这就为我们提供了平面-到-平面单应映射。标定是3×3矩阵的形式，表示相机的视场和挡风玻璃100之间的比例值、旋转、倾斜和平移。

平面-到-平面单应映射能够校正视角对捕获图像的影响，以及从像素维度到毫米的转换使得裂纹分析模块112能够推导出裂缝102的维度参数。

平面-到-平面单应映射为矩阵，该矩阵将相机106的二维图像平面映射到表示挡风玻璃的平面上。

平面-到-平面单应映射的输出提供了以毫米为单位的正射纠正掩模，指示裂缝的位置和形状。

响应于平面-到-平面单应映射的这个输出，如所理解的将为平面-到-平面单应映射例程的输出，裂纹分析模块112从例程库114调用凸包计算例程。挡风玻璃的视场上以毫米为单位的位置被提供给来自例程库114的凸包计算例程。

总之，凸包是覆盖视场上以毫米为单位的每个位置的空间。凸包计算例程的输出是可以简单地表示为“二进制大对象(blob)”的数据，其与检测到的裂缝102具有相同的尺寸。这使得能够使用blob对检测到的裂缝102执行分析。

然后，裂纹分析模块112从例程库114调用最小圆例程，其实现了凸包的最小圆问题的数值解，该凸包由凸包计算例程输出。该模块输出了包围凸包中的每个点的最小圆，因此为裂缝102提供最小半径。

表示凸包的数据，表示凸包的最小圆问题的解的数据和计算的裂缝半径分别由裂纹分析模块112存储在存储器中，该存储器是处理模块108的本地存储器，或远离处理模块108。

也就是说，裂纹分析模块112已经使用移动电话104的几何参数和相机106的参数来生成裂缝102的半径。

然后，在步骤s814中，可以使用显示模块118显示由最小圆例程输出的参数和圆。

可以由显示模块118提供的示例图像在图11中示出。在这种情况下，最小圆的直径表示为16mm，即半径为8mm。在这种情况下估计的最大裂缝直径为16毫米。这里的效果是估计出了裂缝的最小尺寸并且裂缝的最小尺寸可以用于确定更换挡风玻璃的必要性。

可以通过裂纹分析模块112将估计的半径与更换/修复阈值进行比较，以确定裂缝102是否需要更换或者仅作修复是否足够。

移动电话104上存在的壳体可能会在测量参数中引入错误，因为这将增加移动电话104的长度，但误差通常约为3％。在裂纹分析模块112的计算中置入3％的误差容限并由显示模块118提供在显示器上。

移动电话104的底部和相机106之间的距离也可能无法从设备参数数据库获取。在这种情况下，我们可以估计参数以改善所描述方法的稳健性。

当捕获裂缝102的图像时，我们可以使用内置于移动电话104中的倾角计来获得移动电话的角度。这可用于计算高度h，基于以下等式：

h＝l*sin(θ)

其中l是移动电话104的长度，角θ是从倾角计获得的角度。

类似地，可以使用视角，图像分辨率和传感器尺寸来估计电话角度。

如第一技术所描述的，在优选实施例中，图像处理模块108实现为下载到移动电话的软件构建。这可以通过将软件构建下载为app来实现。软件构建能够处理来自相机的图像数据，以确定裂缝是否易于修复，或者需要更换玻璃窗板是否为优选的解决方案。

裂纹分析模块112可以是下载到移动电话的软件构建，优选地作为与图像处理模块108结合的单个下载。单个下载的软件构建i优选地用于处理来自电话相机的图像数据，并使用在软件中执行的一个或多个算法来分析裂纹。

裂纹分析模块112可运行以基于最小圆的半径在显示器上提供关于是否需要完全更换挡风玻璃的警报。如果裂缝102高于指定阈值，则裂纹分析模块112将指示需要更换挡风玻璃，否则不需更换。警报可以显示在移动电话104的显示器上。

应该注意，上述实施例说明而不是限制本发明，并且本领域技术人员将能够设计许多替代实施例而不脱离所附权利要求所限定的本发明的范围。在权利要求中，括号中的任何附图标记都不应解释对权利要求的限制。词语“包括”(comprising)和“包括”(comprises)等不排除除了在任何权利要求或说明书的整体中列出的元件或步骤之外的元件或步骤的存在。在本说明书中，“包括”(comprises)意指“包括(includes)或由......组成”并且包括“包括”(comprising)意指“包括”(including)或“由……组成”。元素的单数引用不排除这些元素的复数引用，反之亦然。本发明可以借助于包括若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的设备权利要求中，这些装置中的若干个可以由同一个硬件项来体现。单纯在相互不同的从属权利要求中叙述某些方法并不代表这些方法的组合不能用于获益。

参考：

[1]http://www.mathworks.com/tagteam/6419991822v00eddinsfinal.pdf

[2]https://en.wikipedia.org/wiki/closing(morphology)

[3]josepho’rourke“computationalgeometryinc”，cambridgeuniversitypress，2012

[4]richardhartleyandandrewzisserman“multipleviewgeometryincomputervision”，cambridgeuniversitypress，2011.