1.基于图像边缘检测信号相关性的室内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，且保证标记位置点的个数少于子区域个数；

步骤十二：在各标记位置点处采集NR个来自不同AP的信号强度矢量，并将其均值矢量作为各标记位置点的代表矢量RV；

步骤十三：在每个逻辑图中，计算与RV中每个元素Mrssnc与逻辑图中各类心的欧式距离 其中， 表示Mrssnc与第zj个逻辑图中的第f3个类心之间的欧氏距离， 为第zj个逻辑图中的第f3个类心；选择与Mrssnc欧式距离最小的类心所对应的逻辑节点为该RV中元素Mrssnc所对应子区域cpnc存在的映射关系，剔除包含与2个或以上不同子区域存在映射关系的逻辑节点所对应的信号逻辑图；

步骤十四：根据步骤十三得到所有未剔除信号逻辑图，利用映射准则，得到所有未剔除的信号逻辑图，以及相应的与物理环境图的映射关系；

步骤十五：选择对于所有标记点具有最高平均定位精度的信号逻辑图作为最优信号逻辑图；其中，每个标记点的定位精度定义为：在该标记点上采集的正确映射到所属子区域的信号强度矢量个数与信号强度矢量总数的比值， 其中， 表示第x个逻辑图的映射关系中第nc个标记点上的定位精度， 为第nc个标记点上采集的NR个接收信号强度矢量中能够映射到cpnc的接收信号强度矢量的个数；所有标记点的平均定位精度步骤十六：根据步骤十五，选择 最大的逻辑图作为最优信号逻辑图，根据步骤十四所确定的映射准则，得到最优信号逻辑图到物理环境图的映射关系cet cetma∈{1,…,Nod })； 表示第cet个逻辑图，Nod 表示第cet个逻辑图中的逻辑节点个数；

步骤十七：利用Graph Drawing正交算法，对最优信号逻辑图及物理环境图进行绘制；

步骤十八：根据终端新采集的信号强度矢量，在最优信号逻辑图中，计算得到与新采集信号强度矢量具有最小欧式距离的逻辑节点；

步骤十九：根据最优信号逻辑图到物理环境图的映射关系，估计终端所在的子区域。

2.根据权利要求1所述的基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法，其特征在于：所述步骤四包括以下步骤：

1)对于Sj，构建邻接矩阵Mwa，如下式：

其中，wpq为Sj中 与 的相似度，

||·||表示二范数计算；

2)在Mwa中设定相似度门限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(f＝1,…,c)为Mr中第d行的第f个元素，将M个重构矢量进行K均值聚类。

3.根据权利要求2所述的基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法，其特征在于：所述步骤6)将M个重构矢量进行K均值聚类包括以下步骤：(1)确定所需聚类的K均值聚类个数K(K<

(2)逐个计算剩余的M-K个样本与所有初始类心的欧式距离，并将样本分配到与其欧氏距离最小的初始类心所对应的聚类中；其中M-K个样本中的每个样本与每个初始类心的欧氏距离计算为(3)通过以下公式重新计算K个聚类的聚类中心，并将新的K个聚类的聚类中心作为初始类心，其中，Cb表示第b个聚类的所有样本；X表示Cb中的每个样本；Nub表示Cb的样本个数；

Zb表示第b个聚类的新类心；

(4)通过以下公式计算出所有聚类类心偏离度J，

(5)重复步骤(2)步至步骤(4)，直到J达到最小值，Sab(b＝1,…,K)为各聚类类心；

将 从小到大排序，得到各类的类号，同类中的样本的类号相同。

4.根据权利要求1所述的基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法，其特征在于：所述步骤八包括以下步骤：

1)将采集到的RSS序列排列为

2)分别计算Γ中不同接收信号强度矢量的欧式距离，进而得到距离矩阵Mdis，其中，MΓ为序列Γ中的接收信号强度矢量数目；

为Γ中第di1个接收信号强度矢量与第di2个接收信号强度矢

量的欧氏距离； 与 分别表示Γ中第di1与第di2个接收信号强度矢量中来自第di3个AP的信号强度值；

3)设定门限std，当 时，令 当 时，令

4)对矩阵Mdis进行中值滤波处理，包括以下步骤：

(1)选择大小为3×3的窗口win作为中值滤波的基本单元，如下所示，其中，win中共有9个元素，winwin1(win1＝1,2,…,9)为win中的第win1个元素；

(2)将矩阵Mdis做补0处理得到Mdis_n，如下所示，

(3)将win中的win5对齐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，其中，Gx_dis_a＝Mdis((dis_a)+1,(dis_b)+1)+2Mdis((dis_a),(dis_b)+1)+Mdis((dis_a)-1,(dis_b)+1)-[Mdis((dis_a)+1,(dis_b)-1)+2Mdis((dis_a),(dis_b)-1)+Mdis((dis_a)-1,(dis_b)-1)]Gy_dis_b＝Mdis((dis_a)+1,(dis_b)+1)+2Mdis((dis_a)+1,(dis_b))；

+Mdis((dis_a)+1,(dis_b)-1)-[Mdis((dis_a)-1,(dis_b)+1)+2Mdis((dis_a)-1,(dis_b))+Mdis((dis_a)-1,(dis_b)-1)](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

7)对矩阵Mdis的边缘信息进行提取，得到序列Γ中不同接收信号强度矢量距离小于std的序列块的长宽，以及其在Γ中的位置信息；

8)根据步骤7)中得到的co_4，将离散的信号逻辑图进行拼接。

