一种利用物体移动的标签二维定序方法

一种利用物体移动的标签二维定序方法
【专利摘要】本发明公开了一种利用物体移动的标签二维定序方法，目的在于，在不移动标签和阅读器的情况下对物品的二维次序进行探测，方便用户管理应用，所采用的技术方案为：利用移动物体通过影响区域重叠率的方法对物体移动过程中所影响的标签信号两两进行比较，确定处于同一列的标签，通过有向图邻接矩阵的方式对标签集合按列进行分组，再通过放置参考标签，确定分组后的标签集合在空间中的实际位置，得到水平方向的相对顺序；通过无物体移动的静止状态，获取标签阵列相位值，并结合标签集合按列分组结果，对同一列标签构建双曲线，筛选可能点，并对可能点进行线性规划，得到竖直方向的相对顺序，完成标签二维定序。
【专利说明】
-种利用物体移动的标签二维定序方法
技术领域
[0001] 本发明设及物联网无线感知领域，具体设及一种利用物体移动的标签二维定序方 法。
【背景技术】
[0002] 在现今的物流、仓储领域，射频识别技术已经逐步广泛应用起来。相比于之前的二 维码、条形码等技术，射频识别技术的优点在于：（1)除了能够利用标签携带一定的信息，还 可W方便的利用阅读器对信息进行更改。（2)不拘泥于将标签放在指定方位，只要将标签至 于阅读器的读取范围，便可进行读写操作。
[0003] 在常见的物流、仓储、超市等系统中，往往需要确定标签在系统中的位置信息，来 更好的对商品和货物进行管理和存储。现有的标签定位算法精度在几厘米到几十厘米不 等，但运些定位算法都各自具有各自的优缺点。比如，有的定位算法需要部署大量的参考标 签;有的定位算法需要部署大量阅读器天线;有的定位算法需要标签和阅读器天线处于相 对运动的状态。但是在实际的环境中，有时并不需要对标签的准确位置进行定位。比如，在 机场的传送带上，标签贴在行李上通过传送带传送到指定的位置;在超市的货架上，标签贴 在商品上，整齐的摆在货架上。在运些情况中，获取标签的顺序信息就足W得知标签的位置 信息。并且，相比于标签的定位，标签的定序算法降低了对标签位置信息的完整性的要求， 一些既得信息可根据实际环境直接获取或不影响系统的使用，因此，如果能够进行准确的 定序，那么从系统复杂度、硬件需求等方面都优于标签的定位方法。
[0004] 在物流、仓储和超市的环境内，经常需要对货品进行清点和管理，货品的位置顺序 至关重要。并且，在图书馆，超市运类人流量大的环境中，还会出现用户不断改变货品位置 的情况。在运样的情况下，利用人力逐个清点非常费时费力。现有的物品二维定序的算法需 要阅读器天线与标签有相对运动。如在机场的行李传送带上固定一个阅读器天线，对传送 带上运动的每个贴标签的行李进行排序。或在图书馆中，图书管理员将阅读器及阅读器天 线放在手推车上，走过书架前方，在运动的过程中确定标签的顺序。但是运些方法需要阅读 器天线和标签具有相对运动，运无疑加大了对系统部署环境和人力参与的要求。因此，利用 工作环境中任意移动的物体来确定环境内物品的二维顺序具有十分重要的意义和作用。

【发明内容】

