CN107749073B - 一种支持时序特征探查的统计地图绘制方法 - Google Patents

一种支持时序特征探查的统计地图绘制方法 Download PDF

Info

Publication number
CN107749073B
CN107749073B CN201610753059.7A CN201610753059A CN107749073B CN 107749073 B CN107749073 B CN 107749073B CN 201610753059 A CN201610753059 A CN 201610753059A CN 107749073 B CN107749073 B CN 107749073B
Authority
CN
China
Prior art keywords
data
width
time sequence
statistical
geographic area
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
Application number
CN201610753059.7A
Other languages
English (en)
Other versions
CN107749073A (zh
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.)
Computer Network Information Center of CAS
Original Assignee
Computer Network Information Center of CAS
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 Computer Network Information Center of CAS filed Critical Computer Network Information Center of CAS
Priority to CN201610753059.7A priority Critical patent/CN107749073B/zh
Publication of CN107749073A publication Critical patent/CN107749073A/zh
Application granted granted Critical
Publication of CN107749073B publication Critical patent/CN107749073B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/20—Drawing from basic elements, e.g. lines or circles
    • G06T11/203—Drawing of straight lines or curves
    • 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
    • G—PHYSICS
    • G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
    • G09B29/00—Maps; Plans; Charts; Diagrams, e.g. route diagram

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Probability & Statistics with Applications (AREA)
  • Software Systems (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Algebra (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Databases & Information Systems (AREA)
  • Operations Research (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • Educational Technology (AREA)
  • Processing Or Creating Images (AREA)
  • Instructional Devices (AREA)

Abstract

本发明涉及一种支持时序特征探查的统计地图绘制方法，该方法根据时序统计数据的地理区域信息以及时序特征，将时序统计数据绘制到地理区域内部，得到统计地图。进一步还包括将统计指标数据绘制到地理区域边缘内侧的步骤：首先根据时序统计数据，计算需要进行分析和呈现的统计指标；然后根据时序统计数据的地理区域信息，将统计指标数据绘制到地理区域边缘内侧。本发明能够在最大限度的利用空间的情况下，对统计数据的整体特征以及时序特征进行探查，能够实现对时序统计数据整体和时序特征的可视化。

Description

一种支持时序特征探查的统计地图绘制方法
技术领域
本发明属于制图学和信息可视化领域，涉及一种支持时序特征探查的统计地图绘制方法，与传统统计地图相比，该方法能够支持在最大限度的利用空间的情况下，对统计数据的整体特征以及时序特征进行探查。
背景技术
统计数据是对某个具体指标进行统计分析后得出的数据，通过对统计数据的可视化，能够帮助分析者发现数据背后的知识，从而辅助进行问题发现或决策制定。统计地图是运用统计数据反映制图对象数量特征的一种可视化方法，它能够将带有空间信息的统计数据在地图中进行呈现，在辅助分析上，具有非常积极的意义。
然而，对于时序统计数据，除需要了解其整体特征外，还需要了解时序变化特点。由于目前大部分统计地图主要通过颜色填充和区域变形(Cano R G,Buchin K,CastermansT,et al.Mosaic drawings and cartograms[C]//Computer Graphics Forum.2015,34(3):361-370.)两种方法来实现对统计数据的呈现，无法实现对时序数据在统计地图中的显示。有方法使用在地理区域中利用传统的时序数据可视化方法进行呈现(Andrienko N,Andrienko G.Interactive visual tools to explore spatio-temporal variation[C]//Proceedings of the working conference on Advanced visual interfaces.ACM,2004:417-420.)，但一方面这类方法无法对整体信息进行呈现，另一方面由于这些传统时序数据可视化方法需要较大的显示空间才能进行有效的表示，对于时序数据较多或区域较小时的呈现效果不好。有方法是用时间切片(small multiple)，将不同时间的统计数据依次排列在显示区域中(Fabrikant S I,Rebich-Hespanha S,Andrienko N,et al.Novelmethod to measure inference affordance in static small-multiple map displaysrepresenting dynamic processes[J].The Cartographic Journal,2013.)，但该方法降低了一个地图的呈现区域，而且该方法将不同时间点的数据分布于不同的显示区域，使得使用者难以发现数据背后的时序特征。有方法使用动画(animation)，在一个空间区域中，将不同时间对应的统计数据按时间顺序动态呈现(Boyandin I,Bertini E,Lalanne D.AQualitative Study on the Exploration of Temporal Changes in Flow Maps withAnimation and Small‐Multiples[C]//Computer Graphics Forum.BlackwellPublishing Ltd,2012,31(3pt2):1005-1014.)。该方法虽然能够保持一个地图的呈现区域的大小，但由于地图的动态变化，容易导致使用者视觉疲劳，并且也不易发现一个地图内不同区域的时序变化规律。
发明内容
本发明充分利用地理区域中的边缘区域和中心区域，结合基于像素的填充技术，给出一种统计地图绘制方法，实现对时序统计数据整体和时序特征的可视化。
本发明采用的技术方案如下：
一种支持时序特征探查的统计地图绘制方法，根据时序统计数据的地理区域信息以及时序特征，将时序统计数据绘制到地理区域内部，得到统计地图。
进一步地，还包括将统计指标数据绘制到地理区域边缘内侧的步骤：首先根据时序统计数据，计算需要进行分析和呈现的统计指标，所述统计指标包括均值、方差、最值中的一种或多种；然后根据时序统计数据的地理区域信息，将统计指标数据绘制到地理区域边缘内侧。
进一步地，所述将时序统计数据绘制到地理区域内部的实现方法为：
a)按时间顺序对各个地理区域内的统计数据进行排序，并获得各地理区域的时序数据T；
b)计算各地理区域的包围盒bbox(x,y,width,height)。其中x,y,width,height分别表示某地理区域的包围盒的左上角x坐标、左上角y坐标、宽度、高度；
c)对于任一地理区域i，其对应的时序数据为Ti，对应的包围盒为bbox(xi,yi,widthi,heighti)，在该区域中，使用不同宽度的矩形条代表Ti的每一个时间点的数据，对该区域进行填充。其中，矩形条的宽度可以根据所对应的时间点的数据重要程度、数值大小等进行定义。
进一步地，所述将时序统计数据绘制到地理区域内部的实现方法为：
a)按时间顺序对各个地理区域内的统计数据进行排序，并获得各地理区域的时序数据T；
b)计算各地理区域的包围盒bbox(x,y,width,height)。其中x,y,width,height分别表示某地理区域的包围盒的左上角x坐标、左上角y坐标、宽度、高度；
c)对于任一地理区域i，其对应的时序数据为Ti，对应的包围盒为bbox(xi,yi,widthi,heighti)，在该区域中，以等宽度为
Figure BDA0001097475740000021
的矩形条呈现时序数据从早到晚的数据；其中，Num(Ti)表示Ti的数据条目数目。
进一步地，所述将时序统计数据绘制到地理区域内部的实现方法为：
a)按时间顺序对各个地理区域内的统计数据进行排序，并获得各地理区域的时序数据T；
b)计算该地理区域的包围盒bbox(x,y,width,height)。其中x,y,width,height分别表示某地理区域的包围盒的左上角x坐标、左上角y坐标、宽度、高度；
c)对于任一地理区域i，其对应的时序数据为Ti，对应的包围盒为bbox(xi,yi,widthi,heighti)，在该区域中，以等面积为
Figure BDA0001097475740000031
的矩形条呈现时序数据从早到晚的数据；其中，Num(Ti)表示Ti的数据条目数目，每个矩形条的宽度为面积
Figure BDA0001097475740000032
对应的宽度。
进一步地，对边缘区域和中心区域使用相同或不同的颜色映射，来表示边缘区域或中心区域表征的数值由高到低。
进一步地，对不同地理区域的中心区域的时序表征使用不同的方向表示时序数据从早到晚，包括：从左至右的方向表示时序数据从早到晚，从上到下的方向表示时序数据从早到晚。
进一步地，对于地图上不同的地理区域，同时使用多种不同的方向来表示时序数据。
进一步地，所有相邻的地理区域使用不同的方向来表示时序数据。
进一步地，在确定每个地理区域的方向时，结合等面积呈现方法，以各个时序数据对应的等面积区域所占据的矩形宽度方差最小的方向作为最终方向。
本发明的有益效果如下：
本发明提供了一种支持时序特征探查的统计地图绘制方法，与传统统计地图相比，该方法能够在最大限度的利用空间的情况下，对统计数据的整体特征以及时序特征进行探查，能够实现对时序统计数据整体和时序特征的可视化。
附图说明
图1是本发明方法的步骤流程图。
图2是实施例中将统计指标数据绘制到地理区域边缘内侧的示意图。
图3是实施例中将时序数据绘制到地理区域内部的示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面通过具体实施例和附图，对本发明做进一步说明。
一个典型的时序统计数据如下面表1所示：
表1.时序统计数据
Figure BDA0001097475740000033
Figure BDA0001097475740000041
以上面的时序统计数据为例，本实施例的支持时序特征探查的统计地图绘制方法，其流程如图1所示，包括以下步骤：
1.根据时序统计数据，计算需要进行分析和呈现的统计指标。如均值、方差、最值等。
2.根据统计数据的地理区域信息，将统计指标数据绘制到地理区域边缘内侧。由于地理区域之间存在边缘重叠的现象，将统计指标数据绘制到地理区域边缘内侧，可以有效防止边缘重叠导致的视觉错误。以中国省份统计数据(统计指标数据为以省进行区分的数据)为例，其呈现效果如图2所示。
3.根据地理区域信息，以及时序特征，将时序数据绘制到地理区域内部。该步骤的操作步骤如下：
a.按时间顺序对各个地理区域内的统计数据进行排序，并获得各地理区域的时序数据T。
b.计算该地理区域的包围盒bbox(x,y,width,height)。其中x,y,width,height分别表示某地理区域的包围盒的左上角x坐标、左上角y坐标、宽度、高度。
c.对于任一地理区域i，其对应的时序数据为Ti，对应的包围盒为bbox(xi,yi,widthi,heighti)，在该区域中，以等宽度
Figure BDA0001097475740000042
的矩形条呈现时序数据从早到晚的数据。其中，Num(Ti)表示Ti的数据条目数目。以中国省份统计数据(统计指标数据为以省进行区分的数据)为例，其呈现效果如图3所示。另外可以使用不同宽度的矩形条代表Ti的每一个时间点的数据，对该区域进行填充，其中矩形条的宽度可以根据所对应的时间点的数据重要程度、数值大小等进行定义。
对于步骤3，与之相似的另一种实现方式为：
a.按时间顺序对各个地理区域内的统计数据进行排序，并获得各地理区域的时序数据T。
b.计算该地理区域的包围盒bbox(x,y,width,height)。其中x,y,width,height分别表示某地理区域的包围盒的左上角x坐标、左上角y坐标、宽度、高度。
c.对于任一地理区域i，其对应的时序数据为Ti，对应的包围盒为bbox(xi,yi,widthi,heighti)，在该区域中，以等面积
Figure BDA0001097475740000043
的矩形条呈现时序数据从早到晚的数据。其中，Num(Ti)表示Ti的数据条目数目。每个矩形条的宽度为面积
Figure BDA0001097475740000051
对应的宽度。
通过步骤1～3，可以绘制出利用地理区域中的边缘区域和中心区域的时序统计数据统计地图。基于该方法，可以进一步实现个性化的定制，包括：
1.对边缘区域和中心区域使用不同的颜色映射方法。在本发明步骤给出的方法中，使用了相同的颜色映射，比如红色到黄色(对应图1中的深色到浅色)表示某数值由高到低。在某些场景中，可以对两者使用不同的颜色映射，比如红色到黄色表示边缘区域表征的数值由高到低，绿色到蓝色表示中心区域表征的时序数值由高到低。
2.对中心区域的时序表征使用不同的方向。在本发明步骤给出的方法中，以从左至右的方向表示时序数据从早到晚。在某些场景中可以：
a.使用其它的方向。例如，使用从上到下的方向表示时序数据从早到晚。
b.使用多种不同的方向来表示时序数据。为便于区分，所有相邻的地理区域使用不同的方向。如A、B、C三个地理区域相互临接，但使用不同的方向来表示时序数据从早到晚。在确定每个区域的方向时，结合步骤3中的等面积呈现方法，以各个时序数据对应的等面积区域所占据的矩形宽度方差最小的方向(角度k所对应的方向)作为最终的方向。其计算公式如下：
Figure BDA0001097475740000052
其中，Varik表示地理区域i在角度k时的方差，widthik与heightik表示地理区域i在角度k时的包围盒宽度与高度，widthijk表示地理区域i的第j个时序序列在角度k时所占有的宽度，Num(Ti)表示地理区域Ti的时序数据条目数目。
对于地理区域i，通过该公式计算各个角度的方差，并选择方差最小的值对应的方向。
除上述时实施方式外，本发明也可不采用步骤1～2，即不将统计指标数据绘制到地理区域边缘内侧，而直接根据地理区域信息以及时序特征，将时序数据绘制到地理区域内部。
以上实施例仅用以说明本发明的技术方案而非对其进行限制，本领域的普通技术人员可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明的精神和范围，本发明的保护范围应以权利要求书所述为准。

Claims (8)

1.一种支持时序特征探查的统计地图绘制方法，其特征在于，根据时序统计数据的地理区域信息以及时序特征，将时序统计数据绘制到地理区域内部，得到统计地图；
所述将时序统计数据绘制到地理区域内部的实现方法为：
a)按时间顺序对各个地理区域内的统计数据进行排序，并获得各地理区域的时序数据T；
b)计算各地理区域的包围盒bbox(x,y,width,height)，其中x,y,width,height分别表示某地理区域的包围盒的左上角x坐标、左上角y坐标、宽度、高度；
c)对于任一地理区域i，其对应的时序数据为Ti，对应的包围盒为bbox(xi,yi,widthi,heighti)，在该区域中，使用不同宽度的矩形条代表Ti的每一个时间点的数据，对该区域进行填充；
其中，矩形条的宽度根据所对应的时间点的数据重要程度或数值大小进行定义；
或者，所述将时序统计数据绘制到地理区域内部的实现方法为：
a)按时间顺序对各个地理区域内的统计数据进行排序，并获得各地理区域的时序数据T；
b)计算各地理区域的包围盒bbox(x,y,width,height)，其中x,y,width,height分别表示某地理区域的包围盒的左上角x坐标、左上角y坐标、宽度、高度；
c)对于任一地理区域i，其对应的时序数据为Ti，对应的包围盒为bbox(xi,yi,widthi,heighti)，在该区域中，以等宽度为
Figure FDA0002814668920000011
的矩形条呈现时序数据从早到晚的数据；其中，Num(Ti)表示Ti的数据条目数目；
或者，所述将时序统计数据绘制到地理区域内部的实现方法为：
a)按时间顺序对各个地理区域内的统计数据进行排序，并获得各地理区域的时序数据T；
b)计算该地理区域的包围盒bbox(x,y,width,height)，其中x,y,width,height分别表示某地理区域的包围盒的左上角x坐标、左上角y坐标、宽度、高度；
c)对于任一地理区域i，其对应的时序数据为Ti，对应的包围盒为bbox(xi,yi,widthi,heighti)，在该区域中，以等面积为
Figure FDA0002814668920000012
的矩形条呈现时序数据从早到晚的数据；其中，Num(Ti)表示Ti的数据条目数目，每个矩形条的宽度为面积
Figure FDA0002814668920000013
对应的宽度。
2.如权利要求1所述的方法，其特征在于，还包括将统计指标数据绘制到地理区域边缘内侧的步骤：首先根据时序统计数据，计算需要进行分析和呈现的统计指标；然后根据时序统计数据的地理区域信息，将统计指标数据绘制到地理区域边缘内侧。
3.如权利要求1所述的方法，其特征在于，对边缘区域和中心区域使用相同或不同的颜色映射，来表示边缘区域或中心区域表征的数值由高到低。
4.如权利要求1所述的方法，其特征在于，对不同地理区域的中心区域的时序表征使用不同的方向表示时序数据从早到晚，包括：从左至右的方向表示时序数据从早到晚，从上到下的方向表示时序数据从早到晚。
5.如权利要求4所述的方法，其特征在于，对于地图上不同的地理区域，同时使用多种不同的方向来表示时序数据。
6.如权利要求5所述的方法，其特征在于，所有相邻的地理区域使用不同的方向来表示时序数据。
7.如权利要求5所述的方法，其特征在于，在确定每个地理区域的方向时，结合等面积呈现方法，以各个时序数据对应的等面积区域所占据的矩形宽度方差最小的方向作为最终方向。
8.如权利要求7所述的方法，其特征在于，通过如下公式计算各个角度的方差，并选择方差最小的值对应的方向作为所述最终方向：
Figure FDA0002814668920000021
其中，Varik表示地理区域i在角度k时的方差，widthik与heightik表示地理区域i在角度k时的包围盒宽度与高度，widthijk表示地理区域i的第j个时序序列在角度k时所占有的宽度，Num(Ti)表示地理区域Ti的时序数据条目数目。
CN201610753059.7A 2016-08-29 2016-08-29 一种支持时序特征探查的统计地图绘制方法 Active CN107749073B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610753059.7A CN107749073B (zh) 2016-08-29 2016-08-29 一种支持时序特征探查的统计地图绘制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610753059.7A CN107749073B (zh) 2016-08-29 2016-08-29 一种支持时序特征探查的统计地图绘制方法

Publications (2)

Publication Number Publication Date
CN107749073A CN107749073A (zh) 2018-03-02
CN107749073B true CN107749073B (zh) 2021-05-11

Family

ID=61254485

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610753059.7A Active CN107749073B (zh) 2016-08-29 2016-08-29 一种支持时序特征探查的统计地图绘制方法

Country Status (1)

Country Link
CN (1) CN107749073B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001045441A2 (en) * 1999-12-17 2001-06-21 Telefonaktiebolaget Lm Ericsson (Publ) System and method for enabling a user of a mobile radio terminal to influence their own radio quality
KR20140141761A (ko) * 2013-05-30 2014-12-11 전북대학교산학협력단 연결형 마이크로맵 도표를 이용한 통계정보 서비스 방법

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1959730A (zh) * 2006-10-11 2007-05-09 王云飞 区域招商投资系统中的资源交通状况评价方法
CN101593209B (zh) * 2008-12-30 2011-03-30 中国科学院地理科学与资源研究所 一种交互式在线统计地图生成与发布装置及方法
CN101741992B (zh) * 2009-12-22 2012-01-11 北京神州泰岳软件股份有限公司 基于gis的话务量分布呈现方法及系统
JP5881197B2 (ja) * 2012-08-14 2016-03-09 Ｎｅｃソリューションイノベータ株式会社 グラフ描画装置及びグラフ描画方法
CN104408137B (zh) * 2014-11-28 2018-11-13 武汉大学 一种网络统计地图可视化数据制备方法
CN105069020B (zh) * 2015-07-14 2018-09-21 国家信息中心 自然资源数据三维可视化方法和系统
CN105893588B (zh) * 2016-04-05 2019-04-30 东南大学 一种大尺度城市空间能耗数字地图制作及显示方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001045441A2 (en) * 1999-12-17 2001-06-21 Telefonaktiebolaget Lm Ericsson (Publ) System and method for enabling a user of a mobile radio terminal to influence their own radio quality
KR20140141761A (ko) * 2013-05-30 2014-12-11 전북대학교산학협력단 연결형 마이크로맵 도표를 이용한 통계정보 서비스 방법

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Community health assessment using self-organizing maps and geographic information systems;Heather G Basara 等;《International Journal of Health Geographics volume》;20081230;第1-8页 *

