CN110958600B - 基于轨迹相似性的区域人口中一机多卡用户数量判定方法 - Google Patents
基于轨迹相似性的区域人口中一机多卡用户数量判定方法 Download PDFInfo
- Publication number
- CN110958600B CN110958600B CN201811126910.9A CN201811126910A CN110958600B CN 110958600 B CN110958600 B CN 110958600B CN 201811126910 A CN201811126910 A CN 201811126910A CN 110958600 B CN110958600 B CN 110958600B
- Authority
- CN
- China
- Prior art keywords
- imsi
- base stations
- imsis
- users
- machine
- 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
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W8/00—Network data management
- H04W8/18—Processing of user or subscriber data, e.g. subscribed services, user preferences or user profiles; Transfer of user or subscriber data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Life Sciences & Earth Sciences (AREA)
- Algebra (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提供一种基于轨迹相似性的区域人口中一机多卡用户数量判定方法,在居住地在统计区域的常驻人口随机抽取一定数量的IMSI,根据抽样数据中每个IMSI在统计时间段内访问过的基站数和在每个基站的驻留时间,确定每个IMSI的移动轨迹;逐一比对一个IMSI与其他IMSI的移动轨迹,筛选出候选一机多卡用户;采用超几何分布对候选一机多卡用户进行检验,并对检验结果进行多重校正,判定出抽样数据中的一机多卡用户。本发明用于通过移动大数据进行区域人口统计,用判定出的统计区域内一机多卡用户对人口统计结果进行校正,可提高人口统计的准确性。
Description
技术领域
本发明属于移动大数据统计分析应用技术领域,涉及运用移动大数据进行人口统计的方法,具体涉及一种基于轨迹相似性的区域人口中一机多卡用户数量的判定方法。
背景技术
采用移动通信大数据对人口进行监测和统计,可以有效地估计人口规模、掌握人口流向、及时对区域人口的发展趋势进行预警。运用大数据开展人口监测和统计工作中所使用的移动通信大数据是基于移动通信终端设备进行采集和分析的,其前提是假设一个手机用户对应一个手机号码。但现实中存在较多的一机多卡现象,即一个手机用户可能同时携带多个手机号码。一机多卡问题已经严重影响了数据统计的准确性和可靠性,但目前并没有任何有效的解决手段。
发明内容
本发明的目的在于针对现实中存在的一机多卡现象,提供一种基于轨迹相似性的一机多卡的判别方法,对一个用户拥有多个移动号码的情况进行有效识别,以提高基于移动通信大数据进行人口进行监测和统计的准确性。
为实现上述目的,本发明提供以下技术方案:
一种基于轨迹相似性的区域人口中一机多卡用户数量判定方法,其特征在于,包括以下步骤:
(1)基于移动运营商数据,在指定城市的常驻人口中,随机抽取一定数量的居住地在该统计区域的IMSI作为抽样数据;
(2)根据抽样数据中每个IMSI在统计时间段内访问过的基站数和在每个基站的驻留时间,确定每个IMSI的移动轨迹;统计出所有抽样IMSI访问过的基站总数,以及每个基站的IMSI访问数量;
(3)数据预处理:根据每个IMSI的移动轨迹,将抽样数据中在一天内只访问过一个基站的IMSI删除;在上述基站总数中将每天IMSI访问数量超过一定值的基站删除;
(4)对预处理后的IMSI数据,将每个IMSI的移动轨迹与其他IMSI的移动轨迹进行逐一比对,筛选出候选一机多卡用户;筛选方法如下:
a.一个IMSI与另一IMSI一个月内共同访问的基站数量超过一定数量;
b.一个IMSI与另一IMSI每天在共同访问的基站的驻留重叠时长累计超过一定值;
满足上述条件,则将这两个IMSI作为一候选IMSI对;
(5)基于超几何分布对每个候选IMSI对是否属于一机多卡用户进行判定:以经过数据预处理后的基站总数作为N,两个IMSI分别访问的基站数作为m和n,被两个IMSI共同访问的基站数目为X,在两个IMSI访问每个基站概率独立相等的假设下,按以下公式计算被两个IMSI共同访问的基站数目X≥k的概率P:
当计算所得结果小于预设的判定阈值时,则判定该候选IMSI对为一机多卡用户;
(6)按上述方法判定出抽样数据中所有一机多卡用户的数量,并计算出其在抽样数据中的比例,并可据此推算出统计区域内常驻人口中的一机多卡用户数量。
进一步,上述基于轨迹相似性的区域人口中一机多卡用户数量判定方法,所述步骤(5)中,采用bonferroni校正法对所述的判定阈值进行校正:将所述判定阈值除以抽样数据总数,所得的值作为校正阈值;当经过超几何分布计算所得的P值小于所述校正阈值时,则判定该IMSI对为一机多卡用户。
本发明通过把轨迹数据转化成IMSI对数据,并统计每个IMSI对每天来自相同基站的时长及每个月来自相同基站的个数,采用超几何检验方法,统计出每个IMSI对是否属于同一个用户(一机多卡)的概率值,对达到设定概率阈值的IMSI对,判定为一机多卡用户;进一步,还可采用多重检验使判断精度进一步提高。本发明能有效判断出两个IMSI是否属于一机多卡用户,可对传统的假设一个手机用户对应一个IMSI的人口统计方法进行优化,从而提高通过移动大数据进行人口统计的准确性。
附图说明
图1是本发明的流程图;
图2是以抽样IMSI在一天中访问的基站数目为基准,对用户的活动范围进行统计图。
具体实施方式
本发明所采用的数据来自移动运营商的信令数据,包括:用户手机号码-IMSI(国际移动用户识别码:International Mobile Subscriber Identification Number);位置区识别码lac:用于标识不同的位置区;基站编号ci:与位置区识别码(lac)结合,用于识别网络中覆盖的小区;IMSI进入基站的时间、离开该基站的时间。上述数据来自同一移动运营商。一机多卡属于不同运营商的情况不在本发明的研究范围内。通过上述数据,可刻画出IMSI的移动轨迹。
本发明的基本思路是:基于移动信令数据,根据两个IMSI在一定时间内共同访问的基站数量,及在共同访问基站的驻留时间重叠时长,判断两个IMSI的轨迹相似度,若达到一定相似度,则认为两个IMSI有可能是一机多卡用户,然后采用超几何分布进行检验,检验结果小于设定的阈值,则判定这两个IMSI是一机多卡用户。
为提高判定结果的准确度,还可采用Bonferroni校正对判定进一步结果筛选。
本发明的实现主要基于以下理论:
1、超几何分布。
超几何分布(hypergeometric distribution)是统计学上一种离散概率分布,它描述了由有限个物件中抽出n个物件,成功抽出指定种类的物件的次数(不归还)。
在统计学的假设检验分析中,先假设样本数据来自零假设的总体,假设检验方法即根据零假设下检测统计量的分布计算出零假设的拒绝域,当样本统计量落在拒绝域即认为样本拒绝零假设,即零假设不真。而来自零假设的样本统计落入拒绝域的概率被称为显著性水平,习惯上设为5%。
本发明将超几何分布理论用于对两个IMSI共同访问的基站数目进行检验,零假设即为两个待分析的IMSI是完全独立的,其同时访问某个基站的现象是随机产生的。以两个IMSI同时访问的基站次数作为统计量,可以认为这个统计量是服从超几何分布的,基于超几何分布对IMSI对共同访问的基站数目X进行检验:以全部基站数作为N,两个IMSI分别访问的基站数作为m和n,在IMSI访问每个基站概率独立相等的假设下,被两个IMSI共同访问的基站数目为X,则两个IMSI共同访问基站数目X≥k的概率为:
将计算结果所得的P值与预先设定的判定阈值比较,对两个IMSI是否是一机多卡用户进行判定。该判定阈值可取超几何分布惯常的显著性水平,即5%。计算结果小于5%,即判定两个IMSI是一机多卡用户。
2、Bonferroni校正。
Bonferroni校正是一种较为严格的多重检验校正方法,即对于同一个数据集有n次(n>=2)假设检验时,那么用于每一假设的统计显著水平,应为仅检验一个假设时的显著水平的1/n。
bonferroni校正后阈值会极大提高,因此可以使用bonferroni校正后的阈值来对上述判定结果进行进一步筛选。
图1是本发明的实现方法流程图。本发明的具体实施方式如下:
(1)基于移动运营商数据,在指定城市的常驻人口中,随机抽取一定数量的居住地在该统计区域的IMSI作为抽样数据。
本发明要将每个IMSI与其余IMSI一一比对,对于人口数量较大的城市,如果将全部用户作为分析对象,计算量过大,考虑到计算时间问题,本发明对城市各个区域分别进行统计。统计区域可以是一个市辖区,也可以是一个街道,一个小区,或一个按特定要求指定的片区。采用随机抽样的方法在统计区域抽取一定数量的IMSI作为分析对象,以抽样代替总体,可大大提高运算速度。在本发明的一个具体实施例中,在城市某小区随机抽取了3000个IMSI作为研究对象。在具体应用时,可以考虑多次抽样以减小误差。
本发明主要应用于人口统计领域。城市人口统计中一般将常驻人口作为主要统计指标,因此本发明的统计对象由常驻人口中抽取。本发明对城市常驻人口的定义,是根据根据移动用户在指定城市中的居住时间而定,指一个月内,在该城市驻留超过15天且每天驻留时长超过10小时的移动用户。
本发明需要采集IMSI用户在一定时间段内的数据,该时间段一般至少为一个月。为保证数据的连续性,需要采集居住地在统计区域的用户。本发明对移动用户的居住地如下定义:移动用户在满足城市常驻人口的前提下,其每天晚上21点至次日凌晨7点时段,在某个基站驻留时间最长,则把此基站所在的区域定义为用户该日的居住地;按月统计用户在各区域的居住天数,居住天数最多的区域为用户本月的居住地。
上述移动用户在城市及某个基站的驻留时间可通过移动运营商的轨迹数据匹配出。
(2)根据抽样数据中每个IMSI在一定时间段内每天访问过的基站数和在每个基站的驻留时间,确定每个IMSI的移动轨迹;统计出所有抽样IMSI访问过的基站总数,以及每个基站每天的IMSI访问数量。上述时间段具体实施时可限定为一个月。
(3)数据预处理:根据每个IMSI的移动轨迹,将抽样数据中在一天内只访问过一个基站的IMSI删除;在上述基站总数中将每天IMSI访问数量超过一定值的基站删除。
不同的IMSI具有不同的行动规律。IMSI持有者活动范围较小,则在一天中该IMSI连接的基站数也较少,而IMSI持有者的活动范围较大,则其一天中该IMSI访问的基站数目也会比较多。这样,IMSI在一天中访问的基站数目,可以作为反应其持有者活动范围的一项指标。
图2所示,是本发明的一个具体实施例中,以抽样IMSI在一天中访问的基站数目为基准,对用户的活动范围进行统计图。
图2中纵轴为IMSI出现的频数,横轴为IMSI访问的基站数。从图2中可以看部分IMSI每天访问的基站数量较少,这种活动范围较小的用户,其轨迹数据对区分一机多卡问题意义较小,数据预处理过程中将这类一天之中只访问过一个基站的IMSI数据删除。
图2显示大部分的基站每天来访问的IMSI数量较小,但是少数基站每天有大量IMSI访问。造成该情况的最主要原因,是所选数据来自于同一小区。由于该类基站富集大量的IMSI信息,所以来自这些基站所预测为一机双卡的IMSI对的准确率较低,在数据预处理过程中将此类数据删除。具体方法是按所有基站一天之中访问IMSI数量从多到小排序,将IMSI访问量最大的前1%的基站删除。
(4)对预处理后的IMSI数据,将每个IMSI的移动轨迹与其他IMSI的移动轨迹进行逐一比对,筛选出候选一机多卡用户;候选一机多卡用户满足以下条件:
a.一个IMSI与另一IMSI一个月内共同访问的基站数量超过一定数量;该数量具体实施时设定为至少5个;
b.出一个IMSI与另一IMSI每天在共同访问的基站的驻留重叠时长超过一定值;具体实施时,重叠时长可设定为至少5小时。
满足上述条件,则将这两个IMSI作为一候选IMSI对。
(5)基于超几何分布对每个候选IMSI对是否属于一机多卡用户进行判定:以经过数据预处理后的基站总数作为N,两个IMSI分别访问的基站数作为m和n,被两个IMSI共同访问的基站数目为X,在两个IMSI访问每个基站概率独立相等的假设下,按以下公式计算被两个IMSI共同访问的基站数目X≥k的概率P:
预先设定一个判定阈值,当计算所得P值小于预设的判定阈值时,则判定该候选IMSI对为一机多卡用户。
上述判定阈值根据超几何分布的惯常阈值,一般设定为5%。即当计算所得P值小于5%时,则判定该候选IMSI对为一机多卡用户。
本发明具体实施时,为提高判定结果的准确性,可采用bonferroni校正法对所述的判定阈值进行校正,将所述判定阈值除以抽样数据总数,所得的值作为校正阈值。
假设在上述抽样数据中找到的候选一机多卡的IMSI对为N、判定阈值为5%的情况下,校正阈值为0.05/N。当经过超几何分布计算所得的P值小于0.05/N时,则判定该IMSI对为一机多卡用户。
(6)按上述方法判定出抽样数据中所有一机多卡的用户数量,并计算出其在抽样数据中的比值,得出抽样数据中一机多卡用户的比例,并可据此推算出统计区域内常驻人口中的一机多卡用户数量。
Claims (7)
1.基于轨迹相似性的区域人口中一机多卡用户数量的判定方法,其特征在于,包括以下步骤:
(1)基于移动运营商数据,在指定城市的常驻人口中,随机抽取一定数量的居住地在该统计区域的IMSI作为抽样数据;
(2)根据抽样数据中每个IMSI在统计时间段内访问过的基站数和每天在每个基站的驻留时间,确定每个IMSI的移动轨迹;统计出所有抽样IMSI访问过的基站总数,以及每个基站的IMSI访问数量;
(3)数据预处理:根据每个IMSI的移动轨迹,将抽样数据中在一天内只访问过一个基站的IMSI删除;在上述基站总数中将每天IMSI访问数量超过一定值的基站删除;
(4)对预处理后的IMSI数据,将每个IMSI的移动轨迹与其他IMSI的移动轨迹进行逐一比对,筛选出候选一机多卡用户,筛选方法如下:
a.一个IMSI与另一IMSI一个月内共同访问的基站数量超过一定数量;
b.一个IMSI与另一IMSI每天在共同访问的基站的驻留重叠时长累计超过一定值;
满足上述条件,则将这两个IMSI作为一候选IMSI对;
(5)基于超几何分布对每个候选IMSI对是否属于一机多卡用户进行判定:以经过数据预处理后的基站总数作为N,两个IMSI分别访问的基站数作为m和n,被两个IMSI共同访问的基站数目为X,在两个IMSI访问每个基站概率独立相等的假设下,按以下公式计算被两个IMSI共同访问的基站数目X≥k的概率P:
当计算所得结果小于预设的判定阈值时,则判定该候选IMSI对为一机多卡用户;
(6)按上述方法判定出抽样数据中所有一机多卡的用户数量,并计算出其在抽样数据中的比例,并可据此推算出统计区域内常驻人口中的一机多卡用户数量。
2.根据权利要求1所述的基于轨迹相似性的区域人口中一机多卡用户数量的判定方法,其特征在于:步骤(1)所述城市常驻人口,是指一个月内,在该城市驻留超过15天且每天驻留时长超过10小时的移动用户;所述的居住地,是指移动用户在满足城市常驻人口的前提下,其每天晚上21点至次日凌晨7点时段,在某个基站驻留时间最长,则把此基站所在的区域定义为用户该日的居住地;按月统计用户在各区域的居住天数,居住天数最多的区域为用户本月的居住地;移动用户在基站的驻留时间通过移动运营商的轨迹数据匹配出。
3.根据权利要求1所述的基于轨迹相似性的区域人口中一机多卡用户数量的判定方法,其特征在于:步骤(2)所述的统计时间段是一个月。
4.根据权利要求1所述的基于轨迹相似性的区域人口中一机多卡用户数量的判定方法,其特征在于:步骤(3)中,数据预处理时,将基站总数中每天IMSI访问量最大的前1%的基站删除。
5.根据权利要求1所述的基于轨迹相似性的区域人口中一机多卡用户数量的判定方法,其特征在于:步骤(4)中,将两个IMSI作为候选一机多卡用户,具体应满足以下条件:
a.一个IMSI与另一IMSI一个月内共同访问的基站数量超过5个;
b.一个IMSI与另一IMSI每天在共同访问的基站的驻留重叠时长累计超过5小时。
6.根据权利要求1所述的基于轨迹相似性的区域人口中一机多卡用户数量的判定方法,其特征在于:步骤(5)中所述的判定阈值为5%。
7.根据权利要求1所述的基于轨迹相似性的区域人口中一机多卡用户数量的判定方法,其特征在于:所述步骤(5)中,采用bonferroni校正法对所述的判定阈值进一步校正:将所述判定阈值除以抽样数据总数,所得的值作为校正阈值;当经过超几何分布计算所得的P值小于所述校正阈值时,则判定该IMSI对为一机多卡用户。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811126910.9A CN110958600B (zh) | 2018-09-26 | 2018-09-26 | 基于轨迹相似性的区域人口中一机多卡用户数量判定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811126910.9A CN110958600B (zh) | 2018-09-26 | 2018-09-26 | 基于轨迹相似性的区域人口中一机多卡用户数量判定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110958600A CN110958600A (zh) | 2020-04-03 |
CN110958600B true CN110958600B (zh) | 2022-06-28 |
Family
ID=69966243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811126910.9A Active CN110958600B (zh) | 2018-09-26 | 2018-09-26 | 基于轨迹相似性的区域人口中一机多卡用户数量判定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110958600B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111970685B (zh) * | 2020-10-23 | 2021-01-15 | 上海世脉信息科技有限公司 | 一种大数据环境下一人多卡识别方法 |
CN117150319A (zh) * | 2023-10-30 | 2023-12-01 | 北京艾瑞数智科技有限公司 | 一种一人多号的识别方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331404A (zh) * | 2013-07-22 | 2015-02-04 | 中国科学院深圳先进技术研究院 | 一种基于用户手机上网数据的用户行为预测方法和装置 |
CN104753728A (zh) * | 2013-12-27 | 2015-07-01 | 中国移动通信集团山东有限公司 | 一种用户流量监测方法及系统 |
CN106231555A (zh) * | 2016-08-03 | 2016-12-14 | 深圳市数字城市工程研究中心 | 一种人口空间分布的模拟方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101450776B1 (ko) * | 2007-12-04 | 2014-10-14 | 삼성전자주식회사 | Imsi 정보를 이용해 캠핑을 수행하는 단말기 및 그제어방법 |
-
2018
- 2018-09-26 CN CN201811126910.9A patent/CN110958600B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331404A (zh) * | 2013-07-22 | 2015-02-04 | 中国科学院深圳先进技术研究院 | 一种基于用户手机上网数据的用户行为预测方法和装置 |
CN104753728A (zh) * | 2013-12-27 | 2015-07-01 | 中国移动通信集团山东有限公司 | 一种用户流量监测方法及系统 |
CN106231555A (zh) * | 2016-08-03 | 2016-12-14 | 深圳市数字城市工程研究中心 | 一种人口空间分布的模拟方法及装置 |
Non-Patent Citations (1)
Title |
---|
王瑜等.基于网络活跃用户终端信息采集与分析的4G用户提升方案.《信息通信技术》.2017,(第05期), * |
Also Published As
Publication number | Publication date |
---|---|
CN110958600A (zh) | 2020-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110958599B (zh) | 一种基于轨迹相似性的一机多卡用户判别方法 | |
CN108156626B (zh) | 轨道交通无线网络质量评估方法、装置及介质 | |
CN110958600B (zh) | 基于轨迹相似性的区域人口中一机多卡用户数量判定方法 | |
CN109996258B (zh) | 无线网络利用率评估方法、装置、计算设备及存储介质 | |
CN109992605B (zh) | 基于手机信令数据的人口识别方法及系统 | |
CN104427549A (zh) | 一种网络问题分析方法及系统 | |
CN114173356A (zh) | 网络质量的检测方法、装置、设备及存储介质 | |
CN112001829A (zh) | 一种基于手机信令数据的人口分布判断方法 | |
CN111866847B (zh) | 移动通信网络数据分析方法、设备及计算机存储介质 | |
CN115665665A (zh) | 移动路径识别方法、识别装置、电子设备和可读存储介质 | |
CN111459702B (zh) | 基于mdt数据的室内分布系统故障监控方法和装置 | |
CN109996272B (zh) | 室分外泄信号的检测方法、装置、设备及介质 | |
CN110012436B (zh) | 用户位置确定方法、装置、设备及计算机可读存储介质 | |
CN110933601B (zh) | 目标区域确定方法、装置、设备和介质 | |
CN111372194A (zh) | 一种手机换卡用户智能化识别方法 | |
CN108573265B (zh) | 人流量统计方法及统计系统 | |
CN109982373B (zh) | Lte频繁回落点分析方法、装置、计算设备及存储介质 | |
CN102209336B (zh) | 一种基于邻区测量的邻区优化方法 | |
CN114828055A (zh) | 用户业务感知评估方法、装置、设备、介质和程序产品 | |
CN110621027B (zh) | 短驻留用户发现和短驻留原因分析方法 | |
CN109714771B (zh) | 一种基于信令数据的基站分析工作中的基站预处理方法 | |
CN115310735A (zh) | 商圈边界识别的方法、装置、设备及存储介质 | |
CN114205820A (zh) | 携带伪基站的可疑用户检测方法、装置、计算机设备 | |
CN114095947B (zh) | 一种网络覆盖评估方法、装置、电子设备及存储介质 | |
KR100687408B1 (ko) | 비동기식 wcdma망에서의 핸드오버 처리를 위한 네이버리스트 자동 최적화 방법 및 장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |