CN110348685B - 城市工业用地空间干扰度测评方法、系统、设备及介质 - Google Patents
城市工业用地空间干扰度测评方法、系统、设备及介质 Download PDFInfo
- Publication number
- CN110348685B CN110348685B CN201910494995.4A CN201910494995A CN110348685B CN 110348685 B CN110348685 B CN 110348685B CN 201910494995 A CN201910494995 A CN 201910494995A CN 110348685 B CN110348685 B CN 110348685B
- Authority
- CN
- China
- Prior art keywords
- land
- industrial
- residential
- interference
- various
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 28
- 238000012545 processing Methods 0.000 claims abstract description 24
- 238000013439 planning Methods 0.000 claims description 41
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000013278 delphi method Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 9
- 238000004140 cleaning Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 12
- 238000004458 analytical method Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 7
- 238000009826 distribution Methods 0.000 description 6
- 238000004590 computer program Methods 0.000 description 5
- 230000004927 fusion Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000013523 data management Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000002093 peripheral effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 239000008239 natural water Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000013077 scoring method Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种城市工业用地空间干扰度测评方法、系统、设备及介质，所述方法包括：获取各类用地数据；对各类用地数据进行处理，生成相应的专题地图；根据所述专题地图，对各类用地进行数据统计，得到居住用地和自然空间用地受工业用地干扰区域的面积，以及居住用地和自然空间用地受工业用地干扰边界的长度；确定工业用地对居住用地和自然空间用地的干扰影响权重；根据所述长度、所述面积和所述权重，计算得到各类用地的形态紧凑度、所述面积占自身用地面积的比例以及所述长度占自身用地周长的比例。本发明可以减轻规划师对于方案差异判读的工作量，提升协调工作的效率，同时在一定程度上为各类用地合理性的测评提供新的思路。
Description
技术领域
本发明涉及一种空间干扰度测评方法，尤其是一种城市工业用地空间干扰度测评方法、系统、设备及介质，属于属于城乡规划领域的技术应用。
背景技术
城市是人类经济和社会活动最重要的载体。随着城市的快速发展，城市用地规模的急剧扩张使得用地布局的成为城市规划中非常重要的环节。然而，我国目前对于城市用地布局合理性的研究尚待完善与深化。
在城市规划中，工业生产作为城市经济活动的重要组成部分，工业用地的布局也占据着重要的地位。城市工业用地可分为三类：一类工业用地、二类工业用地、三类工业用地。一类工业用地是对居住和公共设施等环境基本无干扰和污染的工业用地；二类工业用地是对居住和公共设施等环境具有一定干扰和污染的；三类工业用地则对周围的居住和公共设施具有严重的干扰和污染影响。因此，不同的工业用地布局方案对周围用地空间的干扰程度也不同，目前，还没有一种以定量的形式清晰的描述工业用地空间干扰度的方法。
目前，在现有的城市规划用地布局方案制定过程中，对用地布局合理性的判断仍处于感性思维，定性的判断，缺乏定量的清晰测度方法，对工业用地布局对其他用地的干扰影响也没有一套明确的测度其影响程度的指标和方法。
发明内容
本发明的第一个目的是为了解决上述现有技术的缺陷，提供了一种城市工业用地空间干扰度测评方法，该方法可以实现对城市工业用地空间干扰度的定量测评，优化工业用地布局的合理性判断，减轻规划师对于方案差异判读的工作量，提升协调工作的效率，同时在一定程度上为各类用地合理性的测评提供新的思路。
本发明的第二个目的在于提供一种城市工业用地空间干扰度测评系统。
本发明的第三个目的在于提供一种计算机设备。
本发明的第四个目的在于提供一种存储介质。
本发明的第一个目的可以通过采取如下技术方案达到：
一种城市工业用地空间干扰度测评方法，所述方法包括：
获取各类用地数据；其中，所述各类用地包括工业用地、居住用地和自然空间用地；
对各类用地数据进行处理，生成相应的专题地图；
根据所述专题地图，对各类用地进行数据统计，得到居住用地和自然空间用地受工业用地干扰区域的面积，以及居住用地和自然空间用地受工业用地干扰边界的长度；
确定工业用地对居住用地和自然空间用地的干扰影响权重；
根据所述长度、所述面积和所述权重，计算得到各类用地的形态紧凑度、所述面积占自身用地面积的比例以及所述长度占自身用地周长的比例。
进一步的，所述获取各类用地数据，具体包括：
获取用地规划CAD图纸；
对用地规划CAD图纸进行数据清理和矢量化处理，得到各类用地数据。
进一步的，所述对各类用地数据进行处理，生成各类用地的地图，具体包括：
将各类用地的线数据通过拓扑生成各类用地的面；
设定工业缓冲区的范围，制作工业用地对周围用地干扰区域缓冲区和干扰边界缓冲区，进行相交处理得到与居住用地、自然空间用地重合的区域和边界，并生成相应的专题地图。
进一步的，所述工业用地对居住用地和自然空间用地的干扰影响权重根据模糊德尔菲法确定，且所述工业用地对居住用地和自然空间用地的干扰影响权重包括一类工业用地对居住用地和自然空间用地的干扰影响权重、二类工业用地对居住用地和自然空间用地的干扰影响权重以及三类工业用地对居住用地和自然空间用地的干扰影响权重。
进一步的，所述各类用地的形态紧凑度，如下式：
其中，CF表示用地空间的分形维数；P表示各类用地周长，A表示各类用地面积。
进一步的，所述方法还包括：
将各类用地的形态紧凑度、所述面积占自身用地面积的比例以及所述长度占自身用地周长的比例生成柱状图。
本发明的第二个目的可以通过采取如下技术方案达到：
一种城市工业用地空间干扰度测评系统，所述系统包括：
获取模块，用于获取各类用地数据；其中，所述各类用地包括工业用地、居住用地和自然空间用地；
生成模块，用于对各类用地数据进行处理，生成相应的专题地图；
统计模块，用于根据所述专题地图，对各类用地进行数据统计，得到居住用地和自然空间用地受工业用地干扰区域的面积，以及居住用地和自然空间用地受工业用地干扰边界的长度；
确定模块，用于确定工业用地对居住用地和自然空间用地的干扰影响权重；
计算模块，用于根据所述长度、所述面积和所述权重，计算得到各类用地的形态紧凑度、所述面积占自身用地面积的比例以及所述长度占自身用地周长的比例。
进一步的，所述获取模块，具体包括：
图纸获取单元，用于获取用地规划CAD图纸；
矢量化处理单元，对用地规划CAD图纸进行数据清理和矢量化处理，得到各类用地数据。
进一步的，所述生成模块，具体包括：
面生成单元，用于将各类用地的线数据通过拓扑生成各类用地的面；
专题地图生成单元，用于设定工业缓冲区的范围，制作工业用地对周围用地干扰区域缓冲区和干扰边界缓冲区，进行相交处理得到与居住用地、自然空间用地重合的区域和边界，并生成相应的专题地图。
本发明的第三个目的可以通过采取如下技术方案达到：
一种计算机设备，包括处理器以及用于存储处理器可执行程序的存储器，其特征在于，所述处理器执行存储器存储的程序时，实现上述的城市工业用地空间干扰度测评方法。
本发明的第四个目的可以通过采取如下技术方案达到：
一种存储介质，存储有程序，所述程序被处理器执行时，实现上述的城市工业用地空间干扰度测评方法。
本发明相对于现有技术具有如下的有益效果：
本发明能够提高空间规划中对工业用地布局影响的判断效率，帮助选择更合适的工业用地空间区位与布局形式以减小对周围用地负面影响，应用于规划方案比选中可提供新的评选方法和技术支撑；此外，本发明不仅可以应用于工业用地的空间干扰测评，也可拓展应用于生态水系、绿地空间的外部性影响测评。
附图说明
图1为本发明实施例1的城市工业用地空间干扰度测评方法流程图。
图2a为本发明实施例1的方案一的工业用地、居住用地、自然空间用地布局图。
图2b为本发明实施例1的方案二的工业用地、居住用地、自然空间用地布局图。
图3a为本发明实施例1的方案一的工业用地对周围用地干扰区域缓冲区分布图。
图3b为本发明实施例1的方案二的工业用地对周围用地干扰区域缓冲区分布图。
图4a为本发明实施例1的方案一的工业用地对周围用地干扰边界缓冲区分布图。
图4b为本发明实施例1的方案二的工业用地对周围用地干扰边界缓冲区分布图。
图5a为本发明实施例1的方案一的受工业用地干扰的区域图。
图5b为本发明实施例1的方案二的受工业用地干扰的区域图。
图6a为本发明实施例1的方案一的受工业用地干扰的边界图。
图6b为本发明实施例1的方案二的受工业用地干扰的边界图。
图7为本发明实施例1的干扰程度语义变数与隶属度的曲线图。
图8为本发明实施例1的两个用地规划方案的各类用地的形态紧凑度对比曲线图。
图9为本发明实施例1的两个用地规划方案的一类、二类工业用地对其他用地干扰区域面积占比的曲线图。
图10为本发明实施例1的两个用地规划方案的一类、二类工业用地对其他用地干扰边界长度占比的曲线图。
图11为本发明实施例1的加权计算后两个用地规划方案受工业用地干扰区域占自身比例的曲线图。
图12为本发明实施例1的加权计算后两个用地规划方案受工业用地干扰边界占自身比例的曲线图。
图13为本发明实施例2的城市工业用地空间干扰度测评系统的结构框图。
图14为本发明实施例2的获取模块的结构框图。
图15为本发明实施例2的第一生成模块的结构框图。
图16为本发明实施例3的计算机设备的结构框图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。
实施例1：
如图1所示，本实施例提供了一种城市工业用地空间干扰度测评方法，该方法应用于两个用地规划方案(两个用地规划方案分别即为方案一和方案二)的比选，包括以下步骤：
S1、获取各类用地数据。
进一步地，该步骤S1具体包括：
S101、获取用地规划CAD图纸。
用地规划CAD图纸可以直接获取现成的用地规划CAD图纸(方案一和方案二的CAD图纸)，如果没有现成的用地规划CAD图纸，可以通过“湘源控规”软件手动描图绘制。
S102、对用地规划CAD图纸进行数据清理和矢量化处理，得到各类用地数据。
各类用地包括工业用地、居住用地和自然空间用地，其中工业用地包括一类工作用地M1、二类工作用地M2和三类工作用地M3，居住用地包括一类居住用地R1和二类居住用地R2，自然空间用地包括公园绿地(G1)、防护绿地(G2)、自然水域(E1)、农林用地(E2)。
S1021、通过“湘源控规”软件打开用地规划CAD图纸，对其进行数据清理，打开图层管理器，仅保留所需的工业用地、居住用地和自然空间用地的线图层和用地图层，关闭其他图层，另存为“用地CAD文件”。
S1022、使用“湘源控规“打开“用地CAD文件”提取CAD线要素，输入OL命令—孤立图层，点选居住用地(R1)图层线框，单击空格确定，即单独显示居住用地(R1)的线框，然后另存为“居住R1.dwg”文件；同样的方法得到“工业M1.dwg”、“工业M2.dwg”、“工业M3.dwg”、“居住R1.dwg”、“居住R2.dwg”、“公园绿地G1.dwg”、“防护绿地G2.dwg”、“自然水域E1.dwg”、“农林用地E2.dwg”。
S1023、对各类用地的CAD文件进行矢量化处理，将CAD文件转为ArcGIS可用的shapfile文件，该shapfile文件即为各类用地数据。
1)用ArcGIS软件，新建个人地理数据库：打开ArcCatalog，点击链接到文件夹，选择文件存储的位置。单击“内容”列表空白处，选择，新建—个人地理数据库，方案一和方案二的个人地理数据库可以分别命名为“方案一的数据库”和“方案二的数据库”。
2)新建地图文档，设置默认地理数据库：打开ArcMap，点击菜单栏，文件—地图文档属性—默认地理数据库，选择个人地理数据库，勾选“存储数据源的相对路径”。
3)设置图层坐标系：点击菜单栏——内容列表，显示内容列表，右键点击图层，选择属性—坐标系—地理坐标系—World—WGS1984，将其设置为WGS1984坐标系。
4)链接到工作目录文件夹，添加各类用地的CAD文件。
5)在ArcMap中，将CAD文件中的Polyline右键—数据—导出数据，分别导出为shapfile文件，并分别命名为M1、M2、R1、R2、G1、G2、E1、E2，并添加到图层中。
S1024、为矢量数据定义投影坐标系，将文件单位由经纬度转化为以米为单位。
打开目录列表中的系统工具箱—数据管理(Data Management Tools)—投影和变换，输入用地的shapfile文件，输出为投影坐标系—UTM—WGS1984，将各个用地分别定义为UTM下WGS1984的投影坐标系，即转换为了以米为单位的地图文档。
S2、对各类用地数据进行处理，生成相应的专题地图。
进一步地，该步骤S2具体包括：
S201、将各类用地的线数据通过拓扑生成各类用地的面。
S2011、左键单击编辑器—开始编辑，打开自定义—工具条—高级编辑，选择某一类用地的图层，如M1，选择M1中闭合的线段—点击拓扑工具条中的构造面，生成所有面后完成拓扑，得到M1面即工业用地的区域范围；同样的方法得到M2面、R1面、R2面、G1面、G2面、E1面、E2面。
S2012、点击ArcMap页面左下角的布局视图，将数据视图切换为布局视图，点击工具栏的放大与缩小工具调整合适的图像显示范围与大小，打开图层管理器，打开拓扑生成的各类工业用地的面，调整图例颜色，点击文件-导出地图，图2a示出了方案一的工业用地、居住用地和自然空间用地布局，图2b示出了方案二的工业用地、居住用地和自然空间用地布局。
S202、设定工业缓冲区的范围，制作工业用地对周围用地干扰区域缓冲区和干扰边界缓冲区，进行相交处理得到与居住用地、自然空间用地重合的区域和边界，并生成相应的专题地图。
S2021、制作工业用地对周围用地干扰区域缓冲区和干扰边界缓冲区，导出生成相应的专题地图。
工业用地对周边用地的干扰区域以缓冲区的形式表现，不同类型的工业用地对周边用地的影响范围不同，影响范围与周边用地重合的区域即为受工业用地干扰的区域。
1)在居住区中步行可达的范围为500m，以500m为二类工业(M2)用地影响的缓冲区，以300m为一类工业用地(M1)的影响范围：在ArcMap中，使用系统工具箱—分析工具(Analysis Tools)—邻域分析—多环缓冲区，输入一类工业用地(M1面)，设置缓冲距离300m，输出为“M1面buf300”即一类工业用地300m的缓冲范围，同理得到“M2面buf500”，即二类工业用地500m的缓冲范围，导出生成相应的专题地图，该专题地图即为工业用地对周围用地干扰区域缓冲区分布图，方案一的工业用地对周围用地干扰区域缓冲区分布图如图3a所示，方案二的工业用地对周围用地干扰区域缓冲区分布图如图3b所示。
2)在城市总体规划中，城市道路(包括主干道、次干道、支路)基本不超过60m，故选择60m作为工业对周围其他用地影响的边界范围。用S2021中的方法得到一类工业用地和二类工业用地60m的影响范围，“M1面buf60”、“M2面buf60”，导出生成相应的专题地图，该专题地图即为工业用地对周围用地干扰边界缓冲区分布图，方案一的工业用地对周围用地干扰边界缓冲区分布图如图4a所示，方案二的工业用地对周围用地干扰边界缓冲区分布图如图4b所示。
S2022、获取居住用地和自然空间用地受工业用地干扰的区域，导出生成相应的专题地图。
1)分别计算工业缓冲区与居住用地和自然空间用地的重叠区域，即“M1面buf300”、“M2面buf500”与R2面、G1面、G2面、E1面、E2面重合的区域。打开系统工具箱—分析工具(Analysis Tools)—提取分析—裁剪，分别输入工业缓冲区和居住用地及自然空间用地作为输入要素，分别输入M1面、M2面作为裁剪要素，点击确定，输出得到M1R2面、M1G1面、M1G2面、M1E1面、M1E2面、M2R2面、M2G1面、M2G2面、M2E1面、M2E2面，即为各类用地分别受M1、M2干扰的区域。
2)点击ArcMap页面左下角的布局视图，将数据视图切换为布局视图，点击工具栏的放大与缩小工具调整合适的图像显示范围与大小，打开图层管理器，打开受工业干扰区域图层，调整图例颜色，点击文件-导出地图，导出相应的专题地图，该专题地图即为受工业用地干扰的区域图，方案一的受工业用地干扰的区域图如图5a所示，方案二的受工业用地干扰的区域图如图5b所示。
S2023、获取居住用地和自然空间用地受工业用地干扰的边界，导出生成相应的专题地图。
1)工业用地与其他用地相邻会对周围用地的边界产生干扰，对不同地块所产生的干扰长度不同。以60m为工业用地对周围用地边界干扰的范围所得到的“M1面buf60”，“M2面buf60”，打开系统工具箱—分析工具(Analysis Tools)—叠加分析—相交，输入工业缓冲区和居住用地及自然空间用地的线要素，即R2、G1、G2、E1、E2，点击确定，输出得到M1R2线、M1G1线、M1G2线、M1E1线、M1E2线、M2R2线、M2G1线、M2G2线、M2E1线、M2E2线，即为各类用地分别受M1、M2干扰的边界。
2)根据要显示的地图内容有选择地进行图层的开关，叠加不同图层后导出生成相应的专题地图，该专题地图即为受工业用地干扰的边界图，方案一的受工业用地干扰的边界图如图6a所示，方案二的受工业用地干扰的边界图如图6b所示。
S3、根据所述专题地图，对各类用地进行数据统计，得到居住用地和自然空间用地受工业用地干扰区域的面积，以及居住用地和自然空间用地受工业用地干扰边界的长度。
在ArcGIS中，打开各专题地图的shapefile文件，将各类用地融合得到统一的要素，在属性表中对各类用地进行数据统计，得到居住用地和自然空间用地受工业用地干扰区域的面积，以及居住用地和自然空间用地受工业用地干扰边界的长度，数据汇总粘贴到Excel。
进一步地，该步骤S3具体包括：
S301、在ArcMap中打开各专题地图的shapefile文件，通过融合工具将各类用地融合，并在其属性表中，对工业用地、居住用地和自然空间用地的长度和面积进行统计。
S3011、将各专题地图的shapefile文件分别进行融合得到统一的要素，打开系统工具箱—数据管理(Data Management Tools)—制图综合—融合，将要素M1、M2、R2、G1、G2、E1、E2分别作为输入要素，选择融合字段，输出M1_Dissolve、M2_Dissolve、R2_Dissolve、G1_Dissolve、G2_Dissolve、E1_Dissolve、E2_Dissolve。
S3012、右键单击融合得到的图层，打开属性表，属性表中SHAPE_Length即为要素的长度，SHAPE_Area即为要素面积，右键单击复制。新建一个Excel表格，将复制得到的长度(P)和面积(A)数据粘贴到表格中。用此方法依次得到方案中M1、M2、R2、G1、G2、E1、E2各用地的周长和面积。
S302、对居住用地和自然空间用地受工业用地干扰区域的面积进行统计。
S3021、打开居住用地和自然空间用地受工业用地干扰的区域，即M1R2面、M1G1面、M1G2面、M1E1面、M1E2面、M2R2面、M2G1面、M2G2面、M2E1面、M2E2面。右键单击图层，打开属性表，右键单击SHAPE_Area，点击统计，将显示的列表下的“总和”(Ai)数据复制粘贴到Excel表中。
S3022、依次对每个图层进行同样的操作，即可分别得到受干扰的面积统计数据。
S303、对居住用地和自然空间用地受工业用地干扰边界的长度进行统计。
打开居住用地和自然空间用地受工业用地干扰的边界，即M1R2线、M1G1线、M1G2线、M1E1线、M1E2线、M2R2线、M2G1线、M2G2线、M2E1线、M2E2线。右键单击图层，打开属性表，右键单击SHAPE_Length，点击统计，将显示的列表下的“总和”(Pi)数据复制粘贴到Excel表中。
S4、确定工业用地对居住用地和自然空间用地的干扰影响权重。
工业用地对居住用地和自然空间用地的干扰影响权重包括一类工业用地对居住用地和自然空间用地的干扰影响权重、二类工业用地对居住用地和自然空间用地的干扰影响权重以及三类工业用地对居住用地和自然空间用地的干扰影响权重。
在对一类、二类、三类工业用地对其他用地干扰影响程度判断上，为追求干扰影响权重设定的客观性，因此邀请相关领域专家对三种工业类型进行打分，综合应用了模糊德尔菲法(Fuzzy Delphi)。
进一步地，该步骤S4具体包括：
S401、共选取专家36名，专家库的组成如下：星桥国际私人有限公司22位专家(有新加坡纬壹科技城、苏州工业园、天津生态城的开发管理经验)；雅思柏RSP设计事务所和中国城市规划设计研究院10位专家(中新广州知识城项目的主创设计方和咨询团队，在规划编制过程中，对智慧新区有深入了解)；高等教育学府中对用地空间形态有丰硕研究成果的4位专家(作为此领域的专家，其建议对本实施例的研究具有重要意义)。
以语意方式收集专家表达评价的信息，即以“干扰程度”语义变数为界面，将语意尺度均匀划分为三角模数，以此量化计算隶属度，如图7所示；然后将语意表达的词汇转化为模糊数(Fuzzy number)，再将模糊集合转化为明确数(Crisp score)。
S402、根据模糊德尔菲法，确定工业用地对居住用地和自然空间用地的干扰影响权重。
模糊德尔菲法的具体步骤如下：
S4021、评估因子的平均值。
先以式(1)计算各评估因素的综合三角模糊数的平均数，获得单一因素数值，并作为模糊化分析之用。
Vm＝(1/N)[nm1(0,0,1/(L-1))+nm2(0,1/(L-1)),2/(L-1)+···+nm2((K-2)/(L-1),(K-1)/(L-1)),K/(L-1))+···+nmL((K-2)/(L-1),1,1)] (1)
其中，Vm为第m个评估因素的平均值；nml为第m个评估因素的第L个语意尺度被选择的次数；N为nmL的加总；L为语意尺度的划分数目；k为第k个语意尺度。
S4022、模糊数的反模糊化分析。
反模糊化分析的目的在于将三角形模糊数转化为明确的数值，以进行因素之间的优劣势排序比较与数值相关分析，操作方法与步骤如下：定义语义尺度集合的最大集和最小集。
其中，μmax与μmin分别与Vs的右界与左界产生交集，已知Vs＝(a,b,c)代表三个点坐标(a，0)、(b，1)、(c，0)，则其三角模糊数所建构的模糊线性函数:
y＝(x-a)/(b-a)
y＝(x-c)/(b-c)
由μmax(x)最大隶属函数与Vs模糊函数所形成的函数关系如式(4)所示：
μmax(m)＝supmin[μmax(x),μVs(x)] (4)
其中，右界值的两个焦点坐标为：(a/(1+a-b),a/(1+a-b))与(c/(1+c-b),c/(1+c-b))。
取y坐标值(隶属值)较大者的y值代表μR(x)。
同理，由μmin(x)最小集隶属函数与Vs模糊函数所形成的函数关系如公式(5)所示：
μmin(m)＝supmin[μmin(x),μVs(x)] (5)
其中，左界值的两个焦点坐标为：(b/(1+b-a),(1-a)/(1+b-a))与(b/(1+b-c),(1-c),(1+b-c))
取y坐标值(隶属值)较大者的y值代表μL(x)。
最后，以反模糊化运算模糊数的反模糊得：
μT(m)＝0.5[μR(m)+1-μL(m)] (6)
由此可知，三角形模糊越趋向右边界(满意程度越高)，其反模糊比值μT(m)越大。因此，利用此方法建立的模糊语意尺度，可以避免传统的评分法进行模糊评估，并使受访者可以以日常口语化方式表达评估值。
根据模糊德尔菲法，得到各类工业用地对居住用地和自然空间用地的干扰影响权重，如下表所示。
表1各类工业用地对居住用地和自然空间用地的干扰影响权重
S5、根据所述长度、所述面积和所述权重，计算得到各类用地的形态紧凑度、所述面积占自身用地面积的比例以及所述长度占自身用地周长的比例。
对汇总的Excel数据进行处理，根据各类工业用地对居住用地和自然空间用地的干扰影响权重进行加权计算，计算各类用地的形态紧凑度、居住用地和自然空间用地受工业用地干扰区域的面积占自身用地面积的比例，以及居住用地和自然空间用地受工业用地干扰边界的长度占自身用地周长的比例。
进一步地，该步骤S5具体包括：
S501、计算各类用地的形态紧凑度。
形态紧凑度反映用地边界形状的紧凑程度；用地形态越紧凑，用地布局越紧凑，从而导致形态紧凑度升高，形态紧凑度的测定如下式：
其中，CF表示用地空间的分形维数；P表示各类用地周长，A表示各类用地面积。
S502、计算居住用地和自然空间用地受工业用地干扰区域的面积占自身用地面积的比例，如下式：
其中，A表示各类用地自身的面积，Ai表示各类用地分别受工业用地干扰的面积。
S503、计算居住用地和自然空间用地受工业用地干扰边界的长度占自身用地周长的比例，如下式：
其中，P表示各类用地自身的周长，Pi表示各类用地分别受工业用地干扰的长度。
S504、加权计算，根据模糊德尔菲法得到的各类工业用地对居住用地和自然空间用地的干扰影响权重(一类工业为0.1、二类工业为0.3、三类工业为0.6)，按照比重统计一类工业和二类工业的影响，分别汇总计算工业用地总体对周围的影响。
为了方便对两个用地规划方案进行对比，本实施例的城市工业用地空间干扰度测评方法，还可包括：
S6、将各类用地的形态紧凑度、所述面积占自身用地面积的比例以及所述长度占自身用地周长的比例生成柱状图。
对各类用地的形态紧凑度、受工业用地干扰区域的面积占自身用地面积的比例以及受工业用地干扰边界的长度占自身用地周长的比例这三个指标的计算结果表格用Excel生成柱状图进行可视化，具体将两个用地规划方案的三个指标的计算结果统计结果放到同一个Excel表格下进行对比，插入柱状图，添加数据标签得到可视化结果，图8为两个用地规划方案的各类用地的形态紧凑度对比曲线图，图9为两个用地规划方案的一类、二类工业用地对其他用地干扰区域所占比例的曲线图，图10为两个用地规划方案的一类、二类工业用地对其他用地干扰边界所占比例的曲线图；加权计算后，两个用地规划方案受工业用地干扰区域占自身比例的曲线图如图11所示，两个用地规划方案受工业用地干扰边界占自身比例的曲线图如图12所示。
综合对比，可以看出方案一的工业形态紧凑度较低，一类和二类工业的形态紧凑度均低于方案二，工业对周围用地产生的干扰影响也最弱，因此从自然保护和良好的居住环境的角度出发，选择方案一的用地布局模式较为合适。
因此，本实施例的城市工业用地空间干扰度测评方法可以总结为：工业用地对周边的消极影响与其自身的形态紧凑度有关，形态紧凑度越高，空间布局越集中，对周围的用地干扰影响越小，可推广至其他类用地的消极或积极影响判定。
应当注意，尽管以特定顺序描述了上述实施例的方法操作，但是这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。相反，描绘的步骤可以改变执行顺序。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。
实施例2：
如图13所示，本实施例提供了一种城市工业用地空间干扰度测评系统，该系统包括获取模块1301、第一生成模块1302、统计模块1303、确定模块1304、计算模块1305和第二生成模块1306，各个模块的具体功能如下：
所述获取模块1301，用于获取各类用地数据；其中，所述各类用地包括工业用地、居住用地和自然空间用地。
所述生成模块1302，用于对各类用地数据进行处理，生成相应的专题地图。
所述统计模块1303，用于根据所述专题地图，对各类用地进行数据统计，得到居住用地和自然空间用地受工业用地干扰区域的面积，以及居住用地和自然空间用地受工业用地干扰边界的长度。
所述确定模块1304，用于确定工业用地对居住用地和自然空间用地的干扰影响权重。
所述计算模块1305，用于根据所述长度、所述面积和所述权重，计算得到各类用地的形态紧凑度、所述面积占自身用地面积的比例以及所述长度占自身用地周长的比例。
所述第二生成模块1306，用于将各类用地的形态紧凑度、所述面积占自身用地面积的比例以及所述长度占自身用地周长的比例生成柱状图。
如图14所示，所述获取模块1301，具体包括：
图纸获取单元13011，用于获取用地规划CAD图纸。
矢量化处理单元13012，对用地规划CAD图纸进行数据清理和矢量化处理，得到各类用地数据。
如图15所示，所述第一生成模块1302，具体包括：
拓扑生成单元13021，用于将各类用地的线数据通过拓扑生成各类用地的面。
制作处理单元13022，用于设定工业缓冲区的范围，制作工业用地对周围用地干扰区域缓冲区和干扰边界缓冲区，进行相交处理得到与居住用地、自然空间用地重合的区域和边界，并生成相应的专题地图。
需要说明的是，本实施例提供的系统仅以上述各功能模块的划分进行举例说明，在实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。
实施例3：
本实施例提供了一种计算机设备，该计算机设备可以是计算机，如图16所示，其包括通过系统总线1601连接的处理器1602、存储器、输入装置1603、显示器1604和网络接口1605，该处理器用于提供计算和控制能力，该存储器包括非易失性存储介质1606和内存储器1607，该非易失性存储介质1606存储有操作系统、计算机程序和数据库，该内存储器1607为非易失性存储介质中的操作系统和计算机程序的运行提供环境，处理器1602执行存储器存储的计算机程序时，实现上述实施例1的城市工业用地空间干扰度测评方法，如下：
获取各类用地数据；其中，所述各类用地包括工业用地、居住用地和自然空间用地；
对各类用地数据进行处理，生成相应的专题地图；
根据所述专题地图，对各类用地进行数据统计，得到居住用地和自然空间用地受工业用地干扰区域的面积，以及居住用地和自然空间用地受工业用地干扰边界的长度；
确定工业用地对居住用地和自然空间用地的干扰影响权重；
根据所述长度、所述面积和所述权重，计算得到各类用地的形态紧凑度、所述面积占自身用地面积的比例以及所述长度占自身用地周长的比例。
实施例4：
本实施例提供了一种存储介质，该存储介质为计算机可读存储介质，其存储有计算机程序，计算机程序被处理器执行时，实现上述实施例1的城市工业用地空间干扰度测评方法，如下：
获取各类用地数据；其中，所述各类用地包括工业用地、居住用地和自然空间用地；
对各类用地数据进行处理，生成相应的专题地图；
根据所述专题地图，对各类用地进行数据统计，得到居住用地和自然空间用地受工业用地干扰区域的面积，以及居住用地和自然空间用地受工业用地干扰边界的长度；
确定工业用地对居住用地和自然空间用地的干扰影响权重；
根据所述长度、所述面积和所述权重，计算得到各类用地的形态紧凑度、所述面积占自身用地面积的比例以及所述长度占自身用地周长的比例。
本实施例中的存储介质可以是磁盘、光盘、计算机存储器、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、U盘、移动硬盘等介质。
综上所述，本发明能够提高空间规划中对工业用地布局影响的判断效率，帮助选择更合适的工业用地空间区位与布局形式以减小对周围用地负面影响，应用于规划方案比选中可提供新的评选方法和技术支撑；此外，本发明不仅可以应用于工业用地的空间干扰测评，也可拓展应用于生态水系、绿地空间的外部性影响测评。
以上所述，仅为本发明专利较佳的实施例，但本发明专利的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明专利所公开的范围内，根据本发明专利的技术方案及其发明专利构思加以等同替换或改变，都属于本发明专利的保护范围。
Claims (8)
1.一种城市工业用地空间干扰度测评方法，其特征在于，所述方法包括：
获取各类用地数据；其中，所述各类用地包括工业用地、居住用地和自然空间用地；
对各类用地数据进行处理，生成相应的专题地图；
根据所述专题地图，对各类用地进行数据统计，得到居住用地和自然空间用地受工业用地干扰区域的面积，以及居住用地和自然空间用地受工业用地干扰边界的长度；
确定工业用地对居住用地和自然空间用地的干扰影响权重；其中，所述工业用地对居住用地和自然空间用地的干扰影响权重根据模糊德尔菲法确定，且所述工业用地对居住用地和自然空间用地的干扰影响权重包括一类工业用地对居住用地和自然空间用地的干扰影响权重、二类工业用地对居住用地和自然空间用地的干扰影响权重以及三类工业用地对居住用地和自然空间用地的干扰影响权重；
计算各类用地的形态紧凑度；计算居住用地和自然空间用地受工业用地干扰区域的面积占自身用地面积的比例；计算居住用地和自然空间用地受工业用地干扰边界的长度占自身用地周长的比例；加权计算，根据模糊德尔菲法得到的各类工业用地对居住用地和自然空间用地的干扰影响权重，按照比重统计一类工业和二类工业的影响，分别汇总计算工业用地总体对周围的影响；
所述形态紧凑度的计算如下式：
其中，CF表示用地空间的分形维数；P表示各类用地周长，A表示各类用地面积；
所述面积占自身用地面积的比例，如下式：
其中，A表示各类用地自身的面积，Ai表示各类用地分别受工业用地干扰的面积；
所述长度占自身用地周长的比例，如下式：
其中，P表示各类用地自身的周长，Pi表示各类用地分别受工业用地干扰的长度。
2.根据权利要求1所述的城市工业用地空间干扰度测评方法，其特征在于，所述获取各类用地数据，具体包括：
获取用地规划CAD图纸；
对用地规划CAD图纸进行数据清理和矢量化处理，得到各类用地数据。
3.根据权利要求1所述的城市工业用地空间干扰度测评方法，其特征在于，所述对各类用地数据进行处理，生成各类用地的地图，具体包括：
将各类用地的线数据通过拓扑生成各类用地的面；
设定工业缓冲区的范围，制作工业用地对周围用地干扰区域缓冲区和干扰边界缓冲区，进行相交处理得到与居住用地、自然空间用地重合的区域和边界，并生成相应的专题地图。
4.根据权利要求1-3任一项所述的城市工业用地空间干扰度测评方法，其特征在于，所述方法还包括：
将各类用地的形态紧凑度、所述面积占自身用地面积的比例以及所述长度占自身用地周长的比例生成柱状图。
5.一种城市工业用地空间干扰度测评系统，其特征在于，所述系统包括：
获取模块，用于获取各类用地数据；其中，所述各类用地包括工业用地、居住用地和自然空间用地；
生成模块，用于对各类用地数据进行处理，生成相应的专题地图；
统计模块，用于根据所述专题地图，对各类用地进行数据统计，得到居住用地和自然空间用地受工业用地干扰区域的面积，以及居住用地和自然空间用地受工业用地干扰边界的长度；
确定模块，用于确定工业用地对居住用地和自然空间用地的干扰影响权重；其中，所述工业用地对居住用地和自然空间用地的干扰影响权重根据模糊德尔菲法确定，且所述工业用地对居住用地和自然空间用地的干扰影响权重包括一类工业用地对居住用地和自然空间用地的干扰影响权重、二类工业用地对居住用地和自然空间用地的干扰影响权重以及三类工业用地对居住用地和自然空间用地的干扰影响权重；
计算模块，用于计算各类用地的形态紧凑度；计算居住用地和自然空间用地受工业用地干扰区域的面积占自身用地面积的比例；计算居住用地和自然空间用地受工业用地干扰边界的长度占自身用地周长的比例；加权计算，根据模糊德尔菲法得到的各类工业用地对居住用地和自然空间用地的干扰影响权重，按照比重统计一类工业和二类工业的影响，分别汇总计算工业用地总体对周围的影响；
所述形态紧凑度的计算如下式：
其中，CF表示用地空间的分形维数；P表示各类用地周长，A表示各类用地面积；
所述面积占自身用地面积的比例的计算，如下式：
其中，A表示各类用地自身的面积，Ai表示各类用地分别受工业用地干扰的面积；
所述长度占自身用地周长的比例的计算，如下式：
其中，P表示各类用地自身的周长，Pi表示各类用地分别受工业用地干扰的长度。
6.根据权利要求5所述的城市工业用地空间干扰度测评系统，其特征在于，所述生成模块，具体包括：
面生成单元，用于将各类用地的线数据通过拓扑生成各类用地的面；
专题地图生成单元，用于设定工业缓冲区的范围，制作工业用地对周围用地干扰区域缓冲区和干扰边界缓冲区，进行相交处理得到与居住用地、自然空间用地重合的区域和边界，并生成相应的专题地图。
7.一种计算机设备，包括处理器以及用于存储处理器可执行程序的存储器，其特征在于，所述处理器执行存储器存储的程序时，实现权利要求1-4任一项所述的城市工业用地空间干扰度测评方法。
8.一种存储介质，存储有程序，其特征在于，所述程序被处理器执行时，实现权利要求1-4任一项所述的城市工业用地空间干扰度测评方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910494995.4A CN110348685B (zh) | 2019-06-10 | 2019-06-10 | 城市工业用地空间干扰度测评方法、系统、设备及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910494995.4A CN110348685B (zh) | 2019-06-10 | 2019-06-10 | 城市工业用地空间干扰度测评方法、系统、设备及介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110348685A CN110348685A (zh) | 2019-10-18 |
| CN110348685B true CN110348685B (zh) | 2021-12-21 |
Family
ID=68181668
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910494995.4A Active CN110348685B (zh) | 2019-06-10 | 2019-06-10 | 城市工业用地空间干扰度测评方法、系统、设备及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110348685B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111611335B (zh) * | 2020-05-09 | 2023-05-26 | 浙江数治空间规划设计有限公司 | 一种国土空间适用性评价方法、系统及存储介质 |
| CN115168972B (zh) * | 2022-08-04 | 2023-04-11 | 上海市农业科学院 | 基于有限空间识别的面源污染生态拦截工程优化方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101510273A (zh) * | 2008-08-28 | 2009-08-19 | 南京大学 | 快速统计城市规划用地面积的方法和计算机程序 |
| CN102147890A (zh) * | 2011-04-11 | 2011-08-10 | 复旦大学 | 城市土地使用与交通一体化规划决策支持方法及系统 |
| CN103714244A (zh) * | 2013-12-08 | 2014-04-09 | 同济大学 | 一种基于用地兼容性的土地利用混杂度评价方法 |
| CN105930865A (zh) * | 2016-04-18 | 2016-09-07 | 重庆市勘测院 | 城市建设用地分类提取及评估方法 |
| CN109598056A (zh) * | 2018-11-30 | 2019-04-09 | 华南理工大学 | 城市建设用地形态紧凑度的测度方法、系统及存储介质 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9208612B2 (en) * | 2010-02-12 | 2015-12-08 | The University Of North Carolina At Chapel Hill | Systems and methods that generate height map models for efficient three dimensional reconstruction from depth information |
-
2019
- 2019-06-10 CN CN201910494995.4A patent/CN110348685B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101510273A (zh) * | 2008-08-28 | 2009-08-19 | 南京大学 | 快速统计城市规划用地面积的方法和计算机程序 |
| CN102147890A (zh) * | 2011-04-11 | 2011-08-10 | 复旦大学 | 城市土地使用与交通一体化规划决策支持方法及系统 |
| CN103714244A (zh) * | 2013-12-08 | 2014-04-09 | 同济大学 | 一种基于用地兼容性的土地利用混杂度评价方法 |
| CN105930865A (zh) * | 2016-04-18 | 2016-09-07 | 重庆市勘测院 | 城市建设用地分类提取及评估方法 |
| CN109598056A (zh) * | 2018-11-30 | 2019-04-09 | 华南理工大学 | 城市建设用地形态紧凑度的测度方法、系统及存储介质 |
Non-Patent Citations (1)
| Title |
|---|
| Evaluation of urban land-use structures with a view to sustainable development;Nguyen Xuan Thinh et al.;《Environment Impact Assessment Review》;20021031;第1-21页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110348685A (zh) | 2019-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Torkayesh et al. | Landfill location selection for healthcare waste of urban areas using hybrid BWM-grey MARCOS model based on GIS | |
| Malczewski | Integrating multicriteria analysis and geographic information systems: the ordered weighted averaging (OWA) approach | |
| Paudel et al. | Assessing landscape changes and dynamics using patch analysis and GIS modeling | |
| Bornmann et al. | Mapping excellence in the geography of science: An approach based on Scopus data | |
| Malczewski et al. | Exploring multicriteria decision strategies in GIS with linguistic quantifiers: A case study of residential quality evaluation | |
| Carr et al. | Smart land-use analysis: the LUCIS model land-use conflict identification strategy | |
| CN110569580B (zh) | 城市街道空间活力模拟方法、系统、计算机设备及介质 | |
| Dujmović et al. | LSP suitability maps | |
| CN107038263B (zh) | 一种基于数据图谱、信息图谱和知识图谱的搜索优化方法 | |
| CN110348685B (zh) | 城市工业用地空间干扰度测评方法、系统、设备及介质 | |
| Pert et al. | Participatory development of a new interactive tool for capturing social and ecological dynamism in conservation prioritization | |
| Tsang et al. | An approach for aggregating upstream catchment information to support research and management of fluvial systems across large landscapes | |
| Richter et al. | Meshing Centroidal Voronoi Tessellation with spatial statistics to optimize waste management regions | |
| Argent et al. | A modelling shell for participatory assessment and management of natural resources | |
| Alexander | Searching for a robust operation of Lake Mead | |
| Li et al. | Urban Morphology and Historical Urban Landscape Conservation and Management | |
| Wang | Environmental Informatics for Environmental Planning and Management. | |
| Torres et al. | A prioritization tool for SUDS planning in large cities by coupling an urban drainage model with mixed integer linear programming | |
| Mokallaf Sarband et al. | Adaptation of a compromise programming approach for evaluating the localized impacts of water allocation | |
| Yan et al. | Multilevel Robustness for 2D Vector Field Feature Tracking, Selection and Comparison | |
| Qiao et al. | Research on environmental planning method based on neural network and artificial intelligence technology | |
| Qin | Local ideal point method for GIS-based multicriteria analysis: A case study in London, Ontario | |
| Geng et al. | Toward digital soil mapping in Canada: Existing soil survey data and related expert knowledge | |
| Germanaitė et al. | Case study of spatial pattern description, identification and application methodology | |
| Cai et al. | A schema of ecological environment sensitivity evaluation based on GIS |
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 |