CN111681313B - 一种基于数字地形图的空间视域分析方法及电子设备 - Google Patents

一种基于数字地形图的空间视域分析方法及电子设备 Download PDF

Info

Publication number
CN111681313B
CN111681313B CN202010418855.1A CN202010418855A CN111681313B CN 111681313 B CN111681313 B CN 111681313B CN 202010418855 A CN202010418855 A CN 202010418855A CN 111681313 B CN111681313 B CN 111681313B
Authority
CN
China
Prior art keywords
elevation
line segment
viewpoint
slope
line
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
CN202010418855.1A
Other languages
English (en)
Other versions
CN111681313A (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.)
Yunnan Xunzie Technology Co ltd
Original Assignee
Yunnan Xunzie Technology Co ltd
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 Yunnan Xunzie Technology Co ltd filed Critical Yunnan Xunzie Technology Co ltd
Priority to CN202010418855.1A priority Critical patent/CN111681313B/zh
Publication of CN111681313A publication Critical patent/CN111681313A/zh
Application granted granted Critical
Publication of CN111681313B publication Critical patent/CN111681313B/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
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models
    • 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
    • G09B29/003—Maps
    • G09B29/005—Map projections or methods associated specifically therewith

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Educational Technology (AREA)
  • Educational Administration (AREA)
  • Business, Economics & Management (AREA)
  • Remote Sensing (AREA)
  • Mathematical Physics (AREA)
  • Computer Graphics (AREA)
  • Image Generation (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明的实施例公开一种基于数字地形图的空间视域分析方法及电子设备，涉及机场导航台选址建设、城市规划、日照分析、地形分析、公共安全等领域。所述方法包括：在具有高程信息的目标数字地形图中，以待分析的目标点为视点，真北方向为基准方向，顺时针按照预定角度间隔作射线；在第每条射线方向上，以所述视点为起点，构造首尾相连的多条有向线段；对于每条射线方向上的每条有向线段，筛选出其所处范围内高程值大于上一有向线段的终点对应的最大虚拟高程值的高程线；计算每条有向线段与筛选出的高程线的交点，计算交点与所述视点连线的斜率，确定所述视点至有向线段的终点之间的最大遮挡视线斜率。本发明无需插值计算，能够提高计算效率。

Description

一种基于数字地形图的空间视域分析方法及电子设备
技术领域
本发明涉及机场导航台选址建设、城市规划、日照分析、地形分析、公共安全技术领域，尤其涉及一种基于数字地形图的空间视域分析方法及电子设备。
背景技术
随着计算机技术、网络技术的发展，为大规模、大范围的空间分析提供了可能。空间视域分析是空间分析的一种重要技术，在机场导航台选址建设、城市规划、日照分析、公共安全等领域极其重要。以机场建设为例，导航台的雷达发射信号受到高度超过导航台或天线的障碍物(地物和地貌)的影响，因此准确提供导航台遮蔽角对机场的设计工作至关重要；城市规划建设中，以人所在位置为视点，分析其可视范围作为建构筑物的规划依据；公共安全领域，以防控人员为视点，可分析其安全射击范围；此外，在防空领域，也需要对雷达等探测器进行雷达视域计算和分析。
空间视域分析技术，其实质是以视点为中心，计算每个方向的遮蔽角。遮蔽角为视点和视点所在水平面向上算起视线被地形地物遮挡的最大垂直张角。遮蔽角计算的本质是LOS算法，其核心是以视点为中心，向目标点引一条射线，判断该射线是否中途被遮挡，如果被遮挡则两点不通视，反之则两点通视。图1为现有的目标位置遮蔽角计算示意图，如图1中所示，O点为视点中心点位置，O’为视点在地面的投影点，O点某一视线方向有障碍物P1、P2、P3三点，则该视线方向与障碍物顶端的连线和地平线之间的夹角为遮蔽角，角值范围0～90°(一般只计算仰角)。图1中，OP3视线方向同地平线之间的夹角最大，一般将此最大夹角称为视点方向的最大遮蔽角。假设O点的海拔高程值为Ho，O点不同视线方向障碍物的最大仰角点位假设为Pj点，Pj点的海拔高程值为Hj，OPj两点水平距离为s，则不考虑地球曲率时O点在不同方向的遮蔽角计算如下式所示：
αi＝atan[(Hj-Ho)/s]   (1)
计算视点位置的遮蔽角，一般以视点O为中心，自真北方向点起沿360°方位，按照一定间隔(1°～3°)进行计算，间隔越小，计算出的可视域越准确。将全部方向障碍物最大仰角点连线，则可以得到视点O在一定区域内的可视区域范围，即遮蔽角覆盖范围图。
目前，遮蔽角的测量方法有实测法和DEM分析法。实测法中，障碍物的高度可采用全站仪悬高法、三角高程法、GNSS RTK等方式测量，障碍物的水平位置可采用全站仪极坐标法或GNSS RTK施测，外业工作量大，效率较低。DEM分析法中，地形的高程通过DEM网格的高程插值计算获得。
DEM分析法主要基于格网DEM数据，采用LOS(Light of Sight)方法进行通视分析。近年来，为应对不同数据源及提高计算效率，国内外学者提出了不同处理方法，如黄太山提出的GIS空间叠纸分析计算法，应申等人提出的双增量地形可视计算法，刘骝等人提出的基于LOS法改进的关键斜率法，张斌等人提出的基于多分辨率DEM的改进LOS法，以及剖面图差积法、JANUS LOS法、DYNTACS LOS法、ModSAF LOS法、Bresenham LOS等其他改进方法。这些方法均基于格网DEM计算遮蔽角，取得了很好的计算效果，但是，DEM分析法对于已有成果数据依赖较高，DEM分析法无法直接使用未经格网化的地形图数据来计算遮蔽角，插值过程比较费时，在一定程度上影响计算速度。
发明内容
有鉴于此，本发明实施例提供一种基于数字地形图的空间视域分析方法及电子设备，能够直接使用未经网格化的地形图数据计算遮蔽角，无需插值计算，能够有效提高空间视域分析的计算速度和计算效率。
第一方面，本发明实施例提供一种基于数字地形图的空间视域分析方法，包括：
在具有高程信息的目标数字地形图中，以待分析的目标点为视点，真北方向为基准方向，顺时针按照预定角度间隔作射线，得到m条射线；其中，m为大于等于3的正整数；
在第i条射线方向上，以所述视点为起点，构造首尾相连的n条有向线段Li,1、Li,2、…Li,n；其中，i＝1,2,…,m；
对于第i条射线方向上的第j条有向线段Li,j，筛选出其所处范围内高程值大于上一有向线段Li,j-1的终点对应的最大虚拟高程值的高程线；其中，j＝1,2,…,n；所述上一有向线段Li,j-1的终点对应的最大虚拟高程值由所述视点至所述上一有向线段Li,j-1的终点之间的最大遮挡视线斜率确定；
计算有向线段Li,j与筛选出的高程线的交点，确定所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率。
结合第一方面，在第一方面的第一种可实施方式中，在所述对于第i条射线方向上的第j条有向线段Li,j，筛选出其所处范围内高程值大于上一有向线段Li,j-1的终点对应的最大虚拟高程值的高程线之前，还包括：
将有向线段Li,j在预定宽度范围占所述目标数字地形图的区域作为有向线段Li,j对应的缓冲区；
其中，所述对于第i条射线方向上的第j条有向线段Li,j，筛选出其所处范围内高程值大于上一有向线段Li,j-1的终点对应的最大虚拟高程值的高程线，包括：
筛选出有向线段Li,j对应的缓冲区内高程值大于Hi,j-1的高程线；其中，Hi,j-1为有向线段Li,j-1的终点对应的最大虚拟高程值，Hi,0为视点的高程值。
结合第一方面的第一种可实施方式，在第一方面的第二种可实施方式中，所述计算有向线段Li,j与筛选出的高程线的交点，确定所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率，包括：
计算有向线段Li,j在其对应缓冲区内与筛选出的各高程线Gl的交点Ai,j,l；
计算各交点Ai,j,l至所述视点的斜率Ki,j,l；
依次判断所述各交点Ai,j,l至所述视点的斜率Ki,j,l是否大于第i条射线方向上的最大斜率Ki,max；
若交点Ai,j,l至所述视点的斜率Ki,j,l大于第i条射线方向上的最大斜率Ki,max，则令Ki,max＝Ki,j,l，并记录对应交点Ai,j,l的位置坐标；其中，Ki,max的初始值为0；
当遍历完全部交点后，将当前所述第i条射线方向上的最大斜率Ki,max确定为所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率；
根据所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率，计算有向线段Li,j的终点对应的最大虚拟高程值Hi,j。
结合第一方面的第一种可实施方式，在第一方面的第三种可实施方式中，所述筛选出有向线段Li,j对应的缓冲区内高程值大于Hi,j-1的高程线，包括：
筛选出有向线段Li,j对应的缓冲区内高程值大于Hi,j-1且属于预设的第一高程值集合的高程线；其中，所述第一高程值集合中的若干个高程值按照大小顺序排列后，任意两个相邻高程值的高程差等于第一预定高程差。
结合第一方面的第三种可实施方式，在第一方面的第四种可实施方式中，在第i条射线方向上，在确定出所述视点至最后一条有向线段Li,n的终点之间的最大遮挡视线斜率后，所述方法还包括：
获取所述第i条射线方向上的最大斜率Ki,max对应交点的高程值；
建立第二高程值集合；所述第二高程值集合中的若干个高程值按照大小顺序排列后，任意两个相邻高程值的高程差等于第二预定高程差，所述第二高程值集合中的最小一个高程值等于上一步骤获取的高程值与第二预定高程差之和，所述第二预定高程差小于所述第一预定高程差；
对于第i条射线方向上的第j条有向线段Li,j，筛选出有向线段Li,j对应的缓冲区内高程值大于Hi,j-1且属于所述第二高程值集合的高程线，随后返回执行所述计算有向线段Li,j与筛选出的高程线的交点，确定所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率的步骤。
结合第一方面，在第一方面的第五种可实施方式中，所述在第i条射线方向上，以所述视点为起点，构造首尾相连的n条有向线段Li,1、Li,2、...Li,n，包括：
在第i条射线方向上，以所述视点为起点，构造首尾相连且长度相等的n条有向线段Li,1、Li,2、...Li,n。
结合第一方面至第一方面的第五种可实施方式中的任一种可实施方式，在第一方面的第六种可实施方式中，在所述在具有高程信息的目标数字地形图中，以待分析的目标点为视点，真北方向为基准方向，顺时针按照预定角度间隔作射线，得到m条射线之前，还包括步骤：
对具有高程信息的初始数字地形图中的高程线进行抽稀，得到具有抽稀后的高程线的目标数字地形图。
结合第一方面的第六种可实施方式，在第一方面的第七种可实施方式中，所述对具有高程信息的初始数字地形图中的高程线进行抽稀，包括：
采用道格拉斯-普克Douglas-Peuker抽稀算法对具有高程信息的初始数字地形图中的高程线进行抽稀。
结合第一方面的第七种可实施方式，在第一方面的第八种可实施方式中，根据所述初始数字地形图等高距和比例尺，以及预定遮蔽角计算精度，确定采用道格拉斯-普克Douglas-Peuker抽稀算法时使用的阈值。
第二方面，本发明实施例提供一种电子设备，所述电子设备包括：壳体、处理器、存储器、电路板和电源电路，其中，电路板安置在壳体围成的空间内部，处理器和存储器设置在电路板上；电源电路，用于为上述电子设备的各个电路或器件供电；存储器用于存储可执行程序代码；处理器通过读取存储器中存储的可执行程序代码来运行与可执行程序代码对应的程序，用于执行前述任一实施方式所述的基于数字地形图的空间视域分析方法。
本发明实施例提供的基于数字地形图的空间视域分析方法及电子设备，在具有高程信息的目标数字地形图中，以待分析的目标点为视点，真北方向为基准方向，顺时针按照预定角度间隔作射线，得到m条射线，对于每条射线方向，通过以所述视点为起点，构造首尾相连的n条有向线段，在确定了视点到每段有向线段的终点之间的最大遮挡视线斜率后，即可计算出当该段有向线段的终点与视点之间的斜率等于该最大遮挡视线斜率时该有向线段的终点对应的最大虚拟高程值，在计算下一有向线段对应的区域时，只需要筛出高于所述下一有向线段的起点对应的最大虚拟高程值的高程线来与该有向线段求交点。该方法在从视点往外方向计算的过程中，根据已经计算出的交点和视点连线的斜率，确定其中的最大遮挡视线斜率，随后根据该最大遮挡视线斜率不断提高离视点更远区域高程值的筛选阈值，该方法不需要使用对现有的地形图数据进行插值计算格网，能够直接通过计算所作射线和高程线的交点，来确定最大遮蔽角，处理的数据对象不同于现有的DEM分析法，方法步骤也完全不同，能够显著提高视域分析过程的计算速度，提高计算效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。
图1为现有的目标位置遮蔽角计算示意图；
图2为本发明实施例一提供的一种基于数字地形图的空间视域分析方法流程示意图；
图3为通过第i条射线的一个竖直平面内的有向线段的终点对应的最大虚拟高程值的确定示意图；
图4为本发明实施例二提供的一种基于数字地形图的空间视域分析方法流程示意图；
图5为目标数字地形图中缓冲区域的划分示意图；
图6为步骤S205的优选实施方法流程图；
图7为采用图6所示方法确定视点至有向线段的终点之间的最大遮挡视线斜率和所述终点的最大虚拟高程值的原理图；
图8为本发明实施例二提供的另一种基于数字地形图的空间视域分析方法流程示意图；
图9为本发明提供的一种电子设备实施例的结构示意图。
具体实施方式
为解决现有的DEM分析法对于已有成果数据依赖较高，无法直接使用未经格网化的地形图数据来计算遮蔽角，插值过程比较费时，在一定程度上影响计算速度的问题，本发明提供一种新的基于具有等高线、建构筑物结构线等的数字地形图的空间视域分析算法，该方法基于LOS原理，但针对高程线的特点进行了有效利用，基本思想是按照视线方向分段最大斜率比较的方法，按照已经计算出的最大斜率快速筛选每段符合求交条件(已经计算出的最大斜率)的高程线，达到提高计算速率的目的。下面结合附图对本发明实施例进行详细描述。
应当明确，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
实施例一
图2为本发明实施例一提供的一种基于数字地形图的空间视域分析方法流程示意图。参看图2，本发明提供的方法可以包括如下步骤S101-S104：
S101、在具有高程信息的目标数字地形图中，以待分析的目标点为视点，真北方向为基准方向，顺时针按照预定角度间隔作射线，得到m条射线。
其中，m为大于等于3的正整数。本步骤中，在具有高程信息的目标数字地形图中，以视点为中心，向外做等距射线模仿视线，并在后续步骤中计算每个射线方向的高点与视点之间连线所成斜率。
S102、在第i条射线方向上，以所述视点为起点，构造首尾相连的n条有向线段Li,1、Li,2、...Li,n。
其中，i＝1,2,…,m。
在一优选实施例中，步骤S201中，在第i条射线方向上，以所述视点为起点，构造首尾相连且长度相等的n条有向线段Li,1、Li,2、…Li,n，即：有向线段Li,j的终点和有向线段Li,j+1的起点重合，且有向线段Li,j和有向线段Li,j+1共线，j＝1,2,…,n。
S103、对于第i条射线方向上的第j条有向线段Li,j，筛选出其所处范围内高程值大于上一有向线段Li,j-1的终点对应的最大虚拟高程值的高程线；
其中，j＝1,2,…,n；所述上一有向线段Li,j-1的终点对应的最大虚拟高程值由所述视点至所述上一有向线段Li,j-1的终点之间的最大遮挡视线斜率确定。具体地，当j＝1时，将视点的高程值Hi,0作为Li,1的上一有向线段的终点对应的最大虚拟高程值，直接筛选出Li,1所处范围内高程值大于视点的高程值Hi,0的高程线；当j≥2时，筛选出有向线段Li,j所处范围内高程值大于上一有向线段Li,j-1的终点对应的最大虚拟高程值Hi,j-1的高程线。
优选地，j≥2时，有向线段Li,j-1的终点对应的最大虚拟高程值Hi,j-1可根据下式(2)计算：
Hi,j-1＝ai,j-1×Ki,max   (2)
其中，Hi,0为视点的高程值，ai,j-1为目标数字地形图中有向线段Li,j-1的终点Pi,j-1与视点O之间的距离，Ki,max为视点O至有向线段Li,j-1的终点之间的最大遮挡视线斜率。值得说明的是，由于步骤S101是在目标数字地形图所处平面(地平面)中做的射线，因此，有向线段Li,j-1的终点Pi,j-1与视点O之间的距离ai,j-1实际上也可以理解为在地平面中，视点O和有向线段Li,j-1的终点在地平面上的投影之间的距离，若通过第i条射线做一竖直(垂直于地平面方向)平面，则本发明中所述的第i条射线可以理解为该竖直平面内以视点O在地平面内的投影O’为起点，以地平方向为延伸方向的一条射线。例如图3所示为通过第i条射线的一个竖直平面，其中O为视点，O’为视点在地平面内的投影点，线段O’O的长度为视点O的高程值Hi,0；图3中以Pi,1、Pi,2、…、Pi,j-1、Pi,j分别表示有向线段Li,1、Li,2、…、Li,j-1、Li,j的终点，则图4中线段O’Pi,1的长度即为有向线段Li,1的长度，线段Pi,1Pi,2的长度即为有向线段Li,2的长度，线段O’Pi,j-1的长度即为ai,j-1；图3中所示的高低起伏的曲线为第i条射线方向上的遮挡物顶端连线，在O’和有向线段Li,j-1的终点Pi,j-1这一段范围内的最高遮挡物顶端为Bi,max，线段OBi,max相对于地平面的斜率为Ki,max，将线段OBi,max延长至点Pi,j-1的正上方的点P’i,j-1，将Pi,j-1P’i,j-1的间距作为有向线段Li,j-1的终点Pi,j-1对应的最大虚拟高程值Hi,j-1。
S104、计算有向线段Li,j与筛选出的高程线的交点，确定所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率。
优选地，步骤S104中，在确定出视点至有向线段Li,j的终点之间的最大遮挡视线斜率时，还记录视点至有向线段Li,j的终点之间的最大遮挡视线斜率对应的交点位置信息，及计算出该斜率的有向线段和相应高程线的交点。
此步骤中，若上一步骤中有向线段Li,j所处范围内有高程值大于上一有向线段Li,j-1的终点对应的最大虚拟高程值的高程线，则此步骤中计算出有向线段Li,j与筛选出的高程线的交点，并计算各交点至所述视点的(对地)斜率，并将本次计算出的斜率中的最大值确定为所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率；若上一步骤中有向线段Li,j所处范围内没有高程值大于上一有向线段Li,j-1的终点对应的最大虚拟高程值的高程线，则，此步骤中将视点至上一有向线段Li,j-1的终点之间的最大遮挡视线斜率确定为所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率。例如图3中所示，由于有向线段Li,j所处范围内的障碍物的高程值均未大于点Pi,j-1对应的最大虚拟高程值Hi,j-1，则视点至有向线段Li,j的终点之间的最大遮挡视线斜率仍为线段OBi,max相对于水平面的斜率Ki,max。
本实施例提供的基于数字地形图的空间视域分析方法，在具有高程信息的目标数字地形图中，以待分析的目标点为视点，真北方向为基准方向，顺时针按照预定角度间隔作射线，得到m条射线，对于每条射线方向，通过以所述视点为起点，构造首尾相连的n条有向线段，在确定了视点到每段有向线段的终点之间的最大遮挡视线斜率后，即可计算出当该段有向线段的终点与视点之间的斜率等于该最大遮挡视线斜率时该有向线段的终点对应的最大虚拟高程值，在计算下一有向线段对应的区域时，只需要筛出高于所述下一有向线段的起点对应的最大虚拟高程值的高程线来与该有向线段求交点。该方法在从视点往外方向计算的过程中，根据已经计算出的交点和视点连线的斜率，确定其中的最大遮挡视线斜率，随后根据该最大遮挡视线斜率不断提高离视点更远区域高程值的筛选阈值，该方法不需要使用对现有的地形图数据进行插值计算格网，能够直接通过计算所作射线和高程线的交点，来确定最大遮蔽角，处理的数据对象不同于现有的DEM分析法，方法步骤也完全不同，能够显著提高视域分析过程的计算速度，提高计算效率。
实施例二
图4为本发明实施例二提供的一种基于数字地形图的空间视域分析方法流程示意图。参看图4，本实施例提供的方法可以包括如下步骤S201-S205：
S201、在具有高程信息的目标数字地形图中，以待分析的目标点为视点，真北方向为基准方向，顺时针按照预定角度间隔作射线，得到m条射线。
S202、在第i条射线方向上，以所述视点为起点，构造首尾相连的n条有向线段Li,1、Li,2、…Li,n；
其中，i＝1,2,…,m。
本实施例中，步骤S201-S202的具体实施方式类似于步骤S101-S102，此处不再赘述。
S203、将有向线段Li,j在预定宽度范围占所述目标数字地形图的区域作为有向线段Li,j对应的缓冲区。
本实施例中，若预定宽度为2W，优选地，有向线段Li,j对应的缓冲区为以该有向线段Li,j为对称轴，以2W为宽，以有向线段Li,j的长度为长的矩形区域。例如图5所示的目标数字地形图(图中未示出等高线)中，O点为视点，以真北方向为基准方向，顺时针按照90°间隔作射线，得到分别位于0°、90°、180°、270°方向的4条射线，若取n＝3，在每条射线方向上，构造首尾相连的3条有向线段，例如对于0°方向的第1条射线，构造了有向线段L1,1、L1,2、L1,3，对于L1,2段，图5中的矩形阴影区域为有向线段L1,2对应的缓冲区，图5中，该缓冲区的左侧边界位于0°射线左侧W距离处，该缓冲区的右侧边界位于0°射线右侧W距离处，即图中的线段BP1,2＝P1,2A＝W。例如，W可以根据目标地形图的范围、比例尺、计算精度等取值为5m，10m等。
S204、对于第i条射线方向上的第j条有向线段Li,j，筛选出有向线段Li,j对应的缓冲区内高程值大于上一有向线段Li,j-1的终点对应的最大虚拟高程值Hi,j-1的高程线。
其中，Hi,j-1为有向线段Li,j-1的终点对应的最大虚拟高程值，Hi,0为视点的高程值。
本实施例中，有向线段Li,j-1的终点对应的最大虚拟高程值Hi,j-1与步骤S103中所述计算原理相同，此处不再赘述。值得说明的是，通过设置缓冲区，即可通过空间计算，快速筛选出缓冲区内高程值大于Hi,j-1的高程线，相对于通过一条射线来和数字地形图中多条高程线求交确定交点位置后，再来筛选有向线段Li,j这一段上高程值大于Hi,j-1的高程线的方法，进一步提高了计算效率。
S205、计算有向线段Li,j与筛选出的高程线的交点，确定所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率。
在一优选实施方式中，如图6所示，步骤S205可以包括如下步骤S2050-S2059：
S2050、判断本次筛选出的高程线数量是否大于等于1；若是，则执行步骤S2051，否则，执行步骤S2058；
S2051：计算有向线段Li,j在其对应缓冲区内与筛选出的各高程线Gl的交点；
S2052、判断交点数是否大于0；若是，则执行步骤S2053，否则，执行步骤S2058；
S2053、计算第q个交点Ai,j,l至所述视点的斜率Ki,j,l；其中，q的初始值为1；
S2054、判断第q个交点Ai,j,l至所述视点的斜率Ki,j,l是否大于记录的第i条射线方向上的最大斜率Ki,max；若是，则执行步骤S2055，否则，执行步骤S2056；
其中，Ki,max的初始值为0；
S2055、令Ki,max＝Ki,j,l，并记录对应交点Ai,j,l的位置坐标；
S2056、判断q是否小于Q；若是，则执行步骤S2057，否则，执行步骤S2058；
其中，Q为有向线段Li,j在其对应缓冲区内与筛选出的各高程线的交点总数。
S2057、令q＝q+1，随后返回执行步骤S2053。
S2058、确定所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率等于当前记录的第i条射线方向上的最大斜率Ki,max。
本实施例中，若本轮计算出某个交点Ai,j,l至所述视点的斜率Ki,j,l大于记录的第i条射线方向上的最大斜率Ki,max，则将记录的第i条射线方向上的最大斜率Ki,max更新为刚计算出的更大的斜率值，直至遍历比较完所有交点至所述视点的斜率与当前记录的第i条射线方向上的最大斜率Ki,max的大小。
S2059、根据所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率，计算有向线段Li,j的终点对应的最大虚拟高程值Hi,j。
本实施例中，在第i条射线所处的竖直平面内，以所述视点为起点，以所述视点至有向线段Li,j的终点之间的最大遮挡视线斜率Ki,max为斜率值创建一条高程值射线，计算有向线段Li,j的终点沿竖直方向在当前创建的这条高程值射线上的投影，具体计算方法类似上述公式(2)，此处不再赘述。
下面采用一个具体的实施例，对图6所示方法实施例的实现原理进行详细说明。
图7所示图纸平面和图3所示平面类似，为通过第i条射线的一个竖直平面，其中O为视点，图7中，视点O在地平面上的投影为O’，视点处的高程值O O’＝Hi,0，Pi,1、Pi,2、Pi,3、Pi,4分别为有向线段Li,1、Li,2、Li,3、Li,4的终点，图7中所示的高低起伏的曲线为第i条射线方向上的遮挡物顶端连线，在视点O和有向线段Li,1的终点Pi,1之间的最高遮挡物顶端为点Bi,max1，根据本发明实施例提供的前述方法，通过OBi,max1连线计算斜率后可确定Pi,1的最大虚拟高程值为Pi,1P’i,1；随后，对于有向线段Li,2，筛选出其缓冲区内高程值大于Pi,1P’i,1的高程线并求交点，计算各交点与视点O的连线斜率，确定出计算出的斜率中交点Bi,max2与视点O的连线斜率大于OBi,max1斜率，因此将当前记录的第i条射线方向上的最大斜率Ki,max由OBi,max1的斜率更新为OBi,max2的斜率；随后，对于有向线段Li,3，由于其缓冲区内没有高程值大于Pi,2P’i,2的高程线，因此有向线段Li,3的终点的最大虚拟高程值仍采用OBi,max2的斜率来计算，得到有向线段Li,3的终点Pi,3的最大虚拟高程值Pi,3P’i,3；随后，对于有向线段Li,4，筛选出其缓冲区内高程值大于Pi,3P’i,3的高程线并求交点，计算各交点与视点O的连线斜率，确定出计算出的斜率中，交点Bi,max4与视点O的连线斜率为各交点中的最大，且OBi,max4的斜率大于OBi,max2的斜率，因此，将当前记录的第i条射线方向上的最大斜率Ki,max由OBi,max2的斜率更新为OBi,max4的斜率。
在一优选实施例中，上述基于数字地形图的空间时域分析方法中，在步骤S103/S204之前，所述方法还可包括建立第一高程值集合的步骤；所述第一高程值集合中的若干个高程值按照大小顺序排列后，任意两个相邻高程值的高程差等于第一预定高程差；本实施例中，若预设有第一高程值集合，则在步骤S103/S204中，对于第i条射线方向上的第j条有向线段Li,j，具体筛选出有向线段Li,j对应的缓冲区内高程值大于Hi,j-1且属于预设的第一高程值集合的高程线。例如：若图5中有向线段的L1,1的终点对应的最大虚拟高程值H1,1为110m，若不使用本实施例提供的方法，在步骤S103/S204中，筛选出有向线段L1,2对应的缓冲区(图5中的阴影区域)内高程值大于110m的高程线有230m，300m，320m共3根满足条件的高程线，但是，若预先设置的第一高程值集合为{100m,200m,300m,400m}，及第一预定高程差为100m，则在步骤S103/S204中筛选出有向线段L1,2对应的缓冲区内高程值大于110m且属于第一高程值集合的高程线有300m这一根，显然，本实施例提供的方法，可以根据视域分析时的分析精度来设置第一预定高程差，在不需要精度很高的情况下，通过增大第一预定高程差，可以进一步减少筛选出的每个缓冲区内的高程线的数量，从而减少下一步中有向线段和高程线求交点的次数，提高计算速度。
图8为本发明实施例二提供的另一种基于数字地形图的空间视域分析方法流程示意图。参看图8，本实施例提供的方法可以包括如下步骤S301-S307：
S301、在具有高程信息的目标数字地形图中，以待分析的目标点为视点，真北方向为基准方向，顺时针按照预定角度间隔作射线，得到m条射线。
S302、在第i条射线方向上，以所述视点为起点，构造首尾相连的n条有向线段Li,1、Li,2、…Li,n；
其中，i＝1,2,…,m。
S303、将有向线段Li,j在预定宽度范围占所述目标数字地形图的区域作为有向线段Li,j对应的缓冲区；
其中，j＝1,2,…,n。
本实施例中，步骤S301-S303的具体实施方式类似于步骤S201-S203，此处不再赘述。
S304、对于第i条射线方向上的第j条有向线段Li,j，筛选出有向线段Li,j对应的缓冲区内高程值大于上一有向线段Li,j-1的终点对应的最大虚拟高程值Hi,j-1且属于预设的第一高程值集合的高程线。
例如：若图5中有向线段的L1,1的终点对应的最大虚拟高程值H1,1为110m，若预先设置的第一高程值集合为{100m,200m,300m,400m}，本步骤中筛选出有向线段L1,2对应的缓冲区内高程值大于110m且属于第一高程值集合的高程线为300m这一根高程线，显然，通过此步骤，可以初步确定有向线段的L1,1对应的缓冲区内大于110m且跨距较大的几根高程线。
S305、计算有向线段Li,j与筛选出的高程线的交点，确定所述视点至最后一条有向线段Li,n的终点之间的最大遮挡视线斜率。
本实施例中，采用步骤S301-S305的方法，在第i条射线方向上，当遍历完j后，即确定出视点至最后一条有向线段Li,n的终点之间的最大遮挡视线斜率后，才执行下一步骤S306。
S306、获取第i条射线方向上的最大斜率Ki,max对应交点的高程值。
S307、建立第二高程值集合；
其中，所述第二高程值集合中的若干个高程值按照大小顺序排列后，任意两个相邻高程值的高程差等于第二预定高程差，所述第二高程值集合中的最小一个高程值等于上一步骤获取的高程值与第二预定高程差之和，且所述第二预定高程差小于所述第一预定高程差。优选地，所述第二高程值集合中的最大一个高程值小于上一步骤获取的高程值与第一预定高程差的整数倍之和，即：若上一步骤获取的高程值为x1，第一预定高程差为Δx，则第二高程值集合中的最大一个高程值x2max<x1+c×Δx，其中，c取值为正整数。例如：若步骤S306中获取到图5所示的第1条射线方向上的最大斜率Ki,max对应交点的高程值为x1＝300m，第一预定高程差为Δx＝100m，则可以设c＝1，第二预定高程差为20m，此步骤中第二高程值集合为{320m,340m,360m,380m}，若设c＝2，第二预定高程差为30m，此步骤中第二高程值集合为{330m,360m,390m,420m,450m,480m}。其中，c和第二预定高程差的取值根据实际需要的计算精度进行设定。
S308、对于第i条射线方向上的第j条有向线段Li,j，筛选出有向线段Li,j对应的缓冲区内高程值大于上一有向线段Li,j-1的终点对应的最大虚拟高程值Hi,j-1且属于所述第二高程值集合的高程线，随后返回执行步骤S305。
本实施例步骤S308中，j的取值又重新从j＝1开始，计算第i条射线方向上的第j条有向线段对应的缓冲区内高程值大于上一有向线段的终点对应的最大虚拟高程值且属于所述第二高程值集合的高程线与有向线段Li,j的交点，进一步筛选出每个有向线段缓冲区内与视点连线斜率可能大于之前步骤S305中确定的视点至最后一条有向线段Li,n的终点之间的最大遮挡视线斜率的交点，直至遍历完第第i条射线方向上的所有有向线段为止。
图8所示实施例中，通过设置第一高程值集合来对每个缓冲区内高于上一有向线段终点对应的最大虚拟高程值的高程线来进行首轮筛选，初略计算出每个射线方向上的最大遮挡视线斜率，在最大程度下去除了不必要的高程线与有向线段的求交计算，随后根据该最大遮挡视线斜率对应交点的高程值，通过设置第二高程值集合，重新使用本发明提供的方法计算，进一步提高计算精度。显然，在本实施例的基础上，根据计算精度的要求，还可以根据第二轮计算的结果，进一步设置第三高程值集合、第四高程值集合等进行类似的第三、第四轮计算，此处不再赘述。
优选地，c取值为1，即所述第二高程值集合中的最大一个高程值小于上一步骤获取的高程值与第一预定高程差之和。例如：若第一预定高程差为100m，在第一次执行步骤S306后获取到图5所示的第1条射线方向上的最大斜率Ki,max对应交点的高程值为300m，若第二预定高程差为20m，则第二高程值集合为{320m,340m,360m,380m}，若第二轮计算后，得到第1条射线方向上的最大斜率Ki,max对应交点的高程值为320m，还需要第三轮计算，则可以取第三预定高程差为5m，建立第三高程值集合为{325m,330m}后，对于第i条射线方向上的第j条有向线段Li,j，筛选出有向线段Li,j对应的缓冲区内高程值大于上一有向线段Li,j-1的终点对应的最大虚拟高程值Hi,j-1且属于所述第三高程值集合的高程线，随后返回执行步骤S305。
实施例三
本发明的上述实施例提供的分析方法中，需要构造射线同高程线进行求交运算，而地形图的等高线为使曲线显示美观，一般都进行了插值，有冗余节点，高程线的节点太多将直接影响求交效率。本实施例提供的基于数字地形图的空间视域分析方法，为减少高程线与射线求交计算的次数，在上述任一实施例的基础上，在所述在具有高程信息的目标数字地形图中，以待分析的目标点为视点，真北方向为基准方向，顺时针按照预定角度间隔作射线，得到m条射线之前，还包括：对具有高程信息的初始数字地形图中的高程线进行抽稀，得到具有抽稀后的高程线的目标数字地形图的步骤。对高程线进行抽稀，可以减少高程线的和所作射线的求交计算数量，进一步提高计算速度。
优选地，采用道格拉斯-普克Douglas-Peuker抽稀算法对具有高程信息的初始数字地形图中的高程线进行抽稀。在一可选实施例中，根据所述初始数字地形图等高距和比例尺，以及预定遮蔽角计算精度，确定采用道格拉斯-普克Douglas-Peuker抽稀算法时使用的阈值。
优选地，在设置了各有向线段对应的缓冲区后，可以对抽稀后的高程线进行分区索引。本实施例中，对高程线进行分区索引，就是对每个缓冲区进行第一索引编号，并对每根高程线按照其所在缓冲区和高程属性值(高程值)进行第二索引编号。这样，在求交计算时，首先根据第一索引编号查找缓冲分区，然后可以根据第二索引编号快速找到需要计算的高程线，显然，在求交计算的时候，直接根据相关索引编号找到需要计算的高程线，减少判断次数。例如：在所述目标数字地形图中，以视点为圆心，真北方向为基准方向，顺时针按照一定的间隔(1°)作射线，然后以该射线一定宽度(根据精度要求调整，如5m或10m)范围做缓冲区。若射线间隔为1°，则有360个缓冲区。为进一步缩小搜索范围，对每个射线方向的缓冲区再划分为n段(对应于将每条射线从圆心出发划分为n段有向线段)，则得到360*n个缓冲区，同建立栅格金字塔类似，对每个缓冲区内的高程线的高程值建立索引，便于快速地选取特定高程值的高程线，用于后续步骤中的求交计算。
本发明实施例还提供一种电子设备。图9为本发明电子设备一个实施例的结构示意图，可以实现本发明图2或图4所示实施例的流程，如图9所示，上述电子设备可以包括：壳体11、处理器12、存储器13、电路板14和电源电路15，其中，电路板14安置在壳体11围成的空间内部，处理器12和存储器13设置在电路板14上；电源电路15，用于为上述电子设备的各个电路或器件供电；存储器13用于存储可执行程序代码；处理器12通过读取存储器13中存储的可执行程序代码来运行与可执行程序代码对应的程序，用于执行前述任一实施例所述的基于数字地形图的空间视域分析方法。
该电子设备以多种形式存在，包括但不限于：
(1)移动通信设备：这类设备的特点是具备移动通信功能，并且以提供话音、数据通信为主要目标。这类终端包括：智能手机(例如iPhone)、多媒体手机、功能性手机，以及低端手机等。
(2)超移动个人计算机设备：这类设备属于个人计算机的范畴，有计算和处理功能，一般也具备移动上网特性。这类终端包括：PDA、MID和UMPC设备等，例如iPad。
(3)便携式娱乐设备：这类设备可以显示和播放多媒体内容。该类设备包括：音频、视频播放模块(例如iPod)，掌上游戏机，电子书，以及智能玩具和便携式车载导航设备。
(4)服务器：提供计算服务的设备，服务器的构成包括处理器、硬盘、内存、系统总线等，服务器和通用的计算机架构类似，但是由于需要提供高可靠的服务，因此在处理能力、稳定性、可靠性、安全性、可扩展性、可管理性等方面要求较高。
(5)其他具有数据交互功能的电子设备。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中的各个实施例均采用相关的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
为了描述的方便，描述以上装置是以功能分为各种单元/模块分别描述。当然，在实施本发明时可以把各单元/模块的功能在同一个或多个软件和/或硬件中实现。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(Random AccessMemory，RAM)等。
以上所述，仅为本发明的具体实施方式，但本发明的防护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的防护范围之内。因此，本发明的防护范围应以权利要求的防护范围为准。

Claims (8)

1.一种基于数字地形图的空间视域分析方法，其特征在于，包括：
在具有高程信息的目标数字地形图中，以待分析的目标点为视点，真北方向为基准方向，顺时针按照预定角度间隔作射线，得到m条射线；其中，m为大于等于3的正整数；
在第i条射线方向上，以所述视点为起点，构造首尾相连的n条有向线段L i,1、 L i,2、…L i,n ；其中，i=1,2,…,m；
对于第i条射线方向上的第j条有向线段L i,j ，筛选出其所处范围内高程值大于上一有向线段L i,j-1的终点对应的最大虚拟高程值的高程线；其中，j=1,2,…,n；所述上一有向线段L i,j-1的终点对应的最大虚拟高程值由所述视点至所述上一有向线段L i,j-1的终点之间的最大遮挡视线斜率确定；
计算有向线段L i,j 与筛选出的高程线的交点，确定所述视点至有向线段L i,j 的终点之间的最大遮挡视线斜率；
在所述对于第i条射线方向上的第j条有向线段L i,j ，筛选出其所处范围内高程值大于上一有向线段L i,j-1的终点对应的最大虚拟高程值的高程线之前，还包括：
将有向线段L i,j 在预定宽度范围占所述目标数字地形图的区域作为有向线段L i,j 对应的缓冲区；
其中，所述对于第i条射线方向上的第j条有向线段L i,j ，筛选出其所处范围内高程值大于上一有向线段L i,j-1的终点对应的最大虚拟高程值的高程线，包括：
筛选出有向线段L i,j 对应的缓冲区内高程值大于H i,j-1的高程线；其中，H i,j-1为有向线段L i,j-1的终点对应的最大虚拟高程值，H i,0为视点的高程值；
所述计算有向线段L i,j 与筛选出的高程线的交点，确定所述视点至有向线段L i,j 的终点之间的最大遮挡视线斜率，包括：
计算有向线段L i,j 在其对应缓冲区内与筛选出的各高程线G l 的交点A i,j,l ；
计算各交点A i,j,l 至所述视点的斜率K i,j,l ；
依次判断所述各交点A i,j,l 至所述视点的斜率K i,j,l 是否大于第i条射线方向上的最大斜率K i,max ；
若交点A i,j,l 至所述视点的斜率K i,j,l 大于第i条射线方向上的最大斜率K i,max ，则令K i,max =K i,j,l ，并记录对应交点A i,j,l 的位置坐标；其中，K i,max 的初始值为0；
当遍历完全部交点后，将当前所述第i条射线方向上的最大斜率K i,max 确定为所述视点至有向线段L i,j 的终点之间的最大遮挡视线斜率；
根据所述视点至有向线段L i,j 的终点之间的最大遮挡视线斜率，计算有向线段L i,j 的终点对应的最大虚拟高程值H i,j 。
2.根据权利要求1所述的基于数字地形图的空间视域分析方法，其特征在于，所述筛选出有向线段L i,j 对应的缓冲区内高程值大于H i,j-1的高程线，包括：
筛选出有向线段L i,j 对应的缓冲区内高程值大于H i,j-1且属于预设的第一高程值集合的高程线；其中，所述第一高程值集合中的若干个高程值按照大小顺序排列后，任意两个相邻高程值的高程差等于第一预定高程差。
3.根据权利要求2所述的基于数字地形图的空间视域分析方法，其特征在于，在第i条射线方向上，在确定出所述视点至最后一条有向线段L i,n 的终点之间的最大遮挡视线斜率后，所述方法还包括：
获取所述第i条射线方向上的最大斜率K i,max 对应交点的高程值；
建立第二高程值集合；所述第二高程值集合中的若干个高程值按照大小顺序排列后，任意两个相邻高程值的高程差等于第二预定高程差，所述第二高程值集合中的最小一个高程值等于上一步骤获取的高程值与第二预定高程差之和，所述第二预定高程差小于所述第一预定高程差；
对于第i条射线方向上的第j条有向线段L i,j ，筛选出有向线段L i,j 对应的缓冲区内高程值大于H i,j-1且属于所述第二高程值集合的高程线，随后返回执行所述计算有向线段L i,j 与筛选出的高程线的交点，确定所述视点至有向线段L i,j 的终点之间的最大遮挡视线斜率的步骤。
4.根据权利要求1所述的基于数字地形图的空间视域分析方法，其特征在于，所述在第i条射线方向上，以所述视点为起点，构造首尾相连的n条有向线段L i,1、 L i,2、…L i,n ，包括：
在第i条射线方向上，以所述视点为起点，构造首尾相连且长度相等的n条有向线段L i,1、 L i,2、…L i,n 。
5.根据权利要求1-4任一项所述的基于数字地形图的空间视域分析方法，其特征在于，在所述在具有高程信息的目标数字地形图中，以待分析的目标点为视点，真北方向为基准方向，顺时针按照预定角度间隔作射线，得到m条射线之前，还包括步骤：
对具有高程信息的初始数字地形图中的高程线进行抽稀，得到具有抽稀后的高程线的目标数字地形图。
6.根据权利要求5所述的基于数字地形图的空间视域分析方法，其特征在于，所述对具有高程信息的初始数字地形图中的高程线进行抽稀，包括：
采用道格拉斯-普克 Douglas-Peuker抽稀算法对具有高程信息的初始数字地形图中的高程线进行抽稀。
7. 根据权利要求6所述的基于数字地形图的空间视域分析方法，其特征在于，根据所述初始数字地形图等高距和比例尺，以及预定遮蔽角计算精度，确定采用道格拉斯-普克Douglas-Peuker抽稀算法时使用的阈值。
8.一种电子设备，其特征在于，所述电子设备包括：壳体、处理器、存储器、电路板和电源电路，其中，电路板安置在壳体围成的空间内部，处理器和存储器设置在电路板上；电源电路，用于为上述电子设备的各个电路或器件供电；存储器用于存储可执行程序代码；处理器通过读取存储器中存储的可执行程序代码来运行与可执行程序代码对应的程序，用于执行前述任一权利要求1-7所述的基于数字地形图的空间视域分析方法。
CN202010418855.1A 2020-05-18 2020-05-18 一种基于数字地形图的空间视域分析方法及电子设备 Active CN111681313B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010418855.1A CN111681313B (zh) 2020-05-18 2020-05-18 一种基于数字地形图的空间视域分析方法及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010418855.1A CN111681313B (zh) 2020-05-18 2020-05-18 一种基于数字地形图的空间视域分析方法及电子设备

Publications (2)

Publication Number Publication Date
CN111681313A CN111681313A (zh) 2020-09-18
CN111681313B true CN111681313B (zh) 2023-04-28

Family

ID=72452031

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010418855.1A Active CN111681313B (zh) 2020-05-18 2020-05-18 一种基于数字地形图的空间视域分析方法及电子设备

Country Status (1)

Country Link
CN (1) CN111681313B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112435337B (zh) * 2020-11-13 2024-09-10 郑亮 一种景观可视域分析方法及系统
CN114494598B (zh) * 2022-01-25 2023-03-21 南京师范大学 一种城市三维可视空间比指标的优化方法
CN114912268B (zh) * 2022-05-09 2023-02-10 中电普信(北京)科技发展有限公司 一种计算机仿真中实时快速通视判断方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102074050A (zh) * 2011-03-01 2011-05-25 哈尔滨工程大学 大规模地形绘制的分形多分辨率简化方法
CN102693328A (zh) * 2011-03-24 2012-09-26 黄丹丹 Arcview gis 3.2技术在山地建筑设计上的应用发明
CN103136393A (zh) * 2011-11-28 2013-06-05 中国电子科技集团公司第五十四研究所 一种基于网格划分的区域覆盖率计算方法
CN103292784A (zh) * 2013-05-14 2013-09-11 哈尔滨工程大学 一种基于声线轨迹的水下地表地形可视性分析方法
CN105740782A (zh) * 2016-01-25 2016-07-06 北京航空航天大学 一种基于单目视觉的驾驶员换道过程量化方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AUPR301601A0 (en) * 2001-02-09 2001-03-08 Commonwealth Scientific And Industrial Research Organisation Lidar system and method
US8368686B2 (en) * 2004-05-26 2013-02-05 Sony Online Entertainment Llc Resource management for rule-based procedural terrain generation
GB2526866A (en) * 2014-06-05 2015-12-09 Univ Bristol Apparatus for and method of inspecting surface topography of a moving object
CN105260523B (zh) * 2015-09-25 2018-06-26 中国科学院南京土壤研究所 一种分布式并行空间可视域分析方法
CN105893675B (zh) * 2016-03-31 2019-01-29 东南大学 一种基于天空可视域评价的开放空间周边建筑形态优化控制方法
CN107036544B (zh) * 2017-05-24 2019-12-06 山东科技大学 一种基于三维激光扫描技术的大型储罐安全检测系统与方法
CN107886570B (zh) * 2017-09-22 2020-03-24 中国矿业大学 一种兼顾精度与效率的可视域计算方法
CN107862743B (zh) * 2017-09-22 2020-08-11 中国矿业大学 一种基于投影的规则网格地形可视域计算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102074050A (zh) * 2011-03-01 2011-05-25 哈尔滨工程大学 大规模地形绘制的分形多分辨率简化方法
CN102693328A (zh) * 2011-03-24 2012-09-26 黄丹丹 Arcview gis 3.2技术在山地建筑设计上的应用发明
CN103136393A (zh) * 2011-11-28 2013-06-05 中国电子科技集团公司第五十四研究所 一种基于网格划分的区域覆盖率计算方法
CN103292784A (zh) * 2013-05-14 2013-09-11 哈尔滨工程大学 一种基于声线轨迹的水下地表地形可视性分析方法
CN105740782A (zh) * 2016-01-25 2016-07-06 北京航空航天大学 一种基于单目视觉的驾驶员换道过程量化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王栋 ; 刘林 ; 张顺发 ; 张云雷 ; .基于SRTM高程数据的变步长海面快速可视域分析方法.指挥控制与仿真.2017,(第04期),全文. *

Also Published As

Publication number Publication date
CN111681313A (zh) 2020-09-18

Similar Documents

Publication Publication Date Title
CN111681313B (zh) 一种基于数字地形图的空间视域分析方法及电子设备
CN111505662B (zh) 一种无人驾驶车辆定位方法及系统
AU2013237642A1 (en) Multiple Viewshed Analysis
CN112305559A (zh) 基于地面定点激光雷达扫描的输电线距离测量方法、装置、系统和电子设备
CN103292784B (zh) 一种基于声线轨迹的水下地表地形可视性分析方法
CN103596263B (zh) 一种室内定位中获取位置的方法及装置
CN107071721B (zh) 基于手机定位数据的停留点提取方法
CN111121785A (zh) 一种基于图搜索的无道路路径规划方法
CN113742437B (zh) 地图更新方法、装置、电子设备和存储介质
Miky et al. A combined contour lines iteration algorithm and Delaunay triangulation for terrain modeling enhancement
CN112002012B (zh) 城市区域的可视性分析方法
CN113936109A (zh) 高精地图点云数据的处理方法、装置、设备以及存储介质
CN108680898A (zh) 室内定位方法、装置、介质及电子设备
CN108268821A (zh) 一种室内外场景识别方法及装置
CN116704024A (zh) 一种联合多类约束的车载激光点云位姿图优化方法及系统
CN116772846A (zh) 一种无人机航迹规划方法、装置、设备及介质
CN116437292A (zh) 一种室内复杂环境下的行人定位匹配方法
Guan et al. Fast approximate viewshed analysis based on the regular-grid digital elevation model: X-type partition proximity-direction-elevation spatial reference line algorithm
CN115147482A (zh) 一种位姿初始化方法
Yang et al. The research on real-time map-matching algorithm
Feld et al. Approximated environment features with application to trajectory annotation
Hao-Nguyen et al. A new algorithm for viewshed computation on raster terrain
CN110602635A (zh) 一种室内地图匹配增强定位方法、设备及存储设备
CN111709432B (zh) 复杂城市环境下InSAR地面点提取方法、装置、服务器及存储介质
CN114003848B (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