CN114543684B - 一种结构位移测量方法 - Google Patents

一种结构位移测量方法 Download PDF

Info

Publication number
CN114543684B
CN114543684B CN202210441171.2A CN202210441171A CN114543684B CN 114543684 B CN114543684 B CN 114543684B CN 202210441171 A CN202210441171 A CN 202210441171A CN 114543684 B CN114543684 B CN 114543684B
Authority
CN
China
Prior art keywords
measured
displacement
smart phone
point
plane
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210441171.2A
Other languages
English (en)
Other versions
CN114543684A (zh
Inventor
袁杨翔
贾哿媛
武斯珩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences Beijing
Original Assignee
China University of Geosciences Beijing
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Geosciences Beijing filed Critical China University of Geosciences Beijing
Priority to CN202210441171.2A priority Critical patent/CN114543684B/zh
Publication of CN114543684A publication Critical patent/CN114543684A/zh
Application granted granted Critical
Publication of CN114543684B publication Critical patent/CN114543684B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/70—Determining position or orientation of objects or cameras

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明提供了一种结构位移测量方法，涉及工程结构检测技术领域，具体包括：准备步骤：将智能手机架设在第一位置，使得被测物体以及与被测物体在同一平面的已知结构尺寸的已知结构体均可在智能手机的视场中显示，将此时的三维角度均置为零；三维角度获取步骤：将智能手机放置到与被测平面平行的第二位置并记录此时的三维角度后将手机复位；拍摄计算步骤：在已知结构体上选取多个已知点，对被测物体变形或位移发生前后进行连续拍摄图像，计算出各个已知点的像素位移与实际位移的比例关系，得到待测点的实际位移。上述结构位移测量方法仅使用智能手机即可实现，具有数据易于获得、操作难度低、不易受外界条件限制、成本小、易于推广等优点。

Description

一种结构位移测量方法
技术领域
本发明涉及工程结构检测技术领域，尤其是涉及一种结构位移测量方法。
背景技术
现有的结构位移测量方法主要分为接触式和非接触式两种。接触式测量方法需要在结构表面安装相关设备，费时费力，是较为传统的测量方式。随着技术的发展，近年发展起来的先进结构测量方式以非接触式为主。基于机器视觉的测量方法以其系统安装简便、成本低廉、精度可控等显著优点，在非接触式方法中最为常见。
基于机器视觉的非接触式测量方法的不足——以工业相机为主体的测量系统，为了满足工程结构测量需要，除了需要进行精确的图像匹配，在图像匹配之前，还需要对相机系统进行标定，从而计算得到图像像素位移到物理位移的关系。在进行工程结构测量时，相机成像的视场范围很大，不可能用传统的拍摄标定板的方法进行标定。因而，基于此类方法的所有设备一般都包含测距机、倾角仪等辅助设备，由测得的辅助参数进行标定。
由此，各种辅助设备大大提高了该方法实现的成本，系统不再紧凑便携，受外界条件限制较大。
发明内容
本发明的目的在于提供一种结构位移测量方法，仅使用具有拍照和陀螺仪功能的智能手机即可实现，具有数据易于获得、操作难度低、减小成本、易于推广等优点。
为实现上述目的，本发明提供以下技术方案：
本发明提供一种结构位移测量方法，包括：
准备步骤：将智能手机架设在第一位置，使得被测物体以及与所述被测物体在同一平面的已知结构尺寸的已知结构体均可在所述智能手机的视场中显示，使用所述智能手机中的陀螺仪将此时的三维角度均置为零；
三维角度获取步骤：将所述智能手机放置到与被测平面平行的第二位置并记录此时的三维角度后将手机复位；
拍摄计算步骤：在所述已知结构体上选取多个已知点，对所述被测物体变形或位移发生前后进行连续拍摄图像，得到各个所述已知点以及待测点的像素坐标，根据第一位置拍摄面与所述被测平面的相对三维角度、各个所述已知点之间的实际尺寸以及各个所述已知点在所述第一位置拍摄面中的像素坐标，计算出各个所述已知点的像素位移与实际位移的比例关系，根据所述比例关系以及所述待测点的像素位移得到所述待测点的实际位移。
进一步地，所述拍摄计算步骤中，计算出所述已知点的像素位移与实际位移的比例关系包括：
确定所述被测平面的空间方程，计算每一个所述已知点在所述被测平面上的位置与在所述第一位置拍摄面上的位置之间的距离；
根据所述距离算得每个所述已知点的像素位移与实际位移的比例关系。
进一步地，设定所述第一位置拍摄面为
Figure M_220425100944750_750496001
，所述被测平面为
Figure M_220425100944766_766157002
，所述
Figure M_220425100944781_781765003
与所述
Figure M_220425100944813_813071004
的相对三维角度为
Figure M_220425100944828_828623005
，所述拍摄计算步骤中，确定所述被测平面的空间方程包括：
建立所述
Figure M_220425100944859_859871001
与所述
Figure M_220425100944875_875500002
的几何关系：以所述
Figure M_220425100944891_891123003
为z=0平面，所述
Figure M_220425100944906_906753004
的中心为坐标原点
Figure M_220425100944939_939427005
（0,0,0），所述智能手机的相机光心
Figure M_220425100944955_955570006
的坐标为（0,0,-
Figure M_220425100944986_986816007
），
Figure M_220425100945002_002447008
为所述智能手机的相机焦距，所述
Figure M_220425100945033_033717009
的方程设为
Figure M_220425100945064_064942010
；
确定方程系数
Figure M_220425100945096_096224001
、
Figure M_220425100945111_111831002
和
Figure M_220425100945127_127444003
：设定所述
Figure M_220425100945147_147485004
的法向量
Figure M_220425100945179_179230005
为（0,0,1），设定
Figure M_220425100945210_210451006
旋转
Figure M_220425100945241_241712007
以后，新的法向量为
Figure M_220425100945257_257327008
，
Figure M_220425100945288_288570009
为所述
Figure M_220425100945304_304226010
的法向量；
Figure M_220425100945319_319842001
；
Figure M_220425100945431_431188001
；
则
Figure M_220425100945462_462400001
。
进一步地，确定所述被测平面的空间方程还包括将各个所述已知点的像素坐标转换为三维直角坐标系内对应坐标：
所述
Figure M_220425100945493_493659001
上各个所述已知点记为：
Figure M_220425100945524_524439002
，i=1，2，
Figure M_220425100945557_557121003
，n；
所述
Figure M_220425100945572_572765001
上各个所述已知点记为：
Figure M_220425100945588_588387002
，i=1，2，
Figure M_220425100945619_619641003
，n，其中相邻的两个所述已知点之间的实际距离已知；
各个所述已知点的像素坐标分别为：
Figure M_220425100945650_650900001
，则设定在三维直角坐标系下对应的坐标为：
Figure M_220425100945682_682187002
，各个所述已知点的像素坐标与三维坐标之间的换算关系为：
Figure M_220425100945728_728998003
，
Figure M_220425100945761_761715004
，其中：
Figure M_220425100945792_792993001
为所述智能手机拍摄所得图像中每个像素的实际尺寸，所述智能手机的分辨率为
Figure M_220425100945808_808595002
。
进一步地，确定所述被测平面的空间方程还包括确定方程系数
Figure M_220425100945839_839866001
：
相机光心
Figure M_220425100945855_855488001
与
Figure M_220425100945871_871148002
的连线
Figure M_220425100945886_886724003
的方程为：
Figure M_220425100945917_917520004
，直线
Figure M_220425100945950_950194005
与
Figure M_220425100945965_965832006
的交点
Figure M_220425100945997_997081007
坐标为
Figure M_220425100946012_012691008
；
其中，
Figure M_220425100946043_043959001
，
Figure M_220425100946075_075193002
为每个
Figure M_220425100946090_090840003
所对应系数D的值，i=1，2，
Figure M_220425100946122_122067004
，n；
Figure M_220425100946138_138654001
与
Figure M_220425100946154_154351002
之间距离的估计值为
Figure M_220425100946170_170410003
；
Figure M_220425100946201_201680001
；
Figure M_220425100946248_248542001
与
Figure M_220425100946279_279791002
之间距离的实际值为
Figure M_220425100946295_295459003
，当有多组约束时，设为方差最小，求解超定方程：
Figure M_220425100946311_311038001
；
Figure M_220425100946359_359371001
；
由此求得方程系数D。
进一步地，所述拍摄计算步骤中，确定所述被测平面的空间方程后，求所述
Figure M_220425100946390_390642001
上的已知点
Figure M_220425100946421_421877002
到所述
Figure M_220425100946437_437522003
上的已知点
Figure M_220425100946453_453155004
的距离
Figure M_220425100946468_468756005
，具体包括：
Figure M_220425100946500_500023001
的空间坐标为
Figure M_220425100946515_515645002
，
Figure M_220425100946547_547101003
坐标为
Figure M_220425100946563_563976004
，其中，
Figure M_220425100946595_595232005
，则：
Figure M_220425100946610_610862001
。
进一步地，在所述准备步骤中，使用三脚架将所述智能手机架设在所述第一位置。
进一步地，在所述准备步骤中，所述已知结构体的已知结构尺寸为桥墩高和桥跨长中至少一者。
进一步地，在所述准备步骤中，所述被测物体至少占所述智能手机的视场的3/4。
进一步地，在所述拍摄计算步骤中，所述已知点的个数不少于三个。
本发明提供的一种结构位移测量方法能产生如下有益效果：
本发明提供的一种结构位移测量方法可直接利用智能手机的拍照功能对被测物体进行连续拍摄，对比计算初始状态照片和后续拍摄的照片，从而计算出待测点的像素位移；再利用智能手机的内置陀螺仪得到第一位置拍摄面与被测平面之间的三维角度关系，根据已知点之间的实际尺寸以及各个已知点在第一位置拍摄面中的像素坐标，计算出已知点的像素位移与实际位移的全场比例关系；最后根据比例关系以及待测点的像素位移得到待测点的实际位移。
相对于现有技术来说，本方法仅使用具有拍照和陀螺仪功能的智能手机即可实现结构位移测量方法，且测量数据不受被测物体尺寸以及与手机之间距离的影响与限制，数据易于获得、操作难度低，且极大程度上减小了成本，减少了测量任务准备工作的工作量，使得该方法更易于推广，尤其是在工程测量教学方面。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图；
图1为本发明实施例提供的智能手机架设在第一位置时的结构示意图；
图2为本发明实施例提供的第一位置拍摄面的示意图；
图3为本发明实施例提供的第一位置拍摄面与被测平面的相对位置示意图。
图标：1-被测平面；2-第一位置拍摄面；3-三脚架；4-待测点。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明，并不用于限制本发明。
本实施例在于提供一种结构位移测量方法，如图1至图3所示，包括：
准备步骤：将智能手机架设在第一位置，使得被测物体以及与被测物体在同一平面的已知结构尺寸的已知结构体均可在智能手机的视场中显示，使用智能手机中的陀螺仪将此时的三维角度均置为零；
三维角度获取步骤：将智能手机放置到与被测平面1平行的第二位置并记录此时的三维角度后将手机复位；
拍摄计算步骤：在已知结构体上选取多个已知点，对被测物体变形或位移发生前后进行连续拍摄图像，得到各个已知点以及待测点4的像素坐标，根据第一位置拍摄面2与被测平面1的相对三维角度、各个已知点之间的实际尺寸以及各个已知点在第一位置拍摄面2中的像素坐标，计算出各个已知点的像素位移与实际位移的比例关系，根据比例关系以及待测点4的像素位移得到待测点4的实际位移。
在传统的结构位移测量方法中，需要使用工业相机配合测距机、倾角仪等辅助设备，因此其具有受外界条件限制大、前期准备工作量大、成本高等缺点。而上述实施例中所公开的结构位移测量方法摆脱了传统工业相机配合测距机、倾角仪等辅助设备进行测量的方式，转为使用具有拍照和陀螺仪功能的智能手机，智能手机具有小巧、易获取以及便于固定等特点，测量人员操作难度低，检测过程不易受外界条件限制且数据易于获取，可极大程度上减少前期准备工作量，降低成本。
具体在操作时，首先进行准备步骤，将智能手机架设在第一位置，如图1所示，具体可以使用三脚架3实现智能手机的架设，使得被测物体以及与被测物体在同一平面的已知结构尺寸的已知结构体均可在智能手机视场中显示，为保证被测物体的清晰度，被测物至少占智能手机视场的3/4，具体可以占3/4、4/5或5/6，随后使用智能手机中的陀螺仪将此时的三维角度均置为零。
其中，以桥梁为例，如图2所示，已知结构体的已知结构尺寸的信息可以为桥跨长L1、拉索高L2、桥墩高L3、桥墩宽L4中的至少一者，当然已知结构体并不限于以上结构，检测人员可根据实际情况选择合适的物体作为已知结构体。
通过上述准备步骤可将智能手机所处第一位置的拍摄面作为基面，建立三维坐标系，便于后续的数据记录以及计算。
随后进行三维角度获取步骤，在此步骤中，将智能手机放置到与被测平面1平行的第二位置，并使用陀螺仪功能记录此时的三维角度后将手机复位至第一位置。
将智能手机放置到与被测平面1平行的第二位置具体可以通过将智能手机贴合在被测物体的被测平面1上实现，从而准确的得到第一位置拍摄面2与被测平面1的相对三维角度。
最后进行拍摄计算步骤，在此步骤中，需要在已知结构体上选取多个已知点，具体可以为两个、三个、四个、五个等；随后对被测物体变形或位移发生前后进行连续拍摄图像，得到各个已知点以及各个待测点4的像素坐标；通过已知结构尺寸的已知结构体得到各个已知点之间的距离；根据第一位置拍摄面2与被测平面1的相对三维角度、各个已知点之间的实际尺寸以及各个已知点在第一位置拍摄面2中的像素坐标，根据三维空间的几何原理，计算出各个已知点的像素位移与实际位移的比例关系；最后根据比例关系以及待测点4的像素位移得到待测点4的实际位移。
在上述步骤中，计算出已知点的像素位移与实际位移的比例关系具体操作可以如下：
首先，确定被测平面1的空间方程，计算每一个已知点在被测平面1上的位置与在第一位置拍摄面2上的位置之间的距离，在上述过程中，可以根据定角旋转公式以及成像模型空间结构关系，确定被测平面1的空间方程，再计算每一个已知点在被测平面1上的位置与在第一位置拍摄面2上的位置之间的距离。
随后，根据上述距离以及经典的小孔成像斜光轴标定算法，可以算得每个已知点的像素位移与实际位移的比例关系，上述过程以及随后的根据比例关系和待测点4的像素位移得到待测点4的实际位移均属于现有技术的计算方法，其中待测点4的像素位移可以通过变形前后的图像经过数字图像相关算法，做二维相关对比匹配得到，上述计算方法具体可以参考相关教科书或者其他相关文献。因此拍摄计算步骤中的改进在于如何仅通过智能手机测得的数据算得上述距离，本实施例中为节省篇幅，根据上述距离算得每个已知点的像素位移与实际位移的比例关系以及最终得到待测点4的实际位移不再详细赘述。
以下对如何计算每一个已知点在被测平面1上的位置与在第一位置拍摄面2上的位置之间的距离进行具体说明：
首先，可以设定第一位置拍摄面2为
Figure M_220425100946657_657728001
，被测平面1为
Figure M_220425100946673_673349002
，
Figure M_220425100946688_688987003
与
Figure M_220425100946720_720235004
的相对三维角度为
Figure M_220425100946737_737771005
，即使用右手坐标系，
Figure M_220425100946753_753538006
通过先绕
Figure M_220425100946785_785167007
轴旋转
Figure M_220425100946800_800793008
角，再绕
Figure M_220425100946832_832055009
轴旋转
Figure M_220425100946847_847675010
角，最后绕
Figure M_220425100946863_863306011
轴旋转
Figure M_220425100946894_894547012
角后，与
Figure M_220425100946925_925799013
平面的方向一致，即与
Figure M_220425100946943_943340014
平行。拍摄计算步骤中，确被测平面1的空间方程包括：
建立
Figure M_220425100946959_959479001
与
Figure M_220425100946975_975099002
的几何关系：如图3所示，以
Figure M_220425100946990_990747003
为z=0平面，
Figure M_220425100947021_021992004
的中心为坐标原点
Figure M_220425100947037_037622005
（0,0,0），智能手机的相机光心
Figure M_220425100947068_068867006
的坐标为（0,0,-
Figure M_220425100947084_084474007
），
Figure M_220425100947100_100107008
为智能手机的相机焦距，
Figure M_220425100947131_131356009
的方程设为
Figure M_220425100947147_147994010
，
Figure M_220425100947163_163637011
、
Figure M_220425100947195_195334012
和
Figure M_220425100947210_210947013
和D为方程的待求参数；
确定方程系数
Figure M_220425100947226_226566001
、
Figure M_220425100947242_242232002
和
Figure M_220425100947273_273454003
：设定
Figure M_220425100947289_289071004
的法向量
Figure M_220425100947304_304687005
为（0,0,1），设定
Figure M_220425100947337_337380006
旋转
Figure M_220425100947353_353528007
以后，新的法向量为
Figure M_220425100947369_369157008
，
Figure M_220425100947400_400435009
为
Figure M_220425100947416_416014010
的法向量，可以算得；
Figure M_220425100947447_447286001
；
Figure M_220425100947534_534149001
；
则
Figure M_220425100947565_565918001
，由此，在未知平面
Figure M_220425100947612_612811002
中仅有一个未知参数D。
在上述实施例的基础上，进一步地，确定被测平面的空间方程还包括将各个已知点的像素坐标转换为三维直角坐标系内对应坐标：
Figure M_220425100947628_628435001
上各个已知点记为：
Figure M_220425100947659_659728002
，i=1，2，
Figure M_220425100947690_690918003
，n；
Figure M_220425100947706_706560001
上各个已知点记为：
Figure M_220425100947738_738746002
，i=1，2，
Figure M_220425100947770_770517003
，n，其中相邻的两个已知点之间的实际距离已知，即各个已知点之间的实际距离是取自已知尺寸结构体的，例如已知桥长，则已知点分布在桥的两端点，此时已知点之间的距离就是桥长；
各个已知点的像素坐标分别为：
Figure M_220425100947786_786144001
，则设定在三维直角坐标系下对应的坐标为：
Figure M_220425100947817_817414002
，各个已知点的像素坐标与三维坐标之间的换算关系如下：
Figure M_220425100947864_864274003
，
Figure M_220425100947895_895514004
，其中：
Figure M_220425100947926_926776001
为智能手机拍摄所得图像中每个像素的实际尺寸，智能手机的分辨率为
Figure M_220425100947943_943339002
，即
Figure M_220425100947975_975107003
为智能手机的水平像素数，
Figure M_220425100947990_990734004
为智能手机的垂直像素数，例如智能手机的分辨率为2532x1170，则
Figure M_220425100948006_006391005
为2532，
Figure M_220425100948021_021983006
为1170。
在上述实施例的基础上，进一步地，确定被测平面的空间方程还包括确定方程系数
Figure M_220425100948053_053249001
：
如图3所示，相机光心
Figure M_220425100948068_068868001
与
Figure M_220425100948084_084467002
的连线
Figure M_220425100948115_115741003
的方程为：
Figure M_220425100948134_134290004
，直线
Figure M_220425100948166_166021005
与
Figure M_220425100948197_197306006
的交点为
Figure M_220425100948212_212920007
，
Figure M_220425100948244_244151008
坐标为
Figure M_220425100948259_259764009
；
其中，
Figure M_220425100948291_291037001
，
Figure M_220425100948322_322306002
为每个
Figure M_220425100948340_340824003
所对应系数D的值，i=1，2，
Figure M_220425100948372_372568004
，n；
Figure M_220425100948388_388210001
与
Figure M_220425100948403_403837002
之间距离的估计值为
Figure M_220425100948419_419459003
；
即
Figure M_220425100948450_450698001
；
Figure M_220425100948513_513223001
与
Figure M_220425100948545_545909002
之间距离的实际值为
Figure M_220425100948562_562022003
，即
Figure M_220425100948593_593295004
的值便为已知结构尺寸，由于
Figure M_220425100948608_608888005
的个数可以有多个，因此有多组已知尺寸的
Figure M_220425100948640_640168006
，当有多组约束时，设为方差最小，求解超定方程：
Figure M_220425100948671_671396001
；
Figure M_220425100948702_702632001
；
由此求得方程系数D的最优解。
在上述实施例的基础上，进一步地，拍摄计算步骤中，确定被测平面的空间方程后，求
Figure M_220425100948751_751476001
上的已知点
Figure M_220425100948767_767088002
到
Figure M_220425100948798_798355003
上的已知点
Figure M_220425100948829_829123004
的距离
Figure M_220425100948845_845241005
，具体包括：
Figure M_220425100948876_876473001
的空间坐标为
Figure M_220425100948892_892120002
坐标为
Figure M_220425100948923_923327003
，其中，
Figure M_220425100948959_959972004
，则：
Figure M_220425100948991_991226001
。
得到
Figure M_220425100949038_038099001
后，可使用上述距离并根据现有技术中的小孔成像斜光轴标定算法以及数字图像相关算法算得待测点4实际位移信息。
综上，上述实施例所提及的结构位移测量方法具有以下优点：
1、智能手机内置陀螺仪的普遍性使得该方法的使用成本大大减小，同时智能手机的便携性也使得测量任务的准备工作更轻松，操作难度也更低。
2、该测量方法仅用到手机中的拍照功能和手机中的内置陀螺仪测角度，这两类数据不受被测物体的尺寸、与手机之间的距离影响与限制，数据易于获得，能精确测量大结构的位移/挠度等相关参数，使得该方法更易于推广，尤其是在工程测量教学方面。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (5)

1.一种结构位移测量方法，其特征在于，包括：
准备步骤：将智能手机架设在第一位置，使得被测物体以及与所述被测物体在同一平面的已知结构尺寸的已知结构体均可在所述智能手机的视场中显示，使用所述智能手机中的陀螺仪将此时的三维角度均置为零；
三维角度获取步骤：将所述智能手机放置到与被测平面平行的第二位置并记录此时的三维角度后将手机复位；
拍摄计算步骤：在所述已知结构体上选取多个已知点，对所述被测物体变形或位移发生前后进行连续拍摄图像，得到各个所述已知点以及待测点的像素坐标，根据第一位置拍摄面与所述被测平面的相对三维角度、各个所述已知点之间的实际尺寸以及各个所述已知点在所述第一位置拍摄面中的像素坐标，计算出各个所述已知点的像素位移与实际位移的比例关系，根据所述比例关系以及所述待测点的像素位移得到所述待测点的实际位移；
所述拍摄计算步骤中，计算出所述已知点的像素位移与实际位移的比例关系包括：
确定所述被测平面的空间方程，计算每一个所述已知点在所述被测平面上的位置与在所述第一位置拍摄面上的位置之间的距离；
根据所述距离算得每个所述已知点的像素位移与实际位移的比例关系；
设定所述第一位置拍摄面为
Figure M_220609134623930_930029001
，所述被测平面为
Figure M_220609134623978_978875002
，所述
Figure M_220609134624010_010105003
与所述
Figure M_220609134624041_041321004
的相对三维角度为
Figure M_220609134624056_056982005
，所述拍摄计算步骤中，确定所述被测平面的空间方程包括：
建立所述
Figure M_220609134624119_119477001
与所述
Figure M_220609134624135_135097002
的几何关系：以所述
Figure M_220609134624166_166364003
为z=0平面，所述
Figure M_220609134624197_197596004
的中心为坐标原点
Figure M_220609134624228_228873005
（0,0,0），所述智能手机的相机光心
Figure M_220609134624259_259689006
的坐标为（0,0,-
Figure M_220609134624291_291355007
），
Figure M_220609134624338_338243008
为所述智能手机的相机焦距，所述
Figure M_220609134624355_355281009
的方程设为
Figure M_220609134624387_387051010
；
确定方程系数
Figure M_220609134624433_433963001
、
Figure M_220609134624465_465233002
和
Figure M_220609134624496_496451003
：设定所述
Figure M_220609134624543_543296004
的法向量
Figure M_220609134624759_759583005
为（0,0,1），设定
Figure M_220609134624822_822593006
旋转
Figure M_220609134624853_853853007
以后，新的法向量为
Figure M_220609134624885_885127008
，
Figure M_220609134624916_916359009
为所述
Figure M_220609134624949_949035010
的法向量；
Figure M_220609134624980_980803001
；
Figure M_220609134625058_058932001
；
则
Figure M_220609134625090_090183001
；
确定所述被测平面的空间方程还包括将各个所述已知点的像素坐标转换为三维直角坐标系内对应坐标：
所述
Figure M_220609134625137_137043001
上各个所述已知点记为：
Figure M_220609134625156_156096002
，i=1，2，
Figure M_220609134625203_203476003
，n；
所述
Figure M_220609134625234_234728001
上各个所述已知点记为：
Figure M_220609134625250_250351002
，i=1，2，
Figure M_220609134625297_297213003
，n，其中相邻的两个所述已知点之间的实际距离已知；
各个所述已知点的像素坐标分别为：
Figure M_220609134625312_312826001
，则设定在三维直角坐标系下对应的坐标为：
Figure M_220609134625363_363648002
，各个所述已知点的像素坐标与三维坐标之间的换算关系为：
Figure M_220609134625394_394869003
，
Figure M_220609134625426_426194004
，其中：
Figure M_220609134625457_457363001
为所述智能手机拍摄所得图像中每个像素的实际尺寸，所述智能手机的分辨率为
Figure M_220609134625488_488623002
；
确定所述被测平面的空间方程还包括确定方程系数
Figure M_220609134625504_504255001
：
相机光心
Figure M_220609134625535_535506001
与
Figure M_220609134625552_552090002
的连线
Figure M_220609134625583_583837003
的方程为：
Figure M_220609134625599_599462004
，直线
Figure M_220609134625630_630700005
与
Figure M_220609134625662_662010006
的交点
Figure M_220609134625693_693208007
坐标为
Figure M_220609134625708_708833008
；
其中，
Figure M_220609134625741_741984001
，
Figure M_220609134625771_771352002
为每个
Figure M_220609134625802_802585003
所对应系数D的值，i=1，2，
Figure M_220609134625818_818227004
，n；
Figure M_220609134625849_849464001
与
Figure M_220609134625865_865081002
之间距离的估计值为
Figure M_220609134625896_896396003
；
Figure M_220609134625944_944653001
；
Figure M_220609134626007_007695001
与
Figure M_220609134626023_023322002
之间距离的实际值为
Figure M_220609134626038_038923003
，当有多组约束时，设为方差最小，求解超定方程：
Figure M_220609134626070_070228001
；
Figure M_220609134626117_117068001
；
由此求得方程系数D；
所述拍摄计算步骤中，确定所述被测平面的空间方程后，求所述
Figure M_220609134626165_165428001
上的已知点
Figure M_220609134626181_181549002
到所述
Figure M_220609134626212_212745003
上的已知点
Figure M_220609134626259_259641004
的距离
Figure M_220609134626306_306479005
，具体包括：
Figure M_220609134626337_337734001
的空间坐标为
Figure M_220609134626359_359196002
，
Figure M_220609134626390_390980003
坐标为
Figure M_220609134626422_422216004
，其中，
Figure M_220609134626453_453476005
，则：
Figure M_220609134626484_484734001
。
2.根据权利要求1所述的结构位移测量方法，其特征在于，在所述准备步骤中，使用三脚架将所述智能手机架设在所述第一位置。
3.根据权利要求1所述的结构位移测量方法，其特征在于，在所述准备步骤中，所述已知结构体的已知结构尺寸为桥墩高和桥跨长中至少一者。
4.根据权利要求1所述的结构位移测量方法，其特征在于，在所述准备步骤中，所述被测物体至少占所述智能手机的视场的3/4。
5.根据权利要求1所述的结构位移测量方法，其特征在于，在所述拍摄计算步骤中，所述已知点的个数不少于三个。
CN202210441171.2A 2022-04-26 2022-04-26 一种结构位移测量方法 Active CN114543684B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210441171.2A CN114543684B (zh) 2022-04-26 2022-04-26 一种结构位移测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210441171.2A CN114543684B (zh) 2022-04-26 2022-04-26 一种结构位移测量方法

Publications (2)

Publication Number Publication Date
CN114543684A CN114543684A (zh) 2022-05-27
CN114543684B true CN114543684B (zh) 2022-07-12

Family

ID=81667497

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210441171.2A Active CN114543684B (zh) 2022-04-26 2022-04-26 一种结构位移测量方法

Country Status (1)

Country Link
CN (1) CN114543684B (zh)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL115971A (en) * 1995-11-14 1997-01-10 Razon Moshe Computer stereo vision system and method
CN106997605B (zh) * 2017-03-22 2019-11-19 浙江大学 一种通过智能手机采集脚型视频和传感器数据获取三维脚型的方法
CN109754429A (zh) * 2018-12-14 2019-05-14 东南大学 一种基于图像的桥梁结构挠度测量方法
CN109460040A (zh) * 2018-12-28 2019-03-12 珠海凯浩电子有限公司 一种通过手机拍摄照片识别地板建立地图系统及方法
CN111091595B (zh) * 2019-12-23 2023-06-02 吉林省广播电视研究所(吉林省广播电视局科技信息中心) 斜视三维测绘方法及测绘系统
CN111947578A (zh) * 2020-07-21 2020-11-17 上海同禾工程科技股份有限公司 一种结构物位移测量系统及其测量方法
CN113822942B (zh) * 2021-09-09 2023-11-17 南京中科逆熵科技有限公司 一种基于二维码的单目摄像头测量物体尺寸的方法

Also Published As

Publication number Publication date
CN114543684A (zh) 2022-05-27

Similar Documents

Publication Publication Date Title
CN110057295B (zh) 一种免像控的单目视觉平面距离测量方法
CN105486235B (zh) 一种球机视频画面中的目标度量方法
CN107588913A (zh) 一种桥梁挠度检测系统及检测方法
CN111220129B (zh) 一种带旋转云台的聚焦测量方法及终端
CN112254663B (zh) 一种基于影像识别的平面变形监测测量方法及其系统
CN110736447B (zh) 一种集成式图像采集设备竖直方向水平位置检校方法
CN102136140B (zh) 一种基于矩形图样的视频图像距离检测方法
CN106447733A (zh) 颈椎活动度及活动轴线位置的确定方法、系统及装置
CN106840010A (zh) 双目立体视觉测量装置
CN102927917A (zh) 多目铁塔视觉测量方法
CN203100724U (zh) 利用相机拍摄的照片进行测距的装置
CN112967312A (zh) 一种面向野外刚体目标的实时鲁棒位移监测方法及系统
CN205300518U (zh) 双目立体视觉测量装置
CN114543684B (zh) 一种结构位移测量方法
CN111649716A (zh) 一种基于全景图像的空间点对点距离测算方法
CN111754584A (zh) 一种远距离大视场相机参数标定系统和方法
CN113240635B (zh) 一种裂缝分辨为基准的结构物检测图像质量测试方法
CN114136544B (zh) 基于高速视频测量的水下振动模拟测试系统及测试方法
CN115290008A (zh) 一种影像测量平行光管的角度标定算法
CN111964604B (zh) 一种基于影像识别的平面变形监测测量方法
CN115018923A (zh) 一种全景环带镜头参数的检测方法及系统
JP4775541B2 (ja) 撮影画像における歪曲収差補正方法
WO2013060017A1 (zh) 臂架振动位移测量方法、系统及包含其的混凝土泵车
CN111387987A (zh) 一种基于图像识别的身高测量方法、装置、设备及存储介质
CN111780683A (zh) 便携式扫描系统及其使用方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CB03 Change of inventor or designer information

Inventor after: Tian Long

Inventor after: Yuan Yangxiang

Inventor after: Jia Geyuan

Inventor after: Wu Sihang

Inventor before: Yuan Yangxiang

Inventor before: Jia Geyuan

Inventor before: Wu Sihang

CB03 Change of inventor or designer information