[0005] 为了解决现有技术中的问题，本发明提出一种利用物体移动的标签二维定序方 法，在不移动标签和阅读器的情况下对物品的二维次序进行探测，方便用户管理应用。
[0006] 为了实现W上目的，本发明所采用的技术方案为:包括W下步骤：
[0007] 1)设置两个竖直放置且紧邻的阅读器天线，货架上的物体贴有标签，阅读器天线 接收标签信号，根据标签信号是否有扰动来判断阅读器天线与标签的通讯范围内有无物体 正在移动；
[000引2)若存在移动物体，通过影响区域重叠率的方法对物体移动过程中所影响的标签 信号两两进行比较，确定处于同一列的标签，通过有向图邻接矩阵的方式对标签集合按列 进行分组;通过在标签集合中放置参考标签，确定分组后的标签集合在空间中的实际位置， 从而得到标签集合在水平方向的相对顺序；
[0009] 3)若不存在移动物体，通过两个阅读器天线接收标签信号，获取标签阵列相位值；
[0010] 4)根据步骤2)中标签集合按列分组结果和步骤3)中的标签阵列相位值，对同一列 标签构建双曲线，筛选可能点，并对可能点进行线性规划，从而得到标签集合在竖直方向的 相对顺序，标签集合在水平方向和竖直方向的相对顺序构成标签的二维序列，完成标签二 维定序。
[0011] 所述步骤2)中两两标签信号的影响区域重叠率的计算公式如下：
[0012] (1)
[0013] i和j分别表示两个标签的编号，S(ti)和S(tj)表示两个标签影响区域归一化后和 时间轴所围成的区域，Rate(i，j)表示了两个标签共同影响区域的面积与总共影响区域的 面积之比；
[0014] 对于每一个标签i，选择与其重叠率最大的Μ个标签，作为与其同一列的候选标签， 运些标签构成的集合称为Φ (i)，在标签关系矩阵Anxn中，将所有符合条件j e Φ (i)的元素 (i，j)置为1，N为标签阵列的标签个数，公式如下：
[0015] |· (2)
[0016] aij为矩阵A的第i行第j列的元素，在标签关系矩阵A中，对于每一行，都有Μ个元素 置为1，其余的元素置为0，置为1的元素表示着编号为其行和列的两个标签处于同一列中。
[0017] 所述标签关系矩阵A与其转置矩阵AT进行矩阵相乘，当矩阵A中第i行第j列的元素 (1^'）为1，且矩阵4了中第1'行第^''列的元素。'^'）为1时，矩阵相乘的结果为1，相乘得到 的矩阵为权重矩阵R，即R = A · AT，权重矩阵R中第P行第q列的元素 rpq代表了编号为P和q的 两个标签所共同认定为同列标签的个数，表达式如下：
[001 引 rpq= I Φ (P) η Φ (q) I (3)
[0019] 元素 rpq越大，那么p和q所在同一列的可能性越高，通过判断rpq与阔值thre的大小 关系，最终判定P和q是否处在同一列。
[0020] 所述步骤4)中构建双曲线的步骤如下：
[0021] 标签T到两个阅读器天线Ai和A2的距离AiT和A2T满足下列限制条件：
[0022] A1T-A2TI <AiA2 (4)
[002引 AiA2的距离已知，AiA2 = d，d为两个阅读器天线之间的垂直距离，AiT和A2T的大小由 标签T在阅读器天线Ai和A2处收到的相位值计算，公式如下：
[0024]
[0025]
(5)
[00%]其中，ki和k2分别是标签Τ到达阅读器天线Ai和A2时所经过的整数周期的个数，θι和 曰2分别是标签Τ在阅读器天线Ai和Α2处收到的相位值，Θ*是标签自身的初始相位值，λ为阅读 器天线的电磁波波长；
[0027]将公式(5)代入公式(4)中得到：
[002引 (6)
[0029] 其中，k = lu-k2，Δ Θ =目广目2，目1，目2满足目1，目2 e (0，231)，则A Θ满足条件Δ Θ e (-化，2 3I)，将Δ目E(-化，23I)代入到公式(6)中，得至化的取值范围为k = 0，±l，±2;
[0030] W阅读器天线Ai，A2为双曲线的焦点，构建双曲线，双曲线上的点到双曲线两个焦 点的差值2a相同，根据k的取值，WA1T-A2T的值作为2a的取值画出双曲线，双曲线与距离阅 读器天线为D的直线的交点，就是标签T的位置，D为货架相对天线的距离。
[0031] 所述标签T位于货架的最底部时，标签T到两个阅读器天线的距离AiT和A2T满足下 列公式：
[0034] 当标签T位于货架的最高处时，标签T到两个天线的距离AiT和A2T满足下列公式：
[0035]
[0036]
[0037] D为货架相对阅读器天线的距离，h为货架的高度，d为两个阅读器天线之间的垂直 距离，ant为阅读器天线的长度，在标签位于货架最底部时A1T-A2T取到最大值，表示为max; 在标签位于货架最高处时A1T-A2T取到最小值，表示为min,符合下列限制条件：
[003引 min《AiT-A2T《max,
(9)。[0040] 所述限制条件公式(9)采用线性规划法用目标函数表示为：
[0039]
[0041]
[0042] 其中，和Μ j)分别是待定标签和错标签在货架上的高度，nil和叫分别是待定标 签和错标签在货架上的层数，错标签是k只有一个确定的取值的标签。
[0043] 与现有技术相比，本发明利用阅读器天线接收标签信号，通过标签信号是否有扰 动来判断是否有物体在移动，利用移动物体通过影响区域重叠率的方法对物体移动过程中 所影响的标签信号两两进行比较，确定处于同一列的标签，通过有向图邻接矩阵的方式对 标签集合按列进行分组，再通过在标签集合中放置参考标签，确定分组后的标签集合在空 间中的实际位置，从而得到标签集合在水平方向的相对顺序;通过无物体移动的静止状态， 获取标签阵列相位值，并结合标签集合按列分组结果，对同一列标签构建双曲线，筛选可能 点，并对可能点进行线性规划，从而得到标签集合在竖直方向的相对顺序，标签集合在水平 方向和竖直方向的相对顺序构成标签的二维序列，完成标签二维定序，本发明利用了移动 的物体对标签阵列的遮挡，与目前的定序方法相比，减少了劳动力投入，可靠性强，并具有 较高的准确率。
[0044] 进一步，当两个标签位于同一列时，物体的移动对其的影响时间是非常接近的；当 两个标签位于不同列时，物体的移动对其影响的时间差距会随着标签之间的水平距离增大 而增大，基于运个现象，将标签受物体影响的信号做归一化后，计算两两标签的影响区域的 重叠率，当两个标签位于同一列时，移动的物体对同一列上的标签影响区域非常接近，因此 重叠率较大;两个当标签位于不同列时，标签影响区域在时间轴上的位置随着标签距离的 增大而增大，因此重叠部分小，重叠率也较同列标签对小，利用计算影响区域重叠率Rate， 能够较准确的计算出两个标签影响区域的相似度。对于每一个标签，选择与其重叠率最大 的Μ个标签，作为与其同一列的候选标签，运些标签构成候选集合，在标签关系矩阵A中，对 于每一行，都有Μ个元素置为1，其余的元素置为0,置为1的元素表示着编号为其行和列的两 个标签，很可能处于同一列中。
[0045] 更进一步，由于在实际环境中容易出现误差，临近的列上的标签可能与该列上的 标签会出现比较大的重叠率，为了避免运些偶然误差对整体结果带来的影响，需要对标签 关系矩阵进行相应的处理。当两个标签位于同一列时，那么运两个标签通常互相都认定对 方为可能位于同一列的标签；当两个标签位于不同的列，出现偶然错误时，那么往往只有一 个标签认定另一个标签为同列标签。
[0046] 进一步，由于移动物体在货架前的移动是任意的、方向不明的，因此，对货架上标 签的影响顺序取决于物体在货架前移动的方向，为了获取物体在货架前移动的实际方向， 我们在货架上水平放置一排参考标签，根据参考标签的影响顺序，确定物体在水平方向移 动的方向，进而确定出每一列标签在实际中的水平顺序。
[0047] 进一步，货架相对阅读器天线的距离，货架的高度等真实环境的部署对k的取值有 一定的限制，根据Ξ角形勾股定理，当标签位于货架的最底部时，标签到两个天线的距离差 取到最大值，当标签位于货架的最高处时，标签到两个天线的距离差取到最小值，引入限制 条件，根据货架和天线的实际部署，将标签到两个阅读器天线的距离差缩小到一个更小的 范围内，能够大大的减少了k的可能值的个数，进一步减少了计算量。
[0048] 更进一步，将只有一个确定的k的取值的标签作为错标签，具有多个可能的k的标 签称之为待定标签，错标签的位置是唯一确定的，每一列标签都均匀分布于货架的每一排 上，且货架的横排之间的间距是已知的，利用线性规划化P)算法，将上述的两个限制条件用 目标函数的形式表示，目标函数是比较计算所得的标签位置与按货架高度估计所得的标签 位置之间的差距，只有当估计的待定标签位置全部正确时，运个差距才是最小的，从而确定 标签的最终位置。
【附图说明】
[0049] 图1是本发明的方法流程图；
[0050] 图2是重叠率计算结果图；
[0051 ]图3是竖直方向定序阅读器天线部署图；
[0化2] 图4是D = 1.2m~2.4m，h = 2m时的A1T-A2T的值域范围图；
[0化3] 图5是D = 1.8m，h = 2m时的k可能解个数图；
[0054] 图6是未加限制条件的不同取值的k构建的双曲线图；
[0055] 图7是增加限制条件的不同取值的k构建的双曲线图；
[0056] 图8通过双曲线+线性规划确定标签纵坐标的结果图。
【具体实施方式】
[0057] 下面结合具体的实施例和说明书附图对本发明作进一步的解释说明。
[005引参见图1，本发明包括W下步骤：
[0059] 1)首先通过标签信号是否有扰动来判断阅读器与标签的通讯范围内有无物体正 在移动；
[0060] 2)对于有物体移动的情况，首先将标签按列为单位分组，按照影响时间的接近程 度确定哪些标签处于同一列，对物体移动过程中所影响的标签信号两两进行比较，通过计 算影响区域重叠率的方法，对每个标签确定与其可能所处同一列的其他标签集合，通过有 向图邻接矩阵的乘法确定标签集合的划分，每个标签集合代表货架中的一列；
[0061] 3)通过放置参考标签，确定标签集合(一列标签)在实际空间中的实际位置，得到 水平方向的相对顺序；
[0062] 4)当无物体移动时，部署两个紧挨的竖直摆放的阅读器天线，获取标签阵列相位 值；
[0063] 5)结合无物体移动时的标签阵列相位值和有物体移动时确定的同列标签分类结 果，确定垂直方向的标签顺序，对同一列标签构建双曲线，筛选可能点，对可能点进行线性 规划，得到同一列标签竖直方向的相对顺序，标签集合在水平方向和竖直方向的相对顺序 构成标签的二维序列，完成标签二维定序。
[0064] 步骤1)中，利用标签信号有无强烈波动来判断工作区域内是否有物体移动。
[0065] 步骤2)中，对于有物体移动的情况，首先确定水平方向标签顺序。对标签阵列中所 有的标签信号进行两两比较。通常来说，当两个标签位于同一列时，物体的移动对其的影响 时间是非常接近的；当两个标签位于不同列时，物体的移动对其影响的时间差距会随着标 签之间的水平距离增大而增大。基于运个现象，将标签受物体影响的信号做归一化后，计算 两两标签的影响区域的重叠率。如公式(1)所示：
[0066]
[0067] 其中，i和j分别表示两个标签的编号，S(i)与S(j)表示两个标签影响区域归一化 后和时间轴所围成的区域。Rate(i，j)表示了两个标签共同影响区域的面积与总共影响区 域的面积之比。根据前面的设想，当标签i和j位于同一列时，移动的物体对同一列上的标签 影响区域非常接近，因此重叠率较大；当标签i和j位于不同列时，标签影响区域在时间轴上 的位置随着标签距离的增大而增大，因此重叠部分小，重叠率也较同列标签对小。利用计算 影响区域重叠率Rate,可W较准确的计算出两个标签影响区域的相似度。具体效果如图2a 和图化所示，其中，经处理的标签信号强度为图中的折线，折线与X轴所围成的部分即为影 响区域，两个标签影响区域交的部分面积的大小为s(ti)ns(tj)(图中阴影区域），两个标签 影响区域并的面积大小为S(ti)US(tj)，图2(a)为同列标签，图2(b)为不同列的标签，当两 个标签处于同一列时，其重叠部分占总的影响区域的面积更大，即Rate更大。
[0068] 对于每一个标签i，选择与其重叠率最大的Μ个标签，作为与其同一列的候选标签， 运些标签构成的集合称为Φα)。在标签关系矩阵Ανχν(Ν为标签阵列的标签个数）中，将所有 符合条件^εφα)的元素(i，j)置为1。如公式(2)所示：
[0069]
[0070] au为矩阵A的第i行第j列的元素。运样，在标签关系矩阵A中，对于每一行，都有Μ个 元素置为1，其余的元素置为0。置为1的元素表示着编号为其行和列的两个标签，很可能处 于同一列中。
[0071] 但是，由于在实际环境中容易出现误差，比如，临近的列上的标签可能与该列上的 标签会出现比较大的重叠率。为了避免运些偶然误差对整体结果带来的影响，需要对标签 关系矩阵A进行相应的处理。通过观察发现，当两个标签位于同一列时，那么运两个标签通 常互相都认定对方为可能位于同一列的标签，即je Φα)且ie ΦΟ);当两个标签位于不 同的列，出现偶然错误时，那么往往只有一个标签认定另一个标签为同列标签，如je Φ (i) 但z'狂抑/)。因此，为了利用运种性质去除一些偶然误差，我们将A与其转置AT进行矩阵相乘， 只有当矩阵A中第i行第j列的元素（i，j)为1，且矩阵AT中第i'行第j'列的元素（i'J'）也为 1时，其相乘的结果才为1。相乘得到的矩阵称为权重矩阵R，即R = A · AT。矩阵R中第P行第q 列的的元素 rpq代表了编号为P和q的两个标签所共同认定为同列标签的个数。例如，标签P认 定1，2,3标签为与其同列的标签，标签q认定2,3,4号标签为与其同列的标签，那么r(p，q) = 2。其表达式如公式(3)所示：
[007^ rpq= I Φ (P) η Φ (q) I (3)
[0073] 直观来说，当标签p和q处在同一列时，其认定为同列的标签重叠的个数会较大;反 之，如果标签P和q不处在同一列，那么其认定为同列的标签重叠个数会较少。因此，矩阵R可 W看做一个权重矩阵，矩阵中的元素 rpq越大，那么P和q所在同一列的可能性越高。通过判断 rpq与阔值thre的大小关系，可W最终判定P和q是否在同一列。rpq最大值为M，按照用户需 求，可W设置合适的阔值thre大小。
[0074] 步骤3)中，由于移动物体在货架前的移动是任意的、方向不明的，因此，对货架上 标签的影响顺序取决于物体在货架前移动的方向。为了获取物体在货架前移动的实际方 向，我们在货架上水平放置一排参考标签，根据参考标签的影响顺序，确定物体在水平方向 移动的方向，进而确定出每一列标签在实际中的水平顺序。例如，经过步骤2)后，我们已经 将整个标签阵列分为不同的集合，每个集合中的标签为同一列的标签，将运些集合随机进 行编号呆 当参考标签检测到移动物体是从货架最右侧移动到最左侧
0 时，假设影响的标签集合依次为3,2,4,6,8-，，那么，可^认为标签集合在实际环境内的部 署从右至左分别为L3，L2，L4，L6，Ls· · ·。
[0075] 步骤4)中，由于根据物体的移动很难确定标签的竖直顺序，因此，采用两个垂直距 离为d的竖直摆放的天线，见附图3,在无物体移动时收集标签阵列的相位值。运样，同一个 标签在两个天线处可采集到两个相位θι和02。
[0076] 步骤5)中，结合步骤3)中将标签分为不同的列的分类结果，利用步骤4)所采集到 的相位值信息，构建空间双曲线，其具体过程如下：
[0077] 由于标签Τ和两个天线Αι，Α2之间构成的Ξ角形Α1Α2Τ需要满足"Ξ角形两边之差小 于第Ξ边"的设定，因此，标签Τ到两个天线的距离AiT和Α2Τ需要满足下面的限制条件：
[007引 IA1T-A2T I <AiA2 (4)
[0079] A1A2的距离已知，即AiA2 = cLAiT，A2T的大小可由标签T在天线Ai，A2处收到的相位粗 略的估计，其表达式如公式(5)所示：
[0082] 其中，kl和k2分别是标签T到达天线Ai，A2时所经过的整数周期的个数，θ?和θ2分别 是标签Τ在天线Αι，Α2处收到的相位大小。0t是标签自身硬件所造成的初始相位值。λ为该应 用频率下电磁波波长。因此，标签Τ到两个天线的距离就可W通过收到的相位值和中间经过 的整周期与波长的乘积估计出来。公式(5)中，共有Ξ个未知数，即ki，k2和Θ*。将公式(5)代 入公式(4)，可得：
[0083]
[0084] 公式(6)中，将未知数Θ*经消减法去除，只留下一个未知数k = ki-k2。通过确定k的 大小，即可获得立角形边关系。由于标签T在天线Ai，A2处收到的相位大小θι，θ2是 观测值，因此满足θι，θ2Ε [0,231)，运样，ΔΘ需满足条件Δθε (-231,231)。将运个限制条件代 入到公式(6)中，可得未知数k的取值范围为k = 0，±l，±2。
[0085] 除此之外，由于标签所标识的物品放置在距离阅读器天线距离为D的货架上，见附 图3。货架相对天线的距离D，货架的高度h等真实环境的部署对未知数k的取值有一定的限 审IJ。根据Ξ角形勾股定理，ant为阅读器天线的长度，阅读器天线采用正方体的天线，当标签 T位于货架的最底部，那么，标签T到两个天线的距离AiT和A2T应符合下面的公式：
[0086]
[0087]
(7)
[0088] 而当标签Τ位于货架的最高处，那么，标签Τ到两个天线的距离AiT和A2T应符合下面 的公式：
[0091]变换标签的位置，可W获得标签到两个天线的距离差A1T-A2T。比如，当货架相对于 天线的距离为D = 1.2m~2.4m，货架的高度h = 2m时，标签到两个天线的距离差A1T-A2T如附 图4所示。我们发现，距离差A1T-A2T呈单调递减的趋势，在标签位于货架最底层时取到最大 值，将其表示为max;在标签位于货架最上层时取到最小值，将其表示为min。因此，标签到两 个天线的距离差A1T-A2T应当符合下面的限制条件：
[0094]因此，根据货架和天线的实际部署，可将两个天线的距离差A1T-A2T缩小到一个更 小的范围内。将公式(9)解出，可大大的减少未知数k可能值的个数。
[00M]为了说明上面的方法，我们举一个具体实例来说明。例如，在系统部署时，已知货 架到天线的距离为1.8米，货架的高度为2米。通过计算，距离差A1T-A2T的取值范围是[- 0.1742,0.0343]，也就是-2.1731《（沈π+ΔΘ)《0.4331，未知数k的可能取值范围可见附图5。 当标签T在两个天线处收到的相位差Δ Θ落在区间（-1.5731，-0.1731) U (0.4331，1.8331)时，只 存在一个可能的k值。而当相位差ΔΘ落在除运个区间之外的范围时，只存在两个可能的k 值。通过运样的方法，即可将公式(6)中5个可能的k值减少的不超过2个。通过双曲线定理， W天线Ai，A2为双曲线的焦点，构建双曲线。由于双曲线上的点到双曲线两个焦点的差值2a 是一样的。因此，每确定一个k的取值时，A1T-A2T就确定下来，将A1T-A2T的值作为2a的取值， 画出一条双曲线，双曲线与距离阅读器天线为D的直线的交点，就是标签的位置。如附图6所 示，未加限制条件时，k有5个可能解（k = 2的情况不满足构建双曲线条件，故未显示在图 中），因此有多条双曲线。而当增加实际部署时的限制条件时，如附图7a和7b所示，k只有1~ 2个可能解，只有1~2个双曲线满足限制条件。我们发现，标签在货架上的实际位置和双曲 线与货架位置的某个交点十分接近。因此，利用构建空间双曲线的方法来判断标签所处的 可能位置是十分有效的。
[0096]由于现在仍然存在不能确定的k的取值，因此标签的位置还没有唯一确定。我们观 察整个货架的标签分布情况发现，每一列的标签在货架上都均匀分布于货架的每一横排 上。且货架横排之间的距离是确定的，并且可W通过提前测量得知的。并且，在上述的举例 中发现，有些标签只有一个确定的k的取值，我们将运样的标签称之为"错标签"；而那些具 有多个可能的k的标签称之为"待定标签"。错标签的位置是确定的。因此，我们可W利用上 面两个限定条件来唯一确定标签位置的可能取值，即（1)错标签的位置时唯一确定的；（2) 每一列标签都均匀分布于货架的每一排上，且货架的横排之间的间距是已知的。利用线性 规划(LP)算法，可将上述的两个限制条件用目标函数的形式表达出来：
[0097]
[0098] 其中，/；(〇和Μ j)分别是待定标签和错标签在货架上的高度，nil和叫分别是待定标 签和错标签在货架上的层数。目标函数是比较计算所得的标签位置与按货架高度估计所得 的标签位置之间的差距。只有当估计的待定标签位置全部正确时，运个差距才是最小的。通 过线性规划所得的结果见附图8,发现利用构建双曲线，筛选可能点的方法具有很好的准确 度。
[0099 ]本发明的方法流程分为水平方向定位和竖直方向定位：
[0100] 水平方向定位模块：
[0101] 对于水平方向定位模块，可分为2个解决步骤，分别为：（1)将标签按列分组；（2)将 分组后的列与实际位置一一对应。在步骤(1)中，将移动的物体对每个标签影响的区域截取 出来，计算两两标签的影响区域重叠率rate,由于同一列的标签具有更接近的影响时间，故 影响区域的重叠面积更大，重叠率rate更大。参见附图2曰，两个同列标签的重叠率rate = 0.974,参见图化，而两个不同列标签的重叠率为rate = 0.688。通过计算重叠率，将整个标 签阵列按列分组。
[0102] 在步骤(2)中，通过部署参考标签，获取移动物体在货架前移动的水平方向，将按 列分组后的标签与实际的相对位置一一对应起来。
[0103] 竖直方向定位模块：
[0104] 对于竖直方向定位模块，也可分为2个解决步骤，分别为：（1)构建空间双曲线，筛 选可能解；（2)利用线性规划方法，确定同列标签的上下顺序。实验部署参见图3,将两个阅 读器天线垂直叠加放置与与货架距离为D的地面上。
[0105] 在步骤（1)中，利用实际环境内货架的高度hW及货架与阅读器天线的距离D对标 签实际可能位置的限制，确定k的取值。参见图4和5,经过上述的限制，标签到两个天线之间 的距离A1T-A2T限定在一定的范围内，且k的取值也减少到不超过2个。参见图7a和7b，增加限 制条件后，根据k值的大小，可构建1~2条双曲线，双曲线与距离为D的货架的交点，就是标 签在货架上可能的位置。
[0106] 在步骤(2)中，利用线性规划的方法，构建目标函数，筛选所有的可能解，当标签顺 序估计正确时，目标函数会取得最优解，实验效果参见附图8,其中a~e为货架的横排编号， 从图8可W看出，利用线性规划，构建目标函数的方法可W有效地排除错误解，找到标签位 置的可能解，具有较好的正确率。
[0107] 本发明出了一种可靠性强，准确度高的标签阵列二维定序方法，与之前的定序方 法相比，利用了移动的物体对标签阵列的遮挡，减少了劳动力投入，并具有较高的准确率， 本发明是一种可靠的且有效的标签阵列二维定序方法。
【主权项】
1. 一种利用物体移动的标签二维定序方法，其特征在于，包括以下步骤： 1) 设置两个竖直放置且紧邻的阅读器天线，货架上的物体贴有标签，阅读器天线接收 标签信号，根据标签信号是否有扰动来判断阅读器天线与标签的通讯范围内有无物体正在 移动； 2) 若存在移动物体，通过影响区域重叠率的方法对物体移动过程中所影响的标签信号 两两进行比较，确定处于同一列的标签，通过有向图邻接矩阵的方式对标签集合按列进行 分组;通过在标签集合中放置参考标签，确定分组后的标签集合在空间中的实际位置，从而 得到标签集合在水平方向的相对顺序； 3) 若不存在移动物体，通过两个阅读器天线接收标签信号，获取标签阵列相位值； 4) 根据步骤2)中标签集合按列分组结果和步骤3)中的标签阵列相位值，对同一列标签 构建双曲线，筛选可能点，并对可能点进行线性规划，从而得到标签集合在竖直方向的相对 顺序，标签集合在水平方向和竖直方向的相对顺序构成标签的二维序列，完成标签二维定 序。2. 根据权利要求1所述的一种利用物体移动的标签二维定序方法，其特征在于，所述步 骤2)中两两标签信号的影响区域重叠率的计算公式如下：i和j分别表示两个标签的编号，S(t〇和S(k)表示两个标签影响区域归一化后和时间 轴所围成的区域，Rate(iJ)表示了两个标签共同影响区域的面积与总共影响区域的面积 之比； 对于每一个标签i，选择与其重叠率最大的Μ个标签，作为与其同一列的候选标签，这些 标签构成的集合称为Φ (i)，在标签关系矩阵ΑΝΧΝ中，将所有符合条件j G φ (i)的元素(i，j) 置为1，N为标签阵列的标签个数，公式如下： ^ j I ./ G φ(?)师I 彐卩·)} y 1〇減〇 叭」(2) 为矩阵A的第i行第j列的元素，在标签关系矩阵A中，对于每一行，都有M个元素置为 1，其余的元素置为0，置为1的元素表示着编号为其行和列的两个标签处于同一列中。3. 根据权利要求2所述的一种利用物体移动的标签二维定序方法，其特征在于，所述标 签关系矩阵Α与其转置矩阵Ατ进行矩阵相乘，当矩阵Α中第i行第j列的元素（i，j)为1，且矩 阵八 1中第i'行第j'列的元素为1时，矩阵相乘的结果为1，相乘得到的矩阵为权重矩 阵R，即R = A · Ατ，权重矩阵!^中第p行第q列的元素rpq代表了编号为 p和q的两个标签所共同 认定为同列标签的个数，表达式如下： rPq= | Φ (ρ) (? φ (q) (3) 元素越大，那么p和q所在同一列的可能性越高，通过判断与阈值thre的大小关系， 最终判定ρ和q是否处在同一列。4. 根据权利要求1所述的一种利用物体移动的标签二维定序方法，其特征在于，所述步 骤4)中构建双曲线的步骤如下： 标签T到两个阅读器天线AdPA2的距离AiT和A2T满足下列限制条件： A1T-A2T | <AiA2 (4) A1A2的距离已知，AiA2 = d，d为两个阅读器天线之间的垂直距离，ΑιΤ和A2T的大小由标签 T在阅读器天线Αι和A2处收到的相位值计算，公式如下：其中，k#Pk2分别是标签T到达阅读器天线六:和知时所经过的整数周期的个数，Θ#ΡΘ2* 别是标签Τ在阅读器天线AdPA2处收到的相位值，0t是标签自身的初始相位值，λ为阅读器天 线的电磁波波长； 将公式(5)代入公式(4)中得到：其中，k = lu-k2, Δ 0 = 满足 0^02^(0,231),则 ΔΘ满足条件 Δθε (-231,231)， 将Δ θε (-231,231)代入到公式(6)中，得到k的取值范围为k=0, ±1，±2; 以阅读器天线^，如为双曲线的焦点，构建双曲线，双曲线上的点到双曲线两个焦点的 差值2a相同，根据k的取值，以值作为2a的取值画出双曲线，双曲线与距离阅读器 天线为D的直线的交点，就是标签T的位置，D为货架相对天线的距离。5. 根据权利要求4所述的一种利用物体移动的标签二维定序方法，其特征在于，所述标 签T位于货架的最底部时，标签T到两个阅读器天线的距离MT和A 2T满足下列公式：当标签T位于货架的最高处时，标签T到两个天线的距离MT和A2T满足下列公式：D为货架相对阅读器天线的距离，h为货架的高度，d为两个阅读器天线之间的垂直距 离，ant为阅读器天线的长度，在标签位于货架最底部时A1T-A2T取到最大值，表示为max;在 标签位于货架最高处时紐-知!'取到最小值，表示为min，符合下列限制条件：6. 根据权利要求5所述的一种利用物体移动的标签二维定序方法，其特征在于，所述限 制条件公式(9)采用线性规划法用目标函数表示为： min；,(i) <Σ ^ ~ mi - + Σ ~ mi' 1" ^ (10) 其中，和h( j)分别是待定标签和锚标签在货架上的高度，mi和mj分别是待定标签和 锚标签在货架上的层数，锚标签是k只有一个确定的取值的标签。
【文档编号】G06K7/00GK106066977SQ201610409005
【公开日】2016年11月2日
【申请日】2016年6月12日 公开号201610409005.9, CN 106066977 A, CN 106066977A, CN 201610409005, CN-A-106066977, CN106066977 A, CN106066977A, CN201610409005, CN201610409005.9
【发明人】韩劲松, 王鸽, 丁菡, 惠维, 赵季中
【申请人】西安交通大学