CN104463929B - 基于图像边缘检测信号相关性的室内wlan信号地图绘制与映射方法 - Google Patents
基于图像边缘检测信号相关性的室内wlan信号地图绘制与映射方法 Download PDFInfo
- Publication number
- CN104463929B CN104463929B CN201410783008.XA CN201410783008A CN104463929B CN 104463929 B CN104463929 B CN 104463929B CN 201410783008 A CN201410783008 A CN 201410783008A CN 104463929 B CN104463929 B CN 104463929B
- Authority
- CN
- China
- Prior art keywords
- dis
- point
- coordinate
- vector
- mapping
- 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
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法,该方法首先利用随机用户采集的接收信号强度RSS序列,通过谱聚类构建各RSS序列的聚类图;其次,利用图像边缘检测法,构建随机用户在定位目标区域内的信号逻辑图;再者,根据相应的映射准则,建立信号逻辑图中的RSS聚类节点到物理环境图中的区域位置节点之间的映射;最终利用信号逻辑图到物理环境图的映射关系,实现对目标用户的位置估计,同时利用绘图技术对信号逻辑图及物理环境图进行绘制,提高了图的可读性,使得信号逻辑图及物理环境图中各节点的连接关系更加明晰。
Description
技术领域
本发明涉及一种信号地图绘制与映射方法,特别涉及一种基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法。
背景技术
随着无线局域网的广泛部署,利用WLAN(Wireless Local Area Network)基础设施进行定位可以大大降低定位设施和系统维护的投入。目前,基于接收信号强度的室内WLAN定位技术受到了广泛的研究,其中典型的基于RSS(Received Signal Strength)位置指纹的定位系统有微软研究院开发的RADAR系统和马里兰大学提出的Horus系统。RADAR系统在2000年由微软公司首次提出,其基本思想是在离线阶段从目标区域内选定一定数目的参考点,并在每个参考点处测量来自不同接入点AP(Access Point)的信号强度值,建立位置指纹数据库;而在在线阶段,则利用欧式距离最小准则,将终端实时测量得到的信号值与位置指纹数据库中已保存的信号数据进行匹配,找到与终端实时测量信号值的RSS欧氏距离较小的若干个位置指纹,并利用这若干个位置指纹所对应的位置坐标,对终端进行位置估计。RADAR定位系统所采用的定位算法简单且易于实现,但定位精度不高,其最大定位误差超过20m。Horus系统沿用了位置指纹的思想,同时又对RADAR定位系统进行了改进,提出了基于位置指纹的概率定位算法,以及利用最强AP信号聚类等优化措施,使得系统的定位性能得到了很大程度的提高,其最大定位误差在5m以内。
位置指纹定位算法需要在离线阶段对所有参考点处的来自不同AP的信号强度进行采集,从而消耗了大量的人力和物力资源,并且当参考点数目较多时,离线阶段所需的人力和时间开销会大大增加,为此,基于SLAM(Simultaneous Localization and Mapping)的WLAN定位系统应运而生,这类系统利用了其他方式获得的传感器信息,如微机电系统MEMS(Micro Electronic Machine System)中的惯性测量单元IMU(Inertial MeasurementUnit)等,将用户当前接收的WLAN信号强度与传感器采集的信息相结合,从而在不需要大量位置指纹信息采集的情况下,实现对终端的位置估计。
由于位置指纹定位算法在离线阶段需要对每个参考点处来自不同AP的信号强度进行采集,因此,在系统设计前期需要投入大量的人力与时间开销,以完成位置指纹数据库的构建。同时,基于SLAM的WLAN定位算法需要利用传感器采集的信息,所以对终端设备有特殊的要求,并且在大多数基于位置的服务LBS(Location-based Service)中,我们无需计算目标的精确物理位置坐标,而仅需获知目标所在的大致区域及其周边物理环境信息。针对这一目标,本发明提出一种基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法。
发明内容
有鉴于此,本发明的目的在于提供一种基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法,该方法首先利用随机用户采集的信号强度序列,通过聚类及图像边缘检测方法,构建用户在定位目标区域内的信号逻辑图;其次,根据相应的映射准则,建立信号逻辑图中的RSS聚类节点到物理环境图中的区域位置节点之间的映射;最终,利用最优信号逻辑图到物理环境图的映射关系,实现对终端的位置估计,同时利用绘图(GraphDrawing)技术,对信号逻辑图及物理环境图进行绘制,提高图的可读性,使得信号逻辑图及物理环境图中各节点的连接关系更加明晰。
为达到上述目的,本发明提供如下技术方案:
一种基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法,该方法包括以下步骤:
步骤一:用户在定位目标区域内随机采集若干条WLAN接收信号强度序列为其中,Mj为第j条序列的序列长度,z为序列条数;RSSjk=(RSSjk1,…,RSSjkn)(k=1,…,Mj)为第j条序列中第k个接收信号强度矢量;n为AP个数;RSSjkr(r=1,…,n)为第j条序列第k个接收信号强度矢量中第r个AP的信号强度值;
步骤二:根据时间戳顺序对Sj中接收信号强度矢量进行升序排列;将第j条序列中第k个接收信号强度矢量重构为新的n+1维矢量
步骤三:对中所包含的时间戳信息和接收信号强度信息进行加权,得到混合矢量其中加权系数为wts和wrss,wts+wrss=1;
步骤四:对混合矢量进行谱聚类,得到Sj中每条混合矢量的聚类号;
步骤五:通过中值滤波方法,修正每条混合矢量的聚类号及相应类心;
步骤六:根据Sj中相邻聚类之间的转移关系,以连接图的形式得到Sj的类转移图;
步骤七:重复步骤2)至步骤6),得到所有接收信号强度序列的类转移图;
步骤八:通过图像边缘检测技术,确定所有类转移图中类间距离小于门限std的类,合并相应的类,对所有离散的类转移图进行拼接,得到待筛选的信号逻辑图,
步骤九:将定位目标区域内的每个物理叉路口作为区域边界进行子区域划分,并对每个子区域进行序号标记;记区域标号为1,…,A,其中A表示所有的子区域的个数;
步骤十:根据各子区域的物理邻接关系,将定位目标区域表示为一副由各子区域连通的物理环境图;
步骤十一:在定位目标区域内选择少量标记位置点CP(Calibration Point),且保证标记位置点的个数少于子区域个数;
步骤十二:在各标记位置点处采集NR个来自不同AP的信号强度矢量,并将其均值矢量作为各标记位置点的代表矢量RV(Representative Vector);
步骤十三:在每个逻辑图中,计算与RV中每个元素Mrssnc与逻辑图中各类心的欧式距离其中,表示Mrssnc与第zj个逻辑图中的第f3个类心之间的欧氏距离,为第zj个逻辑图中的第f3个类心;选择与Mrssnc欧式距离最小的类心所对应的逻辑节点为该RV中元素Mrssnc所对应子区域cpnc存在的映射关系,剔除包含与2个或以上不同子区域存在映射关系的逻辑节点所对应的信号逻辑图;
步骤十四:根据步骤十三得到所有未剔除信号逻辑图,利用映射准则,得到所有未剔除的信号逻辑图,以及相应的与物理环境图的映射关系;
步骤十五:选择对于所有标记点具有最高平均定位精度的信号逻辑图,作为最优信号逻辑图,对于一个信号逻辑图到物理环境图的映射而言,每个标记点的定位精度定义为在该标记点上采集的正确映射到所属子区域的信号强度矢量个数与信号强度矢量总数的比值,其中,表示第x个逻辑图的映射关系中,第nc个标记点上的定位精度,为第nc个标记点上采集的NR个接收信号强度矢量中能够映射到cpnc的接收信号强度矢量的个数;所有标记点的平均定位精度
步骤十六:根据步骤十五,选择最大的逻辑图作为最优信号逻辑图,根据步骤十四所确定的映射准则,可得该最优信号逻辑图到物理环境图的映射关系 表示第cet个逻辑图,也即最优逻辑图中各逻辑节点所对应的物理区域号,Nodcet表示第cet个逻辑图中的逻辑节点个数;
步骤十七:利用Graph Drawing正交算法,对最优信号逻辑图及物理环境图进行绘制;
步骤十八:根据终端新采集的信号强度矢量,在最优信号逻辑图中,计算得到与新采集信号强度矢量具有最小欧式距离的逻辑节点;
步骤十九:根据最优信号逻辑图到物理环境图的映射关系,估计终端所在的子区域。
本发明的有益效果在于:本发明提供一种基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法利用随机用户采集的信号强度序列进行谱聚类、基于图像边缘检测的序列拼接、信号逻辑图筛选、映射构建及利用Graph Drawing对信号逻辑图和物理环境图进行绘制,最终实现对用户的定位。该方法无需在离线阶段采集每个参考点处来自不同AP的信号强度值,只需随机用户采集若干条时间标记的信号强度序列,且无需使用惯性测量单元对传感器信息进行采集,解决了WLAN位置指纹定位方法中离线阶段需要大量人力及时间开销进行位置指纹采集的问题,该方法提高了图的可读性,使得信号逻辑图及物理环境图中各节点的连接关系更加明晰,具有较高的区域定位精度。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步的详细描述,其中:
图1为本发明所述方法的流程图;
图2为接收信号强度矢量进行时间戳标记的过程;
图3为接收信号强度序列经谱聚类、中值滤波及类号重排处理后所得到聚类号和类转移图;
图4(a)为本发明采集的21条接收信号强度序列所形成的类转移图;
图4(b)为利用Graph Drawing重构得到的21条接收信号强度序列所形成的类转移图;
图5(a)为目标区域所对应的物理连接图;
图5(b)为利用Graph Drawing重构得到的目标区域所对应的物理连接图;
图6为目标区域的物理结构及21条测试路径;
图7(a)为原始图像;
图7(b)为图7(a)经中值滤波处理后得到的图像;
图7(c)为图7(b)经腐蚀处理后得到的图像;
图8为利用Graph Drawing得到的最优逻辑图;
图9为21条接收信号强度序列的RSS取值分布情况;
图10为在不同聚类方法条件下的子区域正确定位概率;
图11为在不同聚类方法条件下的子区域平均正确定位概率;
图12为层次聚类后所形成的聚类树。
具体实施方式
下面将结合附图,对本发明的优选实施例进行详细的描述。
本发明提供的一种基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法,如图1所示,该方法包括以下步骤:
步骤一:用户在定位目标区域内随机采集若干条WLAN接收信号强度序列其中,Mj为第j条序列的序列长度,即接收信号强度矢量个数;z为序列条数;RSSjk=(RSSjk1,…,RSSjkn)(k=1,…,Mj)为第j条序列中第k个接收信号强度矢量;n为AP个数;RSSjkr(r=1,…,n)为第j条序列第k个接收信号强度矢量中来自第r个AP的信号强度值;
步骤二:在Sj中,利用时间戳顺序对接收信号强度矢量进行升序排列,其中,将第j条序列中第k个接收信号强度矢量重构为新的n+1维矢量图2为某一条接收信号强度序列中的接收信号强度矢量进行时间戳标记的过程,以重构得到新的n+1维矢量;
步骤三:对中所包含的时间戳信息和接收信号强度信息进行加权,加权系数为wts和wrss,wts+wrss=1,即得到混合矢量
步骤四:对得到的混合矢量进行谱聚类,以得到Sj中每条混合矢量的聚类号;具体步骤为:
步骤四(一):对于Sj,构建邻接矩阵Mwa,如下式:
其中,wpq为Sj中与的相似度,且
其中,||·||表示二范数计算;
步骤四(二):在Mwa中设定相似度门限thw,当wpq<thw时,令wpq=0;当wpq≥thw时,wpq保持不变;在实验中thw=0.2;
步骤四(三):将Mwa中的每一列求和,得到构造M×M的对角矩阵Dm,如下所示:
步骤四(四):计算拉普拉斯算子L:
L=Dm-Mwa (4)
步骤四(五):计算L的特征值及特征向量,即根据
LV=vV (5)
(vE-L)V=0 (6)
det(vE-L)=0 (7)
其中,E为M×M的单位矩阵;v∈R为L的特征值;V表示L的特征向量,是长度为M的列向量;由公式7可得L的a(a≤M)个所有特征值;由公式6可得a个特征值所对应的a(a≤M)个特征向量;取L的c(c≤a)个最小特征值{vh}(h=1,…,c)及相应的c个特征向量Vh=(eh1,…,ehM)T,其中,vh为第h个特征值,Vh为第h个特征值所对应的特征向量,ehd(d=1,…,M)为第h个特征向量中第d个分量。根据公式14将c个特征向量构成一个M×c矩阵Mr;实验中c=7即为所划分的物理子区域个数;
步骤四(六):将Mr中的每一行定义为一个混合矢量的重构矢量,从而得到M个重构矢量Nd=(ned1,…,nedc),其中,nedf(f=1,…,c)为Mr中第d行的第f个元素,将M个重构矢量进行K均值聚类,即
步骤四(六)一:确定所需聚类的K均值聚类个数K(K<<M),实验中K=7;随机在M个重构矢量中选取K个样本Sab=(sab1,sab2,...,sabc)(b=1,…,K)作为初始类心,其中sabf(f=1,…,c)为Sab中第f个元素;
步骤四(六)二:逐个计算剩余的M-K个样本与所有初始类心的欧式距离,并将样本分配到与其欧氏距离(二范数)最小的初始类心所对应的聚类中;其中M-K个样本中的每个样本与每个初始类心的欧氏距离(二范数)计算为
步骤四(六)三:重新计算K个聚类的聚类中心,并将新的K个聚类的聚类中心作为初始类心,每个聚类的聚类中心计算如公式8所示:
其中,Cb表示第b个聚类的所有样本;X表示Cb中的每个样本;Nub表示Cb的样本个数;Zb表示第b个聚类的新类心。
步骤四(六)四:计算出所有聚类类心偏离度J,如公式9所示:
步骤四(六)五:重复步骤四(六)二步至步骤四(六)四,使得J达到最小值,此时的Sab(b=1,…,K)即为各聚类类心;将从小到大排序,将对应的类从1开始按升序排列到K,即可得到各类的类号;同类中的样本的类号相同;
至此,完成了步骤四的全部操作。
步骤五:利用中值滤波,修正每条混合矢量的聚类号及相应类心;即对于Sj相应的类号序列Cnj=(cnj1,…,cnjM),其中cnji∈{1,…,M}(i=1,…,M)为序列j中第i个接收信号强度矢量所对应的类号,用长度为3的滑动窗口(We-1 We We+1)(e=2,…,M-1)中的We对Cnj中的第二个元素到第M-1个元素进行遍历,其中当We遍历到cnje时,(We-1 We We+1)=(cnj(e-1),cnje,cnj(e+1)),取(We-1 We We+1)中各元素的中值替换cnje;当遍历结束后,当Cnj中的不同的类号总数小于M时,对Cnj中的类号进行重新排序,使不同类号从1开始连续递增,具体步骤为:
步骤五(一):从第1个元素开始遍历Cnj中的M个元素;
步骤五(二):记当前遍历元素为cnji,如果cnji为没有被遍历过的第g个新类号,则Ncg=cnji(1≤g≤M),其中,Ncg表示第g个没有被遍历过的新类号;
步骤五(三):重复步骤五(二),直到遍历完Cnj中的M个元素;
步骤五(四):依次取出Ncg,根据K均值聚类的结果,将Ncg所对应的类心记录在Nzg中,其中,Nzg表示Ncg所对应的类心;
步骤五(五):依次取出Ncg,将Ncg所对应的Cnj中的类号更新为g;
至此,完成了步骤五中修正每条混合矢量的聚类号及相应类心的操作。
步骤六:根据Sj中相邻聚类之间的转移关系,以连接图的形式得到Sj的类转移图;即在Sj对应的Cnj中,每一个不同的元素用一个圆圈表示,圆圈中用相应的元素值表示,依次遍历Cnj中的元素,当当前遍历元素与上一时刻遍历元素不同时用直线连接两元素对应的圆圈,当两元素所对应的圆圈存在多次连接时,只保留一条连接线;
步骤七:重复步骤二至步骤六,以得到所有接收信号强度序列的类转移图;
步骤八:利用图像边缘检测技术,确定所有类转移图中类间距离小于门限std的类并且合并相应的类从而对所有离散的类转移图进行拼接,得到待筛选的信号逻辑图,具体为:
步骤八(一):将所有采集到的RSS序列,按升序标号为1,2,...,z,并排列为序列
步骤八(二):分别计算Γ中不同接收信号强度矢量的欧式距离,进而得到距离矩阵Mdis。其中Mdis表示为:
其中,MΓ为序列Γ中的接收信号强度矢量数目;为Γ中第di1个接收信号强度矢量与第di2个接收信号强度矢量的欧氏距离,具体为:
其中,与分别表示Γ中第di1与第di2个接收信号强度矢量中来自第di3个AP的信号强度值;
步骤八(三):设定门限std,当时,令当时,令实验中std=50;
步骤八(四):对矩阵Mdis做中值滤波处理,具体步骤为:
步骤八(四)一:选择大小为3×3的窗口win作为中值滤波的基本单元,即
其中,win中共有9个元素,winwin1(win1=1,2,…,9)为win中的第win1个元素;
步骤八(四)二:将矩阵Mdis做补0处理得到Mdis_n,即
步骤八(四)三:将win中的win5对齐Mdis_n中的元素Mdis_n(dis_n1,dis_n2),其中dis_n1,dis_n2=2,3,...,MΓ+1,Mdis_n(dis_n1,dis_n2)为Mdis_n中第dis_n1行,第dis_n2列的元素;
步骤八(四)四:取出当前win对应于Mdis_n中的元素,将所有元素按从大到小排列并依次标号为1,2,…,9,且将矩阵Mdis中第(dis_n1)-1行第(dis_n2)-1列的元素更新为标号为5的元素;
步骤八(五):对矩阵Mdis做腐蚀处理,具体步骤为:
步骤八(五)一:选择大小为5×5的窗口ero,具体为:
步骤八(五)二:将ero中第三行第三列的元素依次对齐Mdis中的元素Mdis(dis_1,dis_2),其中dis_1=3,4,...,MΓ-2,dis_2=3,4,...,MΓ-2,Mdis(dis_1,dis_2)为Mdis中第dis_1行第dis_2列的元素;
步骤八(五)三:当对齐于ero的所有Mdis中元素全为1时,令Mdis(dis_1,dis_2)=1,而当对齐于ero的所有Mdis中元素不全为1时,令Mdis(dis_1,dis_2)=0;
步骤八(六):对矩阵Mdis做边缘检测处理,具体步骤为:
步骤八(六)一:确定3×3的边缘检测窗口edg,具体为:
步骤八(六)二:将edg中的元素edg5对齐矩阵Mdis中的元素Mdis(dis_a,dis_b),其中dis_a=2,3,...,MΓ-1,dis_b=2,3,...,MΓ-1,Mdis(dis_a,dis_b)为Mdis中第dis_a行第dis_b列的元素;
步骤八(六)三:选定横向及纵向边缘检测的Sobel卷积因子Gx及Gy,具体为:
步骤八(六)四:将edg中对应于矩阵Mdis中的元素分别与Gx及Gy相乘,得到关于当前元素Mdis(dis_a,dis_b)在横向及纵向上的灰度差分值Gx_dis_a及Gy_dis_b,具体为如公式17所示;
步骤八(六)五:得到关于当前元素Mdis(dis_a,dis_b)的灰度差分值Gxy,具体如公式18所示;
Gxy=|Gx_dis_a|+|Gy_dis_b| (18)
步骤八(六)六:设定灰度差分门限Gth,当Gxy≥Gth时,令Mdis(dis_a,dis_b)=1,当Gxy<Gth时,令Mdis(dis_a,dis_b)=0。实验中Gth=5;
步骤八(七):对矩阵Mdis的边缘信息进行提取,得到序列Γ中不同接收信号强度矢量距离小于std的序列块的长宽,以及其在Γ中的位置信息,具体步骤为:
步骤八(七)一:设定空矩阵co_1;
步骤八(七)二:遍历矩阵Mdis中的所有元素Mdis(dis_c,dis_d),其中dis_c=1,2,...,MΓ,dis_d=1,2,...,MΓ,Mdis(dis_c,dis_d)表示矩阵Mdis中第dis_c行第dis_d列的元素;当Mdis(dis_c,dis_d)=1时,将坐标(dis_c,dis_d)记录于co_1中;
步骤八(七)三:设定空矩阵co_2,co_3及co_4;
步骤八(七)四:令遍历号vinu=1;
步骤八(七)五:若vinu>vile则步骤结束,其中,vile为co_1中坐标个数,若vinu≤vile,则搜索co_1中第vinu个坐标co_1(vinu),并令co_1(vinu)为当前处理坐标co_1_(vinu);
步骤八(七)六:将co_1_(vinu)存入co_2中;
步骤八(七)七:若co_1_(vinu,1)≤co_1_(vinu,2),则令vinu=vinu+1并返回步骤八(七)五,其中,co_1_(vinu,1)和co_1_(vinu,2)分别为co_1_(vinu)的横坐标及纵坐标值;
步骤八(七)八:令处理号prnu=1;
步骤八(七)九:若co_1_(vinu,1)>co_1_(vinu,2),则确定co_1_(vinu)的八邻域坐标nef1,nef2,nef3,nef4,nee1,nee2,nee3,nee4,其中,nef1,nef2,nef3,nef4(八邻域中的四邻域坐标)及nee1,nee2,nee3,nee4(八邻域中的对角坐标)分别如公式19所示;
步骤八(七)十:确定坐标nef1,nef2,nef3,nef4在co_1中的个数,若nef1,nef2,nef3,nef4中有fo(fo=0,1,2,3,4)个坐标存在于co_1中,则令四邻域标记量nfo=fo;确定坐标nee1,nee2,nee3,nee4在co_1中的个数,若nee1,nee2,nee3,nee4中有ei(ei=0,1,2,3,4)个坐标存在于co_1中,则令对角标记量nei=ei;
步骤八(七)十一:若prnu=1,则令prnu=prnu+1并转步骤八(七)十二;若prnu≠1,则转步骤八(七)十三;
步骤八(七)十二:若nfo≠0,则将步骤八(七)十中的nfo个坐标中的任意一个设为当前坐标co_1_(vinu),将剩余的nfo-1个四邻域坐标及nei个对角坐标存入co_3并返回步骤八(七)五;若nfo=0,则将步骤八(七)十中的nei个坐标的任意一个设为当前坐标co_1_(vinu),将剩余的nei-1个对角坐标存入co_3并返回步骤八(七)五;
步骤八(七)十三:若nfo≠0,则将步骤八(七)十中的nfo个坐标与co_2中的坐标进行比较,若nfo个坐标不全包含于co_2,则将不包含于co_2中的任意一个坐标设为当前坐标co_1_(vinu),将剩余的不包含于co_2中的四邻域坐标及不包含于co_2中的对角坐标存入co_2并返回步骤八(七)五;若nfo个坐标全部包含于co_2,则将步骤八(七)十中的nei个坐标与co_2中的坐标进行比较,若坐标不全包含于co_2,则将不包含于co_2中的任意一个坐标设为当前坐标co_1_(vinu),将剩余的不包含于co_2中的对角坐标存入co_2并返回步骤八(七)五;若步骤八(七)十中的nei个坐标全部包含于co_2,则将co_2中的坐标与co_3中的坐标进行比较,如果co_2与co_3中有相同的坐标,则取co_2中所有横纵坐标中的最大值和最小值存入co_4中,其中,存储的横纵坐标的最大值和最小值的形式为(xmin,xmax,ymin,ymax),其中(xmin,xmax)及(ymin,ymax)分别为co_2中横纵坐标的最大值和最小值,令vinu=vinu+1,返回步骤八(七)五;如果co_2与co_3中没有相同的坐标,令vinu=vinu+1,返回步骤八(七)五;若nfo=0,则转步骤八(七)十四;
步骤八(七)十四:若nfo=0,则步骤八(七)十中的nei个坐标与co_2所包含的坐标进行比较,若坐标不全部包含于co_2,则将不包含于co_2中的任意一个坐标设为当前坐标co_1_(vinu),将剩余的不包含于co_2中的对角坐标存入co_2并返回步骤八(七)五;若坐标全部包含于co_2,则将co_2中的坐标与co_3中的坐标进行比较,如果co_2与co_3中有相同的坐标,取co_2中所有横纵坐标中的最大值和最小值存入co_4中,其中存储的横纵坐标的最大值和最小值的形式为(xmin,xmax,ymin,ymax),其中,(xmin,xmax)及(ymin,ymax)分别为co_2中横纵坐标的最大值和最小值,令vinu=vinu+1,返回步骤八(七)五;如果co_2与co_3中没有相同的坐标,令vinu=vinu+1,返回步骤八(七)五;
步骤八(八):根据步骤七中得到的co_4,将离散的信号逻辑图进行拼接,即遍历co_4中的所有co4个元素,对于co_4中的第co(co=1,2,...,co4)个元素(xco_min,xco_max,yco_min,yco_max),构造矩阵纵向元素区间[xco_min,xco_max],且由谱聚类所确定的分类结果,可得到此纵向元素区间中每一个接收信号强度矢量所对应的类号及相应的类心。类似地,构造矩阵横向元素区间[yco_min,yco_max],且由谱聚类所确定的分类结果,可得到此横向元素区间中每一个接收信号强度矢量所对应的类号及相应的类心。最后,将所有包含于纵向和横向元素区间中每一个接收信号强度矢量所对应的类心进行合并,即用所有接收信号强度矢量所对应的类心的均值作为合并后的新的类心,此类心代表一个合并的类,并保持该合并的类与其它类之间的连接关系,即将其它类与合并之前所确定的所有类的连接关系修正为与合并之后合并的类的连接关系;在遍历完co_4中的所有元素并进行相应的类合并操作后,就可得到经拼接后的逻辑图,然后,将此拼接后的逻辑图与剩余未拼接的逻辑图共同构成待筛选的逻辑图;
步骤九:将定位目标区域的每个叉路口作为区域边界进行子区域划分,即每个子区域的边界即为目标区域的路口,并对每个子区域进行序号随机标记,记区域标号为1,…,A,其中A表示所有的子区域的个数;实验中A=7;
步骤十:根据各子区域的邻接关系,将定位目标区域表示为各子区域连通的物理环境图,即将每个不同的子区域用一个圆圈表示,圈内的数字是相应的区域号,根据各子区域的物理连接特性,用直线将两个在目标区域中相连的子区域连接起来,将所有的连接关系确定之后,就可得到各子区域连通的物理环境图;
步骤十一:在定位目标区域内选择NC个标记位置点CP(Calibration Points),且保证标记位置点个数少于子区域个数;要求标记位置点在已知的特定的子区域中;CP={cpnc}(nc=1,…,NC;cpnc∈{1,…,A}),其中cpnc表示第nc个标记点所在的区域号;实验中NC=2,分别位于区域1与区域5中,即cpnc∈{1,5};
步骤十二:在各标记位置点处采集NR个来自不同AP的信号强度矢量其中 表示第nc个标记点上采集的第num个接收信号矢量;表示第nc个标记点上所采集的第num个接收信号矢量中来自第r个AP的信号强度,实验中NR=600,并将其均值矢量作为各标记位置点的代表矢量RV(Representative Vector);且RV={Mrssnc},其中,Mrssnc表示第nc个标记点的RV,其中Mrssnc如公式20所示;
步骤十三:在每个逻辑图中,计算与RV中每个元素Mrssnc与逻辑图中各类心的欧式距离其中,表示Mrssnc与第zj个逻辑图中的第f3个类心之间的欧氏距离(二范数),为第zj个逻辑图中的第f3个类心;选择与Mrssnc欧式距离最小的类心所对应的逻辑节点为该RV中元素Mrssnc所对应子区域cpnc存在的映射关系,剔除包含与2个或以上不同子区域存在映射关系的逻辑节点所对应的信号逻辑图,即在CP中,若存在2个或以上的元素映射到某一信号逻辑图中的同一个逻辑节点,则剔除该信号逻辑图;
步骤十四:根据步骤十三得到所有未剔除信号逻辑图,利用映射准则,得到相应的未剔除信号逻辑图与物理环境图的映射关系;每一个信号逻辑图到物理环境图的映射准则具体步骤为:
步骤十四(一):计算物理环境图中各子区域的邻接度AD(Adjacent Degree),各子区域的AD定义为该子区域和其邻接子区域所邻接的子区域个数总和,一个子区域与另一个子区域邻接表示这两个区域能够不经过其他区域到达彼区,此外,在所有子区域的AD中,得到最大ADAmag和最小ADAmig;
步骤十四(二):对于第x(x=1,…,y)个未剔除信号逻辑图Gx,其中,y为未剔除信号逻辑图个数,计算Gx中除标记点所在子区域{cpnc}对应的逻辑节点外的逻辑节点的AD,各逻辑节点的AD定义为该逻辑节点和其邻接逻辑节点所邻接的逻辑节点个数总和,一个逻辑节点与另一个逻辑节点邻接表示这两个逻辑节点能够不经过其他逻辑节点到达彼节点;
步骤十四(三):计算得到Gx中的最大ADAmal和最小ADAmil,对于除标记点所在子区域{cpnc}对应的逻辑节点外的逻辑节点的ADVADl,利用公式21,将其修正为VADg;
步骤十四(四):选择物理环境图中与VADg距离最小的AD所对应的子区域,作为该逻辑节点的初始映射子区域,即此逻辑节点与该子区域具有初始映射关系,此处距离定义为二者之差的绝对值;用表示第x个未剔除信号逻辑图中的第ma个逻辑节点所对应的映射子区域,Nodx表示第x个未剔除信号逻辑图的节点数;标记点所在子区域{cpnc}对应的逻辑节点的映射子区域即为相应的标记点子区域;
步骤十四(五):重复步骤十四(二)至步骤十四(四),可得所有未剔除信号逻辑图与物理环境图的初始映射关系。
对于上述得到的初始映射关系,将利用中心点进行校正,以得到最终的映射关系。具体操作如下:
步骤十四(六):在物理环境图中,寻找每个子区域到其它子区域的Floyd最短路径,对于一个给定的子区域,定义其所对应的中心区域为该子区域到其它子区域的Floyd最短路径的公共区域,当不存在此公共区域时,则此区域的中心区域不存在;其中,一个区域到另一个区域的Floyd最短路径确定方法如下:
步骤十四(六)一:定义区域Phpi(pi=1,…,A)(Phpi∈{1,…,A})到区域Phpj(pj=1,…,A)(Phpj∈{1,…,A})的距离为d(Phpi,Phpj),此处的距离表示Phpi到Phpj经过的区域中不包括源点的区域个数;定义区域Phpi到Phpj之间的路径中插入的所有区域的集合为r(Phpi,Phpj);定义w(Phpi,Phpj)为区域Phpi与Phpj之间的邻接关系,当区域Phpi与Phpj邻接时,w(Phpi,Phpj)=1,否则w(Phpi,Phpj)=0;
步骤十四(六)二:对于所有的Phpi与Phpj,设定初值d(Phpi,Phpj)=w(Phpi,Phpj);r(Phpi,Phpj)={Phpj};kin=1,其中kin为区域Phpi到Phpj之间的路径中插入的区域号;
步骤十四(六)三:若d(Phpi,kin)+d(kin,Phpj)<d(Phpi,Phpj),则d(Phpi,Phpj)=d(Phpi,kin)+d(kin,Phpj),并将kin记入r(Phpi,Phpj)中;否则保持所有参数不变;
步骤十四(六)四:若kin=A,则停止操作,否则kin=kin+1并返回步骤十四(六)三;
步骤十四(七):在每个未剔除信号逻辑图中,根据步骤十四(六)一至十四(六)三,寻找每个逻辑节点到其它逻辑节点的Floyd最短路径,对于一个给定的逻辑节点,定义其所对应的中心逻辑节点为该逻辑节点到其它逻辑节点的Floyd最短路径的公共逻辑节点;当不存在此公共逻辑节点时,则此逻辑节点的中心逻辑节点不存在;
步骤十四(八):对于逻辑图中映射到物理环境图中非公共区域的公共逻辑节点,将其所映射的子区域修正为与其VADg距离最小的AD所对应的公共区域;此处的距离为两者之差的绝对值;
步骤十四(九):重复步骤十四(七)和步骤十四(八),完成对所有初始映射关系的校正。
至此完成了未剔除信号逻辑图与物理环境图的映射关系;
步骤十五:选择对于所有标记点具有最高平均定位精度的信号逻辑图,作为最优信号逻辑图,对于一个信号逻辑图到物理环境图的映射而言,每个标记点的定位精度定义为在该标记点上采集的正确映射到所属子区域的信号强度矢量个数与信号强度矢量总数的比值如公式22所示;所有标记点的平均定位精度定义如公式22所示;
其中,表示第x个逻辑图的映射关系中,第nc个标记点上的定位精度,为第nc个标记点上采集的NR个接收信号强度矢量中能够映射到cpnc的接收信号强度矢量的个数;所有标记点的平均定位精度定义为:
步骤十六:根据步骤十五,选择最大的逻辑图作为最优信号逻辑图,根据步骤十四所确定的映射准则,可得该最优信号逻辑图到物理环境图的映射关系 表示第cet个逻辑图,也即最优逻辑图中各逻辑节点所对应的物理区域号,Nodcet表示第cet个逻辑图中的逻辑节点个数;
步骤十七:利用Graph Drawing正交算法对最优信号逻辑图及物理环境图进行绘制;具体为:
步骤十七(一):对于初始图(最优信号逻辑图或物理环境图)G=(V,E),首先随机选择两个点,一个点标记为s,另一个点标记为t,其中,V是G中所有点的集合,表示为V={vei}(i=1,2,...,nve);E是G中所有边的集合,表示为E={edj}(j=1,2,...,ned);vei表示G中的第i个的点,edj表示G中的第j条边,nve及ned分别表示G中点及边的数目;
步骤十七(二):找出包含s的所有边及所有边包含的点,以s为源点,以其余所有边包含的点为终点作有向边,并在所有边包含的点中,选择所有不为t的终点为新的起点;
步骤十七(三):对于每一个新的起点,找出包含该点的所有边及所有边包含的点,以该新的起点为源点,以其余所有边包含的点为终点作有向边,当某条有向边的两个点都为新的起点时,规定该边的绘制方向为从左往右或从上到下;
步骤十七(四):重复步骤十七(三),直到对所有新的起点都进行了从该点到相应终点的有向边的绘制,然后令步骤十七(三)中的所有不为t的终点为新的起点;
步骤十七(五):重复步骤十七(三)至十七(四),直到遍历完G中所有的点,从而构成了有向图G';
步骤十七(六):在图G'中将点s的值赋为0;并令其为起点;
步骤十七(七):令初始赋值inum=0;赋值间隔iit=1;
步骤十七(八):从起点开始,对图G'中的每条有向边的所关联的两个点赋值,所满足的条件为该边的终点的赋值大于起点的赋值。具体赋值规则为:
步骤十七(八)一:确定包含起点的所有边,再找到所有边所对应的终点;
步骤十七(八)二:将去掉点t的所有终点构成终点集,将所有终点赋值为inum+iit;
步骤十七(八)三:若终点集中存在连接两个终点的有向边,则将该边所对应的终点的值赋为起点值加1,并将步骤十七(八)二所确定的所有终点令为新的起点;
步骤十七(八)四:重复步骤十七(八)一至十七(八)三,直到遍历完除t点以外的所有点;
步骤十七(八)五:确定与点t相邻的点,并将所有与点t相邻的点的赋值中的最大值加1赋给点t,此时,图G'中每个点的赋值为该点在图中的纵坐标;
步骤十七(九):构建另一有向图DΠ,具体步骤为:
步骤十七(九)一:确定有向图DΠ的三部分点集;第一部分点集为图G'中的边所构成的所有最小区域,其中,最小区域定义为图G'中不存在其它边再把该区域进行分割的区域;第二部分点集为图G'以外空间所形成的两个区域,这两个区域满足:区域的并集为图G'以外的所有空间,而交集为空,此外,点s及t在两个区域的边界上,且定义两个区域中的左边区域(或上边区域)为s*,该区域的边界包含三段:第一段为连接s到t的最左边(最上边)的无向路径Π',其中,最左边(最上边)的无向路径为不存另一条连接s到t的无向路径Π使得Π在Π'的左边(上边)或者Π中的某条无向边在Π'的某条无向边的左边(上边);第二段为连接s及s'的无向边Π″,其中,s'为包含于图G'以外空间的某一点,Π″包含于图G'以外空间且s'与s的距离趋于无穷;第三段为连接t及t'的无向边Π″',其中,t'为包含于图G'以外空间的某一点,Π″'包含于图G'以外空间且t'与t的距离趋于无穷。同时,定义两个区域中的右边区域(或下边区域)为为t*(即图G'以外空间去除s*的剩余区域);第三部分点集为图G'中路径所构成的点集,路径的确定步骤为:
步骤十七(九)一(一):确定图G'中所有非s及t点的所有点,构成处理点集;
步骤十七(九)一(二):遍历处理点集中的所有点,对于每一个点,当入边数为奇数时,规定其入边为所有入边的最中间一条边,而当入边数为偶数时,规定其入边为所有入边的最中间两条边中偏左的一条边。其中,一个点的入边为以该点为终点的所有边。对于每一个点,当出边数为奇数时,规定其出边为所有出边的最中间一条边,而当出边数为偶数时,规定其出边为所有入边的最中间两条边中偏右的一条边。其中,一个点的出边为以该点为起点的所有边;
步骤十七(九)一(三):连接所有点所确定的边,构成所有确定的路径,并将图G'中剩余的边作为单独的路径,所有的路径就构成了第三部分点集;
步骤十七(十):确定有向图DΠ的有向边集,具体为:
步骤十七(十)一:对于图DΠ的点集中的每一个代表区域的点(即区域代表点),在图G'中遍历由步骤十七(九)一所确定的所有tnu条路径,对于图G'中第tnu1(tnu1=1,2,...,tnu)条路径若其所包含的某一条边包含于该区域代表点在图G'中对应区域的边界,则在有向图DΠ中,以该区域代表点为起点,以所代表的点为终点作有向边;
步骤十七(十)二:对于图DΠ的点集中的每一个代表路径的点(即路径代表点),在图G'中遍历由步骤十七(九)一所确定的所有anu个区域,对于图G'中第anu1(anu1=1,2,...,anu)个区域若其边界包含该路径代表点在图G'中对应路径的某一条边,则在有向图DΠ中,以该路径代表点为起点,以所代表的点为终点作有向边;
步骤十七(十一):在图DΠ中,将点s*的值赋为-0.5;并令其为起点;
步骤十七(十二):令初始赋值inum=-0.5;赋值间隔iit=0.5;
步骤十七(十三):执行步骤十七(八),其中,在步骤十七(八)中所有对t的操作修改为对t*的操作,且对图DΠ中所有点的赋值修改为图G'中所有边的横坐标;
步骤十七(十四):遍历图G'中的所有点,对于每一个点,其纵坐标为此点的赋值,横坐标有两个,分别为包含该点的路径在有向图DΠ中赋值的最大值与最小值,连接该两个点,即得关于该点的坐标表示线;
步骤十七(十五):遍历图G'中的所有边,对于每一条边,其横坐标为此边所在路径在有向图DΠ中的赋值,纵坐标有两个,分别为该边在有向图G'中的起点和终点的赋值,连接该两个点,即得关于该边的坐标表示线;
步骤十七(十六):在得到G'中所有点和所有边的坐标表示线后,在所有点的坐标表示线上画出一个实点,且实点的位置规则为:对于每一条坐标表示线,实点的位置在与此点的坐标表示线相交的边的坐标表示线的交点上。
步骤十八:根据用户新采集的信号强度矢量在最优信号逻辑图中计算出与新采集信号强度矢量具有最小欧式距离的类心所对应的逻辑节点maus;
步骤十九:根据最优信号逻辑图到物理环境图的映射关系估计用户所处的子区域为实验中的测试数据分别为7个不同区域中的各600个信号;
图3为将某一条接收信号强度序列进行谱聚类、中值滤波及类号重排处理后所得到的每条混合矢量的聚类号,以及由相邻聚类之间的转移关系得到的类转移图;
图4(a)为本发明采集的21条接收信号强度序列所形成的类转移图,图4(b)为利用Graph Drawing重构得到的21条接收信号强度序列所形成的类转移图,从图中可以看出所有的类转移图都以较好的可读性;
图5(a)为目标区域所对应的物理连接图,图5(b)为利用Graph Drawing重构得到的目标区域所对应的物理连接图,从图中也可看出物理连接图展现出了较好的可读性;
图6为目标区域的物理结构及21条测试路径;
图7(a)为原始图像,图像为一个二值图像,像素所对应的坐标值(xpi,ypi)分别为步骤八中序列Γ中的第xpi个接收信号强度矢量与第ypi个接收信号强度矢量,当像素值为1时,则表示第xpi个接收信号强度矢量与第ypi个接收信号强度矢量的欧式距离小于步骤八所设定的阈值std,当像素值为0时,则表示第xpi个接收信号强度矢量与第ypi个接收信号强度矢量的欧式距离大于或等于步骤八所设定的阈值std;图7(b)为图7(a)经中值滤波处理后得到的图像;图7(c)为图7(b)经腐蚀处理后得到的图像,从图中可以看出相关性较大的接收信号强度序列块在经过腐蚀处理后能有效凸显出来;
图8为利用Graph Drawing得到的最优逻辑图;图9给出了21条接收信号强度序列的RSS取值分布情况。
此外,实验中还对Sj中的混合矢量进行了层次聚类,以与本发明采用的谱聚类方法进行对比。层次聚类的基本步骤描述如下:
步骤一:对于Sj中的所有矢量,计算其两两间的欧式距离(二范数运算),以得到距离矩阵MD,如公式24所示:
其中,为第h1个混合RSS和第h2个混合RSS之间的欧氏距离(二范数);
步骤二:找出MD中除对角线以外的最小的元素将Sj第h3个混合RSS和第h4个混合RSS的均值作为一个新的混合矢量,且作为Sj的最后一个元素;同时在Sj中去掉第h3个混合RSS和第h4个混合RSS;
步骤三:重复步骤一和步骤二,直到Sj中只有一个元素位置,至此可得到一棵聚类树,如图12所示;
步骤四:选定聚类数目,若要求聚类数目为2,则根据图12中的L2划分,可得两类为{9,23,3,2,1,20,6,15}和{4,5,7,8,10,11,12,13,14,16,17,18,19,21,22,24,25,26,27,28,29,30}。本发明中把聚类数目设置为所划子区域数,即为7。
图10为在不同聚类方法条件下的子区域正确定位概率;
图11为在不同聚类方法条件下的子区域平均正确定位概率;
图12为层次聚类后所形成的聚类树。根据图12中的L2划分,可得两类为{9,23,3,2,1,20,6,15}和{4,5,7,8,10,11,12,13,14,16,17,18,19,21,22,24,25,26,27,28,29,30},其中,1到30的数字代表某一条接收信号强度序列中不同接收信号强度矢量的时间戳,可见,利用层次聚类也可对接收信号强度矢量进行聚类。表1给出了本发明所采用的谱聚类与传统层次聚类和K均值聚类所得到的物理映射区域对比结果。
表1
最后说明的是,以上优选实施例仅用以说明本发明的技术方案而非限制,尽管通过上述优选实施例已经对本发明进行了详细的描述,但本领域技术人员应当理解,可以在形式上和细节上对其作出各种各样的改变,而不偏离本发明权利要求书所限定的范围。
Claims (7)
1.基于图像边缘检测信号相关性的室内无线局域网信号地图绘制与映射方法,其特征在于:该方法包括以下步骤:
步骤一:用户在定位目标区域内随机采集若干条无线局域网接收的信号强度序列为的信号,其中,Mj为第j条序列的序列长度,z为序列条数;RSSjk=(RSSjk1,…,RSSjkn)为第j条序列中第k个接收信号强度矢量;n为接入点个数;RSSjkr为第j条序列第k个接收信号强度矢量中第r个接入点的信号强度值,j=1,2,...,z,k=1,...,Mj,r=1,...,n;
步骤二:根据时间戳顺序对Sj中接收信号强度矢量进行升序排列;将第j条序列中第k个接收信号强度矢量重构为新的n+1维矢量
步骤三:对中所包含的时间戳信息和接收信号强度信息进行加权,得到混合矢量其中加权系数为wts和wrss,wts+wrss=1;
步骤四:对混合矢量进行谱聚类,得到Sj中每条混合矢量的聚类号;
步骤五:通过中值滤波方法,修正每条混合矢量的聚类号及相应类心;
步骤六:根据Sj中相邻聚类之间的转移关系,以连接图的形式得到Sj的类转移图;
步骤七:重复步骤二至步骤六,得到所有接收信号强度序列的类转移图;
步骤八:通过图像边缘检测技术,确定所有类转移图中类间距离小于门限std的类,合并相应的类,对所有离散的类转移图进行拼接,得到待筛选的信号逻辑图;
步骤九:将定位目标区域内的每个物理叉路口作为区域边界进行子区域划分,并对每个子区域进行序号标记;记区域标号为1,…,A,其中A表示所有的子区域的个数;
步骤十:根据各子区域的物理邻接关系,将定位目标区域表示为一副由各子区域连通的物理环境图;
步骤十一:在定位目标区域内选择少量标记位置点CP,且保证标记位置点的个数少于子区域个数;
步骤十二:在各标记位置点处采集NR个来自不同接入点的信号强度矢量,并将其均值矢量作为各标记位置点的代表矢量RV;
步骤十三:在每个逻辑图中,计算RV中每个元素Mrssnc与逻辑图中各类心的欧氏距离其中,表示Mrssnc与第zj个逻辑图中的第f3个类心之间的欧氏距离,为第zj个逻辑图中的第f3个类心;选择与Mrssnc欧氏距离最小的类心所对应的逻辑节点为该RV中元素Mrssnc所对应子区域cpnc的映射关系,剔除包含与2个或以上不同子区域存在映射关系的逻辑节点所对应的信号逻辑图;
步骤十四:根据步骤十三得到所有未剔除信号逻辑图,利用映射准则,得到所有未剔除的信号逻辑图,以及相应的与物理环境图的映射关系;
步骤十五:选择对于所有标记点具有最高平均定位精度的信号逻辑图作为最优信号逻辑图;其中,每个标记点的定位精度定义为:在该标记点上采集的正确映射到所属子区域的信号强度矢量个数与信号强度矢量总数的比值,其中,表示第x个逻辑图的映射关系中第nc个标记点上的定位精度,为第nc个标记点上采集的NR个接收信号强度矢量中能够映射到cpnc的接收信号强度矢量的个数;所有标记点的平均定位精度 NC表示标记点个数;
步骤十六:根据步骤十五,选择最大的逻辑图作为最优信号逻辑图,根据步骤十四所确定的映射准则,得到最优信号逻辑图到物理环境图的映射关系ma∈{1,…,Nodcet});表示第cet个逻辑图,Nodcet表示第cet个逻辑图中的逻辑节点个数,y为逻辑图个数;
步骤十七:利用Graph Drawing正交算法,对最优信号逻辑图及物理环境图进行绘制;
步骤十八:根据终端新采集的信号强度矢量,在最优信号逻辑图中,计算得到与新采集信号强度矢量具有最小欧氏距离的逻辑节点;
步骤十九:根据最优信号逻辑图到物理环境图的映射关系,估计终端所在的子区域。
2.根据权利要求1所述的基于图像边缘检测信号相关性的室内无线局域网信号地图绘制与映射方法,其特征在于:所述步骤四包括以下步骤:
1)对于Sj,构建邻接矩阵Mwa,如下式:
其中,wpq为Sj中与的相似度,||·||表示二范数计算;p=1,…,M,q=1,…,M;
2)在Mwa中设定相似度门限thw,当wpq<thw时,令wpq=0;当wpq≥thw时,wpq保持不变;
3)将Mwa中的每一列求和,得到构造M×M的对角矩阵Dm,如下所示:
4)通过以下公式计算拉普拉斯算子L:
L=Dm-Mwa;
5)通过以下公式计算L的特征值及特征向量,
LV=vV;
(vE-L)V=0;
det(vE-L)=0;
其中,E为M×M的单位矩阵;v∈R为L的特征值;V表示L的特征向量,长度为M的列向量;
6)将Mr中的每一行定义为一个混合矢量的重构矢量,从而得到M个重构矢量Nd=(ned1,…,nedc),其中,nedf为Mr中第d行的第f个元素,将M个重构矢量进行K均值聚类,f=1,…,c,Mr为c个特征向量构成的维度。
3.根据权利要求2所述的基于图像边缘检测信号相关性的室内无线局域网信号地图绘制与映射方法,其特征在于:所述步骤6)将M个重构矢量进行K均值聚类包括以下步骤:
(1)确定所需聚类的K均值聚类个数K,K<<M;随机在M个重构矢量中选取K个样本Sab=(sab1,sab2,...,sabc)作为初始类心,其中sabf为Sab中第f个元素,b=1,…,K;
(2)逐个计算剩余的M-K个样本与所有初始类心的欧氏距离,并将样本分配到与其欧氏距离最小的初始类心所对应的聚类中;其中M-K个样本中的每个样本与每个初始类心的欧氏距离计算为
(3)通过以下公式重新计算K个聚类的聚类中心,并将新的K个聚类的聚类中心作为初始类心,
其中,Cb表示第b个聚类的所有样本;X表示Cb中的每个样本;Nub表示Cb的样本个数;Zb表示第b个聚类的新类心;
(4)通过以下公式计算出所有聚类类心偏离度J,
(5)重复步骤(2)步至步骤(4),直到J达到最小值,Sab为各聚类类心;将从小到大排序,得到各类的类号,同类中的样本的类号相同。
4.根据权利要求1所述的基于图像边缘检测信号相关性的室内无线局域网信号地图绘制与映射方法,其特征在于:所述步骤八包括以下步骤:
1)将采集到的RSS序列排列为
2)分别计算Γ中不同接收信号强度矢量的欧氏距离,进而得到距离矩阵Mdis,
其中,MΓ为序列Γ中的接收信号强度矢量数目; 为Γ中第di1个接收信号强度矢量与第di2个接收信号强度矢量的欧氏距离,di1,di2=1,2,...,MΓ;与 分别表示Γ中第di1与第di2个接收信号强度矢量中来自第di3个接入点的信号强度值;
3)设定门限std,当时,令当时,令
4)对矩阵Mdis进行中值滤波处理,包括以下步骤:
(1)选择大小为3×3的窗口win作为中值滤波的基本单元,如下所示,
其中,win中共有9个元素,winwin1为win中的第win1个元素,win1=1,2,…,9
(2)将矩阵Mdis做补0处理得到Mdis_n,如下所示,
(3)将win中的win5对齐Mdis_n Mdis_n中的元素Mdis_n(dis_n1,dis_n2),其中dis_n1=2,3,...,MΓ+1,dis_n2=2,3,...,MΓ+1,Mdis_n(dis_n1,dis_n2)为Mdis_n中第dis_n1行第dis_n2列的元素;
(4)取出当前win对应于Mdis_n中的元素,将所有元素按从大到小排列并依次标号为1,2,…,9,且将矩阵Mdis中第(dis_n1)-1行第(dis_n2)-1列的元素更新为标号为5的元素;
5)对矩阵Mdis进行腐蚀处理,包括以下步骤:
(1)选择大小为5×5的窗口ero,如下所示,
(2)将ero中第三行第三列的元素依次对齐Mdis中的元素Mdis(dis_1,dis_2),其中dis_1=3,4,...,MΓ-2,dis_2=3,4,...,MΓ-2,Mdis(dis_1,dis_2)为Mdis中第dis_1行第dis_2列的元素;
(3)当对齐于ero的所有Mdis中元素全为1时,令Mdis(dis_1,dis_2)=1,而当对齐于ero的所有Mdis中元素不全为1时,令Mdis(dis_1,dis_2)=0;
6)对矩阵Mdis进行边缘检测处理,包括以下步骤:
(1)选择大小为3×3的边缘检测窗口edg,如下所示:
(2)将edg中的元素edg5对齐矩阵Mdis中的元素Mdis(dis_a,dis_b),其中dis_a=2,3,...,MΓ-1,dis_b=2,3,...,MΓ-1,Mdis(dis_a,dis_b)为Mdis中第dis_a行第dis_b列的元素;
(3)选定横向及纵向边缘检测的Sobel卷积因子Gx及Gy,其中
(4)将对齐于edg的Mdis中元素分别与Gx及Gy相乘,得到关于当前元素Mdis(dis_a,dis_b)在横向及纵向上的灰度差分值Gx_dis_a及Gy_dis_b,其中,
(5)得到关于当前元素Mdis(dis_a,dis_b)的灰度差分值Gxy,Gxy=|Gx_dis_a|+|Gy_dis_b|;
(6)设定灰度差分门限Gth,当Gxy≥Gth时,令Mdis(dis_a,dis_b)=1,当Gxy<Gth时,令Mdis(dis_a,dis_b)=0;
7)对矩阵Mdis的边缘信息进行提取,得到序列Γ中不同接收信号强度矢量距离小于std的序列块的长宽,以及其在Γ中的位置信息;
8)根据步骤7)中得到的co_4,将离散的信号逻辑图进行拼接,co_4为横纵坐标中的最大值和最小值集合。
5.根据权利要求4所述的基于图像边缘检测信号相关性的室内无线局域网信号地图绘制与映射方法,其特征在于:所述步骤7)包括以下步骤:
S1)设定空矩阵co_1;
S2)遍历矩阵Mdis中的所有元素Mdis(dis_c,dis_d),其中dis_c=1,2,...,MΓ,dis_d=1,2,...,MΓ,Mdis(dis_c,dis_d)表示矩阵Mdis中第dis_c行第dis_d列的元素;当 Mdis(dis_c,dis_d)=1时,将坐标(dis_c,dis_d)记录于co_1中;
S3)设定空矩阵co_2,co_3及co_4;
S4)令遍历号vinu=1;
S5)若vinu>vile则步骤结束,其中,vile为co_1中坐标个数,若vinu≤vile,则搜索co_1中第vinu个坐标co_1(vinu),并令co_1(vinu)为当前处理坐标co_1_(vinu);
S6)将co_1_(vinu)存入co_2中;
S7)若co_1_(vinu,1)≤co_1_(vniu,2),则令vinu=vinu+1并返回步骤S5),其中,co_1_(vinu,1)和co_1_(vinu,2)分别为co_1_(vinu)的横坐标及纵坐标值;
S8)令处理号prnu=1;
S9)若co_1_(vinu,1)>co_1_(vinu,2),则确定co_1_(vinu)的八邻域坐标nef1,nef2,nef3,nef4,nee1,nee2,nee3,nee4,其中,八邻域中的四邻域坐标nef1,nef2,nef3,nef4及八邻域中的对角坐标nee1,nee2,nee3,nee4为:
S10)确定坐标nef1,nef2,nef3,nef4在co_1中的个数,若nef1,nef2,nef3,nef4中有fo个坐标存在于co_1中,则令四邻域标记量nfo=fo;确定坐标nee1,nee2,nee3,nee4在co_1中的个数,若nee1,nee2,nee3,nee4中有ei个坐标存在于co_1中,则令对角标记量nei=ei;fo=0,1,2,3,4,ei=0,1,2,3,4,
S11)若prnu=1,则令prnu=prnu+1并转步骤S12);若prnu≠1,则转步骤S13);
S12)若nfo≠0,则将步骤S10)中的nfo个坐标中的任意一个设为当前坐标co_1_(vinu),将剩余的nfo-1个四邻域坐标及nei个对角坐标存入co_3并返回步骤S5);若nfo=0,则将步骤S10) 中的nei个坐标的任意一个设为当前坐标co_1_(vinu),将剩余的nei-1个对角坐标存入co_3并返回步骤S5);
S13)若nfo≠0,则将步骤S10)中的nfo个坐标与co_2中的坐标进行比较,若nfo个坐标不全包含于co_2,则将不包含于co_2中的任意一个坐标设为当前坐标co_1_(vinu),将剩余的不包含于co_2中的四邻域坐标及不包含于co_2中的对角坐标存入co_2并返回步骤S5);若nfo个坐标全部包含于co_2,则将步骤S10)中的nei个坐标与co_2中的坐标进行比较,若坐标不全包含于co_2,则将不包含于co_2中的任意一个坐标设为当前坐标co_1_(vinu),将剩余的不包含于co_2中的对角坐标存入co_2并返回步骤S5);若步骤S10)中的nei个坐标全部包含于co_2,则将co_2中的坐标与co_3中的坐标进行比较,如果co_2与co_3中有相同的坐标,则取co_2中所有横纵坐标中的最大值和最小值存入co_4中,其中,存储的横纵坐标的最大值和最小值的形式为(xmin,xmax,ymin,ymax),其中(xmin,xmax)及(ymin,ymax)分别为co_2中横纵坐标的最大值和最小值,令vinu=vinu+1,返回步骤S5);如果co_2与co_3中没有相同的坐标,令vinu=vinu+1,返回步骤S5);若nfo=0,则转步骤S14);
S14)若nfo=0,则步骤S10)中的nei个坐标与co_2所包含的坐标进行比较,若坐标不全部包含于co_2,则将不包含于co_2中的任意一个坐标设为当前坐标co_1_(vinu),将剩余的不包含于co_2中的对角坐标存入co_2并返回步骤S5);若坐标全部包含于co_2,则将co_2中的坐标与co_3中的坐标进行比较,如果co_2与co_3中有相同的坐标,取co_2中所有横纵坐标中的最大值和最小值存入co_4中,其中存储的横纵坐标的最大值和最小值的形式为(xmin,xmax,ymin,ymax),其中,(xmin,xmax)及(ymin,ymax)分别为co_2中横纵坐标的最大值和最小值,令vinu=vinu+1,返回步骤S5);如果co_2与co_3中没有相同的坐标,令vinu=vinu+1,返回步骤S5)。
6.根据权利要求1所述的基于图像边缘检测信号相关性的室内无线局域网信号地图绘制与映射方法,其特征在于:所述步骤十四中每一个信号逻辑图到物理环境图的映射准则具体步骤为:
1)计算物理环境图中各子区域的邻接度AD,各子区域的AD定义为该子区域和其邻接子区域所邻接的子区域个数总和,一个子区域与另一个子区域邻接表示这两个区域能够不经过其他区域到达彼区,此外,在所有子区域的AD中,得到最大和最小
2)对于第x个未剔除信号逻辑图Gx,其中,y为未剔除信号逻辑图个数,计算Gx中除标记点 所在子区域{cpnc}对应的逻辑节点外的逻辑节点的AD,各逻辑节点的AD定义为该逻辑节点和其邻接逻辑节点所邻接的逻辑节点个数总和,一个逻辑节点与另一个逻辑节点邻接表示这两个逻辑节点能够不经过其他逻辑节点到达彼节点;x=1,2,….,y;
3)计算得到Gx中的最大和最小对于除标记点所在子区域{cpnc}对应的逻辑节点外的逻辑节点的通过以下公式将其修正为VADg,
4)选择物理环境图中与VADg距离最小的AD所对应的子区域,作为该逻辑节点的初始映射子区域,即此逻辑节点与该子区域具有初始映射关系,此处距离定义为二者之差的绝对值;用 表示第x个未剔除信号逻辑图中的第ma个逻辑节点所对应的映射子区域,Nodx表示第x个未剔除信号逻辑图的节点数;标记点所在子区域{cpnc}对应的逻辑节点的映射子区域即为相应的标记点子区域;
5)重复步骤2)至步骤4),可得所有未剔除信号逻辑图与物理环境图的初始映射关系,利用中心点对上述得到的初始映射关系进行校正,得到最终的映射关系。
7.根据权利要求1所述的基于图像边缘检测信号相关性的室内无线局域网信号地图绘制与映射方法,其特征在于:所述步骤十七包括以下步骤:
1)初始图G=(V,E),随机选择两个点,一个点标记为s,另一个点标记为t,其中,V是G中所有点的集合,表示为V={vei};E是G中所有边的集合,表示为E{edj};vei表示G中的第i个的点,edj表示G中的第j条边,nve及ned分别表示G中点及边的数目;i=1,2,...,nve,j=1,2,...,ned;
2)找出包含s的所有边及所有边包含的点,以s为源点,以其余所有边包含的点为终点作有向边,并在所有边包含的点中,选择所有不为t的终点为新的起点;
3)对于每一个新的起点,找出包含该点的所有边及所有边包含的点,以该新的起点为源点,以其余所有边包含的点为终点作有向边,当某条有向边的两个点都为新的起点时,规定该边的绘制方向为从左往右或从上到下;
4)重复步骤3),直到对所有新的起点都进行了从该点到相应终点的有向边的绘制,然后令步骤3)中的所有不为t的终点为新的起点;
5)重复步骤3)至步骤4),直到遍历完G中所有的点,从而构成了有向图G';
6)在图G'中将点s的值赋为0;并令其为起点;
7)令初始赋值inum=0;赋值间隔iit=1;
8)从起点开始,对图G'中的每条有向边的所关联的两个点赋值,所满足的条件为该边的终点的赋值大于起点的赋值,具体赋值规则为:
(1)确定包含起点的所有边,再找到所有边所对应的终点;
(2)将去掉点t的所有终点构成终点集,将所有终点赋值为inum+iit;
(3)若终点集中存在连接两个终点的有向边,则将该边所对应的终点的值赋为起点值加1,并将步骤2)所确定的所有终点令为新的起点;
(4)重复步骤(1)至步骤(3),直到遍历完除t点以外的所有点;
(5)确定与点t相邻的点,并将所有与点t相邻的点的赋值中的最大值加1赋给点t,此时,图G'中每个点的赋值为该点在图中的纵坐标;
9)构建另一有向图DΠ,具体步骤为:
确定有向图DΠ的三部分点集;第一部分点集为图G'中的边所构成的所有最小区域,其中,最小区域定义为图G'中不存在其它边再把该区域进行分割的区域;第二部分点集为图G'以外空间所形成的两个区域,这两个区域满足:区域的并集为图G'以外的所有空间,而交集为空,此外,点s及t在两个区域的边界上,且定义两个区域中的左边区域为s*,该区域的边界包含三段:第一段为连接s到t的最左边的无向路径Π',其中,最左边的无向路径为不存在另一条连接s到t的无向路径Π使得Π在Π'的左边或者Π中的某条无向边在Π'的某条无向边的左边;第二段为连接s及s'的无向边Π”,其中,s'为包含于图G'以外空间的某一点,Π”包含于图G'以外空间且s'与s的距离趋于无穷;第三段为连接t及t'的无向边Π’’’,其中,t'为包含于图G'以外空间的某一点,Π’’’包含于图G'以外空间且t'与t的距离趋于无穷,同时,定义两个区域中的右边区域为为t*;第三部分点集为图G'中路径所构成的点集,路径的确定步骤为:
(1)确定图G'中所有非s及t点的所有点,构成处理点集;
(2)遍历处理点集中的所有点,对于每一个点,当入边数为奇数时,规定其入边为所有入边的最中间一条边,而当入边数为偶数时,规定其入边为所有入边的最中间两条边中偏左的一条边;其中,一个点的入边为以该点为终点的所有边,对于每一个点,当出边数为奇数时,规定其出边为所有出边的最中间一条边,而当出边数为偶数时,规定其出边为所有入边的最中间两条边中偏右的一条边,其中,一个点的出边为以该点为起点的所有边;
(3)连接所有点所确定的边,构成所有确定的路径,并将图G'中剩余的边作为单独的路径, 所有的路径就构成了第三部分点集;
10)确定有向图DΠ的有向边集,具体为:
(1)对于图DΠ的点集中的每一个代表区域的点,在图G'中遍历由步骤9)所确定的所有tnu条路径,对于图G'中第tnu1条路径若其所包含的某一条边包含于该区域代表点在图G'中对应区域的边界,则在有向图DΠ中,以该区域代表点为起点,以所代表的点为终点作有向边;
(2)对于图DΠ的点集中的每一个代表路径的点,在图G'中遍历由步骤9)所确定的所有anu个区域,对于图G'中第anu1个区域若其边界包含该路径代表点在图G'中对应路径的某一条边,则在有向图DΠ中,以该路径代表点为起点,以所代表的点为终点作有向边;
11)在图DΠ中,将点s*的值赋为-0.5;并令其为起点;
12)令初始赋值inum=-0.5;赋值间隔iit=0.5;
13)执行步骤8),其中,在步骤8)中所有对t的操作修改为对t*的操作,且对图DΠ中所有点的赋值修改为图G'中所有边的横坐标;
14)遍历图G'中的所有点,对于每一个点,其纵坐标为此点的赋值,横坐标有两个,分别为包含该点的路径在有向图DΠ中赋值的最大值与最小值,连接最大值与最小值两个点,即得关于该点的坐标表示线;
15)遍历图G'中的所有边,对于每一条边,其横坐标为此边所在路径在有向图DΠ中的赋值,纵坐标有两个,分别为该边在有向图G'中的起点和终点的赋值,连接起点和终点的赋值两个点,即得关于该边的坐标表示线;
16)在得到G'中所有点和所有边的坐标表示线后,在所有点的坐标表示线上画出一个实点,且实点的位置规则为:对于每一条坐标表示线,实点的位置在与此点的坐标表示线相交的边的坐标表示线的交点上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410783008.XA CN104463929B (zh) | 2014-12-16 | 2014-12-16 | 基于图像边缘检测信号相关性的室内wlan信号地图绘制与映射方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410783008.XA CN104463929B (zh) | 2014-12-16 | 2014-12-16 | 基于图像边缘检测信号相关性的室内wlan信号地图绘制与映射方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104463929A CN104463929A (zh) | 2015-03-25 |
CN104463929B true CN104463929B (zh) | 2017-07-18 |
Family
ID=52909911
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410783008.XA Active CN104463929B (zh) | 2014-12-16 | 2014-12-16 | 基于图像边缘检测信号相关性的室内wlan信号地图绘制与映射方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104463929B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105188035B (zh) * | 2015-08-11 | 2018-06-15 | 重庆邮电大学 | 基于转移概率热点映射的室内wlan增广流形对齐定位方法 |
CN105205114B (zh) * | 2015-09-06 | 2018-12-04 | 重庆邮电大学 | 一种基于图像处理的Wi-Fi定位指纹数据库构建方法 |
CN105589064B (zh) * | 2016-01-08 | 2018-03-23 | 重庆邮电大学 | Wlan位置指纹数据库快速建立和动态更新系统及方法 |
CN105792138B (zh) * | 2016-02-26 | 2019-03-05 | 重庆邮电大学 | 基于线段密度的室内wlan用户定位与行为分析方法 |
CN105682206A (zh) * | 2016-03-17 | 2016-06-15 | 清华大学 | 一种大规模无线访问接入点的关停方法及装置 |
CN107180066B (zh) * | 2017-01-31 | 2017-12-15 | 张军民 | 基于三维空间编码的三维警用地理信息平台及系统架构 |
US10397739B2 (en) * | 2017-03-03 | 2019-08-27 | Here Global B.V. | Supporting the creation of a radio map |
US9820157B1 (en) * | 2017-03-06 | 2017-11-14 | Wipro Limited | Method and system for localizing spatially separated wireless transmitters |
CN109121143A (zh) * | 2017-06-23 | 2019-01-01 | 联芯科技有限公司 | 一种位置标记方法、终端与计算机可读存储介质 |
CN110075378B (zh) * | 2019-05-08 | 2021-11-19 | 南京市中西医结合医院 | 一种血液透析数据信息监测系统 |
CN111079787B (zh) * | 2019-11-15 | 2023-06-06 | 上海航天控制技术研究所 | 一种位置指纹识别方法 |
CN111447553B (zh) * | 2020-03-26 | 2021-10-15 | 云南电网有限责任公司电力科学研究院 | 一种基于wifi的增强视觉slam方法及装置 |
CN112762928B (zh) * | 2020-12-23 | 2022-07-15 | 重庆邮电大学 | 含有激光slam的odom与dm地标组合移动机器人及导航方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103888979A (zh) * | 2014-03-17 | 2014-06-25 | 南京邮电大学 | 一种基于无线局域网的室内定位方法 |
CN104185275A (zh) * | 2014-09-10 | 2014-12-03 | 北京航空航天大学 | 一种基于wlan的室内定位方法 |
-
2014
- 2014-12-16 CN CN201410783008.XA patent/CN104463929B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103888979A (zh) * | 2014-03-17 | 2014-06-25 | 南京邮电大学 | 一种基于无线局域网的室内定位方法 |
CN104185275A (zh) * | 2014-09-10 | 2014-12-03 | 北京航空航天大学 | 一种基于wlan的室内定位方法 |
Non-Patent Citations (2)
Title |
---|
《A Novel Algorithm for Multipath Fingerprinting in Indoor WLAN Environments》;Shih-Hau Fang et al;《IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS》;20080930;第7卷(第9期);3579-3588 * |
《Principal Component Localization in Indoor WLAN Environments》;Shih-Hau Fang et al;《IEEE TRANSACTIONS ON MOBILE COMPUTING 》;20120131;第11卷(第1期);100-110 * |
Also Published As
Publication number | Publication date |
---|---|
CN104463929A (zh) | 2015-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104463929B (zh) | 基于图像边缘检测信号相关性的室内wlan信号地图绘制与映射方法 | |
CN111299815B (zh) | 一种用于低灰度橡胶垫视觉检测与激光切割轨迹规划方法 | |
CN111797697B (zh) | 基于改进CenterNet的有角度高分遥感图像目标检测方法 | |
CN106646338B (zh) | 一种快速精确的室内定位方法 | |
WO2021143778A1 (zh) | 一种基于激光雷达的定位方法 | |
CN109325418A (zh) | 基于改进YOLOv3的道路交通环境下行人识别方法 | |
CN106413083B (zh) | 基于粗精双步相关图像特征提取的室内wlan定位方法 | |
CN109541634A (zh) | 一种路径规划方法、装置和移动设备 | |
CN105975929A (zh) | 一种基于聚合通道特征的快速行人检测方法 | |
CN107784657A (zh) | 一种基于颜色空间分类的无人机遥感图像分割方法 | |
CN109035276A (zh) | 一种图像边缘提取方法、装置及自动驾驶系统 | |
CN109410238A (zh) | 一种基于PointNet++网络的枸杞识别计数方法 | |
CN106710208A (zh) | 交通状态的获取方法及装置 | |
CN104282026A (zh) | 基于分水岭算法和最小生成树的分布均匀性评估方法 | |
CN104408463B (zh) | 一种高分辨率建设用地图斑识别方法 | |
CN104602341B (zh) | 基于随机用户信号逻辑图映射的室内wlan定位方法 | |
CN102496034A (zh) | 基于直线单词的高空间分辨率遥感图像词袋分类方法 | |
CN105956632A (zh) | 一种检测目标的方法和装置 | |
CN110991418A (zh) | 合成孔径雷达目标图像识别方法及系统 | |
CN106017476A (zh) | 一种生成室内定位导航图模型的方法 | |
AU2022258555B2 (en) | Generating mappings of physical spaces from point cloud data | |
CN104700420A (zh) | 基于Hough变换的椭圆检测方法、系统和虫卵识别方法 | |
CN108241819A (zh) | 路面标记的识别方法及装置 | |
CN107545726A (zh) | 一种公交车行驶速度确定方法及装置 | |
CN110084222A (zh) | 一种基于多目标角点池化神经网络的车辆检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |