CN112700137A - 一种基于地理信息数据的流域单元用水量划分方法 - Google Patents

一种基于地理信息数据的流域单元用水量划分方法 Download PDF

Info

Publication number
CN112700137A
CN112700137A CN202011627663.8A CN202011627663A CN112700137A CN 112700137 A CN112700137 A CN 112700137A CN 202011627663 A CN202011627663 A CN 202011627663A CN 112700137 A CN112700137 A CN 112700137A
Authority
CN
China
Prior art keywords
water consumption
basin
data
population
administrative
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
CN202011627663.8A
Other languages
English (en)
Other versions
CN112700137B (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.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen University
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 Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN202011627663.8A priority Critical patent/CN112700137B/zh
Publication of CN112700137A publication Critical patent/CN112700137A/zh
Application granted granted Critical
Publication of CN112700137B publication Critical patent/CN112700137B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • 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/06—Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种基于地理信息数据的流域单元用水量划分方法，包括以下步骤：S1：输入地理信息数据；S2：裁剪栅格数据：提取每个区域的Polygon数据，并对土地利用栅格数据和人口密度栅格数据进行掩膜裁剪操作；分别计算行政区内各流域的人口比例和各类土地利用面积比例；S3：输入统计数据：输入行政区的城镇与农村人口统计数据以及各项用水量数据；S4：划分用水量：计算各流域的当前城镇人口和农村人口数量；计算行政区内各流域的各项用水量，从而将用水量按照流域单元划分。本发明提供一种基于地理信息数据的流域单元用水量划分方法，解决了目前测算行政区域在流域内的用水量的主要手段是经验性的方法，还未形成一种通用的方法的问题。

Description

一种基于地理信息数据的流域单元用水量划分方法
技术领域
本发明涉及水资源管理技术领域，更具体的，涉及一种基于地理信息数据的流域单元用水量划分方法。
背景技术
水量分配能有效地促进提高水资源利用效率，在水量分配方案的制定中，需要综合考虑行政区域在流域内的用水当前情况。目前，我国仅对大流域的用水量进行统计，并公布相关流域的水资源公报，而小流域的各行政单元的流域内用水量并未进行公布。因此，在制定小流域的水量分配方案时，需要对行政区域在流域内的用水量进行合理测算，作为水量分配方案的依据。目前测算行政区域在流域内的用水量的主要手段是经验性的方法，还未形成一种通用的方法。
现有技术中，如2018-04-06公开的中国专利，一种城镇化下平原区分布式产流计算方法，公开号为CN107886444A，相对于传统技术方法，可直接通过空间数据(土地利用、子流域边界、行政区边界)实现产流量的计算，但不适用于测算行政区域在流域内的用水量。
发明内容
本发明为克服目前测算行政区域在流域内的用水量的主要手段是经验性的方法，还未形成一种通用的方法的技术缺陷，提供一种基于地理信息数据的流域单元用水量划分方法。
为解决上述技术问题，本发明的技术方案如下：
一种基于地理信息数据的流域单元用水量划分方法，包括以下步骤：
S1：输入地理信息数据；
所述地理信息数据包括行政区流域边界数据、土地利用栅格数据和人口密度栅格数据；
S2：裁剪栅格数据：
提取行政区流域边界数据中每个区域的Polygon数据，并根据每个区域的Polygon数据对土地利用栅格数据和人口密度栅格数据进行掩膜裁剪操作；
根据经过掩膜裁剪操作之后的土地利用栅格数据和人口密度栅格数据分别计算行政区内各流域的人口比例和各类土地利用面积比例；
S3：输入统计数据：
输入行政区的城镇与农村人口统计数据以及各项用水量数据；
S4：划分用水量：
根据行政区内各流域的人口比例，计算各流域的当前城镇人口和农村人口数量；
根据行政区内各流域的各类土地利用面积比例以及当前城镇人口和农村人口数量，计算行政区内各流域的各项用水量，从而将用水量按照流域单元划分。
优选的，在步骤S2中，通过以下公式进行掩膜裁剪操作：
Figure BDA0002875316080000021
其中，Rraster为栅格数据，
Figure BDA0002875316080000022
为行政区流域边界矢量图形的第i个对象，i＝1,2,...,N，N为行政区流域边界矢量图形中的对象总数，Mi为Rraster使用
Figure BDA0002875316080000023
进行掩膜裁剪操作的区域栅格数据，若没有重叠区域，则Mi为空。
优选的，在步骤S2中，行政区内各流域的各类土地利用面积比例包括行政区内各流域的耕地面积比例、林地面积比例和建设用地面积比例。
优选的，在步骤S4中，通过以下公式计算各流域的当前城镇人口和农村人口数量：
Pba sin＝kba sin×Pdistrict
其中，Pba sin为行政区内流域的当前城镇人口或农村人口数，kba sin为行政区内流域的人口数比例，Pdistrict为行政区的城镇人口或农村人口数。
优选的，在步骤S4中，根据以下公式计算行政区内各流域的各项用水量：
Figure BDA0002875316080000024
其中，
Figure BDA0002875316080000025
为行政区内流域i的各项用水量，
Figure BDA0002875316080000026
为行政区内流域i相应用水量划分比例，Wdistrict为行政区的各项用水量。
优选的，各项用水量包括农田灌溉用水量、林牧渔业用水量、工业用水量、城镇公共用水量、居民生活用水量和生态环境用水量。
优选的，
农田灌溉用水量根据耕地面积比例计算；
林牧渔业用水量根据林地面积比例计算；
工业用水量根据城镇建设用地面积比例计算；
城镇公共用水量根据城镇人口比例计算；
居民生活用水量根据总人口比例计算；
生态环境用水量根据城镇人口比例计算。
优选的，还包括根据经过掩膜裁剪操作之后的土地利用栅格数据进行地图要素可视化。
优选的，进行地图要素可视化包括以下步骤：
S5.1：构建土地利用图的画布；
S5.2：在画布上依次绘制每个区域的土地利用图；
S5.3：在画布上绘制行政区流域边界；
S5.4：构建指向标的矢量图形，并将指向标绘制在土地利用图上，以及添加图例，得到绘制好的土地利用图；
S5.5：保存绘制好的土地利用图。
优选的，在步骤S5.3中，还包括：
根据以下公式设置土地利用图的图框x轴和y轴的范围，使得行政区流域边界的图框内两侧留有空白：
Xmax＝max(lon)+[max(lon)-min(lon)]/8
Xmin＝min(lon)-[max(lon)-min(lon)]/8
Ymax＝max(lat)+[max(lat)-min(lat)]/8
Ymin＝min(lat)-[max(lat)-min(lat)/8
其中，Xmax为x轴范围的最大值，Xmin为x轴范围的最小值，Ymax为y轴范围的最大值，Ymin为y轴范围的最小值，lon为行政区边界的经度，lat为行政区边界的纬度。
与现有技术相比，本发明技术方案的有益效果是：
本发明提供了一种基于地理信息数据的流域单元用水量划分方法，根据行政区的各项用水量，利用各流域的各类土地利用面积比例以及当前城镇人口和农村人口数量计算各流域的各项用水量，规避了经验性的测算方法。
附图说明
图1为本发明的技术方案实施步骤流程图；
图2为本发明中输入的数据文件的示意图；
图3为本发明中一实施的潮州市县级行政区示意图；
图4为本发明中一实施的潮州市内流域划分示意图；
图5为本发明中一实施例的潮州市潮安区用水量；
图6为本发明中一实施例的潮州市县级行政区2018年人口数量；
图7为本发明中一实施例的潮州市潮安区用水量划分结果；
图8为本发明中一实施例的潮州市土地利用图。
具体实施方式
附图仅用于示例性说明，不能理解为对本专利的限制；
为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；
对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
如图1所示，一种基于地理信息数据的流域单元用水量划分方法，包括以下步骤：
S1：输入地理信息数据；
所述地理信息数据包括行政区流域边界数据、土地利用栅格数据和人口密度栅格数据；
其中，通过Python第三方库geopandas的GeoDataFrame模块中的from_file函数根据行政区流域边界的矢量图形文件的存储路径读取行政区流域边界数据；通过rasterio中的open函数，根据土地利用和人口密度的栅格数据文件的存储路径读取栅格数据，实现地理信息数据的输入；
行政区流域边界的矢量图形文件的文件格式为.shp，土地利用的栅格数据文件的文件格式为.tif，人口密度的栅格数据文件为Grid文件；
S2：裁剪栅格数据：
提取行政区流域边界数据中每个区域的Polygon数据，并根据每个区域的Polygon数据通过Python第三方库rasterio中的mask.mask函数对土地利用栅格数据和人口密度栅格数据进行掩膜裁剪操作，将所需要的栅格数据提取出来；
根据经过掩膜裁剪操作之后的土地利用栅格数据和人口密度栅格数据分别计算行政区内各流域的人口比例和各类土地利用面积比例；
其中，通过rasterio的open函数依次读取经过掩膜裁剪操作之后的区域人口密度栅格数据，采用rasterio中的read函数读取每个栅格的变量值，计算每个区域的人口数，最后以行政区为整体，计算行政区域内的各流域的人口数比例；
通过rasterio的open函数依次读取经过掩膜裁剪操作之后的区域土地利用栅格数据，采用rasterio中的read函数读取每个栅格的变量值，变量值与土地利用类型是一一对应的，根据变量值确定土地利用类型，并以此计算每个区域的耕地面积、林地面积和城镇建设用地面积，最后以行政区为整体，计算行政区域内各流域的耕地面积、林地面积和建设用地面积比例；
S3：输入统计数据：
通过Python第三方库pandas的read_csv函数，根据行政区的城镇与农村人口统计数据文件以及各项用水量数据文件的存储路径读取数据，实现行政区的城镇与农村人口统计数据以及各项用水量数据的输入；
S4：划分用水量：
根据行政区内各流域的人口比例，计算各流域的当前城镇人口和农村人口数量；
根据行政区内各流域的各类土地利用面积比例以及当前城镇人口和农村人口数量，计算行政区内各流域的各项用水量，从而将用水量按照流域单元划分。
实施例2
更具体的，在步骤S2中，通过以下公式进行掩膜裁剪操作：
Figure BDA0002875316080000051
其中，Rraster为栅格数据，
Figure BDA0002875316080000052
为行政区流域边界矢量图形的第i个对象，i＝1,2,...,N，N为行政区流域边界矢量图形中的对象总数，Mi为Rraster使用
Figure BDA0002875316080000053
进行掩膜裁剪操作的区域栅格数据，若没有重叠区域，则Mi为空。
更具体的，在步骤S2中，行政区内各流域的各类土地利用面积比例包括行政区内各流域的耕地面积比例、林地面积比例和建设用地面积比例。
更具体的，在步骤S4中，通过以下公式计算各流域的当前城镇人口和农村人口数量：
Pba sin＝kba sin×Pdistrict
其中，Pba sin为行政区内流域的当前城镇人口或农村人口数，kba sin为行政区内流域的人口数比例，Pdistrict为行政区的城镇人口或农村人口数。
更具体的，在步骤S4中，根据以下公式计算行政区内各流域的各项用水量：
Figure BDA0002875316080000061
其中，
Figure BDA0002875316080000062
为行政区内流域i的各项用水量，
Figure BDA0002875316080000063
为行政区内流域i相应用水量划分比例，Wdistrict为行政区的各项用水量。
更具体的，各项用水量包括农田灌溉用水量、林牧渔业用水量、工业用水量、城镇公共用水量、居民生活用水量和生态环境用水量。
更具体的，
农田灌溉用水量根据耕地面积比例计算；
林牧渔业用水量根据林地面积比例计算；
工业用水量根据城镇建设用地面积比例计算；
城镇公共用水量根据城镇人口比例计算；
居民生活用水量根据总人口比例计算；
生态环境用水量根据城镇人口比例计算。
更具体的，还包括根据经过掩膜裁剪操作之后的土地利用栅格数据进行地图要素可视化。
在具体实施过程中，通过Python的第三方库geopandas和matplotlib进行地图要素可视化。
更具体的，进行地图要素可视化包括以下步骤：
S5.1：采用matplotlib中pyplot.subplots函数构建土地利用图的画布；
S5.2：采用matplotlib中pyplot.pcolormesh函数在画布上依次绘制每个区域的土地利用图；
S5.3：采用geopandas中的plot函数在画布上绘制行政区流域边界；
S5.4：采用matplotlib.patches中的Polygon函数构建指向标的矢量图形，并通过matplotlib中的add_patch函数将指向标绘制在土地利用图上，以及通过legend函数添加图例，得到绘制好的土地利用图；
S5.5：通过matplotlib中的savefig函数保存绘制好的土地利用图到指定路径。
更具体的，在步骤S5.3中，还包括：
根据以下公式设置土地利用图的图框x轴和y轴的范围，使得行政区流域边界的图框内两侧留有空白：
Xmax＝max(lon)+[max(lon)-min(lon)]/8
Xmin＝min(lon)-[max(lon)-min(lon)]/8
Ymax＝max(lat)+[max(lat)-min(lat)]/8
Ymin＝min(lat)-[max(lat)-min(lat)/8
其中，Xmax为x轴范围的最大值，Xmin为x轴范围的最小值，Ymax为y轴范围的最大值，Ymin为y轴范围的最小值，lon为行政区边界的经度，lat为行政区边界的纬度。
在具体实施过程中，通过Python编程实现步骤S2和步骤S4中的数据处理以及地图要素可视化过程，并且用class()和def()的方式封装成类函数，以便于在其他程序进行调用。
实施例3
在本实施例中，采用所述的一种基于地理信息数据的流域单元用水量划分方法对中国粤东地区的潮州市各县级行政区的用水量进行划分：
S1：将行政区各流域用水量划分所需要的潮州市行政区流域边界的矢量图形文件、土地利用栅格数据和中国人口密度栅格数据存储在文件夹中，如图2-4所示，通过geopandas的GeoDataFrame模块中的from_file函数读取矢量图形文件，rasterio中的open函数读取栅格数据文件；
S2：通过Python第三方库rasterio进行栅格数据掩膜裁剪操作，并计算各流域的人口比例和各类土地利用面积比例：
使用for循环，依次读取行政区流域边界的每个图形对象，使用rasterio中的mask.mask对栅格数据进行掩膜裁剪操作；
使用rasterio中的open函数依次新建栅格数据，并依次将掩膜裁剪操作后的栅格数据写入新建的栅格数据文件中；
使用for循环，依次使用rasterio中的open函数读取掩膜裁剪操作后的每个区域的栅格数据，并计算每个行政区的各流域的耕地面积比例、林地面积比例、城镇建设用地面积比例、城镇人口比例和总人口比例；
S3.采用Python第三方库pandas的read_csv函数，根据潮州市行政区用水量和人口数据文件的存储路径读取数据，数据排列格式要求如图5和图6所示，读取的数据格式为pandas.DataFrame，并根据每个区域的人口比例计算每个区域现有的城镇人口和农村人口数量；
S4.根据行政区域内各流域的耕地面积比例、林地面积比例、城镇建设用地面积比例、城镇人口比例和总人口比例，采用for循环，分别对潮安区、潮州市区和饶平县的各流域用水量进行测算：行政区域内流域的农田灌溉用水量根据耕地面积比例计算，林牧渔业用水量根据林地面积比例计算，工业用水量根据城镇建设用地面积比例计算，城镇公共用水量根据城镇人口比例计算，居民生活用水量根据总人口比例计算，生态环境用水量根据城镇人口比例计算；其中，潮安区的流域用水量划分结果如图7所示；并对潮州市的土地利用进行可视化，主要利用Python第三方库geopandas和matplotlib绘制各个地图要素：
S5.1：采用matplotlib中pyplot.subplots函数构建土地利用图的画布；
S5.2：土地利用图采用matplotlib中pyplot.pcolormesh函数依次进行绘制，并通过matplotlib中的colors.ListedColormap函数设置颜色带，并设置颜色带的划分范围；
S5.3：采用geopandas中的plot函数，绘制行政区流域边界，将facecolor的参数设置为none(表示无颜色)，使得填充颜色为无填充，edgecolor的颜色设置为k(表示黑色)，使得边缘的颜色为黑色；
为了使得行政区流域边界的图框内两侧留有空白，增加x轴与y轴的范围；
S5.4：采用matplotlib.patches中的Polygon函数构建指向标的矢量图形，并通过matplotlib中的add_patch将指向标可视化，其中color参数设置为k，表示黑色；通过legend函数添加图例；
S5.5：采用matplotlib中的savefig函数保存潮州市土地利用图到指定路径，并设置dpi参数为1200，保证有较高的分辨率，绘制结果如图8所示。
显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。

Claims (10)

1.一种基于地理信息数据的流域单元用水量划分方法，其特征在于，包括以下步骤：
S1：输入地理信息数据；
所述地理信息数据包括行政区流域边界数据、土地利用栅格数据和人口密度栅格数据；
S2：裁剪栅格数据：
提取行政区流域边界数据中每个区域的Polygon数据，并根据每个区域的Polygon数据对土地利用栅格数据和人口密度栅格数据进行掩膜裁剪操作；
根据经过掩膜裁剪操作之后的土地利用栅格数据和人口密度栅格数据分别计算行政区内各流域的人口比例和各类土地利用面积比例；
S3：输入统计数据：
输入行政区的城镇与农村人口统计数据以及各项用水量数据；
S4：划分用水量：
根据行政区内各流域的人口比例，计算各流域的当前城镇人口和农村人口数量；
根据行政区内各流域的各类土地利用面积比例以及当前城镇人口和农村人口数量，计算行政区内各流域的各项用水量，从而将用水量按照流域单元划分。
2.根据权利要求1所述的一种基于地理信息数据的流域单元用水量划分方法，其特征在于，在步骤S2中，通过以下公式进行掩膜裁剪操作：
Figure FDA0002875316070000011
其中，Rraster为栅格数据，
Figure FDA0002875316070000012
为行政区流域边界矢量图形的第i个对象，i＝1,2,...,N，N为行政区流域边界矢量图形中的对象总数，Mi为Rraster使用
Figure FDA0002875316070000013
进行掩膜裁剪操作的区域栅格数据，若没有重叠区域，则Mi为空。
3.根据权利要求1所述的一种基于地理信息数据的流域单元用水量划分方法，其特征在于，在步骤S2中，行政区内各流域的各类土地利用面积比例包括行政区内各流域的耕地面积比例、林地面积比例和建设用地面积比例。
4.根据权利要求1所述的一种基于地理信息数据的流域单元用水量划分方法，其特征在于，在步骤S4中，通过以下公式计算各流域的当前城镇人口和农村人口数量：
Pbasin＝kbasin×Pdistrict
其中，Pbasin为行政区内流域的当前城镇人口或农村人口数，kbasin为行政区内流域的人口数比例，Pdistrict为行政区的城镇人口或农村人口数。
5.根据权利要求1所述的一种基于地理信息数据的流域单元用水量划分方法，其特征在于，在步骤S4中，根据以下公式计算行政区内各流域的各项用水量：
Figure FDA0002875316070000021
其中，
Figure FDA0002875316070000022
为行政区内流域i的各项用水量，
Figure FDA0002875316070000023
为行政区内流域i相应用水量划分比例，Wdistrict为行政区的各项用水量。
6.根据权利要求5所述的一种基于地理信息数据的流域单元用水量划分方法，其特征在于，各项用水量包括农田灌溉用水量、林牧渔业用水量、工业用水量、城镇公共用水量、居民生活用水量和生态环境用水量。
7.根据权利要求6所述的一种基于地理信息数据的流域单元用水量划分方法，其特征在于，
农田灌溉用水量根据耕地面积比例计算；
林牧渔业用水量根据林地面积比例计算；
工业用水量根据城镇建设用地面积比例计算；
城镇公共用水量根据城镇人口比例计算；
居民生活用水量根据总人口比例计算；
生态环境用水量根据城镇人口比例计算。
8.根据权利要求1所述的一种基于地理信息数据的流域单元用水量划分方法，其特征在于，还包括根据经过掩膜裁剪操作之后的土地利用栅格数据进行地图要素可视化。
9.根据权利要求8所述的一种基于地理信息数据的流域单元用水量划分方法，其特征在于，进行地图要素可视化包括以下步骤：
S5.1：构建土地利用图的画布；
S5.2：在画布上依次绘制每个区域的土地利用图；
S5.3：在画布上绘制行政区流域边界；
S5.4：构建指向标的矢量图形，并将指向标绘制在土地利用图上，以及添加图例，得到绘制好的土地利用图；
S5.5：保存绘制好的土地利用图。
10.根据权利要求9所述的一种基于地理信息数据的流域单元用水量划分方法，其特征在于，在步骤S5.3中，还包括：
根据以下公式设置土地利用图的图框x轴和y轴的范围，使得行政区流域边界的图框内两侧留有空白：
Xmax＝max(lon)+[max(lon)-min(lon)]/8
Xmin＝min(lon)-[max(lon)-min(lon)]/8
Ymax＝max(lat)+[max(lat)-min(lat)]/8
Ymin＝min(lat)-[max(lat)-min(lat)/8
其中，Xmax为x轴范围的最大值，Xmin为x轴范围的最小值，Ymax为y轴范围的最大值，Ymin为y轴范围的最小值，lon为行政区边界的经度，lat为行政区边界的纬度。
CN202011627663.8A 2020-12-30 2020-12-30 一种基于地理信息数据的流域单元用水量划分方法 Active CN112700137B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011627663.8A CN112700137B (zh) 2020-12-30 2020-12-30 一种基于地理信息数据的流域单元用水量划分方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011627663.8A CN112700137B (zh) 2020-12-30 2020-12-30 一种基于地理信息数据的流域单元用水量划分方法

Publications (2)

Publication Number Publication Date
CN112700137A true CN112700137A (zh) 2021-04-23
CN112700137B CN112700137B (zh) 2024-06-25

Family

ID=75513366

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011627663.8A Active CN112700137B (zh) 2020-12-30 2020-12-30 一种基于地理信息数据的流域单元用水量划分方法

Country Status (1)

Country Link
CN (1) CN112700137B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113918572A (zh) * 2021-10-29 2022-01-11 中国农业科学院农业资源与农业区划研究所 基于行政区划边界的农业统计数据自动匹配与校正的方法
CN114328782A (zh) * 2021-12-24 2022-04-12 中科三清科技有限公司 地理矢量数据处理方法、电子设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101305864B1 (ko) * 2013-04-26 2013-09-09 (주)웹솔루스 지자체 행정정보를 활용한 물사용량 평가방법 및 그 시스템
CN106503232A (zh) * 2016-10-28 2017-03-15 中国水利水电科学研究院 一种社会经济用水空间展布方法
CN107886444A (zh) * 2017-10-11 2018-04-06 南京大学 一种城镇化下平原区分布式产流计算方法
CN111475592A (zh) * 2020-03-04 2020-07-31 中山大学 一种面向矢量图形文件的流域地图制作方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101305864B1 (ko) * 2013-04-26 2013-09-09 (주)웹솔루스 지자체 행정정보를 활용한 물사용량 평가방법 및 그 시스템
CN106503232A (zh) * 2016-10-28 2017-03-15 中国水利水电科学研究院 一种社会经济用水空间展布方法
CN107886444A (zh) * 2017-10-11 2018-04-06 南京大学 一种城镇化下平原区分布式产流计算方法
CN111475592A (zh) * 2020-03-04 2020-07-31 中山大学 一种面向矢量图形文件的流域地图制作方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113918572A (zh) * 2021-10-29 2022-01-11 中国农业科学院农业资源与农业区划研究所 基于行政区划边界的农业统计数据自动匹配与校正的方法
CN113918572B (zh) * 2021-10-29 2022-06-07 中国农业科学院农业资源与农业区划研究所 基于行政区划边界的农业统计数据自动匹配与校正的方法
CN114328782A (zh) * 2021-12-24 2022-04-12 中科三清科技有限公司 地理矢量数据处理方法、电子设备及存储介质
CN114328782B (zh) * 2021-12-24 2022-11-01 中科三清科技有限公司 地理矢量数据处理方法、电子设备及存储介质

Also Published As

Publication number Publication date
CN112700137B (zh) 2024-06-25

Similar Documents

Publication Publication Date Title
Gharbia et al. Land use scenarios and projections simulation using an integrated GIS cellular automata algorithms
CN112700137A (zh) 一种基于地理信息数据的流域单元用水量划分方法
Bishop‐Taylor et al. Surface‐water dynamics and land use influence landscape connectivity across a major dryland region
CN108255864A (zh) 基于分布式存储和分布式计算的影像地图服务发布方法
He et al. Offshore Island Connection Line: A new perspective of coastal urban development boundary simulation and multi-scenario prediction
Hamidi Atmospheric investigation of frontal dust storms in Southwest Asia
CN106991330A (zh) 一种三维地图上数据访问权限控制的方法
McCormack et al. Feature-based derivation of drainage networks
Bisselink et al. Impact of a changing climate, land use, and water usage on water resources in the Danube river basin
CN115690286B (zh) 三维地形生成方法、终端设备及计算机可读存储介质
Li et al. Synergistic ecological network approach for sustainable development of highly urbanized area in the Bay Bottom region: A study in Chengyang District, Qingdao
CN113066176A (zh) 一种地图数据处理方法、装置、设备及存储介质
CN106908797A (zh) 基于gis的x波段面雨量雷达数据提取方法
CN106326492B (zh) 空间矢量数据生成方法及装置
Hsiang et al. Countries without borders: reterritorializing the Maldives
CN110992134B (zh) 一种海洋预报产品的定制方法及终端
CN113066150B (zh) 一种地图数据处理方法、装置、设备及存储介质
WARDHANI et al. COMPARATIVE STUDY OF THE PATTERNS AND CHARACTERISTICS URBAN MORPHOLOGY IN THE OLD CITY, CASE STUDY ROAD AND BLOCK PATTERNS BENGKULU AND SINGAPORE.
Haase et al. Urban land use in the global context
Li et al. [Retracted] Analysis of the Spatial Distribution Pattern of the Urban Landscape in the Central Plains under the Influence of Multiscale and Multilevel Morphological Geomorphology
Guo et al. Study of ecosystem service functions in typical receiving areas of the South-to-North Water Diversion Central Route based on a set of long time series
CN113066179A (zh) 一种地图数据处理方法、装置、设备及存储介质
CN106503232A (zh) 一种社会经济用水空间展布方法
Sopsop et al. Spatio-temporal pattern of landscape change due to urbanization: a case of Batangas City
CN115062564B (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