WO2016095682A1 - 一种处理数字地面模型接缝的方法及装置 - Google Patents
一种处理数字地面模型接缝的方法及装置 Download PDFInfo
- Publication number
- WO2016095682A1 WO2016095682A1 PCT/CN2015/095642 CN2015095642W WO2016095682A1 WO 2016095682 A1 WO2016095682 A1 WO 2016095682A1 CN 2015095642 W CN2015095642 W CN 2015095642W WO 2016095682 A1 WO2016095682 A1 WO 2016095682A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- adjacent
- boundary
- boundary points
- horizontal plane
- triangle
- Prior art date
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
Definitions
- the present invention relates to the field of digital terrestrial models, and more particularly to a method and apparatus for processing digital ground model joints.
- the Digital Terrain Model is a very important part of navigation maps and is increasingly being used in a variety of navigation applications. Limited by the DTM data acquisition method, the DTM is typically stored using a pixmap plus a profile, where the pixmap is used to store height information for each pixel and the profile is used to store planar location information for the pixel.
- the DTM is rendered in the form of a triangulation.
- a triangulation is a polygonal mesh composed of multiple triangles.
- Triangulation is widely used in graphics and modeling to simulate the surface of complex objects. Triangulation can be subdivided into regular triangulation and irregular triangulation.
- the triangle in the regular triangulation mesh may be a right-angled triangle with equal side lengths.
- the irregular triangulation is formed by a discrete set of points after thinning and triangulation by a regular triangulation.
- the irregular triangulation can reduce the storage space of the digital ground model and greatly improve the rendering efficiency of the navigation application.
- the irregular triangulation obtained by each pixmap may not be the same, so each irregular triangulation There will be gaps in the edges. For example, there is a gap between two adjacent triangular tiles til1 and tile2 as shown in FIG.
- an object of the present invention is to provide a method and a device for processing a digital ground model joint. Set so that the seams between the triangles are tight and coherent without reducing the overall rendering efficiency.
- a method of processing a digital ground model seam may include: passing a boundary point on an adjacent boundary of two adjacent triangular nets in a digital ground model The coordinates on the horizontal plane are compared to distinguish the boundary points at which the two adjacent triangular meshes project on the horizontal plane, and the non-coincident boundary points are projected on the horizontal plane; the boundary points at which the projections coincide with each other on the horizontal plane are The height coordinates in the axial direction of the vertical horizontal plane are updated to the same value; the boundary points on which the projections on the adjacent borders of the two adjacent triangular nets are not coincident are respectively added to the set of boundary points of the opposite party; based on the added boundary points , forming a new triangle in the triangle mesh.
- an apparatus for processing a digital ground model seam may include: a comparing unit for contiguous by two adjacent triangular nets in a digital ground model Comparing the coordinates of the boundary points on the boundary on the horizontal plane, distinguishing the boundary points at which the two adjacent triangular meshes project on the horizontal plane, and projecting the non-coincident boundary points on the horizontal plane; the coincidence boundary point processing unit, For updating the height coordinates of the boundary points which are coincident with each other on the horizontal plane in the vertical horizontal axis to the same value; the non-coincident boundary point processing unit for placing the adjacent boundaries of the two adjacent triangular nets The boundary points of the projections that are not coincident are respectively added to the set of boundary points of the opposite side; the new triangular division unit is used to form a new triangle in the triangulation based on the added boundary points.
- the height coordinate of the boundary point of the projection coincidence on the adjacent boundary of two adjacent triangular nets in the vertical horizontal plane is updated to the same value, the height difference of the same position between the two triangular mesh boundaries is caused by the embodiment of the present invention.
- the gap can be eliminated, since the boundary points that do not coincide with the projections are complementary in the adjacent boundaries of the two adjacent triangulations, a new triangle is formed in the triangulation net for the newly added boundary points in the triangulation, thus The gap between the two triangular mesh boundaries can be filled by forming a new triangle in the triangulation for the newly added boundary points, so that the adjacent boundaries of the two adjacent triangular meshes are coherent and the gaps are filled. More natural. As shown in FIG.
- two adjacent triangular meshes having slits are used to fill the gaps by applying the method provided by the embodiment of the present invention, and the gap between two adjacent triangular meshes is eliminated, and the two sides of the adjacent boundary are continuously connected. Natural, closer to the real world terrain. Moreover, the formation of a new triangle for newly added boundary points does not add a large area of the plane, and does not reduce the overall rendering efficiency.
- Figure 1 is a simplified schematic view of the gap between two adjacent triangular meshes
- FIG. 2 is a schematic flow chart of a method for processing a digital ground model joint according to an embodiment of the present invention
- FIG. 3 is a simplified schematic diagram of two adjacent triangulation networks after applying a method according to an embodiment of the present invention
- FIG. 4 is a second schematic flowchart of a method for processing a digital ground model joint according to an embodiment of the present invention
- Figure 5 is a three-dimensional schematic view of a gap between two adjacent triangular meshes
- FIG. 6 is a three-dimensional schematic diagram of two adjacent triangular meshes after the slot is processed by applying the method provided by the embodiment of the present invention
- FIG. 7 is a schematic structural diagram of an apparatus for processing a digital ground model joint according to an embodiment of the present invention.
- FIG. 8 is a second schematic structural diagram of an apparatus for processing a digital ground model joint according to an embodiment of the present invention.
- FIG. 2 a flow chart of a method for processing a digital ground model seam is provided by the present invention. As shown in FIG. 2, the method may include:
- the digital terrestrial model may be configured with a configuration file for storing horizontal position information of each pixel.
- the horizontal position information can be represented by geographical latitude and longitude coordinates. It is assumed that (x, y, z) represents the three-dimensional coordinates of a pixel point in the embodiment of the present invention, wherein (x, y) is used to represent the coordinates of the pixel point on the horizontal plane (such as geographical latitude and longitude coordinates), z Indicates the height of the pixel point in the axial direction of the vertical horizontal plane.
- the coordinates (x, y) of the boundary points on the adjacent boundary of any two adjacent triangulations on the horizontal plane can be read from the configuration file of the digital ground model, and then the phases of any two adjacent triangulations are The boundary points on the adjacent boundary are compared in coordinates (x, y) on the horizontal plane.
- the boundary points corresponding to the adjacent boundaries adjacent to the triangulation tile 2 are ⁇ p1, p2, p3 ⁇ .
- the boundary point set corresponding to the adjacent boundary adjacent to the triangulation tile 1 is ⁇ p4, p5 ⁇ .
- the x and y coordinates of the five boundary points p1, p2, p3, p4, and p5 can be read from the configuration file of the digital ground model.
- the criterion for judging whether the boundary points coincide on the horizontal plane can be set according to actual needs. For example, if the gap filling accuracy is required to be high, the criterion for determining the boundary points of the coincidence projections may be set to be exactly the same on the horizontal plane. If the gap filling accuracy is not high, a certain error may be allowed, and the horizontal plane may be allowed. The two boundary points whose coordinate difference is smaller than the error are determined as the boundary points of the projection coincidence, so that the number of projection non-coincident boundary points can be reduced to some extent in the subsequent steps, the correlation calculation of the re-divided triangle is reduced, and the efficiency of filling the gap is improved. .
- step S220 and the step S230 there is no specific execution order between the step S220 and the step S230, and may be performed in parallel or sequentially.
- the adjacent boundary between the triangulations and the set of boundary points corresponding to the adjacent boundaries may be determined according to the data in the digital terrestrial model for describing the triangular structure and the triangular mesh structure composed of the triangles. Therefore, boundary points whose projections do not coincide can be added to the boundary point set corresponding to the adjacent boundary of the opposite triangle network.
- the two adjacent triangulations are respectively a first triangulation mesh and a second triangulation mesh
- all boundary points on the adjacent boundary of the first triangulation mesh that do not coincide with the adjacent boundary of the second triangulation mesh may be added to
- all boundary points on the adjacent boundary of the second triangular network that do not coincide with the first triangular mesh are added to the boundary point set corresponding to the adjacent boundary of the first triangular network.
- p2 is a boundary point where the adjacent boundaries of tile1 and tile2 are not projected on the horizontal plane, and p2 may be added to the set of boundary points corresponding to the adjacent boundary of tile2. After adding p2, the boundary points corresponding to the adjacent boundaries of tile2 are ⁇ p4, p2, p5 ⁇ .
- a triangle mesh is divided into triangles, and each boundary point in the triangle mesh falls into a triangle.
- the division of triangles in the triangulation can be set according to actual needs.
- the present invention is not limited to the specific implementation of how the newly added boundary points form a new triangle in the triangulation.
- FIG. 4 it is a schematic flowchart of a method for processing a digital ground model seam provided by the present invention. As shown in FIG. 4, the method may include:
- the boundary points on the adjacent boundaries of two adjacent triangular meshes may be pre-ordered according to the coordinates on the same latitude of the horizontal plane, such as the y-axis coordinates, and the coordinates of the sorted boundary points are performed. Comparison.
- S420 Calculate an average value of height coordinates of the boundary points of the mutually coincident boundary points on the horizontal plane in the vertical horizontal plane; and update the height coordinates of the boundary points which overlap each other on the horizontal plane in the vertical horizontal plane to the average value. It should be noted that, in order to project the updated heights of the boundary points that overlap each other on the horizontal plane, other calculation methods may be adopted in the embodiment of the present invention, and calculating the average value is only one possible implementation manner of the present invention.
- p1 and p4 are boundary points having the same three-dimensional position
- p3 and p5 are boundary points having the same three-dimensional position.
- the boundary points on the adjacent boundaries of two adjacent triangular meshes may be pre-ordered according to the coordinates on the same latitude of the horizontal plane, such as the y-axis coordinates. Therefore, based on the sort result Quickly find out which two existing boundary points the newly added boundary points fall into. Query the triangle where the two existing boundary points fall, and connect the other vertex of the triangle to the newly added boundary point to divide the triangle into new triangles.
- a triangle associated with the original boundary points p4, p5 of tile 2, that is, p4, p5 may be found.
- FIG. 5 and FIG. 6 are three-dimensional schematic views of two adjacent triangular nets, respectively. After the gap between two adjacent triangulations as shown in FIG. 5 is processed by the method provided by the embodiment of the present invention, two adjacent triangulations after filling the gap shown in FIG. 6 can be obtained.
- the embodiment of the invention can be used to process the gaps that occur when multiple triangle nets are integrated, and can also be used to process the gaps that occur when the pixel map of the grid type digital ground model is triangulated.
- two adjacent triangulation networks in the embodiment of the present invention may be each two adjacent irregular triangulation networks in all the irregular triangulation networks in the digital terrestrial model, or the embodiment of the present invention
- the two adjacent triangulations may be every two adjacent irregular triangulations in the partially irregular triangulation network determined according to the integration requirement in the digital terrestrial model.
- the gap between every two adjacent irregular triangulations in all irregular triangulations in the digital terrestrial model is processed, or each two adjacent ones in the partial irregular triangulation network are processed.
- the irregular triangulation it is possible to first determine the combination of all the irregular triangulations between all the irregular triangulations in the digital terrestrial model, or between some irregular triangulations. For each combination, the gap between two adjacent irregular triangulations in the group is processed using the method provided by the embodiment of the present invention.
- the triangulation may be after the thinning in the digital ground model.
- Regular triangle network may further include the step of obtaining the irregular triangulation by converting the digital terrestrial model into a regular triangulation, thinning and triangulating the regular triangulation.
- an embodiment of the present invention also provides an apparatus for processing a digital ground model seam.
- FIG. 7 is a schematic structural diagram of an apparatus for processing a digital ground model joint according to an embodiment of the present invention.
- the apparatus can include:
- the comparing unit 710 can be configured to distinguish the coordinates of the two adjacent triangular meshes on the horizontal plane by comparing the coordinates of the boundary points on the adjacent edges of two adjacent triangular meshes in the digital ground model. Boundary points, as well as projections of non-coincident boundary points on a horizontal plane.
- the coincidence boundary point processing unit 720 can be used to update the height coordinates of the boundary points that overlap each other on the horizontal plane in the vertical horizontal axis to the same value.
- the non-coincident boundary point processing unit 730 may be configured to add boundary points of projections on adjacent edges of the two adjacent triangular nets to each other's boundary point sets.
- a new triangle dividing unit 740 can be used to form a new triangle within the triangulation based on the added boundary points.
- the apparatus provided by the embodiment of the present invention is configured to update the height coordinates of the boundary points of the projected coincidence points on the adjacent boundaries of two adjacent triangular nets in the vertical horizontal plane to the same value, so that the coincident boundary point processing unit 720
- the gap caused by the height difference of the same position between the triangle mesh boundaries can be eliminated, because the non-coincident boundary point processing unit 730 complements the boundary points where the projections do not coincide in the adjacent boundaries of the two adjacent triangular nets, and new
- the triangle dividing unit 740 forms a new triangle in the triangular mesh for the newly added boundary point in the triangulation network, and therefore, the gap on the horizontal plane between the two triangular mesh boundaries can form a new boundary in the triangulation network for the newly added boundary point.
- Triangle to fill It can be seen that the embodiment of the present invention fills the gap naturally, and forming a new triangle for the newly added boundary point does not additionally increase a large area plane, and does not reduce the overall drawing efficiency.
- the coincident boundary point processing unit 720 may include: an average height calculation sub-unit 721, which may be used to calculate a boundary point that projects on the horizontal plane to coincide with each other in the vertical horizontal axis. The average of the height coordinates on the top.
- the height update sub-unit 722 can be used to project the height coordinates of the boundary points coincident with each other on the horizontal plane in the vertical horizontal axis. New is the average.
- the new triangle dividing unit 740 may include: a query triangle sub-unit 741, which may be used to query the newly added in a triangnet to which a boundary point is newly added.
- the connecting sub-unit 742 can be configured to connect the newly added boundary point to a diagonal vertice of the side where the triangle is located, and connect the newly added boundary point to divide the triangle into a new triangle.
- the device provided by the embodiment of the invention can be used to process the gaps generated when multiple triangulation networks are integrated, and can also be used to process the gaps appearing when the pixel map of the grid digital ground model is triangulated.
- two adjacent triangulation networks in the embodiment of the present invention may be each two adjacent irregular triangulation networks in all the irregular triangulation networks in the digital terrestrial model, or the embodiment of the present invention
- the two adjacent triangulations may be every two adjacent irregular triangulations in the partially irregular triangulation network determined according to the integration requirement in the digital terrestrial model.
- the apparatus provided by the embodiment of the present invention may further obtain the irregular triangulation network by converting the digital terrestrial model into a regular triangulation network, thinning and triangulating the regular triangulation network.
- the updated height of the boundary points of the projected coincidence can be obtained by a simple and efficient average value calculation method, and the associated triangles in the triangular network added by the newly added boundary points are divided into new ones.
- the triangle so that there is no need to perform global triangulation again, the efficiency of filling the gap between the triangulation nets is high, and the seams are tight, without losing the appearance of the real terrain.
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
Abstract
本申请公开了一种处理数字地面模型接缝的方法和装置，该方法可以包括：通过对数字地面模型中两个相邻三角网相邻边界上边界点在水平面上的坐标进行比较，区分出两个相邻三角网在水平面上投影相互重合的边界点，以及，在水平面上投影不重合的边界点；将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更新为相同的值；将两个相邻三角网的相邻边界上的投影不重合的边界点分别添加到对方的边界点集合中；基于添加的边界点，在三角网内形成新的三角形。使相邻三角网相邻边界两边地形连贯、自然，且不会额外增加大面积的平面，不会降低总体绘制效率。
Description
本申请要求于2014年12月16日提交中国专利局、申请号为201410785247.9、发明名称为“一种处理数字地面模型接缝的方法及装置”的中国专利申请的优先权，其全部内容通过引用结合在本申请中。
本发明涉及数字地面模型领域，尤其涉及一种处理数字地面模型接缝的方法及装置。
数字地面模型(DTM，Digital Terrain Model)是导航地图非常重要的组成部分，越来越多地被用于各种导航应用中。受限于DTM数据的采集方法，DTM通常使用像素图加配置文件的方式存储，其中，像素图用于存储每个像素点的高度信息，配置文件用于存储像素点的平面位置信息。
在导航应用中，DTM以三角网的形式渲染。一个DTM中可以有多个三角网，一个三角网是由多个三角形组成的多边形网格。在数字地面模型中，可以按照比例尺的需要，划分各个三角网所覆盖的范围。三角网在图形学和建模中广泛使用，用来模拟复杂物体的表面。三角网可细分为规则三角网和不规则三角网。其中，规则三角网中的三角形可以为边长相等的直角三角形。不规则三角网是由规则三角网抽稀并再三角剖分后，由离散的点集形成。不规则三角网可以缩减数字地面模型的存储空间，且能大幅度提高导航应用的渲染效率。但是，由于一个大范围的DTM往往由多张像素图组成，在对DTM进行抽稀再三角剖分时，每个像素图得出的不规则三角网可能并不相同，所以各个不规则三角网的边缘会有出现缝隙。例如，如图1所示的两个相邻的三角网tile1与tile2之间有缝隙。
为了填补缝隙，目前通常在两个三角网之间填充一个垂直于水平面的平面，从而从视觉上消除缝隙，但是两个三角网之间缝隙处的地形不连贯，填补不自然，且由于增加了面积较大的平面，总体绘制效率会降低。
发明内容
有鉴于此，本发明的目的在于提供一种处理数字地面模型接缝的方法及装
置，以使三角网之间接缝严密、连贯且不降低总体绘制效率的目的。
在本发明实施例的第一个方面，提供了一种处理数字地面模型接缝的方法，例如，该方法可以包括：通过对数字地面模型中两个相邻三角网相邻边界上边界点在水平面上的坐标进行比较，区分出所述两个相邻三角网在水平面上投影相互重合的边界点，以及，在水平面上投影不重合的边界点；将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更新为相同的值；将所述两个相邻三角网的相邻边界上的投影不重合的边界点分别添加到对方的边界点集合中；基于添加的边界点，在三角网内形成新的三角形。
在本发明实施例的第二个方面，提供了一种处理数字地面模型接缝的装置，例如，该装置可以包括：比较单元，用于通过对数字地面模型中两个相邻三角网相邻边界上边界点在水平面上的坐标进行比较，区分出所述两个相邻三角网在水平面上投影相互重合的边界点，以及，在水平面上投影不重合的边界点；重合边界点处理单元，用于将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更新为相同的值；非重合边界点处理单元，用于将所述两个相邻三角网的相邻边界上的投影不重合的边界点分别添加到对方的边界点集合中；新三角形划分单元，用于基于添加的边界点，在三角网内形成新的三角形。
可见，本发明实施例具有如下有益效果：
由于本发明实施例将两个相邻三角网相邻边界上投影重合的边界点在垂直水平面轴向上的高度坐标更新为相同的值，因此，两个三角网边界之间同一位置高度落差造成的缝隙可以被消除，由于将投影不重合的边界点在所述两个相邻三角网的相邻边界中互补，针对三角网内新加入的边界点在该三角网内形成新的三角形，因此，两个三角网边界之间水平面上的缝隙可以通过在三角网中针对新加入的边界点形成新三角形来填补，从而使得两个相邻三角网的相邻边界处地形连贯，缝隙处的填补更加自然。如图6所示的具有缝隙的两个相邻三角网，在应用本发明实施例提供的方法填充缝隙后，两个相邻三角网之间的缝隙消除，而且，相邻边界的两边地形连贯、自然、更加接近现实世界地形外貌。而且，针对新加入的边界点形成新三角形不会额外增加大面积的平面，不会降低总体绘制效率。
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为两个相邻的三角网之间缝隙的简单示意图；
图2为本发明实施例提供的一种处理数字地面模型接缝的方法流程示意图之一；
图3为应用本发明实施例提供的方法处理缝隙后的两个相邻的三角网的简单示意图；
图4为本发明实施例提供的一种处理数字地面模型接缝的方法流程示意图之二；
图5为两个相邻的三角网之间缝隙的三维立体示意图；
图6为应用本发明实施例提供的方法处理缝隙后的两个相邻的三角网的三维立体示意图；
图7为本发明实施例提供的一种处理数字地面模型接缝的装置结构示意图之一；
图8为本发明实施例提供的一种处理数字地面模型接缝的装置结构示意图之二。
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
例如，参见图2，为本发明提供的一种处理数字地面模型接缝的方法流程示意图。如图2所示，该方法可以包括：
S210、通过对数字地面模型中两个相邻三角网相邻边界上边界点在水平面
上的坐标进行比较，区分出所述两个相邻三角网在水平面上投影相互重合的边界点，以及，在水平面上投影不重合的边界点。
例如，在本发明实施例一些可能的实施方式中，所述数字地面模型可以配置有用于保存各个像素点的水平面位置信息的配置文件。其中，水平面位置信息可以用地理上的经纬度坐标来表示。假设本发明实施例中以(x，y，z)表示一个像素点的三维坐标，其中，(x，y)用于表示像素点在水平面上的坐标(如地理上的经纬度坐标)，z用于表示像素点在垂直水平面轴向上的高度。则可以从数字地面模型的配置文件中读取出任意两个相邻三角网的相邻边界上的边界点在水平面上的坐标(x，y)，进而将任意两个相邻三角网的相邻边界上的边界点在水平面上的坐标(x，y)进行比较。
例如，参见图1，为两个相邻的三角网tile1与tile2的简单示意图。在三角网tile1中，与三角网tile2相邻的相邻边界对应的边界点集合为{p1、p2、p3}。在三角网tile2中，与三角网tile1相邻的相邻边界对应的边界点集合为{p4、p5}。则可以从数字地面模型的配置文件中读取出p1、p2、p3、p4、p5这五个边界点的x、y坐标。通过将这五个边界点的x、y坐标进行比较，可以从在这五个边界点中区分出在水平面上投影重合的边界点，以及，在水平面上投影不重合的边界点。
需要说明的是，在实际应用中，可以根据实际需要设置边界点在水平面上是否重合的判断标准。例如，如果对缝隙填补精度要求较高，则可以将投影重合的边界点的判断标准设置为在水平面上的坐标完全相同，如果对缝隙填补精度要求不高，则可以允许一定误差，将水平面上坐标之差小于误差的两个边界点判定为投影重合的边界点，从而可以在后续的步骤中一定程度上减少投影不重合边界点的数量，减少重新划分三角形的相关计算，提高填补缝隙的效率。
S220、将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更新为相同的值。
S230、将所述两个相邻三角网的相邻边界上的投影不重合的边界点分别添加到对方的边界点集合中。
需要说明的是，所述步骤S220与步骤S230之间并无特定的执行顺序，可以并行执行，也可以先后执行。
针对在水平面上投影重合的边界点，无论采取哪种计算更新高度的方法，只要保证将水平面上投影重合的任意两个边界点在垂直水平面轴向上的高度坐标更新为相同的值即可达到填补这两个边界点之间缝隙的目的。因此，本发明在实际实施时，可以根据实际需要设置任意可能的计算方法，本发明对此并不进行限制。
需要说明的是，可以根据数字地面模型中用于描述三角形结构以及三角形组成的三角网结构的数据确定三角网之间的相邻边界，以及，所述相邻边界对应的边界点集合。因此，可以将投影不重合的边界点添加到对方三角网相邻边界对应的边界点集合中。例如，假设所述两个相邻三角网分别为第一三角网以及第二三角网，则可以将第一三角网相邻边界上所有与第二三角网相邻边界不重合的边界点添加到第二三角网相邻边界对应的边界点集合中，将第二三角网相邻边界上所有与第一三角网不重合的边界点添加到第一三角网相邻边界对应的边界点集合中。
例如，在图1所示的示例中，p2是tile1与tile2的相邻边界在水平面上投影不重合的边界点，则可以将p2加入tile2的相邻边界对应的边界点集合中。加入p2后，tile2的相邻边界对应的边界点集合为{p4，p2，p5}。
S240、针对所述三角网内新加入的边界点，在该三角网内形成新的三角形。
一个三角网会被划分为若干个三角形，三角网中的每一个边界点都会落入到一个三角形中。三角形在三角网中的划分方式可以根据实际需要设置。因此，本发明对针对新加入的边界点在三角网中如何形成新的三角形的具体实施方式并不进行限制。
应用本发明实施例提供的方法，由于将两个相邻三角网相邻边界上投影重合的边界点在垂直水平面轴向上的高度坐标更新为相同的值，因此，两个三角网边界之间同一位置高度落差造成的缝隙可以被消除，由于将投影不重合的边界点在所述两个相邻三角网的相邻边界中互补，针对三角网内新加入的边界点在该三角网内形成新的三角形，因此，两个三角网边界之间水平面上的缝隙可以通过在三角网中针对新加入的边界点形成新三角形来填补。可见，应用本发明实施例填充两个相邻三角网之间的缝隙后，两个相邻三角网之间的地形连贯、填补自然，且针对新加入的边界点形成新三角形不会额外增加大面积的平
面，不会降低总体绘制效率。
下面，对本发明实施例一些可能的实施方式进行详细说明。
例如，参见图4，为本发明提供的一种处理数字地面模型接缝的方法流程示意图。如图4所示，该方法可以包括：
S410、通过对数字地面模型中两个相邻三角网相邻边界上边界点在水平面上的坐标进行比较，区分出所述两个相邻三角网在水平面上投影相互重合的边界点，以及，在水平面上投影不重合的边界点。
一些可能的实施方式中，为了便于比较，可以预先将两个相邻三角网相邻边界上边界点按水平面同一纬度上的坐标，如y轴坐标进行排序，将排序后的边界点的坐标进行比较。
S420、计算出在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标的平均值；将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更新为所述平均值。需要说明的是，针对在水平面上投影相互重合的边界点的更新高度，本发明实施例中还可以采取其他计算方法，计算平均值仅为本发明一种可能的实施方式。
例如，假设在图1所示的示例中，p1点与p4点的(x，y)坐标，p3点与p5点的(x，y)坐标相同，则可以从数字地面模型的像素图中读取出p1点、p3、p4、p5的高度坐标z，计算p1点与p4点的高度坐标z的平均值，将p3点与p5点的高度坐标z更新为该平均值，计算p3点与p5点的高度坐标z的平均值，将p3点与p5点的高度坐标z更新为该平均值。更新p1、p4、p3、p5点的高度坐标之后，如图3所示，p1与p4成为三维空间位置相同的边界点、p3与p5成为三维空间位置相同的边界点。
S430、将所述两个相邻三角网的相邻边界上的投影不重合的边界点分别添加到对方的边界点集合中。
S440、在新加入了边界点的三角网内，查询出所述新加入的边界点在该三角网内所在的三角形，将所述新加入的边界点在所述三角形所在边的对角顶点，与所述新加入的边界点连线，使所述三角形划分出新的三角形。
一些可能的实施方式中，可以预先将两个相邻三角网相邻边界上边界点按水平面同一纬度上的坐标，如y轴坐标进行排序。因此，可以根据该排序结果
快速找到新加入的边界点落入哪两个已有的边界点之间。查询出所落入的两个已有的边界点所在的三角形，将该三角形的另一顶点与该新加入的边界点连线，使所述三角形划分出新的三角形。
例如，在图3所示的示例中，将p2添加到tile2的相邻边界对应的边界点集合之后，可以找到与tile 2原有的边界点p4，p5相关联的三角形，即p4，p5、p0三点组成的三角形，进而在p2与该三角形的另一顶点p0间连线，从而将p4，p5、p0三点组成的三角形划分为p2、p4、p0以及p2、p5、p0分别组成的两个三角形。
可见，应用该实施例提供的方法，可以通过简单高效的平均值计算方法得到投影重合的边界点的更新高度、利用新加入的边界点对其加入的三角网内关联的三角形进行划分形成新的三角形，从而无需再次进行全局的三角剖分，填补三角网之间缝隙的效率较高，而且接缝严密，不损失现实地形的外貌特点。例如，图5，图6分别为两个相邻三角网的三维立体示意图。应用本发明实施例提供的方法对如图5所示的两个相邻三角网之间的缝隙进行处理后，可以得到图6所示的填补缝隙后的两个相邻三角网。
本发明实施例可以用来处理多个三角网进行整合时出现的缝隙，也可以用来处理格网式数字地面模型的像素图进行三角剖分抽稀时出现的缝隙。具体地，例如，本发明实施例所述两个相邻三角网可以为，所述数字地面模型中所有不规则三角网中的每两个相邻的不规则三角网，或者，本发明实施例所述两个相邻三角网可以为，所述数字地面模型中根据整合需求所确定的部分不规则三角网中的每两个相邻的不规则三角网。
可以理解的是，在处理数字地面模型中所有不规则三角网中的每两个相邻的不规则三角网之间的缝隙，或者，在处理部分不规则三角网中的每两个相邻的不规则三角网之前，可以先确定数字地面模型中所有不规则三角网之间，或者，部分不规则三角网之间，所有两个相邻的不规则三角网的组合。针对每一种组合，使用本发明实施例提供的方法处理该组中两个相邻的不规则三角网之间的缝隙。
在上面的实施方式中，所述三角网可以是数字地面模型中经过抽稀后的不
规则三角网。因此，本发明实施例还可以包括通过对所述数字地面模型转化为规则三角网、对所述规则三角网进行抽稀以及三角剖分获得所述不规则三角网的步骤。
与上述处理数字地面模型接缝的方法相对应地，本发明实施例还提供了一种处理数字地面模型接缝的装置。
例如，参见图7，为本发明实施例提供的处理数字地面模型接缝的装置结构示意图。如图7所示，该装置可以包括：
比较单元710，可以用于通过对数字地面模型中两个相邻三角网相邻边界上边界点在水平面上的坐标进行比较，区分出所述两个相邻三角网在水平面上投影相互重合的边界点，以及，在水平面上投影不重合的边界点。
重合边界点处理单元720，可以用于将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更新为相同的值。
非重合边界点处理单元730，可以用于将所述两个相邻三角网的相邻边界上的投影不重合的边界点分别添加到对方的边界点集合中。
新三角形划分单元740，可以用于基于添加的边界点，在三角网内形成新的三角形。
配置本发明实施例提供的装置，由于重合边界点处理单元720将两个相邻三角网相邻边界上投影重合的边界点在垂直水平面轴向上的高度坐标更新为相同的值，因此，两个三角网边界之间同一位置高度落差造成的缝隙可以被消除，由于非重合边界点处理单元730将投影不重合的边界点在所述两个相邻三角网的相邻边界中互补，且新三角形划分单元740针对三角网内新加入的边界点在该三角网内形成新的三角形，因此，两个三角网边界之间水平面上的缝隙可以通过在三角网中针对新加入的边界点形成新三角形来填补。可见，本发明实施例对缝隙处的填补自然，且针对新加入的边界点形成新三角形不会额外增加大面积的平面，不会降低总体绘制效率
一些可能的实施方式中，如图8所示，所述重合边界点处理单元720可以包括：平均高度计算子单元721，可以用于计算出在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标的平均值。高度更新子单元722，可以用于将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更
新为所述平均值。
一些可能的实施方式中，如图8所示，所述新三角形划分单元740可以包括：查询三角形子单元741，可以用于在新加入了边界点的三角网内，查询出所述新加入的边界点在该三角网内所在的三角形。连线子单元742，可以用于将所述新加入的边界点在所述三角形所在边的对角顶点，与所述新加入的边界点连线，使所述三角形划分出新的三角形。
本发明实施例提供的装置可以用来处理多个三角网进行整合时出现的缝隙，也可以用来处理格网式数字地面模型的像素图进行三角剖分抽稀时出现的缝隙。具体地，例如，本发明实施例所述两个相邻三角网可以为，所述数字地面模型中所有不规则三角网中的每两个相邻的不规则三角网，或者，本发明实施例所述两个相邻三角网可以为，所述数字地面模型中根据整合需求所确定的部分不规则三角网中的每两个相邻的不规则三角网。
在上面的实施方式中，本发明实施例提供的装置还可以通过对所述数字地面模型转化为规则三角网、对所述规则三角网进行抽稀以及三角剖分获得所述不规则三角网。
可见，配置上述实施例提供的装置，可以由通过简单高效的平均值计算方法得到投影重合的边界点的更新高度、利用新加入的边界点对其加入的三角网内关联的三角形进行划分形成新的三角形，从而无需再次进行全局的三角剖分，填补三角网之间缝隙的效率较高，而且接缝严密，不损失现实地形的外貌特点。
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制。虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明。任何熟悉本领域的技术人员，在不脱离本发明技术方案范围情况下，都可利用上述揭示的方法和技术内容对本发明技术方案做出许多可能的变动和修饰，或修改为等同变化的等效实施例。因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均仍属于本发明技术方案保护的范围内。
Claims (8)
- 一种处理数字地面模型接缝的方法，其特征在于，包括：通过对数字地面模型中两个相邻三角网相邻边界上边界点在水平面上的坐标进行比较，区分出所述两个相邻三角网在水平面上投影相互重合的边界点，以及，在水平面上投影不重合的边界点；将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更新为相同的值；将所述两个相邻三角网的相邻边界上的投影不重合的边界点分别添加到对方的边界点集合中；基于添加的边界点，在三角网内形成新的三角形。
- 根据权利要求1所述的方法，其特征在于，所述将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更新为相同的值包括：计算出在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标的平均值；将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更新为所述平均值。
- 根据权利要求1所述的方法，其特征在于，所述基于添加的边界点，在三角网内形成新的三角形包括：在新加入了边界点的三角网内，查询出所述新加入的边界点在该三角网内所在的三角形；将所述新加入的边界点在所述三角形所在边的对角顶点，与所述新加入的边界点连线，使所述三角形划分出新的三角形。
- 根据权利要求1-3任一项所述的方法，其特征在于，所述两个相邻三角网为，所述数字地面模型中所有不规则三角网中的每两个相邻的不规则三角网；或者，所述两个相邻三角网为，所述数字地面模型中根据整合需求所确定的部分不规则三角网中的每两个相邻的不规则三角网。
- 一种处理数字地面模型接缝的装置，其特征在于，包括：比较单元，用于通过对数字地面模型中两个相邻三角网相邻边界上边界点在水平面上的坐标进行比较，区分出所述两个相邻三角网在水平面上投影相互重合的边界点，以及，在水平面上投影不重合的边界点；重合边界点处理单元，用于将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更新为相同的值；非重合边界点处理单元，用于将所述两个相邻三角网的相邻边界上的投影不重合的边界点分别添加到对方的边界点集合中；新三角形划分单元，用于基于添加的边界点，在三角网内形成新的三角形。
- 根据权利要求5所述的装置，其特征在于，所述重合边界点处理单元包括：平均高度计算子单元，用于计算出在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标的平均值；高度更新子单元，用于将在水平面上投影相互重合的边界点在垂直水平面轴向上的高度坐标更新为所述平均值。
- 根据权利要求5所述的装置，其特征在于，所述新三角形划分单元包括：查询三角形子单元，用于在新加入了边界点的三角网内，查询出所述新加入的边界点在该三角网内所在的三角形；连线子单元，用于将所述新加入的边界点在所述三角形所在边的对角顶点，与所述新加入的边界点连线，使所述三角形划分出新的三角形。
- 根据权利要求5-7任一项所述的装置，其特征在于，所述两个相邻三角网为，所述数字地面模型中所有不规则三角网中的每两个相邻的不规则三角网；或者，所述两个相邻三角网为，所述数字地面模型中根据整合需求所确定的部分不规则三角网中的每两个相邻的不规则三角网。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410785247.9 | 2014-12-16 | ||
| CN201410785247.9A CN105761309B (zh) | 2014-12-16 | 2014-12-16 | 一种处理数字地面模型接缝的方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2016095682A1 true WO2016095682A1 (zh) | 2016-06-23 |
Family
ID=56125874
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/CN2015/095642 WO2016095682A1 (zh) | 2014-12-16 | 2015-11-26 | 一种处理数字地面模型接缝的方法及装置 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN105761309B (zh) |
| WO (1) | WO2016095682A1 (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110750556A (zh) * | 2019-10-26 | 2020-02-04 | 武汉中海庭数据技术有限公司 | 一种地理数据抽稀补缝方法、装置及存储介质 |
| CN111105492A (zh) * | 2019-12-18 | 2020-05-05 | 四川大学 | 一种地形平滑过程中瓦片高程数据快速调度方法 |
| CN111506679A (zh) * | 2019-01-31 | 2020-08-07 | 阿里巴巴集团控股有限公司 | 地图要素数据的生成、显示方法及装置、介质、设备 |
| CN112164137A (zh) * | 2020-10-09 | 2021-01-01 | 山东捷瑞数字科技股份有限公司 | 一种虚拟仿真地形直观绘制及接缝自然过渡的方法 |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106251331A (zh) * | 2016-07-15 | 2016-12-21 | 北京道亨时代科技有限公司 | 一种倾斜测量场景中地物的提取方法 |
| CN108205565A (zh) * | 2016-12-19 | 2018-06-26 | 北京四维图新科技股份有限公司 | 电子地图元素抽稀方法、装置及终端 |
| CN108346170B (zh) * | 2017-01-22 | 2022-08-12 | 腾讯科技（深圳）有限公司 | 电子地图的缝隙填充方法及装置 |
| CN108765535A (zh) * | 2018-05-30 | 2018-11-06 | 链家网（北京）科技有限公司 | 一种三维点云缺失数据区域的标记方法及装置 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060066613A1 (en) * | 2004-09-27 | 2006-03-30 | Elshishiny Hisham E E | Method and system for partitioning the surface of a three dimentional digital object model in order to map a texture |
| KR100609004B1 (ko) * | 2006-02-03 | 2006-08-03 | (자)한진개발공사 | 수치사진측량기술을 이용한 지형복원방법 |
| CN101320488A (zh) * | 2008-06-27 | 2008-12-10 | 于华明 | 全球海洋三角形网格的构建方法 |
| CN102129712A (zh) * | 2011-03-10 | 2011-07-20 | 东南大学 | 基于多地层及三维的土石方数量的三角网模型构建方法 |
| CN103093500A (zh) * | 2013-02-28 | 2013-05-08 | 中铁第四勘察设计院集团有限公司 | 全球层次地形模型开孔建模方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101046893A (zh) * | 2007-01-26 | 2007-10-03 | 北京国泰清济科技发展有限公司 | 快速不规则三角网生成方法 |
| CN102054292A (zh) * | 2009-10-30 | 2011-05-11 | 吴立新 | 基于tin耦合的地上下无缝集成三维空间建模方法 |
| CN102117500B (zh) * | 2009-12-30 | 2015-03-25 | 新奥特（北京）视频技术有限公司 | 一种三维模型构建方法及系统 |
| CN103902343A (zh) * | 2014-04-22 | 2014-07-02 | 江苏省城市规划设计研究院 | 一种基于Delaunay三角网精度控制的瓦片地图下载与拼接方法 |
-
2014
- 2014-12-16 CN CN201410785247.9A patent/CN105761309B/zh active Active
-
2015
- 2015-11-26 WO PCT/CN2015/095642 patent/WO2016095682A1/zh active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060066613A1 (en) * | 2004-09-27 | 2006-03-30 | Elshishiny Hisham E E | Method and system for partitioning the surface of a three dimentional digital object model in order to map a texture |
| KR100609004B1 (ko) * | 2006-02-03 | 2006-08-03 | (자)한진개발공사 | 수치사진측량기술을 이용한 지형복원방법 |
| CN101320488A (zh) * | 2008-06-27 | 2008-12-10 | 于华明 | 全球海洋三角形网格的构建方法 |
| CN102129712A (zh) * | 2011-03-10 | 2011-07-20 | 东南大学 | 基于多地层及三维的土石方数量的三角网模型构建方法 |
| CN103093500A (zh) * | 2013-02-28 | 2013-05-08 | 中铁第四勘察设计院集团有限公司 | 全球层次地形模型开孔建模方法 |
Non-Patent Citations (1)
| Title |
|---|
| LIU, QINGYUAN ET AL.: "Establishment of DTM and Research into Algorithm of Boundary Line", MINE SURVEYING, 31 March 2004 (2004-03-31) * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111506679A (zh) * | 2019-01-31 | 2020-08-07 | 阿里巴巴集团控股有限公司 | 地图要素数据的生成、显示方法及装置、介质、设备 |
| CN111506679B (zh) * | 2019-01-31 | 2023-05-26 | 阿里巴巴集团控股有限公司 | 地图要素数据的生成、显示方法及装置、介质、设备 |
| CN110750556A (zh) * | 2019-10-26 | 2020-02-04 | 武汉中海庭数据技术有限公司 | 一种地理数据抽稀补缝方法、装置及存储介质 |
| CN111105492A (zh) * | 2019-12-18 | 2020-05-05 | 四川大学 | 一种地形平滑过程中瓦片高程数据快速调度方法 |
| CN111105492B (zh) * | 2019-12-18 | 2020-11-17 | 四川大学 | 一种地形平滑过程中瓦片高程数据快速调度方法 |
| CN112164137A (zh) * | 2020-10-09 | 2021-01-01 | 山东捷瑞数字科技股份有限公司 | 一种虚拟仿真地形直观绘制及接缝自然过渡的方法 |
| CN112164137B (zh) * | 2020-10-09 | 2021-05-11 | 山东捷瑞数字科技股份有限公司 | 一种虚拟仿真地形直观绘制及接缝自然过渡的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105761309A (zh) | 2016-07-13 |
| CN105761309B (zh) | 2019-02-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2016095682A1 (zh) | 一种处理数字地面模型接缝的方法及装置 | |
| WO2020192354A1 (zh) | 一种嵌入式城市设计场景仿真方法及系统 | |
| CN109186551B (zh) | 倾斜摄影测量建筑物特征点提取方法、装置与存储介质 | |
| CN103871102B (zh) | 一种基于高程点和道路轮廓面的道路三维精细建模方法 | |
| CN109269472B (zh) | 倾斜摄影测量建筑物特征线提取方法、装置与存储介质 | |
| WO2022227910A1 (zh) | 虚拟场景生成方法、装置、计算机设备和存储介质 | |
| US8600713B2 (en) | Method of online building-model reconstruction using photogrammetric mapping system | |
| CN107622530B (zh) | 一种高效鲁棒的三角网切割方法 | |
| CN110838129A (zh) | 基于倾斜摄影测量的三维建筑物模型轮廓特征线提取方法 | |
| CN104183016A (zh) | 一种快速的2.5维建筑物模型的构建方法 | |
| CN113112603A (zh) | 三维模型优化的方法和装置 | |
| CN114332291A (zh) | 一种倾斜摄影模型建筑物外轮廓规则提取方法 | |
| CN115797592B (zh) | 一种基于倾斜摄影三维模型自动生成建筑物体块的方法和装置 | |
| CN111986326B (zh) | 路面三维模型的构造方法、装置、计算机设备及存储介质 | |
| CN112561788A (zh) | 一种bim模型的二维展开方法及纹理贴图方法、装置 | |
| CN111127649B (zh) | 构建三维体块模型的方法及装置、服务器 | |
| CN110675323B (zh) | 一种三维地图语义处理方法、系统、设备及计算机介质 | |
| CN105354878A (zh) | 一种电子地图数据的制作方法及装置 | |
| CN116086411A (zh) | 数字地形图生成方法、装置、设备和可读存储介质 | |
| CN105184854A (zh) | 针对地下空间扫描点云成果数据的快速建模方法 | |
| WO2020114273A1 (zh) | 商户查找方法、装置、电子设备和存储介质 | |
| CN109064482B (zh) | 一种三维倾斜摄影场景中自动获取房屋轮廓的方法及装置 | |
| JP2002092658A (ja) | ３次元デジタル地図作成装置及び３次元デジタル地図作成のプログラムを記憶した記憶媒体 | |
| CN109163708B (zh) | 无人机大比例尺地形图房屋线的优化方法、装置与介质 | |
| CN114581620A (zh) | 道路的虚拟高程生成方法、装置、计算机设备和存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 15869190 Country of ref document: EP Kind code of ref document: A1 |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 15869190 Country of ref document: EP Kind code of ref document: A1 |