CN109597872B - 一种基于数字地形模型的文显度计算方法 - Google Patents
一种基于数字地形模型的文显度计算方法 Download PDFInfo
- Publication number
- CN109597872B CN109597872B CN201910068170.6A CN201910068170A CN109597872B CN 109597872 B CN109597872 B CN 109597872B CN 201910068170 A CN201910068170 A CN 201910068170A CN 109597872 B CN109597872 B CN 109597872B
- Authority
- CN
- China
- Prior art keywords
- land
- point
- grid
- cultural
- observed
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于数字地形模型的文显度计算方法,包括以下步骤:在GIS中导入文化用地的POI数据以及文化用地所在城市的城市数字高程模型信息,生成数字地形模型图片后分为a*b个格网;将观测点格网与待观测点做连线,将连线穿过的所有格网记为V;将观测点和待观测点中较高者记为M,对V求与M连线的斜率;若V的斜率全部大于待观测点与观测点之间的斜率,则对观测点格网标记为q;文显度=s/S,s为标记为q的格网的总面积,S为a*b个格网的总面积。本发明能够客观定量的评估一个城市文化用地的情况,可为解决城市文化用地评估、规划和监测提供支持和指导。
Description
技术领域
本发明属于城乡规划技术领域,特别涉及一种基于数字地形模型的文显度计算方法。
背景技术
地理信息系统(Geographic Information System,GIS)是一种采集、存储、管理、分析、显示与应用地理信息的计算机技术系统,是分析和处理海量地理数据的通用技术;其主要包括数据输入与编辑、数据管理、数据操作以及数据显示和输出等。地理信息系统作为获取、处理、管理和分析地理空间数据的重要工具、技术和学科,目前得到了广泛关注。
数字地形模型(Digital Terrain Model,DTM)是利用一个任意坐标系中大量选择的已知x、y、z的坐标点对连续地面的一种模拟表示;或者说,DTM就是地形表面形态属性信息的数字表达,是带有空间位置特征和地形属性特征的数字描述,数字高程模型(DigitalElevation Model,DEM)数据加上建筑高程信息可生成DTM。
在快速城镇化的背景下,城市文化逐渐衰落。城市文地是承载城市精神与文化的城市用地,兼具公共文化服务及文化传承创新发展等众多功能。城市文地系统规划旨在解决现代城市规划建设中的文化传承与发展难题。文显度为某一文化用地的显示程度,即能够观测到该文化用地的建设用地面积占城市总建设用地面积之比。文显度是文地系统规划的重要指标,是评估、规划和监测城市精神和文化的重要物质空间可视程度的方法,现亟需一种基于数字地形模型的文显度计算方法。
发明内容
本发明的目的在于提供一种基于数字地形模型的文显度计算方法,以解决上述存在的技术问题。本发明提出了计算文显度的方法,能够客观定量的评估一个城市文化用地的情况,可为解决城市文化用地评估、规划和监测提供支持和指导。
为达到上述目的,本发明采用以下技术方案:
一种基于数字地形模型的文显度计算方法,包括以下步骤:
步骤1,在GIS中导入待计算文显度的文化用地的POI数据;POI数据包括坐标及高度信息;将待计算文显度的文化用地记为P;
步骤2,在GIS中导入带有矢量等高线数据的待计算文显度的文化用地所在城市的城市数字高程模型信息,叠加POI数据中的高度信息得到数字地形模型信息,并生成数字地形模型图片;
步骤3,将生成的数字地形模型图片分为a*b个格网,其中a和b分别为行数和列数;每个格网具有相同高度及相同用地类型;
步骤4,将P定为待观测点;将格网1*1定为观测点格网,将格网1*1高度上升至预设高度作为观测者高度并与P的顶点做连线,将连线穿过的所有格网记为V={V1,V2,V3…Vn},n为连线穿过的格网的总数量;判断比较待观测点与观测点的高度,若待观测点高则将待观测点记为M,观测点记为m,若待观测点低,则待观测点记为m,观测点记为M,按距观测点由近及远的方向依次对V求与M连线的斜率;若V的斜率全部大于待观测点与观测点之间的斜率,则对观测点格网标记为q,否则,则抛弃此观测点格网;
步骤5,重复步骤4,将观测点从格网1*1遍历至格网a*b,最终将Q个格网标记为q;
步骤6,文显度的计算公式为:文显度=s/S;s为Q个标记为q的格网的总面积,S为a*b个格网的总面积。
进一步地,POI数据包括:类别、名称、坐标和高度信息。
进一步地,文化用地包括:
文化精神标识用地,能够促进社会群体提高统一性与文化认同的空间用地;
文化纪念用地,能够供居民日常开展文化纪念仪式活动的空间用地;
宗教文化用地,包括各类宗教信仰所依托的空间用地;
文化遗产用地,具有历史、艺术、科学价值的遗产、遗迹用地;
文化服务设施用地,用于满足居民日常公共文化生活需求的各类空间用地。
进一步地,步骤2具体包括:在locaspace viewer软件的操作选项中选择绘制范围下载,绘制好之后选定下载级数及导出路径信息完成下载;
将下载好的地形文件加载到Arcgis中,继续在操作模块中选择等高线预览,将下载的地形文件导入到locaspace viewer中,调整取点密度和等高线间距至预设数值,在操作模块中选择提取高程进行等高线提取操作,导出grd格式文件;将grd等高线数据生成为shp格式的可读矢量数据;将矢量等高线数据加载到Arcgis中与数字高程模型数据叠加在一起。
进一步地,斜率计算方法为:两点之间的高度差HM-HV及两点之间的绝对距离d的比值,表达式为:K=HM-HV/d。
进一步地,若格网内用地类型为建筑用地,则不考虑作为观测点。
进一步地,酒店和美食城不作为观测点。
进一步地,步骤4中,观测者高度预设值为1.6m。
进一步地,步骤2中,调整取点密度和等高线间距为100m。
与现有技术相比,本发明具有以下有益效果:
本发明创新性的提出了基于数字地形模型的文显度计算方法,可在GIS中根据文化用地POI数据计算出文显度,计算便捷;本发明的基本计算为格网面积的相加运算以及格网累加面积的相除运算,不需借助复杂方程,简单易操作,实现简单。以前没有办法计算文地的文显度,现在通过本发明的方法可对文地的文显度进行快捷有效计算。本发明能够客观定量的评估一个城市文化用地的情况,可为解决城市文化用地评估、规划和监测提供支持和指导。
附图说明
图1是本发明的一种基于数字地形模型的文显度计算方法的流程示意框图;
图2是本发明的实施例中格网划分示意图;
图3是本发明实施例中斜率计算示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细说明。
请参阅图1,本发明的一种基于数字地形模型的文显度计算方法,具体包括以下步骤:
步骤1,在GIS中导入一个待观测的文化用地POI数据,记为P,如中国陕西省西安市大雁塔。
POI数据包括:类别、名称、坐标及楼高信息。文化用地POI数据可以从网络下载。文化用地包括:1)文化精神标识用地,即可促进地区乃至更大范围社会群体提高统一性与文化认同的空间用地;2)文化纪念用地,即可供居民日常开展文化纪念仪式活动,以寄托怀故之思,并抒发情感的空间用地;3)宗教文化用地,即各类宗教信仰所依托的空间用地;4)文化遗产用地,即具有历史、艺术、科学价值的遗产、遗迹用地;5)文化服务设施用地,即满足居民日常公共文化生活需求的各类空间用地。
步骤2,在GIS中导入已下载好的带有矢量等高线数据的城市数字高程模型信息(以下简称DEM),叠加POI的楼高信息得到数字地形模型信息(以下简称DTM),并将DTM生成一张图片。
在locaspace viewer软件的操作选项中选择绘制范围下载,绘制好之后选定下载级数及导出路径等信息即可完成下载。将下载好的谷歌地形文件加载到Arcgis中,继续在操作模块中选择等高线预览,将下载的谷歌地形文件导入到locaspace viewer中,便可直接选择与之前下载文件范围一致的区域进行预览,调整取点密度和等高线间距为100米,在操作模块中选择提取高程进行等高线提取操作,导出grd格式文件。将grd等高线数据生成为shp格式的可读矢量数据。将矢量等高线数据加载到Arcgis中与DEM数据叠加在一起。
步骤3,将生成的DTM图片分为a*b个格网,其中a和b分别为行数和列数,每个格网具有相同高度及相同用地类型。
请参阅图2,将图片划分为6*8的格网,每个格网具有相同高度及相关属性。
步骤4,选择待观测点P,从第1个格网出发,视为观测点。若格网内用地类型为建筑用地,则不考虑作为观测点;否则将格网高度上升1.6米,与P顶点做连线,依次将连线穿过的所有格网记为V。其中,酒店、美食城等建设用地不考虑作为观测点。
例如,将P定为待观测点;将格网1*1定为观测点,将格网1*1高度上升至预设高度作为观测者高度并与P的顶点做连线,将连线穿过的所有格网记为V={V1,V2,V3…Vn},n为连线穿过的格网的总数量;将记为V的格网与P的顶点连线;将观测点m与待观测点P之间的斜率Km与穿过的格网{V1,V2,V3…Vn}与待观测点P之间的斜率{K1,K2,K3…Kn}进行比较;若{K1,K2,K3…Kn}中的值均大于Km,则将格网1*1标记为q。
还可以,区分观测点与待观测点高度,总是将高的记为M。具体的,判断比较待观测点与观测点的高度;若待观测点高则将待观测点记为M,观测点记为m,若待观测点低,则待观测点记为m,观测点记为M。按距观测点由近及远的方向依次对V求与M连线的斜率,记为K,若V的斜率全部大于待观测点与观测点之间的斜率,则对观测点格网进行标记。连线上只要有一个格网的斜率小于待观测点与观测点之间的斜率,则抛弃此观测点,继续第将二个格网作为观测点重复上述操作。斜率计算方法为:V点与M点之间的高度差HM-HV及两点之间的绝对距离d比值,即:K=HM-HV/d。
步骤5,判断比较待观测点与观测点的高度;若待观测点高则将待观测点记为M,观测点记为m,若待观测点低,则待观测点记为m,观测点记为M。按距观测点由近及远的方向依次对V求与M连线的斜率,记为K,若V的斜率全部大于待观测点与观测点之间的斜率,则对观测点格网进行标记。连线上只要有一个格网的斜率小于待观测点与观测点之间的斜率,则抛弃此观测点,继续第将二个格网作为观测点重复上述操作。斜率计算方法为:V点与M点之间的高度差HM-HV及两点之间的绝对距离d比值,即:K=HM-HV/d。
请参阅图3,将大雁塔作为M,人在距大雁塔70米处进行观测,记为m点,中间距大雁塔50处有一中间点V,大雁塔高度为64.517米,中间点为10米,人高度为1.6米,则人与大雁塔之间的斜率为0.89,中间点与大雁塔之间的斜率为1.09。因为中间点与大雁塔之间的斜率大于人与大雁塔之间的斜率,所以抛弃m点,继续考察下一个格网。
步骤6,重复步骤4将所有可作为观测点的格网进行遍历。
步骤7,计算标记的格网面积和si,及除去建设用地格网的总面积和S,进行相除运算,得到文显度,即文显度=si/S。
综上,本发明提供了一种基于数字地形模型(DTM)的文显度计算方法,可为解决城市文化用地评估、规划和监测提供方法支持;本发明能够客观定量的评估一个城市文化用地的情况,同时对城市文化用地规划起到一定的参考作用。
Claims (7)
1.一种基于数字地形模型的文显度计算方法,其特征在于,包括以下步骤:
步骤1,在GIS中导入待计算文显度的文化用地的POI数据;POI数据包括坐标及高度信息;将待计算文显度的文化用地记为P;
步骤2,在GIS中导入带有矢量等高线数据的待计算文显度的文化用地所在城市的城市数字高程模型信息,叠加POI数据中的高度信息得到数字地形模型信息,并生成数字地形模型图片;
步骤3,将生成的数字地形模型图片分为a*b个格网,其中a和b分别为行数和列数;每个格网具有相同高度及相同用地类型;
步骤4,将P定为待观测点;将格网1*1定为观测点格网,将格网1*1高度上升至预设高度作为观测者高度并与P的顶点做连线,将连线穿过的所有格网记为V={V1,V2,V3…Vn},n为连线穿过的格网的总数量;判断比较待观测点与观测点的高度,若待观测点高则将待观测点记为M,观测点记为m,若待观测点低,则待观测点记为m,观测点记为M,按距观测点由近及远的方向依次对V求与M连线的斜率;若V的斜率全部大于待观测点与观测点之间的斜率,则对观测点格网标记为q,否则,则抛弃此观测点格网;
步骤5,重复步骤4,将观测点从格网1*1遍历至格网a*b,最终将Q个格网标记为q;
步骤6,文显度的计算公式为:文显度=s/S;s为Q个标记为q的格网的总面积,S为a*b个格网的总面积。
2.根据权利要求1所述的一种基于数字地形模型的文显度计算方法,其特征在于,POI数据包括:类别、名称、坐标和高度信息。
3.根据权利要求1所述的一种基于数字地形模型的文显度计算方法,其特征在于,文化用地包括:
文化精神标识用地,能够促进社会群体提高统一性与文化认同的空间用地;
文化纪念用地,能够供居民日常开展文化纪念仪式活动的空间用地;
宗教文化用地,包括各类宗教信仰所依托的空间用地;
文化遗产用地,具有历史、艺术、科学价值的遗产、遗迹用地;
文化服务设施用地,用于满足居民日常公共文化生活需求的各类空间用地。
4.根据权利要求1所述的一种基于数字地形模型的文显度计算方法,其特征在于,斜率计算方法为:两点之间的高度差HM-HV及两点之间的绝对距离d的比值,表达式为:
K=(HM-HV)/d。
5.根据权利要求1所述的一种基于数字地形模型的文显度计算方法,其特征在于,若格网内用地类型为建筑用地,则不考虑作为观测点。
6.根据权利要求1所述的一种基于数字地形模型的文显度计算方法,其特征在于,酒店和美食城不作为观测点。
7.根据权利要求1所述的一种基于数字地形模型的文显度计算方法,其特征在于,步骤4中,观测者高度预设值为1.6m。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910068170.6A CN109597872B (zh) | 2019-01-24 | 2019-01-24 | 一种基于数字地形模型的文显度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910068170.6A CN109597872B (zh) | 2019-01-24 | 2019-01-24 | 一种基于数字地形模型的文显度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109597872A CN109597872A (zh) | 2019-04-09 |
CN109597872B true CN109597872B (zh) | 2021-07-09 |
Family
ID=65965359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910068170.6A Active CN109597872B (zh) | 2019-01-24 | 2019-01-24 | 一种基于数字地形模型的文显度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109597872B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100953219B1 (ko) * | 2009-06-30 | 2010-04-15 | 주식회사 일도엔지니어링 | Gis를 이용한 예비 토취장 후보지 선정 시스템 및 그 방법 |
CN102054294A (zh) * | 2011-01-19 | 2011-05-11 | 中国科学院地理科学与资源研究所 | 基于曲面论和优化控制理论的曲面建模方法 |
CN103456041B (zh) * | 2013-08-28 | 2016-06-22 | 中国人民解放军海军大连舰艇学院 | 一种基于s-57电子海图数据的三维地形和雷达地形生成方法 |
KR101750885B1 (ko) * | 2015-09-07 | 2017-06-26 | 민부근 | 행정 구역도를 바탕으로 하는 통합 정보 검색 방법 |
CN109165458A (zh) * | 2018-09-05 | 2019-01-08 | 南京林业大学 | 基于bim的地表建筑太阳辐射估算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10371860B2 (en) * | 2013-05-22 | 2019-08-06 | S.S. Papadopulos & Associates, Inc. | Simultaneous multi-event universal kriging methods for spatio-temporal data analysis and mapping |
-
2019
- 2019-01-24 CN CN201910068170.6A patent/CN109597872B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100953219B1 (ko) * | 2009-06-30 | 2010-04-15 | 주식회사 일도엔지니어링 | Gis를 이용한 예비 토취장 후보지 선정 시스템 및 그 방법 |
CN102054294A (zh) * | 2011-01-19 | 2011-05-11 | 中国科学院地理科学与资源研究所 | 基于曲面论和优化控制理论的曲面建模方法 |
CN103456041B (zh) * | 2013-08-28 | 2016-06-22 | 中国人民解放军海军大连舰艇学院 | 一种基于s-57电子海图数据的三维地形和雷达地形生成方法 |
KR101750885B1 (ko) * | 2015-09-07 | 2017-06-26 | 민부근 | 행정 구역도를 바탕으로 하는 통합 정보 검색 방법 |
CN109165458A (zh) * | 2018-09-05 | 2019-01-08 | 南京林业大学 | 基于bim的地表建筑太阳辐射估算方法 |
Non-Patent Citations (2)
Title |
---|
传统城市空间匾额的语言信息及文化关联性研究--以晋中盆地历史城市为例;李欣鹏 等;《城市文化与城市设计》;20181231;第25卷(第4期);第27-31页 * |
夏门古村落人居环境规划中的"自然智慧"初探;王树声 等;《西安建筑科技大学学报(社会科学版)》;20080930;第27卷(第3期);第50-54页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109597872A (zh) | 2019-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lehner et al. | Digital geoTwin Vienna: Towards a digital twin city as geodata hub | |
Yu et al. | Automated derivation of urban building density information using airborne LiDAR data and object-based method | |
Vivoni et al. | Generation of triangulated irregular networks based on hydrological similarity | |
CN109215470A (zh) | 一种基于bim+gis的电子沙盘系统及其构建方法 | |
CN108921943B (zh) | 一种基于车道级高精度地图的道路三维模型建模方法 | |
Dmowska et al. | High resolution dasymetric model of US demographics with application to spatial distribution of racial diversity | |
CN102708587A (zh) | 一种快速获取三维建筑信息的方法及系统 | |
CN110458333A (zh) | 一种基于POIs数据的人口空间分布预测方法及系统 | |
CN112819340A (zh) | 一种基于多源数据的城市洪涝灾害动态评估方法 | |
CN112819319A (zh) | 一种城市活力与空间社会特征相关性的测度方法及应用 | |
CN107392252A (zh) | 计算机深度学习图像特征并量化感知度的方法 | |
CN113658292A (zh) | 气象数据色斑图生成方法、装置、设备及存储介质 | |
CN106202708B (zh) | 一种将勘察测绘得出的cad地形图快速入库的方法 | |
CN112328929A (zh) | 基于矢量瓦片地图发布技术的人流密度动态可视化方法 | |
Walz et al. | Integrating time and the third spatial dimension in landscape structure analysis | |
Cerreta et al. | A complex values map of marginal urban landscapes: An experiment in Naples (Italy) | |
CN107944089A (zh) | 一种基于现状视廊的地块限高分析系统及其分析方法 | |
CN109597872B (zh) | 一种基于数字地形模型的文显度计算方法 | |
Zhou et al. | Automated extraction of 3D vector topographic feature line from terrain point cloud | |
CN116415499B (zh) | 一种社区舒适感模拟预测方法 | |
CN109376638B (zh) | 一种基于遥感图像和地理信息系统的文地率计算方法 | |
KR100764870B1 (ko) | 수치지도의 등고선 자동입력장치와 그 입력방법 | |
CN114490907A (zh) | 名城管理数据库的构建方法、装置及存储介质 | |
CN111209355B (zh) | 企业亩均效益专题地图可视化制图方法 | |
CN114443787A (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 |