CN113076512B - 土地供应成熟度评价优化方法 - Google Patents

土地供应成熟度评价优化方法 Download PDF

Info

Publication number
CN113076512B
CN113076512B CN202110325947.XA CN202110325947A CN113076512B CN 113076512 B CN113076512 B CN 113076512B CN 202110325947 A CN202110325947 A CN 202110325947A CN 113076512 B CN113076512 B CN 113076512B
Authority
CN
China
Prior art keywords
land
representing
marked
maturity
factor
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
CN202110325947.XA
Other languages
English (en)
Other versions
CN113076512A (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.)
Chongqing Planning And Natural Resources Information Center
Original Assignee
Chongqing Planning And Natural Resources Information Center
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 Chongqing Planning And Natural Resources Information Center filed Critical Chongqing Planning And Natural Resources Information Center
Priority to CN202110325947.XA priority Critical patent/CN113076512B/zh
Publication of CN113076512A publication Critical patent/CN113076512A/zh
Application granted granted Critical
Publication of CN113076512B publication Critical patent/CN113076512B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • 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/08—Construction
    • 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/16—Real estate
    • G06Q50/165—Land development

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Mathematical Physics (AREA)
  • General Business, Economics & Management (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Pure & Applied Mathematics (AREA)
  • Operations Research (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Development Economics (AREA)
  • Primary Health Care (AREA)
  • Computational Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提出了一种土地供应成熟度评价优化方法，包括以下步骤：S1，获取待供应土地相影响的数据因素，得到其土地影响因素数据；S2，对步骤S1中得到的土地影响因素数据进行预处理，得到其土地预处理影响因素数据；S3，对步骤S3中得到的土地预处理影响因素数据形成影响因素矩阵；S4，若计算得到的成熟度系数大于或者等于预设成熟度系数，则对其待供应土地在地图上进行标红处理；若计算得到的成熟度系数小于预设成熟度系数，则对其待供应土地进行标灰处理。本发明能够在地图上以不同色彩度标记其土地成熟度。

Description

土地供应成熟度评价优化方法
技术领域
本发明涉及一种优化处理技术领域，特别是涉及一种土地供应成熟度评价优化方法。
背景技术
现行建设项目生成机制缺乏规划实施时序的统筹安排，造成城市建设不合理、不紧凑，土地利用效率和产出效率不高，土地价值未得到充分发挥。因此要强化规划引领城市发展，统筹土地供应管理，优化建设项目时序安排。调控城市建设的空间布局、总量规模、结构比例和时序节奏，促进城市精明增长、紧凑发展，提高土地利用效率和产出效率，实现土地作为生产要素的整体价值最大化；推动城市开发建设“规划一片、实施一片、完善一片、成熟一片”，推进产城景融合发展，为城市提升“赋能”，推动高质量发展、创造高品质生活。
发明内容
本发明旨在至少解决现有技术中存在的技术问题，特别创新地提出了一种土地供应成熟度评价优化方法。
为了实现本发明的上述目的，本发明提供了一种土地供应成熟度评价优化方法，包括以下步骤：
S1，获取待供应土地相影响的数据因素，得到其土地影响因素数据；
S2，对步骤S1中得到的土地影响因素数据进行预处理，得到其土地预处理影响因素数据；
S3，对步骤S3中得到的土地预处理影响因素数据形成影响因素矩阵；
S4，根据其步骤S3中得到的影响因素矩阵与其所对应的权重系数，得到其成熟度系数；判断其成熟度系数是否大于或者等于预设成熟度系数：
若计算得到的成熟度系数大于或者等于预设成熟度系数，则对其待供应土地在地图上进行标红处理；
若计算得到的成熟度系数小于预设成熟度系数，则对其待供应土地进行标灰处理。
在本发明的一种优选实施方式中，在步骤S1中，与待供应土地相影响的数据因素包括规划人口、底线管控、用地实施情况、土地供应条件、配套设施完善度、交通便捷度、生活便利度和生态宜居程度之一或者任意组合因素；
人口因素包括规划人口指标；
底线管控因素包括是否侵占自然保护地指标、是否侵占永久基本农田指标、是否侵占生态保护红线指标之一或者任意组合；
用地实施情况因素包括周边1km范围(对应15分钟社区生活圈)现状建设用地实施率指标，km表示距离单位千米；
土地供应条件因素包括是否涉及宗地抵押指标、是否完成征地手续指标、是否涉及地灾隐患点指标、是否涉及压覆矿产指标、是否涉及现状林地指标之一或者任意组合。
在本发明的一种优选实施方式中，在步骤S2中，对得到的土地影响因素数据进行预处理，得到其土地预处理影响因素数据的方法包括以下步骤：
S21，对得到的土地影响所有因素数据进行顺序编号，分别为A1、A2、A3、……、Aa，a表示土地影响所有因素数据的总数，且为大于或者等于1的正整数；A1表示土地影响第1因素数据，A2表示土地影响第2因素数据，A3表示土地影响第3因素数据，……，Aa表示土地影响第a因素数据；
S22，对土地影响第b因素数据Ab下的所有指标进行顺序编号，
Figure BDA0002994666870000031
其中，
Figure BDA0002994666870000032
表示空集符号，
Figure BDA0002994666870000035
表示包含符号，{}表示集合符号，分别为Ab,1、Ab,2、Ab,3、……、Ab,b′，b′表示土地影响第b因素数据Ab下的所有指标的总数，且为大于或者等于1的正整数；Ab,1表示土地影响第b因素数据Ab下的第1指标，Ab,2表示土地影响第b因素数据Ab下的第2指标，Ab,3表示土地影响第b因素数据Ab下的第3指标，……，Ab,b′表示土地影响第b因素数据Ab下的第b′指标；c＝1，c″＝1；
S23，依次判断其土地影响第c因素数据Ac下的第c″指标Ac,c″与Ac,c″,max和Ac,c″，min间的关系：
若土地影响第c因素数据Ac下的第c″指标Ac,c″小于或者等于Ac,c″,max，且大于或者等于Ac,c″，min，Ac,c″,max表示土地影响第c因素数据Ac下的第c″指标Ac,c″所对应的预设最大阈值，Ac,c″,min表示土地影响第c因素数据Ac下的第c″指标Ac,c″所对应的预设最小阈值；
Figure BDA0002994666870000033
其中，
Figure BDA0002994666870000034
表示空集符号，
Figure BDA0002994666870000036
表示包含符号，{}表示集合符号；c″表示土地影响第c因素数据Ac下的指标序号；即Ac,c″,min≤Ac,c″≤Ac,c″,max；则将土地影响第c因素数据Ac下的第c″指标Ac,c″保留；
若土地影响第c因素数据Ac下的第c″指标Ac,c″大于Ac,c″,max，或者小于Ac,c″，min，即Ac,c″,min＞Ac,c″，或者Ac,c″,max＜Ac,c″；则将土地影响第c因素数据Ac下的第c″指标Ac,c″筛除；
S24，c″＝c″+1；判断其c″与c′的关系：
若c″≥c′，c′表示土地影响第c因素数据Ac下的所有指标的总数，且为大于或者等于1的正整数；则c＝c+1；执行步骤S25；
若c″＜c′，则返回步骤S23；
S25，c＝c+1；判断c与a的关系：
若c≥a，则预处理完成；
若c＜a，则c″＝1；返回步骤S23。
在本发明的一种优选实施方式中，在步骤S3中，对其得到的土地预处理影响因素数据形成影响因素矩阵的方法包括以下步骤：
S31，获取其土地影响所有因素的总数，设定为M；
S32，获取其土地影响每个因素下的指标数个数，选出其指标数个数最大，设定为N；
S33，构建其矩阵PMN，PMN表示矩阵PMN有M行N列，将其指标值对应依次写入矩阵PMN，若在矩阵PMN中的元素值Pmn处未对应其指标值，Pmn表示在矩阵PMN中第m行第n处的元素值，所述m为小于或者等于M的正整数，所述n为小于或者等于N的正整数，将指标值写为1或0。
在本发明的一种优选实施方式中，在步骤S4中，成熟度系数的计算方法为：
Figure BDA0002994666870000041
其中，Rio表示成熟度系数；
PMN表示影响因素矩阵；
g表示影响因素矩阵PMN所对应的权重系数；
η表示调整系数；
Qij表示供应矩阵；
G表示供应矩阵Qij所对应的权重矩阵；
T表示转置矩阵；
|| ||表示2范数。
在本发明的一种优选实施方式中，在步骤S4中，对其地图标红的处理方法为：
S41，获取地图上已标红所对应的像素值；得到其标红第一像素值和标红第二像素值，标红第一像素值小于标红第二像素值，其标红第一像素值是标红像素值中的最小值，标红第二像素值是标红像素值中的最大值；
S42，判断其标红第一像素值和标红第二像素值所对应的成熟度系数与待标红像素值所对应的成熟度系数关系：
若待标红像素值所对应的成熟度系数大于第二像素值所对应的成熟度系数，则
Figure BDA0002994666870000051
其中，
Figure BDA0002994666870000052
表示待标红像素值，φ1表示第二像素值
Figure BDA0002994666870000053
所对应的成熟度系数，φ2表示第一像素值
Figure BDA0002994666870000054
所对应的成熟度系数，
Figure BDA0002994666870000055
表示第二像素值，
Figure BDA0002994666870000056
表示第一像素值，λ0表示待标红像素值
Figure BDA0002994666870000057
所对应的成熟度系数；int()表示取整函数；
Figure BDA00029946668700000524
为大于或者等于1且小于或者等于10的正整数；
Figure BDA0002994666870000058
μ表示色彩位数；
若待标红像素值所对应的成熟度系数小于第一像素值所对应的成熟度系数，则
Figure BDA0002994666870000059
其中，
Figure BDA00029946668700000510
表示待标红像素值，φ1表示第二像素值
Figure BDA00029946668700000511
所对应的成熟度系数，φ2表示第一像素值
Figure BDA00029946668700000512
所对应的成熟度系数，
Figure BDA00029946668700000513
表示第二像素值，
Figure BDA00029946668700000514
表示第一像素值，λ0表示待标红像素值
Figure BDA00029946668700000515
所对应的成熟度系数；int()表示取整函数；
Figure BDA00029946668700000525
为大于或者等于1且小于或者等于10的正整数；
Figure BDA00029946668700000516
ξ为大于或者等于1且小于或者等于10的正整数；
若待标红像素值所对应的成熟度系数小于或者等于第二像素值所对应的成熟度系数，且大于或者等于第一像素值所对应的成熟度系数，则
Figure BDA00029946668700000517
其中，
Figure BDA00029946668700000518
表示待标红像素值，φ1表示第二像素值
Figure BDA00029946668700000519
所对应的成熟度系数，φ2表示第一像素值
Figure BDA00029946668700000520
所对应的成熟度系数，
Figure BDA00029946668700000521
表示第二像素值，
Figure BDA00029946668700000522
表示第一像素值，λ0表示待标红像素值
Figure BDA00029946668700000523
所对应的成熟度系数；int()表示取整函数；
S43，将待标红像素值标记在待标红地图上。
综上所述，由于采用了上述技术方案，本发明能够在地图上以不同色彩度标记其土地成熟度。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是本发明流程示意框图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
本发明提供了一种土地供应成熟度评价优化方法，如图1所示，包括以下步骤：
S1，获取待供应土地相影响的数据因素，得到其土地影响因素数据；
S2，对步骤S1中得到的土地影响因素数据进行预处理，得到其土地预处理影响因素数据；
S3，对步骤S3中得到的土地预处理影响因素数据形成影响因素矩阵；
S4，根据其步骤S3中得到的影响因素矩阵与其所对应的权重系数，得到其成熟度系数；判断其成熟度系数是否大于或者等于预设成熟度系数：
若计算得到的成熟度系数大于或者等于预设成熟度系数，则对其待供应土地在地图上进行标红处理；
若计算得到的成熟度系数小于预设成熟度系数，则对其待供应土地进行标灰处理。
在本发明的一种优选实施方式中，在步骤S1中，与待供应土地相影响的数据因素包括规划人口、底线管控、用地实施情况、土地供应条件、配套设施完善度、交通便捷度、生活便利度和生态宜居程度之一或者任意组合因素；
人口因素包括规划人口指标；
底线管控因素包括侵占自然保护地指标、侵占永久基本农田指标、侵占生态保护红线指标之一或者任意组合；若为侵占自然保护地，其指标值标记为1，若未侵占自然保护地，其指标值标记为0；若为侵占永久基本农田，其指标值标记为1，若未侵占永久基本农田，其指标值标记为0；若为侵占生态保护红线，其指标值标记为1，若为侵占生态保护红线，其指标值标记为0。
用地实施情况因素包括周边1km范围(对应15分钟社区生活圈)现状在建建设用地实施率指标，km表示距离单位千米。
土地供应条件因素包括涉及宗地抵押指标、完成征地手续指标、涉及地灾隐患点指标、涉及压覆矿产指标、涉及现状林地指标之一或者任意组合。若为涉及宗地抵押，其指标值标记为1，若未涉及宗地抵押，其指标值标记为0；若为完成征地手续，其指标值标记为1，若未完成征地手续，其指标值标记为0；若为涉及地灾隐患点，其指标值标记为1，若未涉及地灾隐患点，其指标值标记为0；若为涉及压覆矿产，其指标值标记为1，若未涉及压覆矿产，其指标值标记为0；若为涉及现状林地，其指标值标记为1，若未涉及现状林地，其指标值标记为0。
在本发明的一种优选实施方式中，在步骤S2中，对得到的土地影响因素数据进行预处理，得到其土地预处理影响因素数据的方法包括以下步骤：
S21，对得到的土地影响所有因素数据进行顺序编号，分别为A1、A2、A3、……、Aa，a表示土地影响所有因素数据的总数，且为大于或者等于1的正整数；A1表示土地影响第1因素数据，A2表示土地影响第2因素数据，A3表示土地影响第3因素数据，……，Aa表示土地影响第a因素数据；
S22，对土地影响第b因素数据Ab下的所有指标进行顺序编号，
Figure BDA0002994666870000081
其中，
Figure BDA0002994666870000082
表示空集符号，
Figure BDA0002994666870000085
表示包含符号，{}表示集合符号，分别为Ab,1、Ab,2、Ab,3、……、Ab,b′，b′表示土地影响第b因素数据Ab下的所有指标的总数，且为大于或者等于1的正整数；Ab,1表示土地影响第b因素数据Ab下的第1指标，Ab,2表示土地影响第b因素数据Ab下的第2指标，Ab,3表示土地影响第b因素数据Ab下的第3指标，……，Ab,b′表示土地影响第b因素数据Ab下的第b′指标；c＝1，c″＝1；
S23，依次判断其土地影响第c因素数据Ac下的第c″指标Ac,c″与Ac,c″,max和Ac,c″，min间的关系：
若土地影响第c因素数据Ac下的第c″指标Ac,c″小于或者等于Ac,c″,max，且大于或者等于Ac,c″，min，Ac,c″,max表示土地影响第c因素数据Ac下的第c″指标Ac,c″所对应的预设最大阈值，Ac,c″,min表示土地影响第c因素数据Ac下的第c″指标Ac,c″所对应的预设最小阈值；
Figure BDA0002994666870000083
其中，
Figure BDA0002994666870000084
表示空集符号，
Figure BDA0002994666870000086
表示包含符号，{}表示集合符号；c″表示土地影响第c因素数据Ac下的指标序号；即Ac,c″,min≤Ac,c″≤Ac,c″,max；则将土地影响第c因素数据Ac下的第c″指标Ac,c″保留；
若土地影响第c因素数据Ac下的第c″指标Ac,c″大于Ac,c″,max，或者小于Ac,c″，min，即Ac,c″,min＞Ac,c″，或者Ac,c″,max＜Ac,c″；则将土地影响第c因素数据Ac下的第c″指标Ac,c″筛除；
S24，c″＝c″+1；判断其c″与c′的关系：
若c″≥c′，c′表示土地影响第c因素数据Ac下的所有指标的总数，且为大于或者等于1的正整数；则c＝c+1；执行步骤S25；
若c″＜c′，则返回步骤S23；
S25，c＝c+1；判断c与a的关系：
若c≥a，则预处理完成；
若c＜a，则c″＝1；返回步骤S23。
在本发明的一种优选实施方式中，在步骤S3中，对其得到的土地预处理影响因素数据形成影响因素矩阵的方法包括以下步骤：
S31，获取其土地影响所有因素的总数，设定为M；
S32，获取其土地影响每个因素下的指标数个数，选出其指标数个数最大，设定为N；
S33，构建其矩阵PMN，PMN表示矩阵PMN有M行N列，将其指标值对应依次写入矩阵PMN，若在矩阵PMN中的元素值Pmn处未对应其指标值，Pmn表示在矩阵PMN中第m行第n处的元素值，所述m为小于或者等于M的正整数，所述n为小于或者等于N的正整数，将指标值写为1或0。
在本发明的一种优选实施方式中，在步骤S4中，成熟度系数的计算方法为：
Figure BDA0002994666870000091
其中，Rio表示成熟度系数；
PMN表示影响因素矩阵；
g表示影响因素矩阵PMN所对应的权重系数；
η表示调整系数；
Qij表示供应矩阵；
G表示供应矩阵Qij所对应的权重矩阵；
T表示转置矩阵；
|| ||表示2范数。
在本发明的一种优选实施方式中，在步骤S4中，对其地图标红的处理方法为：
S41，获取地图上已标红所对应的像素值；得到其标红第一像素值和标红第二像素值，标红第一像素值小于标红第二像素值，其标红第一像素值是标红像素值中的最小值，标红第二像素值是标红像素值中的最大值；
S42，判断其标红第一像素值和标红第二像素值所对应的成熟度系数与待标红像素值所对应的成熟度系数关系：
若待标红像素值所对应的成熟度系数大于第二像素值所对应的成熟度系数，则
Figure BDA0002994666870000101
其中，
Figure BDA0002994666870000102
表示待标红像素值，φ1表示第二像素值
Figure BDA0002994666870000103
所对应的成熟度系数，φ2表示第一像素值
Figure BDA0002994666870000104
所对应的成熟度系数，
Figure BDA0002994666870000105
表示第二像素值，
Figure BDA0002994666870000106
表示第一像素值，λ0表示待标红像素值
Figure BDA0002994666870000107
所对应的成熟度系数；int()表示取整函数；
Figure BDA00029946668700001024
为大于或者等于1且小于或者等于10的正整数；
Figure BDA0002994666870000108
μ表示色彩位数；
若待标红像素值所对应的成熟度系数小于第一像素值所对应的成熟度系数，则
Figure BDA0002994666870000109
其中，
Figure BDA00029946668700001010
表示待标红像素值，φ1表示第二像素值
Figure BDA00029946668700001011
所对应的成熟度系数，φ2表示第一像素值
Figure BDA00029946668700001012
所对应的成熟度系数，
Figure BDA00029946668700001013
表示第二像素值，
Figure BDA00029946668700001014
表示第一像素值，λ0表示待标红像素值
Figure BDA00029946668700001015
所对应的成熟度系数；int()表示取整函数；
Figure BDA00029946668700001025
为大于或者等于1且小于或者等于10的正整数；
Figure BDA00029946668700001016
ξ为大于或者等于1且小于或者等于10的正整数；
若待标红像素值所对应的成熟度系数小于或者等于第二像素值所对应的成熟度系数，且大于或者等于第一像素值所对应的成熟度系数，则
Figure BDA00029946668700001017
其中，
Figure BDA00029946668700001018
表示待标红像素值，φ1表示第二像素值
Figure BDA00029946668700001019
所对应的成熟度系数，φ2表示第一像素值
Figure BDA00029946668700001020
所对应的成熟度系数，
Figure BDA00029946668700001021
表示第二像素值，
Figure BDA00029946668700001022
表示第一像素值，λ0表示待标红像素值
Figure BDA00029946668700001023
所对应的成熟度系数；int()表示取整函数；
S43，将待标红像素值标记在待标红地图上。
在步骤S4中，对其地图标灰的处理方法为：
S41，获取地图上已标灰所对应的像素值；得到其标灰第一像素值和标灰第二像素值，标灰第一像素值小于标灰第二像素值，其标灰第一像素值是标灰像素值中的最小值，标灰第二像素值是标灰像素值中的最大值；
S42，判断其标灰第一像素值和标灰第二像素值所对应的成熟度系数与待标灰像素值所对应的成熟度系数关系：
若待标灰像素值所对应的成熟度系数大于第二像素值所对应的成熟度系数，则
Figure BDA0002994666870000111
其中，
Figure BDA0002994666870000112
表示待标灰像素值，φ1′表示第二像素值
Figure BDA0002994666870000113
所对应的成熟度系数，φ2′表示第一像素值
Figure BDA0002994666870000114
所对应的成熟度系数，
Figure BDA0002994666870000115
表示第二像素值，
Figure BDA0002994666870000116
表示第一像素值，λ0′表示待标灰像素值
Figure BDA0002994666870000117
所对应的成熟度系数；int()表示取整函数；
Figure BDA00029946668700001117
为大于或者等于1且小于或者等于10的正整数；
Figure BDA0002994666870000118
ζ表示亮度或灰度最大值；
若待标灰像素值所对应的成熟度系数小于第一像素值所对应的成熟度系数，则
Figure BDA0002994666870000119
其中，
Figure BDA00029946668700001110
表示待标灰像素值，φ1′表示第二像素值
Figure BDA00029946668700001111
所对应的成熟度系数，φ2′表示第一像素值
Figure BDA00029946668700001112
所对应的成熟度系数，
Figure BDA00029946668700001113
表示第二像素值，
Figure BDA00029946668700001114
表示第一像素值，λ0′表示待标灰像素值
Figure BDA00029946668700001115
所对应的成熟度系数；int()表示取整函数；
Figure BDA00029946668700001118
为大于或者等于1且小于或者等于10的正整数；
Figure BDA00029946668700001116
ξ′为大于或者等于1且小于或者等于10的正整数；
若待标灰像素值所对应的成熟度系数小于或者等于第二像素值所对应的成熟度系数，且大于或者等于第一像素值所对应的成熟度系数，则
Figure BDA0002994666870000121
其中，
Figure BDA0002994666870000122
表示待标灰像素值，φ1′表示第二像素值
Figure BDA0002994666870000123
所对应的成熟度系数，φ2′表示第一像素值
Figure BDA0002994666870000124
所对应的成熟度系数，
Figure BDA0002994666870000125
表示第二像素值，
Figure BDA0002994666870000126
表示第一像素值，λ0′表示待标灰像素值
Figure BDA0002994666870000127
所对应的成熟度系数；int()表示取整函数；
S43，将待标灰像素值标记在待标灰地图上。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。

Claims (5)

1.一种土地供应成熟度评价优化方法，其特征在于，包括以下步骤：
S1，获取待供应土地相影响的数据因素，得到其土地影响因素数据；
S2，对步骤S1中得到的土地影响因素数据进行预处理，得到其土地预处理影响因素数据；
S3，对步骤S3中得到的土地预处理影响因素数据形成影响因素矩阵；
S4，根据其步骤S3中得到的影响因素矩阵与其所对应的权重系数，得到其成熟度系数；判断其成熟度系数是否大于或者等于预设成熟度系数：
若计算得到的成熟度系数大于或者等于预设成熟度系数，则对其待供应土地在地图上进行标红处理；对其地图标红的处理方法为：
S41，获取地图上已标红所对应的像素值；得到其标红第一像素值和标红第二像素值，标红第一像素值小于标红第二像素值，其标红第一像素值是标红像素值中的最小值，标红第二像素值是标红像素值中的最大值；
S42，判断其标红第一像素值和标红第二像素值所对应的成熟度系数与待标红像素值所对应的成熟度系数关系：
若待标红像素值所对应的成熟度系数大于第二像素值所对应的成熟度系数，则
Figure FDA0003262846500000011
其中，
Figure FDA0003262846500000019
表示待标红像素值，φ1表示第二像素值
Figure FDA0003262846500000012
所对应的成熟度系数，φ2表示第一像素值
Figure FDA0003262846500000013
所对应的成熟度系数，
Figure FDA0003262846500000014
表示第二像素值，
Figure FDA0003262846500000015
表示第一像素值，λ0表示待标红像素值
Figure FDA0003262846500000016
所对应的成熟度系数；int()表示取整函数；
Figure FDA0003262846500000017
为大于或者等于1且小于或者等于10的正整数；
Figure FDA0003262846500000018
μ表示色彩位数；
若待标红像素值所对应的成熟度系数小于第一像素值所对应的成熟度系数，则
Figure FDA0003262846500000021
其中，
Figure FDA0003262846500000022
表示待标红像素值，φ1表示第二像素值
Figure FDA0003262846500000023
所对应的成熟度系数，φ2表示第一像素值
Figure FDA0003262846500000024
所对应的成熟度系数，
Figure FDA0003262846500000025
表示第二像素值，
Figure FDA0003262846500000026
表示第一像素值，λ0表示待标红像素值
Figure FDA0003262846500000027
所对应的成熟度系数；int()表示取整函数；
Figure FDA0003262846500000028
为大于或者等于1且小于或者等于10的正整数；
Figure FDA0003262846500000029
ξ为大于或者等于1且小于或者等于10的正整数；
若待标红像素值所对应的成熟度系数小于或者等于第二像素值所对应的成熟度系数，且大于或者等于第一像素值所对应的成熟度系数，则
Figure FDA00032628465000000210
其中，
Figure FDA00032628465000000211
表示待标红像素值，φ1表示第二像素值
Figure FDA00032628465000000212
所对应的成熟度系数，φ2表示第一像素值
Figure FDA00032628465000000213
所对应的成熟度系数，
Figure FDA00032628465000000214
表示第二像素值，
Figure FDA00032628465000000215
表示第一像素值，λ0表示待标红像素值
Figure FDA00032628465000000216
所对应的成熟度系数；int()表示取整函数；
S43，将待标红像素值标记在待标红地图上；
若计算得到的成熟度系数小于预设成熟度系数，则对其待供应土地进行标灰处理。
2.根据权利要求1所述的土地供应成熟度评价优化方法，其特征在于，在步骤S1中，与待供应土地相影响的数据因素包括规划人口、底线管控、用地实施情况、土地供应条件、配套设施完善度、交通便捷度、生活便利度和生态宜居程度之一或者任意组合因素；
人口因素包括规划人口指标；
底线管控因素包括是否侵占自然保护地指标、是否侵占永久基本农田指标、是否侵占生态保护红线指标之一或者任意组合；
用地实施情况因素包括周边1km范围现状在建建设用地实施率指标，km表示距离单位千米；
土地供应条件因素包括是否涉及宗地抵押指标、是否完成征地手续指标、是否涉及地灾隐患点指标、是否涉及压覆矿产指标、是否涉及现状林地指标之一或者任意组合。
3.根据权利要求2所述的土地供应成熟度评价优化方法，其特征在于，在步骤S2中，对得到的土地影响因素数据进行预处理，得到其土地预处理影响因素数据的方法包括以下步骤：
S21，对得到的土地影响所有因素数据进行顺序编号，分别为A1、A2、A3、……、Aa，a表示土地影响所有因素数据的总数，且为大于或者等于1的正整数；A1表示土地影响第1因素数据，A2表示土地影响第2因素数据，A3表示土地影响第3因素数据，……，Aa表示土地影响第a因素数据；
S22，对土地影响第b因素数据Ab下的所有指标进行顺序编号，
Figure FDA0003262846500000031
其中，
Figure FDA0003262846500000032
表示空集符号，
Figure FDA0003262846500000033
表示包含符号，{}表示集合符号，分别为Ab,1、Ab,2、Ab,3、……、Ab,b′，b′表示土地影响第b因素数据Ab下的所有指标的总数，且为大于或者等于1的正整数；Ab,1表示土地影响第b因素数据Ab下的第1指标，Ab,2表示土地影响第b因素数据Ab下的第2指标，Ab,3表示土地影响第b因素数据Ab下的第3指标，……，Ab,b′表示土地影响第b因素数据Ab下的第b′指标；c＝1，c″＝1；
S23，依次判断其土地影响第c因素数据Ac下的第c″指标Ac,c″与Ac,c″,max和Ac,c″，min间的关系：
若土地影响第c因素数据Ac下的第c″指标Ac,c″小于或者等于Ac,c″,max，且大于或者等于Ac,c″，min，Ac,c″,max表示土地影响第c因素数据Ac下的第c″指标Ac,c″所对应的预设最大阈值，Ac,c″,min表示土地影响第c因素数据Ac下的第c″指标Ac,c″所对应的预设最小阈值；
Figure FDA0003262846500000034
其中，
Figure FDA0003262846500000035
表示空集符号，
Figure FDA0003262846500000036
表示包含符号，{}表示集合符号；c″表示土地影响第c因素数据Ac下的指标序号；即Ac,c″,min≤Ac,c″≤Ac,c″,max；则将土地影响第c因素数据Ac下的第c″指标Ac,c″保留；
若土地影响第c因素数据Ac下的第c″指标Ac,c″大于Ac,c″,max，或者小于Ac,c″，min，即Ac,c″,min＞Ac,c″，或者Ac,c″,max＜Ac,c″；则将土地影响第c因素数据Ac下的第c″指标Ac,c″筛除；
S24，c″＝c″+1；判断其c″与c′的关系：
若c″≥c′，c′表示土地影响第c因素数据Ac下的所有指标的总数，且为大于或者等于1的正整数；则c＝c+1；执行步骤S25；
若c″＜c′，则返回步骤S23；
S25，c＝c+1；判断c与a的关系：
若c≥a，则预处理完成；
若c＜a，则c″＝1；返回步骤S23。
4.根据权利要求3所述的土地供应成熟度评价优化方法，其特征在于，在步骤S3中，对其得到的土地预处理影响因素数据形成影响因素矩阵的方法包括以下步骤：
S31，获取其土地影响所有因素的总数，设定为M；
S32，获取其土地影响每个因素下的指标数个数，选出其指标数个数最大，设定为N；
S33，构建其矩阵PMN，PMN表示矩阵PMN有M行N列，将其指标值对应依次写入矩阵PMN，若在矩阵PMN中的元素值Pmn处未对应其指标值，Pmn表示在矩阵PMN中第m行第n处的元素值，所述m为小于或者等于M的正整数，所述n为小于或者等于N的正整数，将指标值写为1或0。
5.根据权利要求3所述的土地供应成熟度评价优化方法，其特征在于，在步骤S4中，成熟度系数的计算方法为：
Figure FDA0003262846500000041
其中，Rio表示成熟度系数；
PMN表示影响因素矩阵；
g表示影响因素矩阵PMN所对应的权重系数；
η表示调整系数；
Qij表示供应矩阵；
G表示供应矩阵Qij所对应的权重矩阵；
T表示转置矩阵；
||||表示2范数。
CN202110325947.XA 2021-03-26 2021-03-26 土地供应成熟度评价优化方法 Active CN113076512B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110325947.XA CN113076512B (zh) 2021-03-26 2021-03-26 土地供应成熟度评价优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110325947.XA CN113076512B (zh) 2021-03-26 2021-03-26 土地供应成熟度评价优化方法

Publications (2)

Publication Number Publication Date
CN113076512A CN113076512A (zh) 2021-07-06
CN113076512B true CN113076512B (zh) 2021-10-26

Family

ID=76610552

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110325947.XA Active CN113076512B (zh) 2021-03-26 2021-03-26 土地供应成熟度评价优化方法

Country Status (1)

Country Link
CN (1) CN113076512B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1554219A (zh) * 2003-12-26 2004-12-15 颜怀玮 生产优质蔬菜饲牧食油花卉观赏植物新种源诸葛菜种子的方法
CN102323221A (zh) * 2011-05-31 2012-01-18 河南农业大学 烟叶成熟度检测方法及检测装置
CN108848892A (zh) * 2018-08-13 2018-11-23 张振花 一种基于无人机的联动式农业收割系统
CN110033044A (zh) * 2019-04-17 2019-07-19 北京建筑大学 基于eagle矩阵的地表覆盖产品整合或融合方法
CN110081867A (zh) * 2019-01-31 2019-08-02 浙江习瑞科技有限责任公司 土地空间规划测量仪

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160108419A1 (en) * 2013-05-08 2016-04-21 Basf Plant Science Company Gmbh Plants having one or more enhanced yield-related traits and a method for making the same
EP3220734A4 (en) * 2014-11-19 2018-06-27 Lamoureux, Normand Method and system for increasing degree days for the cultivation of particular fruit bearing vines in unfavourable climatic regions
CN105868906A (zh) * 2016-03-28 2016-08-17 国网河南省电力公司电力科学研究院 一种优化的区域发展成熟度分析方法
CN106542929A (zh) * 2016-11-02 2017-03-29 广东省农业科学院作物研究所 用于烤烟种植的土壤调理剂及其应用
US10028451B2 (en) * 2016-11-16 2018-07-24 The Climate Corporation Identifying management zones in agricultural fields and generating planting plans for the zones
IL268175B (en) * 2017-01-20 2022-08-01 Greenphyto Pte Ltd A system and method for agriculture
CN108416111A (zh) * 2018-02-08 2018-08-17 广州地理研究所 暴雨洪水关键影响因子与土地利用覆被关系模型构建方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1554219A (zh) * 2003-12-26 2004-12-15 颜怀玮 生产优质蔬菜饲牧食油花卉观赏植物新种源诸葛菜种子的方法
CN102323221A (zh) * 2011-05-31 2012-01-18 河南农业大学 烟叶成熟度检测方法及检测装置
CN108848892A (zh) * 2018-08-13 2018-11-23 张振花 一种基于无人机的联动式农业收割系统
CN110081867A (zh) * 2019-01-31 2019-08-02 浙江习瑞科技有限责任公司 土地空间规划测量仪
CN110033044A (zh) * 2019-04-17 2019-07-19 北京建筑大学 基于eagle矩阵的地表覆盖产品整合或融合方法

Also Published As

Publication number Publication date
CN113076512A (zh) 2021-07-06

Similar Documents

Publication Publication Date Title
Arets et al. Greenhouse gas reporting for the LULUCF sector in the Netherlands: Methodological background, update 2021
CN105046456B (zh) 一种基于gis的土地全生命周期智慧监管方法
Walz et al. Indicator for a monitoring of Germany’s landscape attractiveness
CN111738561A (zh) 一种基于深度学习的区域水土流失动态监测方法
CN103903014B (zh) 一种基于卫星遥感技术的荒漠化判释方法
CN107239756A (zh) 基于高分卫星遥感数据结合地物类型的人口密度分析系统
CN107491758A (zh) 长江流域水体信息提取及其空间编码方法
CN113076512B (zh) 土地供应成熟度评价优化方法
Mackovčin A multi-level ecological network in the Czech Republic: Implementating the Territorial System of Ecological Stability
Gabrovec et al. Land registers as a source of studying long-term land-use changes
Blair et al. Ecosystem services valuation for estuarine and coastal restoration in Florida
Boyacıoğlu et al. Exploring the Applicability of International Village Design Guidelines in Turkey
CN112925842B (zh) 规划实施项目库生成方法
Zhou Use of GIS technology for land resource inventories and modelling for sustainable regional development
Ratner Wetlands management in Vietnam: issues and perspectives
Mellor et al. Remote Sensing Victoria’s Public Land Forests—A Two Tiered Synoptic Approach
Garbisch Information needs in the planning process for wetland creation and restoration
Heron et al. Application of the Climate Vulnerability Index for the ruins of Kilwa Kisiwani and the ruins of Songo Mnara, Tanzania
Xiang et al. " DEPSIR-LM" Model-based Study on Shaping of the Eco-security Pattern for Biodiversity--A Case Study of Guangzhou.
Yakubu et al. Urban expansion and land use changes in Akure, Nigeria using geospatial technique
Sinha Geography and rural development
Hammond et al. A procedure for land capability analysis in Southern Africa, based on computer overlay techniques
Yen Protocol for mapping climate-risks and associated adaptation plans under Asian Mega-Deltas Initiative
Fabos et al. Landscape plan formulation and evaluation
Dammag et al. Optimization Strategies of Sustainable Urban Land Use Planning Using a Hybridized Technique Based on SWOT-AHP, and GIS in Ibb City, Yemen

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