5.根据权利要求4所述的基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法，其特征在于：所述步骤7)包括以下步骤：(1)设定空矩阵co_1；

(2)遍历矩阵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中；

(3)设定空矩阵co_2,co_3及co_4；

(4)令遍历号vinu＝1；

(5)若vinu>vile则步骤结束，其中，vile为co_1中坐标个数，若vinu≤vile，则搜索co_1中第vinu个坐标co_1(vinu)，并令co_1(vinu)为当前处理坐标co_1_(vinu)；

(6)将co_1_(vinu)存入co_2中；

(7)若co_1_(vinu,1)≤co_1_(vinu,2)，则令vinu＝vinu+1并返回步骤(5)，其中，co_1_(vinu,1)和co_1_(vinu,2)分别为co_1_(vinu)的横坐标及纵坐标值；

(8)令处理号prnu＝1；

(9)若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为：(10)确定坐标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；

(11)若prnu＝1，则令prnu＝prnu+1并转步骤(12)；若prnu≠1，则转步骤(13)；

(12)若nfo≠0，则将步骤(10)中的nfo个坐标中的任意一个设为当前坐标co_1_(vinu)，将剩余的nfo-1个四邻域坐标及nei个对角坐标存入co_3并返回步骤(5)；若nfo＝0，则将步骤(10)中的nei个坐标的任意一个设为当前坐标co_1_(vinu)，将剩余的nei-1个对角坐标存入co_3并返回步骤(5)；

(13)若nfo≠0，则将步骤(10)中的nfo个坐标与co_2中的坐标进行比较，若nfo个坐标不全包含于co_2，则将不包含于co_2中的任意一个坐标设为当前坐标co_1_(vinu)，将剩余的不包含于co_2中的四邻域坐标及不包含于co_2中的对角坐标存入co_2并返回步骤(5)；

若nfo个坐标全部包含于co_2，则将步骤(10)中的nei个坐标与co_2中的坐标进行比较，若坐标不全包含于co_2，则将不包含于co_2中的任意一个坐标设为当前坐标co_1_(vinu)，将剩余的不包含于co_2中的对角坐标存入co_2并返回步骤(5)；若步骤(10)中的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，返回步骤(5)；如果co_2与co_3中没有相同的坐标，令vinu＝vinu+1，返回步骤(5)；若nfo＝0，则转步骤(14)；

(14)若nfo＝0，则步骤(10)中的nei个坐标与co_2所包含的坐标进行比较，若坐标不全部包含于co_2，则将不包含于co_2中的任意一个坐标设为当前坐标co_1_(vinu)，将剩余的不包含于co_2中的对角坐标存入co_2并返回步骤(5)；若坐标全部包含于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，返回步骤(5)；如果co_2与co_3中没有相同的坐标，令vinu＝vinu+1，返回步骤(5)。

6.根据权利要求1所述的基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法，其特征在于：所述步骤十四中每一个信号逻辑图到物理环境图的映射准则具体步骤为：

1)计算物理环境图中各子区域的邻接度AD，各子区域的AD定义为该子区域和其邻接子区域所邻接的子区域个数总和，一个子区域与另一个子区域邻接表示这两个区域能够不经过其他区域到达彼区，此外，在所有子区域的AD中，得到最大得到最大ADAmag和最小ADAmig；

2)对于第x(x＝1,…,y)个未剔除信号逻辑图Gx，其中，y为未剔除信号逻辑图个数，计算Gx中除标记点所在子区域{cpnc}对应的逻辑节点外的逻辑节点的AD，各逻辑节点的AD定义为该逻辑节点和其邻接逻辑节点所邻接的逻辑节点个数总和，一个逻辑节点与另一个逻辑节点邻接表示这两个逻辑节点能够不经过其他逻辑节点到达彼节点；

3)计算得到Gx中的最大ADAmal和最小ADAmil，对于除标记点所在子区域{cpnc}对应的逻辑节点外的逻辑节点的ADVADl，通过以下公式将其修正为VADg，

4)选择物理环境图中与VADg距离最小的AD所对应的子区域，作为该逻辑节点的初始映射子区域，即此逻辑节点与该子区域具有初始映射关系，此处距离定义为二者之差的绝对值；用 表示第x个未剔除信号逻辑图中的第ma个逻辑节点所对应x

的映射子区域，Nod表示第x个未剔除信号逻辑图的节点数；标记点所在子区域{cpnc}对应的逻辑节点的映射子区域即为相应的标记点子区域；

5)重复步骤2)至步骤4)，可得所有未剔除信号逻辑图与物理环境图的初始映射关系，利用中心点对上述得到的初始映射关系进行校正，得到最终的映射关系。

7.根据权利要求1所述的基于图像边缘检测信号相关性的室内WLAN信号地图绘制与映射方法，其特征在于：所述步骤十七包括以下步骤：

1)初始图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中点及边的数目；

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(tnu1＝1,2,...,tnu)条路径Πtnu1，若其所包含的某一条边包含于该区域代表点在图G'中对应区域的边界，则在有向图DΠ中，以该区域代表点为起点，以 所代表的点为终点作有向边；

(2)对于图DΠ的点集中的每一个代表路径的点，在图G'中遍历由步骤9)所确定的所有anu个区域，对于图G'中第anu1(anu1＝1,2,...,anu)个区域 若其边界包含该路径代表点在图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'中所有点和所有边的坐标表示线后，在所有点的坐标表示线上画出一个实点，且实点的位置规则为：对于每一条坐标表示线，实点的位置在与此点的坐标表示线相交的边的坐标表示线的交点上。