Also Published As

Publication number Publication date
CN107749073A (zh) 2018-03-02

Similar Documents

Publication Publication Date Title
US20120050285A1 (en) 3d building generalization for digital map applications
Liqiang et al. A spatial cognition-based urban building clustering approach and its applications
CN106354760A (zh) 基于变形统计地图的多视图时空数据可视化方法与应用
Ge et al. Vectorial boundary-based sub-pixel mapping method for remote-sensing imagery
CN107967702A (zh) 逆向渲染绘制的热力图可视化方法
CN102270236A (zh) 一种基于栅格化gis 空间关系判断方法及其系统
CN109902120B (zh) 一种面向栅格文本降水数据的可视化分析方法
CN109859109B (zh) 一种系列比例尺pdf地图无缝组织与显示方法
Shen et al. A new approach to simplifying polygonal and linear features using superpixel segmentation
US20110012917A1 (en) Dynamic generation of images to facilitate information visualization
CN108230433A (zh) 3d地图散点数据展示方法及系统、信息数据处理终端
CN112328929A (zh) 基于矢量瓦片地图发布技术的人流密度动态可视化方法
CN114239508A (zh) 表格还原方法和装置、存储介质及电子设备
CN104537024A (zh) 一种海量位置签到数据的可视化方法
CN108375806A (zh) 气象地图的获取方法和装置
CN113157988B (zh) 由ofd格式文件表示地理信息的方法和装置
Nguyen et al. Towards faithful graph visualizations
Yue et al. A function-based linear map symbol building and rendering method using shader language
CN107749073B (zh) 一种支持时序特征探查的统计地图绘制方法
She et al. An efficient method for rendering linear symbols on 3D terrain using a shader language
CN106873931A (zh) 栅格数据可视化装置、方法及计算机系统
CN112066998A (zh) 用于航线地图的渲染方法及系统
Kardos et al. The visualisation of uncertainty for spatially referenced census data using hierarchical tessellations
Du et al. Banded choropleth map
CN102074004B (zh) 空间实体遮挡类型的判定方法及装置

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