设备位置信息的获得方法、装置、终端设备以及存储介质与流程

1.本发明涉及数字化测量技术领域，特别涉及一种设备位置信息的获得方法、装置、终端设备以及计算机可读存储介质。


背景技术：

2.在现代工业领域的数字化测量活动中，坐标系的配准是首先要解决的基础问题，坐标系的配准是指确定出测量坐标系与全局坐标系之间的转换关系。
3.坐标系的配准精度对于不同站位或测量设备之间坐标转换的精度起着决定性的作用。当前的坐标系配准常用的奇异值分解法、四元数法、最小二乘法等方法，来获得最终的配准参数，并利用获得的配准参数进行坐标系的配准操作。然后利用配准后的坐标系进行设备的位置信息的测量。
4.但是，采用现有的设备位置信息的获得方法，获得的位置信息的准确率较低。


技术实现要素：

5.本发明的主要目的是提供一种设备位置信息的获得方法、装置、终端设备以及计算机可读存储介质，旨在解决现有技术中采用现有的设备位置信息的获得方法，获得的位置信息的准确率较低的技术问题。
6.为实现上述目的，本发明提出一种设备位置信息的获得方法，所述方法包括以下步骤：
7.获取多个测量点在测量坐标系下的多个测量坐标值、多个所述测量点在全局坐标系下的多个名义坐标值和多个所述测量点与测量设备之间的多个测量距离，多个所述测量点为基准设备上的被标记的点；
8.利用多个所述测量距离，构建配准残差修正矩阵；
9.对多所述测量坐标值进行重心化处理，获得多个预处理测量坐标值，并对多个所述名义坐标值进行重心化处理，获得多个预处理名义坐标值；
10.利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，计算配准参数；
11.获取目标设备在所述测量坐标系的目标测量位置信息，并基于所述目标测量位置信息和所述配准参数，获得所述目标设备在所述全局坐标系下的最终测量位置信息。
12.可选的，所述利用多个所述测量距离，构建配准残差修正矩阵的步骤，包括：
13.利用每个所述测量点的测量距离，计算每个所述测量点的最大测量误差；
14.利用多个所述测量点对应的多个最大测量误差，计算多个所述测量点对应的多个坐标转换权重；
15.利用多个所述坐标转换权重，构建所述配准残差修正矩阵。
16.可选的，所述对多所述测量坐标值进行重心化处理，获得多个预处理测量坐标值的步骤，包括：
17.利用多个所述测量坐标值，计算多个所述测量点的测量重心坐标值；
18.利用所述测量重心坐标值，对多个所述测量坐标值进行坐标转换，获得多个所述预处理测量坐标值。
19.可选的，所述利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，计算配准参数的步骤之前，所述方法包括：
20.获取多个所述预处理测量坐标值与多个所述预处理名义坐标值之间的预设转换关系；
21.利用初始残差矩阵对所述预设转换关系进行误差修正，获得修正后的转换关系；
22.利用所述所述配准残差修正矩阵对所述修正后的转换关系进行再次修正，获得结果转换关系；
23.所述利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，计算配准参数的步骤，包括：
24.利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述结果转换关系，计算配准参数。
25.可选的，所述利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述结果转换关系，计算配准参数的步骤，包括：
26.利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述结果转换关系，创建待优化修正残差矩阵；
27.利用优化目标，对所述待优化修正残差矩阵进行优化，获得所述待优化修正残差矩阵对应的优化最小值；
28.基于所述优化最小值，获得所述配准参数。
29.可选的，所述基于所述优化最小值，获得所述配准参数的步骤之前，所述方法还包括：
30.获取与所述优化目标对应的等价优化目标；
31.利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，创建中间过程矩阵；
32.获取所述中间过程矩阵与奇异值分解矩阵的第一映射关系；
33.所述基于所述优化最小值，获得所述配准参数的步骤，包括：
34.利用所述等价优化目标、所述中间过程矩阵、所述优化最小值、所述第一映射关系和所述奇异值分解矩阵，计算所述配准参数。
35.可选的，所述奇异值分解矩阵包括第一正交矩阵、第二正交矩阵和对角矩阵；所述配准参数包括旋转矩阵和平移矩阵；所述利用所述等价优化目标、所述中间过程矩阵、所述优化最小值、所述第一映射关系和所述奇异值分解矩阵，计算所述配准参数的步骤，包括：
36.获取所述第一正交矩阵和所述第二正交矩阵与所述旋转矩阵的第二映射关系；
37.基于所述等价优化目标、所述最小值、所述中间过程矩阵、所述第一映射关系、所述第二映射关系、所述正交矩阵、所述第二正交矩阵和所述对角矩阵，获得所述旋转矩阵；
38.基于所述旋转矩阵、所述预设转换关系、多个所述测量坐标值和多个所述名义坐标值，获得所述平移矩阵。
39.此外，为实现上述目的，本发明还提出了一种设备位置信息的获得装置，所述装置
包括：
40.获取模块，用于获取多个测量点在测量坐标系下的多个测量坐标值、多个所述测量点在全局坐标系下的多个名义坐标值和多个所述测量点与测量设备之间的多个测量距离，多个所述测量点为基准设备上的被标记的点；
41.构建模块，用于利用多个所述测量距离，构建配准残差修正矩阵；
42.获得模块，用于对多所述测量坐标值进行重心化处理，获得多个预处理测量坐标值，并对多个所述名义坐标值进行重心化处理，获得多个预处理名义坐标值；
43.计算模块，用于利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，计算配准参数；
44.测量模块，用于获取目标设备在所述测量坐标系的目标测量位置信息，并基于所述目标测量位置信息和所述配准参数，获得所述目标设备在所述全局坐标系下的最终测量位置信息。
45.此外，为实现上述目的，本发明还提出了一种终端设备，所述终端设备包括：存储器、处理器及存储在所述存储器上并在所述处理器上运行设备位置信息的获得程序，所述设备位置信息的获得程序被所述处理器执行时实现如上述任一项所述的设备位置信息的获得方法的步骤。
46.此外，为实现上述目的，本发明还提出了一种计算机可读存储介质，所述计算机可读存储介质上存储有设备位置信息的获得程序，所述设备位置信息的获得程序被处理器执行时实现如上述任一项所述的设备位置信息的获得方法的步骤。
47.本发明技术方案提出了一种设备位置信息的获得方法，所述方法包括以下步骤：获取多个测量点在测量坐标系下的多个测量坐标值、多个所述测量点在全局坐标系下的多个名义坐标值和多个所述测量点与测量设备之间的多个测量距离，多个所述测量点为基准设备上的被标记的点；利用多个所述测量距离，构建配准残差修正矩阵；对多所述测量坐标值进行重心化处理，获得多个预处理测量坐标值，并对多个所述名义坐标值进行重心化处理，获得多个预处理名义坐标值；利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，计算配准参数；获取目标设备在所述测量坐标系的目标测量位置信息，并基于所述目标测量位置信息和所述配准参数，获得所述目标设备在所述全局坐标系下的最终测量位置信息。
48.现有的坐标系配准方法，获得的最终转换关系未引入测量点与测量设备之间的距离对应的测量误差，从而使得配准参数的准确率较低，导致目标设备的位置信息转换准确率较低，目标设备在全局坐标系下的最终测量位置信息准确率较低。而利用本发明的方法，引入了测量点与测量设备之间的距离对应的测量误差，利用该测量误差构建的配准残差修正矩阵进行最终配准参数的计算，从而提高了配准参数的准确率，使得目标设备的位置信息转换准确率较高，达到了提高目标设备在全局坐标系下的最终测量位置信息准确率的技术效果。
附图说明
49.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本
发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图示出的结构获得其他的附图。
50.图1为本发明实施例方案涉及的硬件运行环境的终端设备结构示意图；
51.图2为本发明设备位置信息的获得方法第一实施例的流程示意图；
52.图3为本发明设备位置信息的获得装置第一实施例的结构框图。
53.本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
54.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
55.参照图1，图1为本发明实施例方案涉及的硬件运行环境的终端设备结构示意图。
56.通常，终端设备包括：至少一个处理器301、存储器302以及存储在所述存储器上并可在所述处理器上运行的设备位置信息的获得程序，所述设备位置信息的获得程序配置为实现如前所述的设备位置信息的获得方法的步骤。
57.处理器301可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器301可以采用dsp(digital signal processing，数字信号处理)、 fpga(field－programmable gate array，现场可编程门阵列)、pla(programmablelogic array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器301也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称cpu(central processingunit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器301 可以在集成有gpu(graphics processing unit，图像处理器)，gpu用于负责显示屏所需要显示的内容的渲染和绘制。处理器301还可以包括ai(artificialintelligence，人工智能)处理器，该ai处理器用于处理有关设备位置信息的获得方法操作，使得设备位置信息的获得方法模型可以自主训练学习，提高效率和准确度。
58.存储器302可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器302还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器302中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器301所执行以实现本技术中方法实施例提供的设备位置信息的获得方法。
59.在一些实施例中，终端还可选包括有：通信接口303和至少一个外围设备。处理器301、存储器302和通信接口303之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与通信接口303相连。具体地，外围设备包括：射频电路304、显示屏305和电源306中的至少一种。
60.通信接口303可被用于将i/o(input/output，输入/输出)相关的至少一个外围设备连接到处理器301和存储器302。在一些实施例中，处理器301、存储器302 和通信接口303被集成在同一芯片或电路板上；在一些其他实施例中，处理器301、存储器302和通信接口303中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。
61.射频电路304用于接收和发射rf(radio frequency，射频)信号，也称电磁信号。射频电路304通过电磁信号与通信网络以及其他通信设备进行通信。射频电路304将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路304包括：天线系统、rf收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路304可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于：城域网、各代移动通信网络(2g、3g、4g及5g)、无线局域网和/或wifi(wireless fidelity，无线保真)网络。在一些实施例中，射频电路304 还可以包括nfc(near field communication，近距离无线通信)有关的电路，本技术对此不加以限定。
62.显示屏305用于显示ui(user interface，用户界面)。该ui可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏305是触摸显示屏时，显示屏305 还具有采集在显示屏305的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器301进行处理。此时，显示屏305还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏305 可以为一个，电子设备的前面板；在另一些实施例中，显示屏305可以为至少两个，分别设置在电子设备的不同表面或呈折叠设计；在再一些实施例中，显示屏305可以是柔性显示屏，设置在电子设备的弯曲表面上或折叠面上。甚至，显示屏305还可以设置成非矩形的不规则图形，也即异形屏。显示屏305可以采用 lcd(liquidcrystal display，液晶显示屏)、oled(organic light-emitting diode, 有机发光二极管)等材质制备。
63.电源306用于为电子设备中的各个组件进行供电。电源306可以是交流电、直流电、一次性电池或可充电电池。当电源306包括可充电电池时，该可充电电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。本领域技术人员可以理解，图1中示出的结构并不构成对终端设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
64.此外，本发明实施例还提出一种计算机可读存储介质，所述计算机可读存储介质上存储有设备位置信息的获得程序，所述设备位置信息的获得程序被处理器执行时实现如上文所述的设备位置信息的获得方法的步骤。因此，这里将不再进行赘述。另外，对采用相同方法的有益效果描述，也不再进行赘述。对于本技术所涉及的计算机可读存储介质实施例中未披露的技术细节，请参照本技术方法实施例的描述。确定为示例，程序指令可被部署为在一个终端设备上执行，或者在位于一个地点的多个终端设备上执行，又或者，在分布在多个地点且通过通信网络互连的多个终端设备备上执行。
65.本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，上述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，上述的计算机可读存储介质可为磁碟、光盘、只读存储记忆体(read-onlymemory，rom)或随机存储记忆体(random accessmemory，ram)等。
66.基于上述硬件结构，提出本发明设备位置信息的获得方法的实施例。
67.参照图2，图2为本发明设备位置信息的获得方法第一实施例的流程示意图，所述方法用于终端设备，所述方法包括以下步骤：
68.步骤s11：获取多个测量点在测量坐标系下的多个测量坐标值、多个所述测量点在全局坐标系下的多个名义坐标值和多个所述测量点与测量设备之间的多个测量距离，多个所述测量点为基准设备上的被标记的点。
69.基准设备作为多个测量点的实体依照物，可以是任何类型的设备，在基准设备上设定多个标记点，一个标记点即为一个待测量点。
70.需要说明的是，本发明的执行主体是终端设备，终端设备安装有设备位置信息的获得程序，终端设备执行设备位置信息的获得程序时，实现本发明的设备位置信息的获得方法的步骤。
71.具体的，在本发明中，用户可以基于需求设定多个测量点，例如n个测量点，n为大于1的自然数，多个测量点构成测量点集。对于测量点集中的每个测量点，利用测量设备测量其在测量坐标系下的坐标值——测量坐标值，并测量其在全局坐标系下的坐标值(也就是理论坐标值)——名义坐标值，同时，测量其与测量设备之间距离——测量距离，一个测量点对应一个测量坐标值、一个名义坐标值和一个测量距离，即，对于n个测量点对应n个测量坐标值、n 个名义坐标值和n个测量距离。
72.其中，测量设备可以是lecia at930激光跟踪仪等。
73.步骤s12：利用多个所述测量距离，构建配准残差修正矩阵。
74.具体的，所述利用多个所述测量距离，构建配准残差修正矩阵的步骤，包括：利用每个所述测量点的测量距离，计算每个所述测量点的最大测量误差；利用多个所述测量点对应的多个最大测量误差，计算多个所述测量点对应的多个坐标转换权重；利用多个所述坐标转换权重，构建所述配准残差修正矩阵。
75.在实际应用中，对于任意一个测量点i的测量坐标值p
′i(测量点i的名义坐标值为pi)，其对应的最大测量误差按照设备制造商提供的经验公式进行计算，经验公式如下：
76.δ
′i＝a+b
·di
77.其中，δ
′i为p
′i对应的最大测量误差，di为测量点i与测量设备之间的测量距离，a与b为设备生产商提供的预设参数，例如测量设备为lecia at930激光跟踪仪时，a＝15μm，b＝6μm/m。
78.按照上述经验公式，求得n个测量点对应的n个最大测量误差，表示如下：
79.δ
′
＝(δ
′1,δ
′2,δ
′3,
…
,δ
′n)
80.其中，δ
′
为n个测量点对应的n个最大测量误差组成的测量误差集。
81.然后，继续基于上述结构，利公式一求解坐标转换权重，一个测量点对应一个坐标转换权重，公式一如下：
[0082][0083]
其中，λi为测量点i的坐标转换权重。然后，利用坐标转换权重构建配准残差修正矩阵，构建方式参照公式二，公式二如下：
[0084][0085]
其中，λ为所述配准残差修正矩阵。
[0086]
在该实施例中，提出了一种配准残差修正矩阵的构建方式，配准残差修正矩阵考虑到了测量点与测试设备之间的测试距离，使得配准残差修正矩阵考虑了测量误差，配准残差修正矩阵的准确率较高。
[0087]
步骤s13：对多所述测量坐标值进行重心化处理，获得多个预处理测量坐标值，并对多个所述名义坐标值进行重心化处理，获得多个预处理名义坐标值。
[0088]
需要对每个测量坐标值和每个名义坐标值进行转换，获得对应的预处理测量坐标值和预处理名义坐标值。
[0089]
具体的，所述对多所述测量坐标值进行重心化处理，获得多个预处理测量坐标值的步骤，包括：利用多个所述测量坐标值，计算多个所述测量点的测量重心坐标值；利用所述测量重心坐标值，对多个所述测量坐标值进行坐标转换，获得多个所述预处理测量坐标值。
[0090]
同理，对多个所述名义坐标值进行重心化处理，获得多个预处理名义坐标值的步骤，包括：利用多个所述名义坐标值，计算多个所述测量点的名义重心坐标值；利用所述名义重心坐标值，对多个所述名义坐标值进行坐标转换，获得多个所述预处理名义坐标值。
[0091]
具体的，利用公式三求得名义重心坐标值和测量重心坐标值，公式三如下：
[0092][0093]
其中，pc为所述名义重心坐标值，p
′c为所述测量重心坐标值。
[0094]
然后继续利用公式四，求得预处理名义坐标值和预处理测量坐标值，公式四如下：
[0095][0096]
其中，为测量点i的预处理名义坐标值，为测量点i的预处理测量坐标值。
[0097]
步骤s14：利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，计算配准参数。
[0098]
其中，配准参数包括求得的已知的旋转矩阵和平移矩阵，对于任意一个测量点，配准参数对具有以下固有关系：
[0099]
pi＝rp
′i+t
[0100]
其中，r为旋转矩阵，t为平移矩阵。本发明的目的即是求得具体的r和 t。
[0101]
但是，实际应用中，存在制造误差和测量误差等误差，上述固有关系具体会被修正为新的固有关系，如下：
[0102]
pi＝rp
′i+t+ε
[0103]
其中，ε为残差矩阵(下文所述初始残差矩阵)，用于修正各种误差。
[0104]
进一步的，所述利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，计算配准参数的步骤之前，所述方法包括：获取多个所述预处理测量坐标值与多个所述预处理名义坐标值之间的预设转换关系；利用初始残差矩阵对所述预设转换关系进行误差修正，获得修正后的转换关系；利用所述所述配准残差修正矩阵对所述修正后的转换关系进行再次修正，获得结果转换关系；所述利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，计算配准参数的步骤，包括：利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述结果转换关系，计算配准参数。
[0105]
对于预处理名义坐标值和预处理测量坐标值存在关系——预设转换关系，预设转换关系表示为公式五，公式五如下：
[0106][0107]
与上述新的固有关系同理，此时，对应的预处理名义坐标值和预处理测量坐标值存在关系表示为公式六，公式六如下：
[0108][0109]
公式六即为利用初始残差矩阵对所述预设转换关系进行误差修正，获得的修正后的转换关系。
[0110]
然后进一步利用公式七进行修正，公式七如下：
[0111][0112][0113]
其中，即为所述结果转换关系。
[0114]
此时，利用这个结果转换关系、多个预处理测量坐标值和多个预处理名义坐标值，计算配准参数。
[0115]
在该实施例中，对初始残差矩阵进行修正，获得结转换关系，使得结果转换关系引入了具体的测量设备与测量点之间距离的误差，结果转换关系的准确度较高。
[0116]
进一步的，所述利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述结果转换关系，计算配准参数的步骤，包括：利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述结果转换关系，创建待优化修正残差矩阵；利用优化目标，对所述待优化修正残差矩阵进行优化，获得所述待优化修正残差矩阵对应的优化最小值；基于所述优化最小值，获得所述配准参数。
[0117]
其中，所述基于所述优化最小值，获得所述配准参数的步骤之前，所述方法还包括：获取与所述优化目标对应的等价优化目标；利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，创建中间过程矩阵；获取所述中间过程矩阵
与奇异值分解矩阵的第一映射关系；所述基于所述优化最小值，获得所述配准参数的步骤，包括：利用所述等价优化目标、所述中间过程矩阵、所述优化最小值、所述第一映射关系和所述奇异值分解矩阵，计算所述配准参数。
[0118]
将上述已知的λ、和带入结果转换关系式，即可获得待优化修正残差矩阵。在本发明中，以ε
′
的2-范数作为优化目标，优化目标如下：
[0119][0120]
其中，ε
′
ix
、ε
′
iy
和ε
′
iz
为待优化修正残差矩阵中测量点i对应的x、y和z三个方向的分量。该优化目标输出的是一个最小值，即优化最小值，然后继续利用优化最小值，获得所述配准参数。
[0121]
同时，优化目标对应的等价优化目标如下：
[0122]
max(trance(rh))
[0123]
其中，trance(rh)为rh的迹，h为中间过程矩阵，中间过程矩阵按照公式八求解，公式八如下：
[0124][0125]
同时，中间过程矩阵与奇异值分解矩阵的第一映射关系如下：
[0126]
h＝q
t
λv
[0127]
其中，奇异值分解矩阵包括第一正交矩阵、第二正交矩阵和对角矩阵，q 为第一正交矩阵，v第二正交矩阵，λ对角矩阵。
[0128]
此时优化目标和等价优化目标的可以如下表示：
[0129][0130]
利用ε
′
的2-范数作为优化目标，使得获得的优化最小值的准确率和精确度均较高，使得利用优化最小值获得的配准参数准确率较高。从而进一步实现了提高配准精准度的技术效果。
[0131]
进一步的，所述利用所述等价优化目标、所述中间过程矩阵、所述优化最小值、所述第一映射关系和所述奇异值分解矩阵，计算所述配准参数的步骤，包括：获取所述第一正交矩阵和所述第二正交矩阵与所述旋转矩阵的第二映射关系；基于所述等价优化目标、所述最小值、所述中间过程矩阵、所述第一映射关系、所述第二映射关系、所述正交矩阵、所述第二正交矩阵和所述对角矩阵，获得所述旋转矩阵；基于所述旋转矩阵、所述预设转换关系、多个所述测量坐标值和多个所述名义坐标值，获得所述平移矩阵。
[0132]
具体的，所述第一正交矩阵和所述第二正交矩阵与所述旋转矩阵的第二映射关系如下：
[0133]
r＝v
tq[0134]
利用优化目标、等价优化目标、公式八、第一映射关系、第二映射关系、 q、v和λ，
进行求解，获得具有具体值的旋转矩阵。
[0135]
此时，继续基于求得的旋转矩阵，利用公式九进行平移矩阵的求解，公式九如下：
[0136]
t＝p-rp
′
[0137]
其中，p为全部测量点的名义坐标值组成的集合，p
′
为全部测量点的测量坐标值组成的集合。
[0138]
至此，已经求得旋转矩阵r和平移矩阵t，也即求得配准参数。
[0139]
步骤s15：利获取目标设备在所述测量坐标系的目标测量位置信息，并基于所述目标测量位置信息和所述配准参数，获得所述目标设备在所述全局坐标系下的最终测量位置信息。
[0140]
此时，利用已经求得旋转矩阵r和平移矩阵t，对所述全局坐标系和所述测量坐标系进行配准操作。然后，利用配准后的坐标系进行数据测量。
[0141]
目标设备即是进行位置信息后获得的设备，例如大型制造机器。目标设备在测量坐标系下，由测量设备测得位置信息——目标测量位置信息，然后利用配准后的全局坐标系和测量坐标系对应的配准参数，对目标测量位置信息进行转换，获得目标设备在全局坐标系下的位置信息——最终测量位置信息。
[0142]
配准操作实际就是利用旋转矩阵r和平移矩阵t构建全局坐标系和测量坐标系的映射关系，以使在实际情况中，直接利用该映射关系进行坐标转换(位置信息的转换)操作。
[0143]
为了验证本发明的方法的优势，进行如下仿真：
[0144]
表1为4个测量点的名义坐标值，表1如下：
[0145]
表1
[0146]
测量点x(mm)y(mm)z(mm)p102000p21000200200p33002000250p4200040003000
[0147]
表2为表1的测量点对应的测量坐标值，表2如下：
[0148]
测量点x(mm)y(mm)z(mm)p1’
278.1712764-450.3486783715.6983118p2’
1003.551311-1126.296289477.1220501p3’
1729.624244663.8269817504.9926376p4’
5262.993078631.94660341907.414945
[0149]
其中，表2中的各个坐标值是对表1中的各个坐标值，按照六自由度参数 ([10
°
,25
°
,-45
°
,150,-600,750])转换获得。
[0150]
在测量过程中使用lecia at960激光跟踪仪在站位点(0,0,0)对以上点进行测量，并进行站位转化。由于激光跟踪仪测量误差的存在，各点的坐标值在实际测量过程中的测量值可以表示为：
[0151]
p
修
＝p+e
[0152]
其中，e服从均值为0，标准差为δ
′
的正态分布：
[0153]
e～n(0,δ
′
)
[0154]
δ
′
＝15μm+6dm/m
[0155]
按照此种方式，对于表1中的名义坐标值和表2中的测量坐标值进行转换，获得修正后的坐标值，如表3：
[0156]
表3
[0157] x(mm)误差范围(mm)y(mm)误差范围(mm)z(mm)误差范围(mm)p10[-16.2μm,16.2μm]200[-16.2μm,16.2μm]0[-16.2μm,16.2μm]p21000[-21.2μm,21.2μm]200[-21.2μm,21.2μm]200[-21.2μm,21.2μm]p3300[-27.2μm,27.2μm]2000[-27.2μm,27.2μm]250[-27.2μm,27.2μm]p42000[-47.3μm,47.3μm]4000[-47.3μm,47.3μm]3000[-47.3μm,47.3μm]p1’
278.1737[-20.3μm,20.3μm]-450.345[-20.3μm,20.3μm]715.6908[-20.3μm,20.3μm]p2’
1003.549[-23.5μm,24.5μm]-1126.28[-23.5μm,24.5μm]477.1154[-23.5μm,24.5μm]p3’
1729.608[-26.5μm,26.5μm]663.8235[-26.5μm,26.5μm]504.9732[-26.5μm,26.5μm]p4’
5262.993[-48.8μm,48.8μm]631.9333[-48.8μm,48.8μm]1907.442[-48.8μm,48.8μm]
[0158]
分别利用修正后的测量坐标值和修正后的名义坐标值，基于误差范围，进行蒙特卡洛仿真，得到随机数据列。基于仿真后的数据分别利用传统的svd方法和本专利方法进行六自由度参数的求解，结果如表4，表4如下：
[0159]
表4
[0160]
方法α(
°
)β(
°
)γ(
°
)t
x
(mm)ty(mm)tz(mm)svd9.996824.9945-45.0102149.9828-599.8926749.7854本发明9.999924.9997-45.0001149.9954-599.9957749.9888
[0161]
将真值与计算结果对比，可得配准误差分别为表5，表5如下：
[0162]
表5
[0163]
方法δα(
°
)δβ(
°
)δγ(
°
)δt
x
(mm)δty(mm)δtz(mm)svd0.00320.00550.01020.0172-0.10740.2146本发明10-4
0.00030.00010.0046-0.00430.0112
[0164]
可见，本发明的方法的准确度和精度均较高。
[0165]
本发明技术方案提出了一种设备位置信息的获得方法，所述方法包括以下步骤：获取多个测量点在测量坐标系下的多个测量坐标值、多个所述测量点在全局坐标系下的多个名义坐标值和多个所述测量点与测量设备之间的多个测量距离；利用多个所述测量距离，构建配准残差修正矩阵；对多所述测量坐标值进行重心化处理，获得多个预处理测量坐标值，并对多个所述名义坐标值进行重心化处理，获得多个预处理名义坐标值；利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，计算配准参数；获取目标设备在所述测量坐标系的目标测量位置信息，并基于所述目标测量位置信息和所述配准参数，获得所述目标设备在所述全局坐标系下的最终测量位置信息。
[0166]
现有的坐标系配准方法，获得的最终转换关系未引入测量点与测量设备之间的距离对应的测量误差，从而使得配准参数的准确率较低，导致目标设备的位置信息转换准确率较低，目标设备在全局坐标系下的最终测量位置信息准确率较低。而利用本发明的方法，引入了测量点与测量设备之间的距离对应的测量误差，利用该测量误差构建的配准残差修正矩阵进行最终配准参数的计算，从而提高了配准参数的准确率，使得目标设备的位置信息转换准确率较高，达到了提高目标设备在全局坐标系下的最终测量位置信息准确率的技
术效果。
[0167]
参照图3，图3为本发明设备位置信息的获得装置第一实施例的结构框图，所述装置用于终端设备，基于与前述实施例相同的发明构思，所述装置包括：
[0168]
获取模块10，用于获取多个测量点在测量坐标系下的多个测量坐标值、多个所述测量点在全局坐标系下的多个名义坐标值和多个所述测量点与测量设备之间的多个测量距离，多个所述测量点为基准设备上的被标记的点；
[0169]
构建模块20，用于利用多个所述测量距离，构建配准残差修正矩阵；
[0170]
获得模块30，用于对多所述测量坐标值进行重心化处理，获得多个预处理测量坐标值，并对多个所述名义坐标值进行重心化处理，获得多个预处理名义坐标值；
[0171]
计算模块40，用于利用多个所述预处理测量坐标值、多个所述预处理名义坐标值和所述配准残差修正矩阵，计算配准参数；
[0172]
测量模块50，用于获取目标设备在所述测量坐标系的目标测量位置信息，并基于所述目标测量位置信息和所述配准参数，获得所述目标设备在所述全局坐标系下的最终测量位置信息。
[0173]
需要说明的是，由于本实施例的装置所执行的步骤与前述方法实施例的步骤相同，其具体的实施方式以及可以达到的技术效果都可参照前述实施例，这里不再赘述。
[0174]
以上所述仅为本发明的可选实施例，并非因此限制本发明的专利范围，凡是在本发明的发明构思下，利用本发明说明书及附图内容所作的等效结构变换，或直接/间接运用在其他相关的技术领域均包括在本发明的专利保护范围内。