CN108919177B - 一种基于虚拟信源估计与轨迹校正的定位地图构建方法 - Google Patents
一种基于虚拟信源估计与轨迹校正的定位地图构建方法 Download PDFInfo
- Publication number
- CN108919177B CN108919177B CN201810781054.4A CN201810781054A CN108919177B CN 108919177 B CN108919177 B CN 108919177B CN 201810781054 A CN201810781054 A CN 201810781054A CN 108919177 B CN108919177 B CN 108919177B
- Authority
- CN
- China
- Prior art keywords
- positioning
- fingerprint
- sample
- signal
- sub
- 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
Links
- 238000012937 correction Methods 0.000 title claims abstract description 17
- 238000010276 construction Methods 0.000 title claims abstract description 10
- 238000000034 method Methods 0.000 claims abstract description 48
- 239000002245 particle Substances 0.000 claims description 27
- 238000005070 sampling Methods 0.000 claims description 19
- 230000001133 acceleration Effects 0.000 claims description 16
- 238000012216 screening Methods 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 230000001902 propagating effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000012549 training Methods 0.000 description 6
- 230000033001 locomotion Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000001186 cumulative effect Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000012952 Resampling Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 241000764238 Isis Species 0.000 description 1
- 244000062793 Sorghum vulgare Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 235000019713 millet Nutrition 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/0252—Radio frequency fingerprinting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/005—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 with correlation of navigation data from several sources, e.g. map or contour matching
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
- G01C21/206—Instruments for performing navigational calculations specially adapted for indoor navigation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/021—Services related to particular areas, e.g. point of interest [POI] services, venue services or geofences
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/025—Services making use of location information using location based information parameters
- H04W4/026—Services making use of location information using location based information parameters using orientation information, e.g. compass
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/025—Services making use of location information using location based information parameters
- H04W4/027—Services making use of location information using location based information parameters using movement velocity, acceleration information
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/029—Location-based management or tracking services
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/30—Services specially adapted for particular environments, situations or purposes
- H04W4/33—Services specially adapted for particular environments, situations or purposes for indoor environments, e.g. buildings
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于虚拟信源估计与轨迹校正的定位地图构建方法,包括:在定位场景中设置N个信号源,划分P个子区域,利用众包样本数据重构用户粗行走轨迹,获得用户粗行走轨迹中的定位样本,计算定位样本中的指纹中每一信号强度的权值;利用筛选后的定位样本建立指纹库。同时已构建的指纹库又可在用户轨迹重建阶段对用户轨迹进行校正,进一步提高众包样本数据的准确性。通过在每个子区域中,针对每一个信号源,通过虚拟信源估计方法来重构信号源在子区域内的虚拟信号传播曲面实现。本发明能够解决现有定位地图构建方法中,众包样本使用不合理,引致指纹库构建不准确,定位精度低下,且构建的指纹库臃肿,占用大量空间的问题。
Description
技术领域
本发明属于通信和无线网络技术领域,更具体地,涉及一种基于虚拟信源估计与轨迹校正的定位地图构建方法。
背景技术
随着信息时代的爆炸发展,人们对于定位信息的需求也日益增加。由于GPS定位技术基于卫星信号定位的技术,其应用场景被仅仅限制于了没有信号阻挡的室外。面对大型的室内环境,如机场、商场、学校,人们有着非常迫切的高精度室内定位需求。近些年来室内定位技术被广泛研究,诞生出了包括WiFi定位、蓝牙定位、超宽带定位、红外定位等的多种室内定位技术。伴随城市的智能化发展,城市内大多场景都实现了WLAN覆盖。基于WiFi技术的室内定位由于无需额外布置定位信标,用户无需额外购置定位硬件设备,受到广泛青睐,并在近年成为了室内定位研究的重要对象。
目前在国内外,常用的WiFi定位技术可分为两种:基于传播模型的定位技术和基于指纹的定位技术。前者通过利用在接收机所接收的到达时间(TOA,Time of Arrival)、到达往返时间(TDOA,Time Difference of Arrival)或者信号强度(RSSI,Received SignalStrength Indicator),从而推算出目标和信号源之间的距离。当信号源存在多个且已知信源物理位置时,即可通过三角定位得到目标的大致位置。但由于室内环境的复杂多变性,使得普适传播模型不可能精确得到,且在室内家具、物体的遮蔽与多径反射作用下导致该技术定位精度不是很高。而后者则基于在室内不同位置上所接受到的信号源信号强度具有唯一性,从而称作信号强度(RSS,Received Signal Strength)指纹定位。该技术可分为两个阶段,一是训练阶段,二是动态定位阶段。在训练阶段,需要在室内预设参考点,通常为均匀分布的栅格点,并在每个点处进行长时间的指纹采样,在采样完毕后存储到指纹数据库中;在动态定位阶段,将当前接收到的指纹与指纹数据库中的指纹进行比较,选择最近的指纹位置得出精确的目标位置。
两种方法对比下,前者通常需要加装特殊的WLAN模块,而后者所需要的硬件设备则相对简单(只需日常所用的智能手机即可),且定位精度高。因此指纹定位方法成为了一种最常用的非基于测距的室内定位方法。但其代价是需要耗费大量人力物力人工采样数据点,并且指纹数据库存储量巨大。因此最近有相关研究提出利用大量用户所收集的WLAN信号强度数据,也称作众包数据,进行更加快捷方便地指纹数据库构建。
这些研究方法通过收集用户的移动轨迹,并实时采集移动轨迹上设备接收到的WLAN信号强度数据。用户移动轨迹可通过实时的WLAN室内定位方法获取,或通过用户手机惯性传感器采集的运动信息计算而得。然而WLAN室内定位结果误差较大,用户使用手机姿态多变,给移动轨迹的精确获取带来了巨大的挑战。同时已有的指纹库构建方法所构建的指纹库臃肿复杂,不便于实时的室内定位。如何能精确获取用户移动轨迹,收集众包数据,并构建合适的指纹库成为了一个很重要的研究课题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种基于虚拟信源估计与轨迹校正的定位地图构建方法,其目的在于通过虚拟信源参数估计构建准确的轻量化指纹库并能据此获取精确用户移动轨迹,由此解决现有定位技术中获取用户轨迹粗糙,构建指纹库臃肿且误差较大的技术问题。
为实现上述目的,本发明提供了一种基于虚拟信源估计与轨迹校正的定位地图构建方法,包括:
(1)在定位场景中设置N个信号源,划分P个子区域,采集定位场景中用户行走轨迹中每一时刻的加速度值、航向角和N个信号源的信号强度指纹,构成众包样本数据;
(2)利用众包样本数据重构用户粗行走轨迹,获得用户粗行走轨迹中每一步的定位样本,定位样本包括位置坐标和指纹,计算定位样本中的指纹中每一信号强度的权值;
(3)在构建指纹库前,利用位置坐标将定位样本分配至子区域,在构建指纹库后,则利用评估的每个信号源在各子区域的定位重要度和定位样本的指纹对定位样本进行区域定位,得到定位样本所在的子区域,利用定位样本所在的子区域的信号源的虚拟信源传播曲面拟合函数、估计的虚拟信源传播参数的确定度和定位样本中的指纹中每一信号强度的权值进行指纹定位,得到定位样本所在子区域内的位置坐标;
(4)利用熵值法评估每个信号源在各子区域的定位重要度,利用各子区域中定位样本的物理空间簇和信号空间簇之间的交集计算定位样本的权值;
(5)筛选出权值大于等于权值阈值的定位样本,通过定位样本的指纹,在各子区域内对每个信号源建立一个虚拟信源传播曲面拟合函数,得到估计的虚拟信源传播参数,计算估计的虚拟信源传播参数的确定度;
(6)将各子区域内所有信号源的虚拟信源传播曲面函数通过等距栅格采样,构建各子区域的区域指纹,利用各子区域内每个信号源的估计的虚拟信源传播参数和各子区域的区域指纹构建定位场景的指纹库;
(7)对于定位场景中除步骤(1)之外的其他用户行走轨迹,采集其他用户行走轨迹中每一时刻的加速度值、航向角和N个信号源的信号强度指纹,更新众包样本数据,然后执行步骤(2)-(6),更新定位场景的指纹库。
进一步地,步骤(2)包括:
利用众包样本数据中的加速度值计算用户行走过的步数T′与步长lt,将步数T′、步长lt与航向角θt结合,通过步行者航位推算算法重构用户粗行走轨迹,获得用户粗行走轨迹中每一步的定位样本定位样本包括位置坐标和指纹计算定位样本中的指纹中每一信号强度的权值。
进一步地,定位样本中的指纹中每一信号强度的权值为:
进一步地,区域定位的具体实现方式为:
进一步地,定位样本所在的子区域的位置坐标为:
其中,为定位样本所在的子区域的位置坐标,φn(PC,β,xi,yi,x0,y0,h0)为定位样本所在的子区域的第n个信号源的虚拟信源传播曲面拟合函数,为定位样本所在的子区域的第n个信号源的估计的虚拟信源传播参数的确定度,θtn为第t步定位样本所在的子区域从第n个信号源收到的信号强度的权值,(x,y)为定位样本所在的子区域的任意坐标。
进一步地,步骤(3)还包括:利用粒子滤波算法对用户粗行走轨迹进行校正。
进一步地,步骤(5)包括:
当定位样本的权值大于等于权值阈值时,通过定位样本的指纹,在各子区域内对每个信号源建立一个虚拟信源传播曲面拟合函数,利用信赖域映射算法,当子区域中第n个信号源的虚拟信源传播曲面拟合函数φn(PC,β,xi,yi,x0,y0,h0)与第i步定位样本的指纹中接收到的第n个信号源的信号强度rin之间的平方差最小时,得到估计的虚拟信源传播参数,即:
其中,ωi为第i步定位样本的权值,θin为第i步定位样本的指纹中接收到的第n个信号源的信号强度的权值,M′j表示定位样本的权值大于等于权值阈值时的所有定位样本。
进一步地,计算估计的虚拟信源传播参数的确定度包括:
令第k个虚拟信源的估计的虚拟信源传播参数PC′,β′,h0′不变,将第k个虚拟信源置于任意位置(x0,y0)时的代价函数Hk(x0,y0)为:
其中,ωi为第i步定位样本的权值,θi为第i步定位样本的指纹中接收到的信号源的信号强度的权值,φn(PC′,β′,xi,yi,x0,y0,h0′)为第k个虚拟信源的虚拟信源传播曲面拟合函数,rin为第i步定位样本的指纹中接收到的第n个信号源的信号强度,S′为采样点集合;
代价函数Hk(x0,y0)在(x0′,y0′)处取得最小值,代价函数Hk(x0,y0)在(x0′,y0′)处的斜率平均值为:
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
(1)本发明利用众包样本数据重构用户粗行走轨迹,获得用户粗行走轨迹中每一步的定位样本,定位样本包括位置坐标和指纹,计算定位样本中的指纹中每一信号强度的权值;将由于多径效应造成分布不均的信号强度,依据其对定位的重要性重新进行了加权分配,减少了大量无用样本对定位结果的影响。
(2)本发明对通过用户众包数据生成的粗行走轨迹进行了校正,减少了众包样本的定位误差,使得在后续步骤中构建的指纹库更加精确;由于利用各子区域中定位样本的物理空间簇和信号空间簇之间的交集计算定位样本的权值,评估了带噪音的众包数据的真实性,选取了较合理且具有代表性的数据赋予较高的权值,同样使得构建的指纹数据库更加准确,提高了定位精度。
(3)本发明在各子区域内对每个信号源建立一个虚拟信源传播曲面拟合函数,得到估计的虚拟信源传播参数,据此参数可以计算出子空间内信号源在任意位置的信号强度的分布,将子空间内所有信号源的虚拟信源传播曲面参数汇总即可得到子空间的指纹库。此指纹库只需少量虚拟信源参数即可构建,且能获得任意位置点的指纹,使得构建的指纹库不仅轻量化,且存储信息丰富。
(4)本发明构建了信号源的虚拟信源传播曲面,更加准确地反映了信号强度在子空间内的分布情况,更加有效地预测子区域中所有位置的指纹信息,极大地提高了指纹库的精确度。本发明评估了估计的虚拟信源参数的确定度,为指纹库中信号源的重要性划分了层次,进一步提高了指纹库的精确度。
(5)本发明通过先前构建好的指纹库进行指纹定位,对后续上传的用户轨迹数据进行粒子滤波校正,能获得高精度的用户轨迹校正结果。精确的用户轨迹结果在很多场合,譬如轨迹记录、用户行为分析等,都具有很强的实用性。
附图说明
图1是本发明实施例提供的一种基于虚拟信源估计与轨迹校正的定位地图构建方法的流程图;
图2是本发明实施例提供的定位样本权值确定算法的流程图;
图3是本发明实施例提供的粒子群算法的流程图;
图4是本发明实施例1提供的的定位场景图;
图5是本发明实施例1提供的三种类似的曲面拟合定位方法栅格曲面定位法、原始曲面定位法、聚类曲面定位法以及本发明选择样本曲面定位方法在图4的A房间内建立的关于信号源3的信号强度拟合曲面;
图6本发明实施例1提供的所有房间使用不同的众包样本数据量时,近邻栅格定位、栅格曲面定位、原始曲面定位、聚类曲面定位及选择样本曲面定位的定位效果对比;
图7本发明实施例1提供的当所有房间的众包样本数据包含不同大小的噪声干扰时,近邻栅格定位、栅格曲面定位、原始曲面定位、聚类曲面定位及选择样本曲面定位的定位效果对比;
图8本发明实施例1提供的所有房间的众包样本数据量和包含的噪声干扰固定时,近邻栅格定位、栅格曲面定位、原始曲面定位、聚类曲面定位及选择样本曲面定位的定位误差累积分布示意图;
图9本发明实施例1提供的在建立指纹库后,用户静态站立收集样本指纹后,利用现有指纹库分别进行原始栅格定位、曲线拟合定位、带权曲线拟合定位、虚拟信源估计定位及带权虚拟信源估计定位后的定位误差累积分布示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示,本发明提供了一种基于虚拟信源估计与轨迹校正的定位地图构建方法,包括:
(1)在定位场景中设置N个信号源,划分P个子区域,采集定位场景中用户行走轨迹中每一时刻的加速度值、航向角和N个信号源的信号强度指纹,构成众包样本数据;
(2)利用众包样本数据重构用户粗行走轨迹,获得用户粗行走轨迹中每一步的定位样本,定位样本包括位置坐标和指纹,计算定位样本中的指纹中每一信号强度的权值;
(3)在构建指纹库前,利用位置坐标将定位样本分配至子区域,在构建指纹库后,则利用评估的每个信号源在各子区域的定位重要度和定位样本的指纹对定位样本进行区域定位,得到定位样本所在的子区域,利用定位样本所在的子区域的信号源的虚拟信源传播曲面拟合函数、估计的虚拟信源传播参数的确定度和定位样本中的指纹中每一信号强度的权值进行指纹定位,得到定位样本所在子区域内的位置坐标;
(4)利用熵值法评估每个信号源在各子区域的定位重要度,利用各子区域中定位样本的物理空间簇和信号空间簇之间的交集计算定位样本的权值;
(5)筛选出权值大于等于权值阈值的定位样本,通过定位样本的指纹,在各子区域内对每个信号源建立一个虚拟信源传播曲面拟合函数,得到估计的虚拟信源传播参数,计算估计的虚拟信源传播参数的确定度;
利用各子区域内每个信号源的估计的虚拟信源传播参数和各子区域的区域指纹构建定位场景的指纹库;
(7)对于定位场景中除步骤(1)之外的其他用户行走轨迹,采集其他用户行走轨迹中每一时刻的加速度值、航向角和N个信号源的信号强度指纹,更新众包样本数据,然后执行步骤(2)-(6),更新定位场景的指纹库。
进一步地,步骤(1)包括:
在定位场景中设置N个信号源,保证定位场景中任何位置都能使终端设备接收到来自至少一个信号源的信号,并将给定的定位场景按照是否有明显墙壁遮挡划分为P个开阔子区域。采集定位场景中用户行走轨迹中每一时刻的加速度值、航向角和N个信号源的信号强度指纹,构成众包样本数据;令用户手持终端设备在定位场景中任意走动,可以随意穿越多个开阔子区域,不限制用户的运动。以用户开始行走时为0时刻,终端设备以一定的采样频率f采集数据。令表示终端设备在时刻t接收到的原始样本,其中分别表示终端设备传感器在时刻t接收到的总加速度值与航向角和接收到的N个信号源的信号强度指纹。rtn表示时刻t接收到的指纹中从第n个信号源接收到的信号强度大小,其中n=1,2,...,N,t=1,2,...,T,T为用户一段完整行走轨迹中收集的数据总数,M,P,N都是自然数。
进一步地,步骤(2)包括:
利用众包样本数据中的加速度值计算用户行走过的步数T′与步长lt,将步数T′、步长lt与航向角θt结合,通过步行者航位推算算法重构用户粗行走轨迹,获得用户粗行走轨迹中每一步的定位样本定位样本包括位置坐标和指纹计算定位样本中的指纹中每一信号强度的权值。
具体地,由于多径效应,接收到的信号强度弱的样本较多,而接收到信号强度强的样本较少,因此样本信号强度越强,其定位贡献度越高,因此其权值应更高。对于中接收到的每个信号源信号强度值,通过下式计算信号强度权值:
进一步地,区域定位的具体实现方式为:
进一步地,定位样本所在的子区域的位置坐标为:
其中,为定位样本所在的子区域的位置坐标,为定位样本所在的子区域的第n个信号源的虚拟信源传播曲面拟合函数,为定位样本所在的子区域的第n个信号源的估计的虚拟信源传播参数的确定度,θtn为第t步定位样本所在的子区域从第n个信号源收到的信号强度的权值,(x,y)为定位样本所在的子区域的任意坐标。
具体地,利用熵值法评估每个信号源在各子区域的定位重要度包含以下子步骤:
用表示子区域Aj中信号源Sk所接收到的所有信号强度值,其中j=1,2,...,P,k=1,2,...,Nj。表示在子区域Aj中从信号源Sk接收到的所有信号强度值的总数,通过下式将信号强度值归一化处理:
利用熵值函数计算该信号源的信号强度分布情况:
依据熵信,计算下列两个权信:
本步骤的优点在于通过熵值法评估了每个区域中每个信号源所提供的信号强度值包含的信息量,从而决定了每个信号源在定位时的重要度,而不是将所有信号源一视同仁。该步骤提高了在区域粗定位阶段时的精度。
如图2所示,随机挑选K个样本作为K个物理空间簇中心点,依据样本之间的物理空间距离将样本划分到这K个中心点所在的簇中;依据样本之间信号空间距离将样本同样K个中心点所在的簇中,最终得到K个物理空间簇与K个信号空间簇
本步骤的优点在于:通过评估带噪音的众包数据的真实性,选取了较合理且具有代表性的数据赋予较高的权值,使得构建的指纹库更加准确,提高了定位精度。
所述根据样本权值进行样本筛选包括:
设定百分比阈值pth与权值阈值ωth,其中pth,ωth∈[0,1]。将所有样本权值从小到大排列成找出在百分比为pth处的权值大小ωk。若ωk≥ωth时不需要筛选样本,否则将中前个样本删除。选地pth=50%,ωth=0.25。
如图3所示,所述粒子滤波算法包含如下子步骤:
在用户初始位置附近生成M个粒了粒子位置为以用户初始位置为均值,方差为σp的高斯分布随机变量,粒子初始行进角度为以用户初始行进角度为均值,方差为σang的高斯分布随机变量,设置为均匀概率优选地,σp=1.5m,σang=2°。
根据行走步长lt与航向角θt,通过下式更新每个粒子所处的位置:
在用户每一步行进过程中,利用获取指纹位置计算上一步位置(xt-1,yt-1)与指纹位置间的距离与角度考虑指纹位置的不可靠性,仅当且才认为与均比较合理且指纹位置可用于粒子滤波,否则计算粒子的质心,作为用户当前步的校正位置(xt,yt):
优选地,lmin=0.3m,lmax=2m,θmin=-120°,θmax=120°。
计算粒子的质心得到计算上一步位置(xt-1,yt-1)与间的距离与角度当时,该校正位置属于用户可达位置,校正结果较可靠,可保留为最终校正位置,否则利用角度与PDR估测的行走步长更新用户位置(xt,yt):
优选地,ρ=0.8。
粒子筛选是在粒子重采样后,对权值较低的样本剔除并重生成的过程。本例中将概率向粒子全部删除,并替换为位置为以用户初始位置为均值,方差为σp的高斯分布随机变量,初始行进角度为以用户初始行进角度为均值,方差为σang的高斯分布随机变量的新粒子。由于概率将会在下次重采样时重新计算,因此不重新赋予新概率。
进一步地,步骤(5)包括:
当定位样本的权值大于等于权值阈值时,通过定位样本的指纹,在各子区域内对每个信号源建立一个虚拟信源传播曲面拟合函数,利用信赖域映射算法,当子区域中第n个信号源的虚拟信源传播曲面拟合函数φn(PC,β,xi,yi,x0,y0,h0)与第i步定位样本的指纹中接收到的第n个信号源的信号强度rin之间的平方差最小时,得到估计的虚拟信源传播参数,即:
其中,ωi为第i步定位样本的权值,θin为第i步定位样本的指纹中接收到的第n个信号源的信号强度的权值,M′j表示定位样本的权值大于等于权值阈值时的所有定位样本。
进一步地,计算估计的虚拟信源传播参数的确定度包括:
令第k个虚拟信源的估计的虚拟信源传播参数PC′,β′,h0′不变,将第k个虚拟信源置于任意位置(x0,y0)时的代价函数Hk(x0,y0)为:
其中,ωi为第i步定位样本的权值,θi为第i步定位样本的指纹中接收到的信号源的信号强度的权值,φn(PC′,β′,xi,yi,x0,y0,h0′)为第k个虚拟信源的虚拟信源传播曲面拟合函数,rin为第i步定位样本的指纹中接收到的第n个信号源的信号强度,S′为采样点集合;
代价函数Hk(x0,y0)在(x0′,y0′)处取得最小值,代价函数Hk(x0,y0)在(x0′,y0′)处的斜率平均值为:
本步骤的优点在于考虑到了众包数据的不可靠性,评估了通过众包数据进行虚拟信源参数估计的确定性,而不是将所得到的所有虚拟信源都认为同等正确,为虚拟信源划分了层次,一定程度上提高了指纹库的精确性。
实施例1
如图4所示,本发明实施例1的场景根据有无明显墙壁遮挡,可划分为3个子区域,包括两端的直走廊与目标区域中心的正方形区域。其中中心正方形区域内为空心栏杆并非实心墙壁,因此可划分为一开阔子区域。其中两端直走廊长29.7m,宽2.4m,中心正方形区域长15.6m,宽18m。在场景内可测得N=492个信号源,该场景中任意位置均可接收到至少一个信号源的信号强度。每个子区域内划分的物理空间簇和信号空间簇的数目分别为49(两端直走廊)与42(中心正方形区域)。
使用小米6智能手机作为终端进行信号强度的测量。采样数据分为两大部分,一类为动态行走采样数据,另一类为静态站立采样数据。动态行走采样数据用于实时构建指纹库,共行走采集了5种用户轨迹,如图4走廊上的行走轨迹所示。其中黑色实线轨迹采集了5条,作为指纹库建立所用训练轨迹,灰色虚线轨迹采集了9条,作为测试指纹库建立效果的测试轨迹。静态站立采样数据用于测试指纹库的准确性,在场景中以1m为间隔均匀站立采样,共采集了602个静态点数据。
本应用实例中,共设置了五种实验方案进行测试。
(a)PDR(Pedestrian Dead Reckoning,步行者航位推算算法)粗行走轨迹。仅根据采集的用户行走时的加速度及航向角重构的用户粗行走轨迹。
(b)PDR+曲线拟合。根据采集的用户行走时的加速度及航向角重构的用户粗行走轨迹,不设定样本权值与信号源重要度,利用曲线拟合算法拟合虚拟信源的信号强度分布与到虚拟信源距离的关系生成指纹库,并用粒子滤波算法校正用户粗行走轨迹。
(c)PDR+带权曲线拟合。根据采集的用户行走时的加速度及航向角重构的用户粗行走轨迹,按照上述实施方式设定权值,利用曲线拟合算法拟合虚拟信源的信号强度分布与到虚拟信源距离的关系生成指纹库,并用粒子滤波算法校正用户粗行走轨迹。
(d)PDR+虚拟信源估计。根据采集的用户行走时的加速度及航向角重构的用户粗行走轨迹,不设定样本权值与信号源重要度,利用上述虚拟信源估计算法估计虚拟信源参数生成指纹库,并用粒子滤波算法校正用户粗行走轨迹。
(e)PDR+带权虚拟信源估计。根据采集的用户行走时的加速度及航向角重构的用户粗行走轨迹,按照上述实施方式设定权值,利用上述虚拟信源估计算法估计虚拟信源参数生成指纹库,并用粒子滤波算法校正用户粗行走轨迹。
图5示出了在使用训练轨迹建立指纹库过程中,18条训练轨迹在通过粒子滤波算法校正后的定位误差结果随训练轨迹序号的变化关系。其中PDR粗行走轨迹的平均定位误差为2.427m,PDR+曲线拟合的平均定位误差为2.144m,PDR+带权曲线拟合的为2.027m,PDR+虚拟信源估计的为1.667m,PDR+带权虚拟信源估计的为1.563m。从图上可以看出PDR粗行走轨迹定位误差基本为最高,而其他四种定位算法的定位误差依次下降,本发明提出的带权虚拟信源估计方法校正后的轨迹误差基本最小。
图6示出在指纹库建立后,利用现有指纹库对测试轨迹进行校正后的平均定位误差随轨迹序号的变化关系。其中PDR粗行走轨迹的平均定位误差为1.437m,PDR+曲线拟合的平均定位误差为1.4663m,PDR+带权曲线拟合的为1.332m,PDR+虚拟信源估计的为1.211m,PDR+带权虚拟信源估计的为1.138m。从图上可以看出本发明提出的带权虚拟信源估计方法校正后的轨迹误差普遍最小。
图7示出了在建立指纹库过程中18条轨迹的PDR粗行走轨迹,及分别经过PDR+曲线拟合、PDR+带权曲线拟合、PDR+虚拟信源估计、PDR+带权虚拟信源估计方法校正后的定位误差结果累积分布示意图。
通过该对比结果可以看出,本发明提出的带权虚拟信源估计方法校正后的轨迹误差最小。
图8示出在建立指纹库后,测试轨迹的PDR粗行走轨迹,及利用现有指纹库分别经过PDR+曲线拟合、PDR+带权曲线拟合、PDR+虚拟信源估计、PDR+带权虚拟信源估计方法进行校正后的累积分布示意图。
通过该对比结果可以看出,本发明提出的指纹库建立方法构建的指纹库对用户粗行走轨迹的校正效果显著,取得了最小的定位误差,也即构建的指纹库更加精准。
图9示出在建立指纹库后,用户静态站立收集样本指纹后,利用现有指纹库分别进行原始栅格定位、曲线拟合定位、带权曲线拟合定位、虚拟信源估计定位及带权虚拟信源估计定位后的定位误差累积分布示意图。
通过该对比结果可以看出,本发明提出的指纹库建立方法构建的指纹库的指纹定位结果取得了最小的定位误差,也即构建的指纹库更加精准。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于虚拟信源估计与轨迹校正的定位地图构建方法,其特征在于,包括:
(1)在定位场景中设置N个信号源,划分P个子区域,采集定位场景中用户行走轨迹中每一时刻的加速度值、航向角和N个信号源的信号强度指纹,构成众包样本数据;
(2)利用众包样本数据重构用户粗行走轨迹,获得用户粗行走轨迹中每一步的定位样本,定位样本包括位置坐标和指纹,计算定位样本中的指纹中每一信号强度的权值;
(3)在构建指纹库前,利用位置坐标将定位样本分配至子区域,在构建指纹库后,则利用评估的每个信号源在各子区域的定位重要度和定位样本的指纹对定位样本进行区域定位,得到定位样本所在的子区域,利用定位样本所在的子区域的信号源的虚拟信源传播曲面拟合函数、估计的虚拟信源传播参数的确定度和定位样本中的指纹中每一信号强度的权值进行指纹定位,得到定位样本在所在子区域内的位置坐标;
(4)利用熵值法评估每个信号源在各子区域的定位重要度,利用各子区域中定位样本的物理空间簇和信号空间簇之间的交集计算定位样本的权值;
(5)筛选出权值大于等于权值阈值的定位样本,通过筛选出的定位样本的指纹,在各子区域内对每个信号源建立一个虚拟信源传播曲面拟合函数,得到估计的虚拟信源传播参数,计算估计的虚拟信源传播参数的确定度;
(6)将各子区域内所有信号源的虚拟信源传播曲面拟合函数通过等距栅格采样,构建各子区域的区域指纹,利用各子区域内每个信号源的估计的虚拟信源传播参数和各子区域的区域指纹构建定位场景的指纹库;
(7)对于定位场景中除步骤(1)之外的其他用户行走轨迹,采集其他用户行走轨迹中每一时刻的加速度值、航向角和N个信号源的信号强度指纹,更新众包样本数据,然后执行步骤(2)-(6),更新定位场景的指纹库。
6.如权利要求1或2所述的一种基于虚拟信源估计与轨迹校正的定位地图构建方法,其特征在于,所述步骤(3)还包括:利用粒子滤波算法对用户粗行走轨迹进行校正。
7.如权利要求1或2所述的一种基于虚拟信源估计与轨迹校正的定位地图构建方法,其特征在于,所述步骤(5)包括:
当定位样本的权值大于等于权值阈值时,通过定位样本的指纹,在各子区域内对每个信号源建立一个虚拟信源传播曲面拟合函数,利用信赖域映射算法,当子区域中第n个信号源的虚拟信源传播曲面拟合函数φn(PC,β,xi,yi,x0,y0,h0)与第i步定位样本的指纹中接收到的第n个信号源的信号强度rin之间的平方差最小时,得到估计的虚拟信源传播参数,即:
其中,ωi为第i步定位样本的权值,θin为第i步定位样本的指纹中接收到的第n个信号源的信号强度的权值,M′j表示定位样本的权值大于等于权值阈值时的所有定位样本。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810781054.4A CN108919177B (zh) | 2018-07-16 | 2018-07-16 | 一种基于虚拟信源估计与轨迹校正的定位地图构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810781054.4A CN108919177B (zh) | 2018-07-16 | 2018-07-16 | 一种基于虚拟信源估计与轨迹校正的定位地图构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108919177A CN108919177A (zh) | 2018-11-30 |
CN108919177B true CN108919177B (zh) | 2020-12-18 |
Family
ID=64414087
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810781054.4A Active CN108919177B (zh) | 2018-07-16 | 2018-07-16 | 一种基于虚拟信源估计与轨迹校正的定位地图构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108919177B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110234072B (zh) * | 2019-06-28 | 2020-12-11 | 徐州工业职业技术学院 | 基于指纹量化的改进加权质心定位方法 |
CN110366100B (zh) * | 2019-07-17 | 2020-11-10 | 京信通信系统(中国)有限公司 | 终端的定位方法、定位装置、可读存储介质及终端设备 |
CN110967667A (zh) * | 2019-11-19 | 2020-04-07 | 重庆大学 | 一种基于众包Wi-Fi指纹定位的室内轨迹采集方法 |
CN113079466B (zh) * | 2020-10-21 | 2022-04-26 | 中移(上海)信息通信科技有限公司 | 指纹数据库的构建方法、装置、设备及计算机存储介质 |
CN112197771B (zh) * | 2020-12-07 | 2021-03-19 | 深圳腾视科技有限公司 | 车辆失效轨迹重构方法、设备以及存储介质 |
CN113810850B (zh) * | 2021-08-11 | 2022-12-13 | 荣耀终端有限公司 | 一种定位数据库的创建方法及电子设备 |
CN115451937A (zh) * | 2022-08-12 | 2022-12-09 | 合肥未来计算机技术开发有限公司 | 复杂环境下安全通道的电子地图生成系统及方法 |
CN115406435B (zh) * | 2022-08-24 | 2024-07-23 | 同济大学 | 基于wlan和mems的室内电子地图构建方法、装置及存储介质 |
CN116405880B (zh) * | 2023-05-31 | 2023-09-12 | 湖北国际贸易集团有限公司 | 一种基于联邦学习的无线电地图构建方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105636201A (zh) * | 2016-03-14 | 2016-06-01 | 中国人民解放军国防科学技术大学 | 基于稀疏信号指纹数据库的室内定位方法 |
CN106792465A (zh) * | 2016-12-26 | 2017-05-31 | 华中科技大学 | 一种基于众包指纹的室内指纹地图构建方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9357346B2 (en) * | 2012-07-16 | 2016-05-31 | Aisle411, Inc. | System and method for indoor location mapping and tracking |
US20160371394A1 (en) * | 2015-06-22 | 2016-12-22 | The Governing Council Of The University Of Toronto | Indoor localization using crowdsourced data |
-
2018
- 2018-07-16 CN CN201810781054.4A patent/CN108919177B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105636201A (zh) * | 2016-03-14 | 2016-06-01 | 中国人民解放军国防科学技术大学 | 基于稀疏信号指纹数据库的室内定位方法 |
CN106792465A (zh) * | 2016-12-26 | 2017-05-31 | 华中科技大学 | 一种基于众包指纹的室内指纹地图构建方法 |
Non-Patent Citations (3)
Title |
---|
FinCCM: Fingerprint Crowdsourcing, Clustering and Matching for Indoor Subarea Localization;Qiuyun Chen 等;《IEEE WIRELESS COMMUNICATIONS LETTERS》;20151231;第4卷(第6期);全文 * |
RMapTAFA:Radio Map Construction Based on Trajectory Adjustment and Fingerprint Amendment;Wensong Li 等;《IEEE Access》;20190123;第7卷;全文 * |
基于指纹众包的室内定位技术研究;陈秋云;《中国优秀硕士学位论文全文数据库 信息科技辑》;20180115(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108919177A (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108919177B (zh) | 一种基于虚拟信源估计与轨迹校正的定位地图构建方法 | |
CN108989976B (zh) | 一种智慧教室内的指纹定位方法及系统 | |
CN108534779B (zh) | 一种基于轨迹纠正和指纹改进的室内定位地图构建方法 | |
CN104796866B (zh) | 室内定位方法和装置 | |
CN103747524B (zh) | 一种基于云平台的Android终端室内定位方法 | |
CN104635203B (zh) | 一种基于粒子滤波算法的无线电干扰源测向定位方法 | |
CN104394588B (zh) | 基于Wi‑Fi指纹和多维尺度分析的室内定位方法 | |
CN110933599B (zh) | 一种融合uwb与wifi指纹的自适应定位方法 | |
CN105717483B (zh) | 一种基于多源定位方式的位置确定方法及装置 | |
CN109379711B (zh) | 一种定位方法 | |
CN103533647A (zh) | 一种基于分簇机制及稳健回归的射频地图自适应定位方法 | |
CN111698774B (zh) | 基于多源信息融合的室内定位方法及装置 | |
CN109059919B (zh) | 一种基于众包样本加权曲面拟合的室内定位方法 | |
CN109348416B (zh) | 基于二分k均值的指纹室内定位方法 | |
CN108450060A (zh) | 基于wi-fi接入点的定位方法、设备 | |
CN114095853B (zh) | 一种生成室内地图的方法和装置 | |
CN105898711B (zh) | 一种基于地磁基准线的定位方法及装置 | |
He et al. | WiFi iLocate: WiFi based indoor localization for smartphone | |
CN115100376A (zh) | 一种基于改进反距离插值法的电磁频谱地图实现方法 | |
CN104101861B (zh) | 测距定位方法及系统 | |
CN105866729B (zh) | 一种基于用户行为特征的室内定位方法及装置 | |
Chua Ching et al. | Mobile indoor positioning using Wi-Fi localization and image processing | |
Ja’afar et al. | Analysis of indoor location and positioning via Wi-Fi signals at FKEKK, UTeM | |
CN105792230B (zh) | 基于rss指纹室内定位系统的接入点和定标点部署方法 | |
Wang et al. | An advanced algorithm for Fingerprint Localization based on Kalman Filter |
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 |