CN104809908A - Method of ZigBee network based vehicle positioning in indoor parking area environment - Google Patents

Method of ZigBee network based vehicle positioning in indoor parking area environment Download PDF

Info

Publication number
CN104809908A
CN104809908A CN201510231092.9A CN201510231092A CN104809908A CN 104809908 A CN104809908 A CN 104809908A CN 201510231092 A CN201510231092 A CN 201510231092A CN 104809908 A CN104809908 A CN 104809908A
Authority
CN
China
Prior art keywords
msub
mrow
positioning
positioning terminal
vehicle
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.)
Granted
Application number
CN201510231092.9A
Other languages
Chinese (zh)
Other versions
CN104809908B (en
Inventor
周鹏
魏英华
张鑫
赵健乐
刘希红
陈艳梅
戴永寿
李建风
吕正阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201510231092.9A priority Critical patent/CN104809908B/en
Publication of CN104809908A publication Critical patent/CN104809908A/en
Application granted granted Critical
Publication of CN104809908B publication Critical patent/CN104809908B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/14Traffic control systems for road vehicles indicating individual free spaces in parking areas
    • G08G1/141Traffic control systems for road vehicles indicating individual free spaces in parking areas with means giving the indication of available parking spaces
    • G08G1/144Traffic control systems for road vehicles indicating individual free spaces in parking areas with means giving the indication of available parking spaces on portable or mobile units, e.g. personal digital assistant [PDA]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W84/00Network topologies
    • H04W84/18Self-organising networks, e.g. ad-hoc networks or sensor networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种室内停车场环境下基于ZigBee网络的车辆定位方法,该方法通过在室内停车场内建立ZigBee无线传感网络,并对该网络进行分时分区域的划分来计算待定位车辆与定位终端结点间的距离,提高了距离的计算精度,进而提高了车辆定位的精度,有助于实现在室内停车场环境下对车辆的快速、准确引导。

The invention discloses a ZigBee network-based vehicle positioning method in an indoor parking lot environment. The method establishes a ZigBee wireless sensor network in the indoor parking lot, and divides the network into time-divided areas to calculate the distance between the vehicle to be positioned and the location. Positioning the distance between terminal nodes improves the calculation accuracy of the distance, thereby improving the accuracy of vehicle positioning, and helps to realize the rapid and accurate guidance of vehicles in the indoor parking lot environment.

Description

室内停车场环境下基于ZigBee网络的车辆定位方法Vehicle positioning method based on ZigBee network in indoor parking lot environment

技术领域technical field

本发明属于车辆定位技术领域,尤其涉及一种室内停车场环境下基于ZigBee网络的车辆定位方法。The invention belongs to the technical field of vehicle positioning, in particular to a ZigBee network-based vehicle positioning method in an indoor parking lot environment.

背景技术Background technique

智能交通在城市道路交通系统中占有非常重要的地位,针对城市交通进行合理的规划、控制、监测和诱导,对提高城市交通路网的效率以及提高城市居民的生活水平,具有重要的意义。车辆定位系统是智能交通的重要支撑技术之一,车辆作为移动结点的正确定位能够直观有效地对城市交通的运行状态进行观测和控制。Intelligent transportation occupies a very important position in the urban road traffic system. Reasonable planning, control, monitoring and guidance for urban traffic is of great significance to improve the efficiency of urban traffic network and improve the living standards of urban residents. The vehicle positioning system is one of the important supporting technologies of intelligent transportation. The correct positioning of the vehicle as a mobile node can directly and effectively observe and control the operation status of urban traffic.

为实现在大型室内停车场环境下将车辆自动引导至事先规划好的目的车位,需首先解决车辆在室内停车场行驶过程中的定位问题。由于在室内停车场环境下,导航仪所接收到的卫星信号通常比较微弱,降低了定位停车的效率和准确性。因此,越来越多的室内停车场选用多点间的无线通信技术来实现室内停车。目前,常用的无线通信技术包括:蓝牙、WIFI、RFID、ZigBee、超宽带通信、GPRS、3G、4G等,综合考虑成本、功耗、作用距离、精度、网络容量等多种因素,一般宜采用基于ZigBee无线传感网络技术的车辆定位方案。基于ZigBee无线传感网络的车辆定位技术是在车辆的智能手机上连接商业化的ZigBee模块,在车辆的移动过程中,该模块通过与车辆周围的ZigBee定位终端结点进行通信实现对车辆位置的定位,然后将位置信息传递给智能手机软件上的导航模块,导航模块利用车辆当前位置信息以及停车场的电子地图,最终实现对目的车位的路径引导。In order to automatically guide the vehicle to the pre-planned destination parking space in a large indoor parking lot environment, it is necessary to first solve the positioning problem of the vehicle in the indoor parking lot. Because in the indoor parking lot environment, the satellite signal received by the navigator is usually relatively weak, which reduces the efficiency and accuracy of positioning parking. Therefore, more and more indoor parking lots use wireless communication technology between multiple points to realize indoor parking. At present, commonly used wireless communication technologies include: Bluetooth, WIFI, RFID, ZigBee, ultra-wideband communication, GPRS, 3G, 4G, etc. Considering various factors such as cost, power consumption, operating distance, accuracy, and network capacity, it is generally appropriate to use Vehicle positioning scheme based on ZigBee wireless sensor network technology. The vehicle positioning technology based on ZigBee wireless sensor network is to connect the commercialized ZigBee module to the smart phone of the vehicle. During the moving process of the vehicle, the module communicates with the ZigBee positioning terminal nodes around the vehicle to realize the location of the vehicle. Positioning, and then transmit the location information to the navigation module on the smartphone software. The navigation module uses the current location information of the vehicle and the electronic map of the parking lot to finally realize the path guidance to the target parking space.

根据ZigBee技术中基于接收信号强度指示RSSI的定位技术理论及无线信号传输领域中的Shadowing模型,接收信号的强度与收发结点间距离的关系为:According to the positioning technology theory based on RSSI in ZigBee technology and the Shadowing model in the field of wireless signal transmission, the relationship between the strength of the received signal and the distance between the transmitting and receiving nodes is:

PP LL (( dd )) == PP LL (( dd 00 )) -- 1010 nlgnlg (( dd dd 00 )) ++ Xx dd

其中,d为发送端与接收端间的距离,单位为m;d0为参考距离,单位为m;PL(d)为接收端的接收信号功率,单位为dBm;PL(d0)为距离d0处的接收信号功率,单位为dBm;Xd是服从高斯分布的随机变量,均值为零,单位为dBm;n为路径衰减指数,它表示RSSI值随着距离增加而衰减的快慢。当进行大量、独立测量时,可令Xd=0(dBm),此时有:其中,为接收信号平均功率,单位为dBm;a为距发送端1m处的接收信号功率的相反数,单位为dBm。例如:如果在距发送端1m处的接收信号平均功率为-30dBm,则参数a的值为30,对上式进行变形后,可得:由此可见,若已知参数a和n的值,经多次测量取平均得到后,即可计算出收发结点间的距离,若具有3个以上的距离测量值,可根据三角形质心定位算法或极大似然估计法确定待定位结点的坐标。a和n的精度是决定测距乃至定位精度的关键因素,这两个参数的取值不仅与无线信号传输时的本地环境有关,还会随时间发生变化,因此给它们的测量带来了很大的困难。比较传统的方法是经验值法,即定位前在定位场合中进行大量的实验测量,通过线性回归的方法确定a和n的取值。具体实施过程为:在待定位区域,首先安装一个固定的发射装置,然后在某位置处摆放接收装置,精确测量发射装置与接收装置间的距离,采集一组RSSI值并求平均;接下来,按事先规划好的移动路线,将接收装置挪到另一位置,重复上述过程;这样的测量过程一直重复下去,直到整个区域被接收装置的移动路线所“覆盖”完毕。对应的系数计算公式为:Among them, d is the distance between the sending end and the receiving end, the unit is m; d 0 is the reference distance, the unit is m; PL (d) is the received signal power of the receiving end, the unit is dBm; P L (d 0 ) is The received signal power at distance d 0 , in dBm; X d is a random variable that obeys Gaussian distribution, with a mean value of zero, in dBm; n is the path attenuation index, which indicates the speed attenuation of the RSSI value as the distance increases. When performing a large number of independent measurements, X d =0(dBm) can be set, at this time: in, is the average power of the received signal, in dBm; a is the opposite number of the received signal power 1m away from the sending end, in dBm. For example: if the average power of the received signal at 1m away from the sending end is -30dBm, then the value of parameter a is 30, after modifying the above formula, we can get: It can be seen that, if the values of the parameters a and n are known, the average value obtained after multiple measurements is After that, the distance between the sending and receiving nodes can be calculated. If there are more than 3 distance measurements, the coordinates of the nodes to be located can be determined according to the triangle centroid positioning algorithm or the maximum likelihood estimation method. The accuracy of a and n is the key factor that determines the distance measurement and even the positioning accuracy. The values of these two parameters are not only related to the local environment when the wireless signal is transmitted, but also change with time, so they bring a lot to their measurement. big difficulty. The more traditional method is the empirical value method, that is, a large number of experimental measurements are carried out in the positioning field before positioning, and the values of a and n are determined by linear regression. The specific implementation process is: in the area to be positioned, first install a fixed transmitting device, then place the receiving device at a certain position, accurately measure the distance between the transmitting device and the receiving device, collect a set of RSSI values and calculate the average; , move the receiving device to another location according to the pre-planned moving route, and repeat the above process; this measurement process is repeated until the entire area is "covered" by the moving route of the receiving device. The corresponding coefficient calculation formula is:

nno ^^ == -- ΣΣ ii == 11 Mm (( lglg dd ii -- lgdlgd ‾‾ )) (( RSSIRSSI ii -- RSSIRSSI ‾‾ )) 1010 ΣΣ ii == 11 Mm (( lglg dd ii -- lgdlgd ‾‾ )) 22 aa ^^ == -- RSSIRSSI ‾‾ -- 1010 nno ^^ ·&Center Dot; lgdlgd ‾‾

式中,M为测量的点数,di为发射机与第i个接收点间的精确距离,为M个距离值取对数后的平均值,RSSIi为在第i个接收点进行的多次测量中获得的RSSI值的平均值,为M个RSSIi值的平均值。从上述测量过程可看出,传统方法的工作量非常大,并且也无法反应a和n随时间和空间的变化。In the formula, M is the number of measured points, d i is the precise distance between the transmitter and the i-th receiving point, is the average value after logarithm of M distance values, RSSI i is the average value of RSSI values obtained in multiple measurements at the i receiving point, is the average value of M RSSI i values. It can be seen from the above measurement process that the workload of the traditional method is very large, and it cannot reflect the changes of a and n with time and space.

由此可见,现有技术有待于进一步的改进和提高。This shows that the prior art needs to be further improved and improved.

发明内容Contents of the invention

本发明为避免上述现有技术存在的不足之处,提供了一种室内停车场环境下基于ZigBee网络的车辆定位方法。In order to avoid the disadvantages of the prior art, the present invention provides a ZigBee network-based vehicle positioning method in an indoor parking lot environment.

本发明所采用的技术方案为:The technical scheme adopted in the present invention is:

室内停车场环境下基于ZigBee网络的车辆定位方法,该方法包括以下步骤:A vehicle positioning method based on ZigBee network under indoor parking lot environment, the method comprises the following steps:

步骤1:在室内停车场内设置ZigBee无线传感网络系统,ZigBee无线传感网络系统,包括一个ZigBee网络,在ZigBee网络中设置一个中心结点、若干个路由器结点以及数量更多的终端结点,终端结点包括用于对待定位车辆进行定位的定位终端结点;待定位车辆的驾乘者持有智能手机,智能手机上配有ZigBee模块,待定位车辆通过智能手机发出定位请求;Step 1: Set up the ZigBee wireless sensor network system in the indoor parking lot. The ZigBee wireless sensor network system includes a ZigBee network, and a central node, several router nodes and more terminal nodes are set in the ZigBee network. The terminal node includes the positioning terminal node for locating the vehicle to be positioned; the driver and passenger of the vehicle to be positioned hold a smart phone, and the smart phone is equipped with a ZigBee module, and the vehicle to be positioned sends a positioning request through the smart phone;

步骤2:按照室内停车场中车位的分布,对上述定位终端结点进行区域划分,使之划分成若干个互不重叠的区域,并对各区域进行编号,每个区域呈三角形或四边形,定位终端结点安装在每个三角形或四边形的顶点上,各区域中,每个定位终端结点和与其临近的两个定位终端结点进行直接通信;Step 2: According to the distribution of parking spaces in the indoor parking lot, the above positioning terminal nodes are divided into several non-overlapping areas, and each area is numbered. Each area is triangular or quadrilateral. Positioning Terminal nodes are installed on the vertices of each triangle or quadrilateral, and in each area, each positioning terminal node communicates directly with two adjacent positioning terminal nodes;

步骤3:对上述每个区域的各定位终端结点的坐标进行精确测量后计算相邻两顶点间的距离,并根据下式推算得出区域中某一个定位终端结点周围环境的两个待定常数ai、niStep 3: After accurately measuring the coordinates of each positioning terminal node in each of the above areas, calculate the distance between two adjacent vertices, and calculate the two undetermined values of the surrounding environment of a certain positioning terminal node in the area according to the following formula: Constants a i , n i :

RSSIRSSI ‾‾ ii 11 [[ dBmdBm ]] == -- (( aa ii ++ 1010 nno ii ·&Center Dot; lglg dd ii 11 )) RSSIRSSI ‾‾ ii 22 [[ dBmdBm ]] == -- (( aa ii ++ 1010 nno ii ·· lglg dd ii 22 ))

其中,di1和di2为区域中第i个定位终端结点与周围第1个和第2个相邻定位终端结点间的距离,为区域中第i个定位终端结点收到周围第1个和第2个相邻定位终端结点所发送信号的平均功率;Among them, d i1 and d i2 are the distances between the i-th positioning terminal node in the area and the surrounding first and second adjacent positioning terminal nodes, and is the average power of the signal sent by the i-th positioning terminal node in the area received by the surrounding first and second adjacent positioning terminal nodes;

步骤4:根据计算出的距离值di1和di2,以及测量值计算出ai和ni,并将区域内所有定位终端结点的ai和ni取平均作为该区域内a和n的取值,得出:Step 4: Based on the calculated distance values d i1 and d i2 , and the measured values and Calculate a i and n i , and take the average of a i and n i of all positioning terminal nodes in the area as the values of a and n in the area, and obtain:

aa ‾‾ == 11 LL ΣΣ ii == 11 LL aa ii nno ‾‾ == 11 LL ΣΣ ii == 11 LL nno ii

其中,为该区域内两个待定常数的平均值,L为该区域内定位终端结点的个数;in, and is the average value of two undetermined constants in this area, and L is the number of positioning terminal nodes in this area;

当在停车场的入口或出口处检测到有车辆进入或离开时,由中心结点通过路由器结点向各定位终端结点发出命令,各定位终端结点与邻近的定位终端结点通信,对ai和ni的值进行更新,并重新计算得出的待定常数待定位车辆提出定位请求后,所有与其处于能直接通信范围内的定位终端结点都会做出响应,向待定位车辆发送定位终端结点所处区域的a和n的平均值,因某些定位终端结点通常被划分至多个区域,所以这类定位终端结点会发送多个a和n的平均值以及所对应区域的编号,待定位车辆的智能手机接收到附近的定位终端结点反馈的信息后,选择编号出现次数最多的区域进行下一步骤;When a vehicle is detected to enter or leave at the entrance or exit of the parking lot, the central node sends commands to each positioning terminal node through the router node, and each positioning terminal node communicates with the adjacent positioning terminal node. The values of a i and ni are updated, and the undetermined constants obtained by recalculation and After the vehicle to be positioned makes a positioning request, all positioning terminal nodes within the direct communication range with it will respond by sending the average value of a and n of the area where the positioning terminal node is located to the vehicle to be positioned. Terminal nodes are usually divided into multiple areas, so this type of positioning terminal node will send the average value of multiple a and n and the number of the corresponding area, and the smart phone of the vehicle to be positioned receives the feedback from the nearby positioning terminal node After the information is displayed, select the area with the highest number of occurrences to proceed to the next step;

步骤5:通过下式计算进入某区域的待定位车辆与选定区域的定位终端结点间的距离d,Step 5: Calculate the distance d between the vehicle to be positioned entering a certain area and the positioning terminal node of the selected area by the following formula,

dd ii == 1010 -- (( aa ‾‾ ++ RSSIRSSI ii ‾‾ )) 1010 nno ‾‾

其中,a和n取步骤4中的平均值,为接收信号平均功率。Among them, a and n take the average value in step 4, is the average power of the received signal.

步骤6:计算出3个以上上述的待定位车辆与定位终端结点间的距离后,通过三角形质心定位算法或极大似然估计法计算出待定位车辆的坐标。Step 6: After calculating the distance between more than 3 above-mentioned vehicles to be positioned and the positioning terminal node, calculate the coordinates of the vehicle to be positioned by the triangle centroid positioning algorithm or the maximum likelihood estimation method.

各所述定位终端结点上均设置有功率放大模块。Each of the positioning terminal nodes is provided with a power amplification module.

所述步骤1中的终端结点还包括用于监测车位状态的监测终端结点,定位终端结点与监测终端结点之间、相邻两监测终端结点之间无通信;所述路由器结点设置在若干个监测终端结点的包围圈的中心位置,路由器结点用于中继通信。The terminal node in described step 1 also comprises the monitoring terminal node that is used to monitor parking space state, between location terminal node and monitoring terminal node, there is no communication between adjacent two monitoring terminal nodes; Described router node The point is set at the center of the encirclement of several monitoring terminal nodes, and the router node is used to relay the communication.

所述步骤2中,中心结点和与之相邻的路由器结点相通信,中心结点安装在停车场的入口处、或出口处、或中心位置处。In the step 2, the central node communicates with its adjacent router node, and the central node is installed at the entrance, exit, or center of the parking lot.

由于采用了上述技术方案,本发明所取得的有益效果为:Owing to adopting above-mentioned technical scheme, the beneficial effect that the present invention obtains is:

本发明提供了一种室内停车场环境下分时分区域的车辆定位算法,有助于将车辆快速、准确地引导到事先规划好的目标车位上,具有定位速度快、精度高、实时性好等显著优点。The present invention provides a vehicle positioning algorithm for time-division and area-division in an indoor parking lot environment, which helps to guide vehicles to a pre-planned target parking space quickly and accurately, and has the advantages of fast positioning speed, high precision, and good real-time performance. Significant advantages.

附图说明Description of drawings

图1为本发明的整体结构示意图。Figure 1 is a schematic diagram of the overall structure of the present invention.

图2为本发明的区域划分示意图。Fig. 2 is a schematic diagram of area division in the present invention.

图3为不考虑a和n随时间和空间变化的情况下,最大定位误差随时间变化的示意图。Fig. 3 is a schematic diagram of the variation of the maximum positioning error with time without considering the variation of a and n with time and space.

图4为不考虑a和n随时间变化,只考虑a和n随空间变化的情况下,最大定位误差随时间变化的示意图。Fig. 4 is a schematic diagram of the variation of the maximum positioning error with time when a and n are not considered to vary with time, but only a and n are considered to vary with space.

图5为既考虑a和n随时间变化,又考虑a和n随空间变化的情况下,最大定位误差随时间变化的示意图。Fig. 5 is a schematic diagram of the variation of the maximum positioning error with time when a and n are considered to vary with time, and a and n are also considered to vary with space.

具体实施方式Detailed ways

下面结合附图和具体的实施例对本发明作进一步的详细说明,但本发明并不限于这些实施例。The present invention will be further described in detail below in conjunction with the accompanying drawings and specific embodiments, but the present invention is not limited to these embodiments.

如图1和图2所示,室内停车场环境下基于ZigBee网络的车辆定位方法,该方法包括以下步骤:As shown in Figure 1 and Figure 2, the vehicle location method based on ZigBee network under indoor parking lot environment, this method comprises the following steps:

步骤1:在室内停车场内设置ZigBee无线传感网络系统,ZigBee无线传感网络系统,包括一个ZigBee网络,在ZigBee网络中设置一个中心结点、若干个路由器结点以及数量更多的终端结点,中心结点和与之相邻的路由器结点相通信,中心结点安装在停车场的入口处、或出口处、或中心位置处,终端结点包括用于对待定位车辆进行定位的定位终端结点;为减少定位终端结点的布设量,各所述定位终端结点上均设置有功率放大模块,从而扩大了定位终端结点的直接通信范围;所述终端结点还包括用于监测车位状态的监测终端结点,定位终端结点与监测终端结点之间、相邻两监测终端结点之间无通信;所述路由器结点设置在若干个监测终端结点的包围圈的中心位置,路由器结点用于中继通信;待定位车辆的驾乘者持有智能手机,智能手机上配有ZigBee模块,待定位车辆通过智能手机发出定位请求;Step 1: Set up the ZigBee wireless sensor network system in the indoor parking lot. The ZigBee wireless sensor network system includes a ZigBee network, and a central node, several router nodes and more terminal nodes are set in the ZigBee network. The central node communicates with the router node adjacent to it. The central node is installed at the entrance, exit, or central position of the parking lot. The terminal node includes a positioning system for locating the vehicle to be positioned. terminal node; in order to reduce the layout amount of the positioning terminal node, each of the positioning terminal nodes is provided with a power amplification module, thereby expanding the direct communication range of the positioning terminal node; the terminal node also includes a Monitoring the monitoring terminal node of the parking space state, there is no communication between the positioning terminal node and the monitoring terminal node, and between two adjacent monitoring terminal nodes; the router node is set in the center of the encirclement circle of several monitoring terminal nodes Location, the router node is used for relay communication; the driver and occupant of the vehicle to be positioned holds a smartphone, which is equipped with a ZigBee module, and the vehicle to be positioned sends a positioning request through the smartphone;

步骤2:按照室内停车场中车位的分布,对上述定位终端结点进行区域划分,使之划分成A、B、C、D……等若干个互不重叠的区域,每个区域呈三角形或四边形,定位终端结点安装在每个三角形或四边形的顶点上,各区域中,每个定位终端结点和与其临近的两个定位终端结点进行直接通信;Step 2: According to the distribution of parking spaces in the indoor parking lot, divide the above-mentioned positioning terminal nodes into several non-overlapping areas such as A, B, C, D, etc., and each area is triangular or Quadrilateral, the positioning terminal node is installed on the vertices of each triangle or quadrilateral, and in each area, each positioning terminal node communicates directly with two adjacent positioning terminal nodes;

步骤3:对上述每个区域的各定位终端结点的坐标进行精确测量后计算相邻两顶点间的距离,并根据下式推算得出区域中某一个定位终端结点周围环境的两个待定常数ai、niStep 3: After accurately measuring the coordinates of each positioning terminal node in each of the above areas, calculate the distance between two adjacent vertices, and calculate the two undetermined values of the surrounding environment of a certain positioning terminal node in the area according to the following formula: Constants a i , n i :

RSSIRSSI ‾‾ ii 11 [[ dBmdBm ]] == -- (( aa ii ++ 1010 nno ii ·· lglg dd ii 11 )) RSSIRSSI ‾‾ ii 22 [[ dBmdBm ]] == -- (( aa ii ++ 1010 nno ii ·· lglg dd ii 22 ))

其中,di1和di2为区域中第i个定位终端结点与周围第1个和第2个相邻定位终端结点间的距离,为区域中第i个定位终端结点收到周围第1个和第2个相邻定位终端结点所发送信号的平均功率;Among them, d i1 and d i2 are the distances between the i-th positioning terminal node in the area and the surrounding first and second adjacent positioning terminal nodes, and is the average power of the signal sent by the i-th positioning terminal node in the area received by the surrounding first and second adjacent positioning terminal nodes;

步骤4:根据计算出的距离值di1和di2,以及测量值计算出ai和ni,各区域的定位终端结点计算出自己的ai和ni后,通过路由器结点传输给中心结点,中心结点计算出各区域a和n的平均值后传输给各定位终端结点,各定位终端结点在自带的数据存储装置中存储自己所属区域(通常有多个)a和n的平均值,得出:Step 4: Based on the calculated distance values d i1 and d i2 , and the measured values and After calculating a i and n i , the positioning terminal nodes in each area calculate their own a i and n i , and transmit them to the central node through the router node, and the central node calculates the average value of a and n in each area It is transmitted to each positioning terminal node, and each positioning terminal node stores the average value of a and n of its own area (usually multiple) in its own data storage device, and draws:

aa ‾‾ == 11 LL ΣΣ ii == 11 LL aa ii nno ‾‾ == 11 LL ΣΣ ii == 11 LL nno ii

其中,为该区域内两个待定常数的平均值,L为该区域内定位终端结点的个数;in, and is the average value of two undetermined constants in this area, and L is the number of positioning terminal nodes in this area;

待定位车辆提出定位请求后,所有与其处于能直接通信范围内的定位终端结点都会做出响应,向待定位车辆发送定位终端结点所处区域的a和n的平均值,因某些定位终端结点通常被划分至多个区域,所以这类定位终端结点会发送多个a和n的平均值以及所对应区域的编号,待定位车辆的智能手机接收到附近的定位终端结点反馈的信息后,选择编号出现次数最多的区域进行下一步骤;After the vehicle to be positioned makes a positioning request, all positioning terminal nodes within the direct communication range with it will respond by sending the average value of a and n of the area where the positioning terminal node is located to the vehicle to be positioned. Terminal nodes are usually divided into multiple areas, so this type of positioning terminal node will send the average value of multiple a and n and the number of the corresponding area, and the smart phone of the vehicle to be positioned receives the feedback from the nearby positioning terminal node After the information is displayed, select the area with the highest number of occurrences to proceed to the next step;

步骤5:通过下式计算进入某区域的待定位车辆与选定区域的定位终端结点间的距离d,Step 5: Calculate the distance d between the vehicle to be positioned entering a certain area and the positioning terminal node of the selected area by the following formula,

dd ii == 1010 -- (( aa ‾‾ ++ RSSIRSSI ii ‾‾ )) 1010 nno ‾‾

其中,a和n取步骤4中的平均值,为接收信号平均功率。Among them, a and n take the average value in step 4, is the average power of the received signal.

步骤6:计算出3个以上上述的待定位车辆与定位终端结点间的距离后,通过三角形质心定位算法或极大似然估计法计算出待定位车辆的坐标;在不同的时间段,随着停车场中车辆的进出等因素,所有区域中的两个待定常数会随时间发生变化,显然,若较频繁地对它们进行更新,可以保证对车辆的定位精度,但将明显增加组成各区域的顶点的功耗。为此,只有当在停车场的入口或出口处检测到有车辆进入或离开时,由中心结点通过路由器结点向各定位终端结点发出命令,各定位终端结点与邻近的定位终端结点通信,对ai和ni的值进行更新,根据重新计算得出的待定常数即可计算得出待定位车辆与定位终端结点间的距离d,进而计算得出待定位车辆的坐标。Step 6: After calculating the distance between more than 3 above-mentioned vehicles to be positioned and the positioning terminal node, calculate the coordinates of the vehicle to be positioned by the triangle centroid positioning algorithm or the maximum likelihood estimation method; in different time periods, follow According to factors such as the entry and exit of vehicles in the parking lot, the two undetermined constants in all areas will change with time. Obviously, if they are updated more frequently, the positioning accuracy of the vehicle can be guaranteed, but it will obviously increase the number of components in each area. The power consumption of the vertices. For this reason, only when a vehicle enters or leaves is detected at the entrance or exit of the parking lot, the central node sends commands to each positioning terminal node through the router node, and each positioning terminal node connects with the adjacent positioning terminal node. point communication, update the values of a i and ni i , according to the undetermined constants obtained by recalculation and The distance d between the vehicle to be positioned and the positioning terminal node can be calculated, and then the coordinates of the vehicle to be positioned can be calculated.

需要说明的是,待定位车辆上智能手机连接的ZigBee模块,其定位请求的广播信息发出后,路由器结点不转发该信息,只有处于和待定位结点能直接通信范围内的定位终端结点会接收该信息并作出响应;对待定位车辆进行定位时,定位终端结点会将自己的精确坐标及该区域内a和n的平均值发送给智能手机的ZigBee模块,智能手机中安装有位置解算软件,计算出待定位车辆的实时位置。之后,智能手机根据内置的停车场的内部地图和计算出的车辆位置,将车辆逐步地引导到目标车位。It should be noted that for the ZigBee module connected to the smartphone on the vehicle to be positioned, after the broadcast information of the positioning request is sent, the router node does not forward the information, and only the positioning terminal node within the direct communication range with the node to be positioned It will receive the information and make a response; when positioning the vehicle to be positioned, the positioning terminal node will send its precise coordinates and the average value of a and n in the area to the ZigBee module of the smartphone, and the location solution is installed in the smartphone Calculation software to calculate the real-time position of the vehicle to be positioned. After that, the smartphone guides the vehicle step by step to the target parking space based on the built-in internal map of the parking lot and the calculated vehicle position.

为验证所设计的算法可提高定位精度,在图2的基础上设定了一些参数,并利用Matlab语言编写了计算机模拟程序。所使用的具体参数如表1所示:In order to verify that the designed algorithm can improve the positioning accuracy, some parameters are set on the basis of Figure 2, and a computer simulation program is written using Matlab language. The specific parameters used are shown in Table 1:

表1 Matlab模拟程序所使用的参数Table 1 Parameters used by Matlab simulation program

除表1中所列的参数以外,设定信噪比为30;每个区域中为计算某定位终端结点与待定位结点间的距离值,所进行的测量次数为100次。In addition to the parameters listed in Table 1, the signal-to-noise ratio is set to 30; in each area to calculate the distance between a positioning terminal node and the node to be located, the number of measurements is 100 times.

图3为不考虑a和n随时间和空间变化的情况下(即:在利用步骤5中的公式计算距离的过程中,各区域的a和n统一使用表1中第3列和第4列的平均值),图3中A至K共11个区域在0至100小时内的最大定位误差随时间的变化。从图中可看出,误差值很大,最大的误差已达35m。Figure 3 does not consider the change of a and n with time and space (that is, in the process of calculating the distance using the formula in step 5, a and n in each area are uniformly used in the third and fourth columns of Table 1 The average value of ), the maximum positioning error of 11 areas from A to K in Fig. 3 changes with time from 0 to 100 hours. It can be seen from the figure that the error value is very large, and the largest error has reached 35m.

图4为只考虑a和n随空间的变化,但不考虑它们随时间变化的情况下(即:在利用步骤5中的公式计算距离的过程中,各区域的a和n,分别取表1中该区域对应行第3列和第4列的值),图3中A至K共11个区域在0至100小时内的最大定位误差随时间的变化。与图4相比,误差有所减小,但最大误差仍达近30m。Figure 4 only considers the changes of a and n with space, but does not consider their changes with time (that is, in the process of calculating the distance using the formula in step 5, the a and n of each area are respectively taken from Table 1 This area corresponds to the values in the third and fourth columns of the row), and the maximum positioning error of the 11 areas from A to K in Figure 3 changes with time from 0 to 100 hours. Compared with Figure 4, the error has been reduced, but the maximum error is still nearly 30m.

图5为既考虑a和n随空间的变化,也考虑它们随时间变化的情况下(即:在利用步骤5中的公式计算距离的过程中,各区域的a和n,分别取当前时刻a和n的即时测量值),图3中A至K共11个区域在0至100小时内的最大定位误差随时间的变化。从图中可看出,误差显著减小,基本控制在2m以内,取得了非常理想的定位效果。Figure 5 is a case where not only the variation of a and n with space, but also their variation with time are considered (that is, in the process of calculating the distance using the formula in step 5, the a and n of each area are respectively taken as the current moment a and n), the maximum positioning error of 11 areas from A to K in Fig. 3 changes with time from 0 to 100 hours. It can be seen from the figure that the error is significantly reduced, basically controlled within 2m, and a very ideal positioning effect has been achieved.

由此可见,本发明中的定位方法显著地提高了室内停车场环境下的车辆定位精度,有助于将车辆快速、准确地引导至事先规划好的目标车位上。It can be seen that the positioning method in the present invention significantly improves the positioning accuracy of the vehicle in the indoor parking lot environment, and helps to guide the vehicle to the target parking space planned in advance quickly and accurately.

本发明中未述及的部分采用或借鉴已有技术即可实现。The parts not mentioned in the present invention can be realized by adopting or referring to the prior art.

本文中所描述的具体实施例仅仅是对本发明的精神所作的举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。The specific embodiments described herein are merely illustrative of the spirit of the invention. Those skilled in the art to which the present invention belongs can make various modifications or supplements to the described specific embodiments or adopt similar methods to replace them, but they will not deviate from the spirit of the present invention or go beyond the definition of the appended claims range.

Claims (4)

1. A vehicle positioning method based on a ZigBee network in an indoor parking lot environment is characterized in that: the method comprises the following steps:
step 1: the method comprises the steps that a ZigBee wireless sensing network system is arranged in an indoor parking lot and comprises a ZigBee network, a central node, a plurality of router nodes and a larger number of terminal nodes are arranged in the ZigBee network, and the terminal nodes comprise positioning terminal nodes used for positioning vehicles to be positioned; a driver of a vehicle to be positioned holds a smart phone, the smart phone is provided with a ZigBee module, and the vehicle to be positioned sends a positioning request through the smart phone;
step 2: according to the distribution of parking spaces in an indoor parking lot, carrying out region division on the positioning terminal nodes to enable the positioning terminal nodes to be divided into a plurality of non-overlapping regions, numbering each region, wherein each region is triangular or quadrilateral, the positioning terminal nodes are installed on the top points of each triangle or quadrilateral, and in each region, each positioning terminal node is in direct communication with two positioning terminal nodes adjacent to the positioning terminal node;
and step 3: accurately measuring the coordinates of each positioning terminal node of each region, calculating the distance between two adjacent vertexes, and calculating two undetermined constants a of the surrounding environment of one positioning terminal node in the region according to the following formulai、ni
<math> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mover> <mi>RSSI</mi> <mo>&OverBar;</mo> </mover> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mo>[</mo> <mi>dBm</mi> <mo>]</mo> <mo>=</mo> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mi>i</mi> </msub> <mo>+</mo> <mn>10</mn> <msub> <mi>n</mi> <mi>i</mi> </msub> <mo>&CenterDot;</mo> <mi>lg</mi> <msub> <mi>d</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>RSSI</mi> <mo>&OverBar;</mo> </mover> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> <mo>[</mo> <mi>dBm</mi> <mo>]</mo> <mo>=</mo> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mi>i</mi> </msub> <mo>+</mo> <mn>10</mn> <msub> <mi>n</mi> <mi>i</mi> </msub> <mo>&CenterDot;</mo> <mi>lg</mi> <msub> <mi>d</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </math>
Wherein d isi1And di2For the distance between the ith positioning terminal node and the 1 st and 2 nd adjacent positioning terminal nodes in the area,andreceiving average power of signals sent by 1 st and 2 nd adjacent positioning terminal nodes for the ith positioning terminal node in the area;
and 4, step 4: according to the calculated distance value di1And di2And measured valueAndcalculate aiAnd niAnd a of all the positioning terminal nodes in the areaiAnd niTaking the average as the value of a and n in the area to obtain:
<math> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mover> <mi>a</mi> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>L</mi> </munderover> <msub> <mi>a</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mover> <mi>n</mi> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>L</mi> </munderover> <msub> <mi>n</mi> <mi>i</mi> </msub> </mtd> </mtr> </mtable> </mfenced> </math>
wherein,andthe average value of two undetermined constants in the region is obtained, and L is the number of positioning terminal nodes in the region;
when detecting that a vehicle enters or leaves at the entrance or exit of the parking lot, the central node sends a command to each positioning terminal node through the router node, and each positioning terminal node communicates with the adjacent positioning terminal node, and the pair aiAnd niIs updated and recalculated to obtain the undetermined constantAndafter a vehicle to be positioned provides a positioning request, all positioning terminal nodes in a direct communication range respond, and send the average values of a and n of the areas where the positioning terminal nodes are located to the vehicle to be positioned;
and 5: calculating the distance d between a vehicle to be positioned entering a certain area and a positioning terminal node of the selected area by the following formula,
<math> <mrow> <msub> <mi>d</mi> <mi>i</mi> </msub> <mo>=</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mfrac> <mrow> <mo>(</mo> <mover> <mi>a</mi> <mo>&OverBar;</mo> </mover> <mo>+</mo> <msub> <mover> <mi>RSSI</mi> <mo>&OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mrow> <mn>10</mn> <mover> <mi>n</mi> <mo>&OverBar;</mo> </mover> </mrow> </mfrac> </mrow> </msup> </mrow> </math>
wherein a and n are the average values in step 4,is the received signal average power.
Step 6: after the distances between more than 3 vehicles to be positioned and the positioning terminal nodes are calculated, the coordinates of the vehicles to be positioned are calculated through a triangular centroid positioning algorithm or a maximum likelihood estimation method.
2. The vehicle positioning method based on the ZigBee network in the indoor parking lot environment according to claim 1, wherein: and each positioning terminal node is provided with a power amplification module.
3. The vehicle positioning method based on the ZigBee network in the indoor parking lot environment according to claim 1, wherein: the terminal nodes in the step 1 further comprise monitoring terminal nodes for monitoring the parking space state, and no communication exists between the positioning terminal nodes and the monitoring terminal nodes and between two adjacent monitoring terminal nodes; the router nodes are arranged in the center positions of the surrounding rings of the monitoring terminal nodes and used for relay communication.
4. The vehicle positioning method based on the ZigBee network in the indoor parking lot environment according to claim 1, wherein: in the step 1, the central node communicates with the adjacent router node, and the central node is installed at an entrance, an exit or a central position of the parking lot.
CN201510231092.9A 2015-05-08 2015-05-08 Vehicle positioning method based on ZigBee network in indoor parking lot environment Expired - Fee Related CN104809908B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510231092.9A CN104809908B (en) 2015-05-08 2015-05-08 Vehicle positioning method based on ZigBee network in indoor parking lot environment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510231092.9A CN104809908B (en) 2015-05-08 2015-05-08 Vehicle positioning method based on ZigBee network in indoor parking lot environment

Publications (2)

Publication Number Publication Date
CN104809908A true CN104809908A (en) 2015-07-29
CN104809908B CN104809908B (en) 2017-01-25

Family

ID=53694703

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510231092.9A Expired - Fee Related CN104809908B (en) 2015-05-08 2015-05-08 Vehicle positioning method based on ZigBee network in indoor parking lot environment

Country Status (1)

Country Link
CN (1) CN104809908B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106023630A (en) * 2015-08-28 2016-10-12 千寻位置网络有限公司 Method and device for automatically judging parking position
CN106683470A (en) * 2016-11-17 2017-05-17 南昌大学 Vehicle positioning method oriented to parking lot
CN106991842A (en) * 2017-05-03 2017-07-28 常州信息职业技术学院 A kind of underground parking parking robot parking carries car localization method
CN107274713A (en) * 2017-07-11 2017-10-20 东华大学 A kind of parking induction method based on cluster
CN110444043A (en) * 2019-09-16 2019-11-12 深圳普智联科机器人技术有限公司 A kind of parking stall cruising inspection system and its method based on location technology
CN111179634A (en) * 2019-12-04 2020-05-19 浙江科技学院 Indoor parking system based on Wi-Fi position fingerprint
CN111818636A (en) * 2020-06-03 2020-10-23 哈尔滨工业大学(威海) A vehicle-mounted bluetooth positioning system and positioning method thereof
CN114141054A (en) * 2021-12-07 2022-03-04 安徽一流云停科技有限公司 Underground parking lot management system based on ZigBee-Bluetooth technology
CN116132953A (en) * 2021-11-15 2023-05-16 大唐移动通信设备有限公司 Positioning method, positioning device and terminal
CN116386378A (en) * 2023-04-03 2023-07-04 惠州市德赛西威汽车电子股份有限公司 A parking space guidance system, method, electronic device and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102883428A (en) * 2012-08-07 2013-01-16 广州星博信息技术有限公司 ZigBee wireless sensor network-based node positioning method
CN103680201A (en) * 2013-12-10 2014-03-26 中国石油大学(华东) Intelligent city parking system
KR101397648B1 (en) * 2013-03-18 2014-05-26 인하대학교 산학협력단 Method for arranging zigbee node in ship and apparatus thereof
CN104244462A (en) * 2013-06-14 2014-12-24 江南大学 Distance measurement positioning method for wireless sensor network based on RSSI

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102883428A (en) * 2012-08-07 2013-01-16 广州星博信息技术有限公司 ZigBee wireless sensor network-based node positioning method
KR101397648B1 (en) * 2013-03-18 2014-05-26 인하대학교 산학협력단 Method for arranging zigbee node in ship and apparatus thereof
CN104244462A (en) * 2013-06-14 2014-12-24 江南大学 Distance measurement positioning method for wireless sensor network based on RSSI
CN103680201A (en) * 2013-12-10 2014-03-26 中国石油大学(华东) Intelligent city parking system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张美燕等: "无线Zigbee传感器RSSI定位技术研究", 《计算机技术与发展》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106023630A (en) * 2015-08-28 2016-10-12 千寻位置网络有限公司 Method and device for automatically judging parking position
CN106023630B (en) * 2015-08-28 2018-11-20 千寻位置网络有限公司 The method and device thereof of automatic discrimination parking position
CN106683470A (en) * 2016-11-17 2017-05-17 南昌大学 Vehicle positioning method oriented to parking lot
CN106683470B (en) * 2016-11-17 2019-11-01 南昌大学 A kind of bus location method towards parking lot
CN106991842A (en) * 2017-05-03 2017-07-28 常州信息职业技术学院 A kind of underground parking parking robot parking carries car localization method
CN107274713A (en) * 2017-07-11 2017-10-20 东华大学 A kind of parking induction method based on cluster
CN110444043A (en) * 2019-09-16 2019-11-12 深圳普智联科机器人技术有限公司 A kind of parking stall cruising inspection system and its method based on location technology
CN111179634A (en) * 2019-12-04 2020-05-19 浙江科技学院 Indoor parking system based on Wi-Fi position fingerprint
CN111179634B (en) * 2019-12-04 2021-04-27 浙江科技学院 Indoor parking system based on Wi-Fi position fingerprint
CN111818636A (en) * 2020-06-03 2020-10-23 哈尔滨工业大学(威海) A vehicle-mounted bluetooth positioning system and positioning method thereof
CN111818636B (en) * 2020-06-03 2022-01-28 哈尔滨工业大学(威海) Vehicle-mounted Bluetooth positioning system and positioning method thereof
CN116132953A (en) * 2021-11-15 2023-05-16 大唐移动通信设备有限公司 Positioning method, positioning device and terminal
CN114141054A (en) * 2021-12-07 2022-03-04 安徽一流云停科技有限公司 Underground parking lot management system based on ZigBee-Bluetooth technology
CN114141054B (en) * 2021-12-07 2022-11-08 安徽一流云停科技有限公司 Underground parking lot management system based on ZigBee-Bluetooth technology
CN116386378A (en) * 2023-04-03 2023-07-04 惠州市德赛西威汽车电子股份有限公司 A parking space guidance system, method, electronic device and storage medium

Also Published As

Publication number Publication date
CN104809908B (en) 2017-01-25

Similar Documents

Publication Publication Date Title
CN104809908B (en) Vehicle positioning method based on ZigBee network in indoor parking lot environment
CN104571112B (en) Pilotless automobile lateral control method based on turning curvature estimation
RU2756439C1 (en) Determination of localisation for operation of a vehicle
US20170307400A1 (en) Device for providing electric-moving-body information and method for providing electric-moving-body information
CN112660148B (en) Method, apparatus, device, and medium for determining turning path of vehicle
CN107044858B (en) Energy efficient navigation route for electric transportation device
CN104391445B (en) Fleet&#39;s collaboration autonomous control method based on observer
SE538458C2 (en) Method, apparatus and system comprising the apparatus for supporting the creation of vehicle trains
CN103256934B (en) Road information provides equipment
CN103366593A (en) Intelligent parking navigation system based on cloud computing
CN103010265B (en) Be applicable to the static train locating method of CBTC system
CN106991842B (en) Parking robot parking and lifting positioning method for underground parking lot
CN204330354U (en) Automatic driving vehicle remote controlled braking Performance Test System
CN107271951A (en) UAV Navigation System and its air navigation aid based on WiFi fingerprint locations
CN103454614B (en) Based on particle filter vehicle between relative positioning method
CN102209379B (en) RSSI-based method for positioning wireless sensor network node
CN104200672A (en) Traffic intersection queuing length detection method and system based on multi-sensor fusion
US20240219507A1 (en) Positioning Calibration Method and System, and Storage Medium
CN103632558A (en) Bionic swarm intelligence-based real-time positioning navigation and motion control method and system for moving vehicle
CN102992196A (en) Implementation method of warehouse gantry crane running positioning system based on ZigBee technology
CN106604236A (en) Indoor positioning system
Zhou et al. RSSI-based mobile robot navigation in grid-pattern wireless sensor network
CN106364526B (en) A kind of tunnel construction new energy rail electric locomotive real-time high-precision automatic station-keeping system
CN104515528A (en) Single-point map matching method based on road section accumulation probability
KR20190056262A (en) Method for estimating position of vehicle and the vehicle

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170125