WO2013057760A1 - Système de génération de modèle numérique d'altitude et procédé de génération d'un modèle numérique d'altitude - Google Patents

Système de génération de modèle numérique d'altitude et procédé de génération d'un modèle numérique d'altitude Download PDF

Info

Publication number
WO2013057760A1
WO2013057760A1 PCT/JP2011/005838 JP2011005838W WO2013057760A1 WO 2013057760 A1 WO2013057760 A1 WO 2013057760A1 JP 2011005838 W JP2011005838 W JP 2011005838W WO 2013057760 A1 WO2013057760 A1 WO 2013057760A1
Authority
WO
WIPO (PCT)
Prior art keywords
flattened
flattened region
sub
contour
curve
Prior art date
Application number
PCT/JP2011/005838
Other languages
English (en)
Inventor
Tao Guo
Original Assignee
Hitachi, 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 Hitachi, Ltd. filed Critical Hitachi, Ltd.
Priority to PCT/JP2011/005838 priority Critical patent/WO2013057760A1/fr
Publication of WO2013057760A1 publication Critical patent/WO2013057760A1/fr

Links

Images

Classifications

    • 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/12—Relief maps

Definitions

  • This invention relates to a technology for generating sub-contour lines from existing contour lines of a digital topographic map by means of interpolation.
  • a digital elevation model is one of the most important fundamental data in a geographic information system (GIS) and is normally acquired from field survey, measurement on aerial photos, 3D matching from stereo pair of remote sensing images, and 3D point cloud of a set of points acquired from laser scanning.
  • GIS geographic information system
  • the DEM is available at a certain scale due to its data source and generation method. Moreover, the scale determines the interval between contour lines.
  • TIN triangulated irregular network
  • the algorithm of the TIN model is to form the triangulated regions from three nodes on contour lines and interpolation is implemented in each triangle. Therefore, the interpolation employing the TIN model is fast and simple, and is thus widely used.
  • the TIN model exhibits a high performance in representing a terrain with a sharp slope.
  • Patent Literatures 1, 2, and 3 various methods are being studied to solve the deformations of terrain in the DEM as disclosed in Patent Literatures 1, 2, and 3.
  • the terrain in the triangle is treated as a horizontal surface.
  • the nodes of a triangle may be placed on a contour line with the same elevation on a terrain such as a saddle region and a gentle slope of a ridge, the interpolation using the TIN model often determines this region as flat, and the terrain of the triangle patch is thus considered as flat, resulting in strangely deforming the terrain, and affecting the generation of the DEM.
  • the weak points of the TIN model such as influence of the above-mentioned flattening, have not been studied for a long time.
  • This invention has an object to propose a method of generating a highly-accurate DEM by gently interpolating the elevation particularly on a saddle region and a ridge, and generating new contour lines.
  • a representative example of the disclosed invention in this application calculates a vector field gradually changing from one boundary to another on a flattened region, smoothly interpolates the elevation, and generates sub-contour lines.
  • a digital elevation model generating system comprising: a storage unit for storing a contour map; a processing unit for carrying out an interpolation process for generating sub-contour lines; a memory for storing a program to be executed by the processing unit; an input unit for receiving an input of the contour map; and an output unit for outputting a generated elevation map.
  • the processing unit includes: a reception module for partitioning the input contour map into a triangulated irregular network; a marking module for marking each of triangle patches as flattened or non-flattened, and merging the triangle patches marked as flattened to form a flattened region; a TIN calculating module for generating curves serving as sub-contour lines by means of interpolation using the triangulated irregular network for the triangle patch marked as non-flattened; a vector field generating module for calculating a vector field in which each vector represents a spatial effect of neighboring contour lines for the flattened region; a sub-contour line generating module for determining a curve on which vectors have the same magnitude; and a sub-contour line connecting module for generating a newly integrated sub-contour line by creating a buffer in a neighborhood of a boundary between a non-flattened region constituted of the triangle patches marked as non-flattened and the flattened region, and by smoothly linking the curve in the non-flattened region
  • a highly-accurate DEM can be generated.
  • Fig. 1 is a block diagram illustrating a configuration of a matching system according to an embodiment of this invention.
  • Fig. 2 is an overall flowchart of processing of generating a digital elevation model according to the embodiment of this invention.
  • Fig. 3 is an explanatory diagram illustrating triangle patches generated in the embodiment of this invention.
  • Fig. 4 is an explanatory diagram illustrating triangle patches generated in the embodiment of this invention.
  • Fig. 5 is an explanatory diagram illustrating generation of effect vectors according to the embodiment of this invention.
  • Fig. 6 is an explanatory diagram illustrating generation of sub-contour lines according to the embodiment of this invention.
  • Fig. 7 is an explanatory diagram illustrating determination of the elevations of the sub-contour lines according to the embodiment of this invention.
  • Fig. 1 is a block diagram illustrating a configuration of a matching system according to an embodiment of this invention.
  • Fig. 2 is an overall flowchart of processing of generating a digital elevation model according to the embodiment of this
  • FIG. 8 is an explanatory diagram illustrating connecting of the sub-contour lines according to the embodiment of this invention.
  • Fig. 9 is an explanatory diagram illustrating an effect vector field generated by the processing of generating the digital elevation model according to the embodiment of this invention.
  • Fig. 10 is an explanatory diagram illustrating a generation of TIN according to the embodiment of this invention.
  • Fig. 1 is a block diagram illustrating a configuration of a digital elevation model (DEM) generating system according to an embodiment of this invention.
  • DEM digital elevation model
  • the DEM generating system of this embodiment is a computer including an input unit 110, a display unit 120, a data storage unit 130, a processing unit 140, an elevation map output unit 150, and a storage unit 160.
  • the input unit 110, the display unit 120, the data storage unit 130, the elevation map output unit 150, and the storage unit 160 are connected with each other through the processing unit 140 (or mutually through a bus).
  • the input unit 110 includes a map input unit 111 and a parameter input unit 112.
  • the map input unit 111 is an interface used for receiving an input of a digital contour map, and is constituted of an optical disc drive, a USB interface, and the like, for example.
  • the parameter input unit 112 is a device used for receiving an input of parameters (such as the number and interval of sub-contour lines to be generated) used for processing of interpolating the sub-contour lines, and is constituted of a keyboard, a mouse, and the like, for example. It should be noted that the parameter input unit 112 may receive an input of individual parameters, or a set of parameters prepared in advance depending on the interval of the sub-contour lines to be generated.
  • the display unit 120 includes a map displaying unit 121 and a parameter displaying unit 122.
  • the map displaying unit 121 is a display device for displaying a map to be processed by this system.
  • the parameter displaying unit 122 is a display device for displaying parameters to be used for the processing of interpolating the sub-contour lines. It should be noted that the map displaying unit 121 and the parameter displaying unit 122 may be constituted of the same display device or independent display devices.
  • the data storage unit 130 is a non-volatile storage device for storing data for the digital contour map to be processed by the system, and is constituted of a magnetic disk device, a non-volatile memory, or the like, for example.
  • the processing unit 140 includes a processor, and carries out processing to be carried out in this system by executing a program loaded on the memory.
  • the elevation map output unit 150 is a device for outputting a result processed by this system as a vector form or an image form, and is constituted of a printer, a plotter, or the like, for example.
  • the storage unit 160 is a storage device for storing the program to be executed by the processing unit 140, and is constituted of a hard disk drive, a non-volatile memory, or the like, for example.
  • the storage unit 160 stores sub-programs 161-167 for installing each of a reception module, a marking module, a TIN calculating module, a vector field generating module, a sub-contour line generating module, a sub-contour line connecting module, and a terrain deviation determination module.
  • DEM generating processing (Fig. 2) is installed by the processing unit 140 executing those programs.
  • the reception sub-program 161 acquires the contour map input to the map input unit 111, and partitions a region to be processed into a triangulated irregular network (TIN).
  • the marking sub-program 162 marks each triangle patch as flattened or non-flattened, and merges the triangle patches marked as flattened to form a flattened region.
  • the TIN calculating sub-program 163 interpolates sub-contour lines by using the TIN in the triangle patches marked as non-flattened (non-flattened region).
  • the vector field generating sub-program 164 generates a vector field constructed by vectors representing a spatial effect of surrounding contour lines for the flattened region.
  • the sub-contour line generating sub-program 165 defines a curve on which a vector at every location has the same magnitude.
  • the sub-contour line connecting sub-program 166 smoothly links the sub-contour lines in a buffer provided in a neighborhood of a boundary between the non-flattened region and the flattened region.
  • the terrain deviation determination sub-program 167 determines whether the sub-contour line is higher or lower than the neighboring contour lines.
  • Fig. 2 is a flowchart of the DEM generating processing according to the embodiment of this invention.
  • the map input unit 111 receives an input of a contour map on which sub-contour lines are to be interpolated (S1).
  • This contour map includes data on a plurality of contour lines.
  • the input contour map is digitized data, and is to be stored in the data storage unit 130.
  • the parameter input unit 112 receives a specification of a region 301 as illustrated in Fig. 3 in which sub-contour lines are to be generated in the input contour map.
  • the processing unit 140 provides points at a predetermined interval on contour lines in the specified region by executing the reception sub-program 161.
  • the interval of the points provided on the contour lines is input to the parameter input unit 112, and a small value is set to the interval if the number of sub-contour lines to be generated is large.
  • the processing unit 140 selects two neighboring points on a certain contour line, and selects a point on a neighboring contour line which has the shortest sum of the distances from the selected two points. It should be noted that the neighboring contour line includes both a contour line at a different elevation and a contour line at the same elevation. Then, the processing unit 140 generates a triangle patch using the selected three points.
  • Fig. 3 illustrates a state in which triangle patches are generated in a specified region 301.
  • the processing unit 140 marks the generated triangle patches to a flattened region or a non-flattened region by executing the marking sub-program 162 (S2). Specifically, this marking compares the elevations of the contour lines on which the points are set, and determines that the triangle patch is a flattened region when the elevations of the three vertices of the triangle patch are equal, and the triangle patch is a non-flattened region when at least one elevation is different.
  • triangle patches AFB, FBG, BGC, GCH, CHD, HDI, and DIE have vertices on contour lines at the same elevation, and are marked as flattened regions.
  • the other triangle patches have at least one vertex on a contour line at a different elevation, and are thus marked as non-flattened regions.
  • triangle patches BFE, BEC, and ECD are formed by points on a contour line at the same elevation as vertices, and are thus marked as flattened regions.
  • triangle patches AHB, BHF, and FHG are formed using at least one point on a contour line at a different elevation as a vertex, and are thus marked as non-flattened regions.
  • a saddle region between contour lines of 200 m illustrated in Fig. 3 a gentle slope of a ridge in a contour line of 100 m extending toward the left in Fig. 4 out of terrains contained in the contour maps are marked as flattened regions.
  • those flattened regions are marked as planes, and thus sub-contour lines cannot be generated.
  • the processing unit 140 merges neighboring triangle patches marked as non-flattened regions to form one wider non-flattened region, and merges neighboring triangle patches marked as flattened regions to form one wider flattened region by executing the marking sub-program 162. It should be noted that it is not necessary to merge triangle patches in a non-flattened region to which the conventional TIN is applied.
  • the processing unit 140 generates sub-contour lines using the TIN model in the merged non-flattened region by executing the TIN calculating sub-program 163 (S3), and generates sub-contour lines according to the method proposed by this specification in the merged flattened region (S4, S5).
  • the method of generating sub-contour lines according to this invention requires a larger calculation amount than the method using the TIN model, but the overall calculation amount can be reduced while accurate sub-contour lines can be generated by applying the method according to this invention only to a flattened region.
  • DEM generation method marks regions in which sub-contour lines are to be generated as non-flattened regions and flattened regions, and applies the method according to this invention only to the flattened regions, but the method according to this invention may be applied to all the regions in which sub-contour lines are to be generated.
  • Step S3 a description is given of the generation of sub-contour lines according to the TIN model applied to the non-flattened regions in Step S3.
  • the TIN model is represented by triangles continuously in contact with each other without overlapping, and a ground surface is represented by a plane in each of the triangles.
  • the triangles can be selected at any location, and the accuracy of the surface model increases even for a location at which the surface shape drastically changes if nodes are properly selected.
  • Delaunay triangulation is a popular triangulation method, which is able to guarantee that a circle drawn through the three nodes of a triangle include no other node.
  • all the vertices of the contour lines can be selected as nodes of the triangles to generate TIN. It is usually efficient to choose widely spaced nodes in regions where there is little variation in elevation whereas in areas of more intense variation in elevation the node density is increased.
  • point set1 ⁇ a1, b1, c1 ⁇ and point set2 ⁇ a2, b2, c2, d2 ⁇ are randomly selected from contour L1 and contour L2 respectively.
  • the Delaunay triangulation can be implemented in an incremental way, that is firstly to form a triangle, for example b1-a2-d2, then add the remains of points into this triangle to make a patch of triangles, then remove the triangles whose bounding circle includes other points. Repeat these steps until all points are included the triangles.
  • the height value at any location p within a triangle can be calculated through an interpolation of height values of the nodes of the triangle. As described above, sub-contour lines whose height values are between contour L1 and L2 can be easily generated.
  • the contour lines are a common data source for digital elevation data used for generating a TIN model.
  • an inflection point of a contour line is selected as a vertex of a triangle constituting the TIN model.
  • Nodes which can partition a wide region in which the elevation hardly changes are usually efficient.
  • density of nodes may be increased in a region in which the elevation changes drastically.
  • Step S4 as illustrated in Fig. 5, the processing unit 140 provides a point A in a flattened region first by executing the vector field generating sub-program 164.
  • Fig. 5 illustrates the one point A, as illustrated in the vector field of Fig. 9, the points A are provided in the form of a mesh at a predetermined interval in the flattened region.
  • contour lines C1 and C2 neighboring the point A are determined, a point Pi is provided on the contour line C1, and a point Qj is provided on the contour line C2.
  • effect vectors Vpij and Vqij are defined at the point A.
  • the directions of the vectors Vpij and Vqij are respectively directions from the point A to the point Pi and to the point Qj, and the magnitudes thereof are represented by the following equations.
  • is the distance between the point A and the point Pi
  • is the distance between the point A and the point Qi.
  • the effect vectors Vpij and Vqij are calculated for all combinations of the point Pi and the point Qj in the flattened region, and the calculated effect vectors Vpij and Vqij are summed to calculate a resultant vector V.
  • the vector V which is the accumulation of all effects from the neighboring two contour lines C1 and C2 is obtained.
  • the resultant vector V is divided by the number of accumulated vectors. In this way, the magnitude of the vector V is 1 on each of the contour lines C1 and C2, and the magnitude of the vector V is 0 on a center line of the contour lines C1 and C2.
  • a vector field of the effect vectors is generated by repeating the processing described above for all the points A.
  • This vector field changes smoothly between the contour lines C1 and C2, and this vector field can be used to interpolate the elevation between existing contour lines.
  • Step S5 a sub-contour line is generated by connecting points at which the magnitudes of calculated effect vectors are the same.
  • the processing unit 140 obtains a center line S of the two contour lines C1 and C2 by executing the sub-contour line generating sub-program 165.
  • the effects from the two contour lines C1 and C2 are the same at a point on this center line, and the magnitude of the effect vector is 0.
  • points A at which the magnitude of the effect vector is 0 or equal to or less than a predetermined small threshold are retrieved in the region between the two contour lines C1 and C2, and a curve connecting the retrieved points A is generated.
  • a cubic spline curve may be used for this curve (center line), for example.
  • the magnitudes of the effect vectors to be selected are determined according to the number of sub-contour lines to be generated. Specifically, when two sub-contour lines are to be generated between the contour lines C1 and C2, two lines are drawn at an equal terrain deviation between the contour lines C1 and C2. In other words, as illustrated in Fig. 6, one sub-contour line is generated each between the contour line C1 and the center line S and between the contour line C2 and the center line S. As described above, the magnitude of the effect vector V is 1 on each of the contour lines C1 and C2, the magnitude of the effect vector V is 0 on the center line S, and the magnitude of the effect vector is thus proportional to the terrain deviation. Therefore, the difference in magnitude of the effect vectors on the sub-contour line is 2/3, which is a value obtained by dividing 2 by (number of sub-contour lines to be generated + 1).
  • one sub-contour line can be generated between the contour line C1 and the center line S.
  • one sub-contour line can be generated between the contour line C2 and the center line S.
  • center line S is designated as a sub-contour line in the case where one sub-contour line is generated between the contour lines C1 and C2.
  • the vectors having the determined magnitude (1/3) are selected in the region between the contour line C1 and the center line S, and points corresponding to the vectors are connected with each other.
  • the selected vectors are the same in magnitude, but are different in direction.
  • the points corresponding to the selected vectors are connected with each other, thereby generating a sub-contour line R.
  • the vectors having the determined magnitude (1/3) are selected in the region between the contour line C2 and the center line S, and points corresponding to the vectors are connected.
  • the selected vectors are the same in magnitude, but are different in direction.
  • the points corresponding to the selected vectors are connected with each other, thereby generating a sub-contour line T.
  • Cubic spline curves may be used for the sub-contour lines R and T, for example.
  • the processing unit 140 determines terrain deviations of the sub-contour lines by executing the terrain deviation determination sub-program 167.
  • An accumulated change in magnitude of the vectors on the contour line C1 to the vectors on the contour line C2 via the center line S is 2.
  • the total difference in magnitude of the selected vectors is 2, and hence the terrain deviation between the sub-contour lines is a value obtained by dividing the terrain deviation between the contour lines C1 and C2 by 2 and multiplying the resultant quotient by the difference in magnitude of the vectors between the sub-contour lines.
  • the terrain deviation (hR-hT) between the newly provided two sub-contour lines R and T is calculated as what is obtained by dividing a terrain deviation 100 m between the contour lines C1 and C2 by 2 and multiplying the resultant quotient by the difference in magnitude 2/3 between the vectors of the sub-contour lines, namely 33.3 (100/3) m.
  • the sub-contour line R represents the elevation 166.7 m
  • the sub-contour line T represents the elevation 133.3 m.
  • a slope of an actual terrain may not be monotonic, and it is necessary to determine whether the interpolated sub-contour line is higher or lower than the existing contour lines C1 and C2.
  • the processing unit 140 calculates the terrain deviation at the boundary of the flattened region, and determines whether the elevation of the sub-contour line is higher or lower than the elevation of the neighboring contour lines according to the sign of the calculated terrain deviation by executing the terrain deviation determination sub-program 167.
  • a case of obtaining the elevation h of the sub-contour line S generated between the contour lines C1 and C2 is considered. If the elevation h1 of the contour line C1 is higher than the elevation h0 of the contour line C0, this case implies that the elevation increases from the contour line C0 toward the contour line C1, and the elevation h of the sub-contour line S is h1+0.5*d. On the other hand, if the elevation h1 of the contour line C1 is lower than the elevation h0 of the contour line C0, this case implies that the elevation decreases from the contour line C0 toward the contour line C1, and the elevation h of the sub-contour line S is h1-0.5*d.
  • the elevation h1 of the contour line C1 is equal to the elevation h0 of the contour line C0, there is a ridge, a valley, or a plain between the contour lines C0 and C1.
  • an operator determines whether the elevation h of the sub-contour line S is h1+0.5*d or h1-0.5*d based on other information.
  • the processing unit 140 connects the generated sub-contour lines to each other by executing the sub-contour line connecting sub-program 166 (S6).
  • S6 sub-contour line connecting sub-program 166
  • this processing of connecting the sub-contour lines is now described.
  • the different methods are used for generating the sub-contour lines between the flattened region and the non-flattened region, and the generated sub-contour lines are not aligned at a boundary between the flattened region and the non-flattened region. It is therefore necessary to smoothly link the sub-contour line in the flattened region and the sub-contour line in the non-flattened region.
  • a rectangular buffer 802 is provided close to the boundary of the flattened region and the non-flattened region.
  • the size d of one side of the buffer 802 in the direction of the extension of the sub-contour lines is preferably determined according to the magnitude of the displacement between the sub-contour line of the flattened region and the sub-contour line of the non-flattened region.
  • the size d of the buffer 802 is preferably increased, and if the displacement between the sub-contour line of the flattened region and the sub-contour line of the non-flattened region is small, the size d of the buffer 802 may be decreased.
  • intersection point P1 between a sub-contour line S1 of the flattened region and the boundary of the buffer 802 is set, a center line 803 in the flattened region in the buffer 802 is set, and an intersection point A1 between the center line 803 and the sub-contour line S1 is set.
  • intersection points P2 and P3 and intersection points A2 and A3 are set on sub-contour lines S2 and S3.
  • intersection point Q1 between the sub-contour line T1 of the non-flattened region and the boundary of the buffer 802 is set, a center line 804 in the flattened region in the buffer 802 is set, and an intersection point B1 between the center line 804 and the sub-contour line T1 is set.
  • intersection points Q2 and Q3 and intersection points B2 and B3 are set on sub-contour lines T2 and T3.
  • a curve U1 connecting the four points, P1, A1, B1, and Q1, with each other is generated.
  • a cubic spline curve may be used for this curve, for example.
  • a newly integrated sub-contour line is generated by using the sub-contour line S1 up to the point P1 in the flattened region, the cubic curve U1 in the buffer 802 (between the point P1 and the point Q1), and the sub-contour line T1 up to the point Q1 in the non-flattened region.
  • sub-contour lines S2 and T2 are connected through a cubic curve U2
  • sub-contour lines S3 and T3 are connected through a cubic curve U3, thereby generating newly integrated sub-contour lines.
  • the DEM may be vector data or image data in the form of a map containing the generated sub-contour lines.
  • This embodiment can generate a highly-accurate DEM.
  • accurate sub-contour lines can be interpolated by eliminating a determination that a saddle region or a gentle slope of a ridge in a contour map is determined as flattened by the TIN model.
  • the method according to this invention can thus be adapted to many applications for 3D terrain visualization.
  • the effect vector field generated by the embodiment of this invention represents the effects of two contour lines as gentle changes between the two contour lines as illustrated in Fig. 9.
  • smooth sub-contour lines can be accurately interpolated.
  • the elevation of any point between two contour lines can be determined without loss of continuity in elevation of surrounding points.
  • sub-contour lines can be calculated in a triangle patch considered as flattened in place of the conventional method employing the TIN model. As a result, a highly-accurate DEM can be generated.
  • this invention can be applied to any region (non-flattened region) partitioned by existing contour lines in addition to the flattened region. This is because the effect vector field contains all effects from surrounding contour lines. Further, the elevation which is interpolated according to this invention provides smoother interpolation at every location than the conventional method employing the TIN model.
  • an optimal balance between efficiency and accuracy is achieved to generate an accurate DEM by simultaneously employing the conventional method employing the TIN model and the method according to this invention on a client terminal.
  • limitation on computer resources can be relieved and a large and accurate DEM can be generated in a batch manner, thereby supplying high quality DEM data to the client terminal.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • Educational Technology (AREA)
  • General Physics & Mathematics (AREA)
  • Image Generation (AREA)

Abstract

L'invention concerne un système de génération de modèle numérique d'altitude, qui comprend : un module de réception destiné à segmenter la carte de contour d'entrée en un réseau irrégulier triangulé ; un module de marquage destiné à marquer chacune des zones triangulaires comme planes ou non planes ; un module de calcul de structure triangulaire TIN destiné à générer des lignes de contour secondaire par interpolation, à l'aide du réseau irrégulier triangulé dans la région non plane ; un module de génération de champ de vecteur destiné à calculer un champ de vecteur, dans lequel chaque vecteur représente un effet spatial de lignes de contour voisines pour la région plane ; un module de génération de ligne de contour secondaire destiné à déterminer une courbe sur laquelle les vecteurs ont la même amplitude ; et un module de connexion de ligne de contour secondaire destiné à générer une ligne de contour secondaire nouvellement intégrée par la création d'un tampon dans un voisinage d'une limite entre une région non plane et la région plane et par liaison lisse de la courbe dans la région non plane et la courbe dans la région plane dans le tampon.
PCT/JP2011/005838 2011-10-19 2011-10-19 Système de génération de modèle numérique d'altitude et procédé de génération d'un modèle numérique d'altitude WO2013057760A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/JP2011/005838 WO2013057760A1 (fr) 2011-10-19 2011-10-19 Système de génération de modèle numérique d'altitude et procédé de génération d'un modèle numérique d'altitude

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2011/005838 WO2013057760A1 (fr) 2011-10-19 2011-10-19 Système de génération de modèle numérique d'altitude et procédé de génération d'un modèle numérique d'altitude

Publications (1)

Publication Number Publication Date
WO2013057760A1 true WO2013057760A1 (fr) 2013-04-25

Family

ID=48140441

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2011/005838 WO2013057760A1 (fr) 2011-10-19 2011-10-19 Système de génération de modèle numérique d'altitude et procédé de génération d'un modèle numérique d'altitude

Country Status (1)

Country Link
WO (1) WO2013057760A1 (fr)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103985249A (zh) * 2014-06-03 2014-08-13 北京航空航天大学 道路的高程信息提取方法及装置
CN104282221A (zh) * 2013-07-01 2015-01-14 北京四维图新科技股份有限公司 一种构建室内地图的道路网络的方法
CN107154070A (zh) * 2016-03-04 2017-09-12 高德软件有限公司 矢量要素与数字地面模型叠加方法及装置
WO2018080969A1 (fr) * 2016-10-24 2018-05-03 Agco International Gmbh Système de cartographie et de guidage terrestres
CN109523630A (zh) * 2016-06-21 2019-03-26 武汉瑞得信息工程有限责任公司 一种基于张力样条插值的等高线生成方法及装置
CN109684424A (zh) * 2018-11-09 2019-04-26 广西壮族自治区遥感信息测绘院 基于离散特征线的地貌数据快速生成与优化方法
CN110955742A (zh) * 2019-11-28 2020-04-03 中国电建集团成都勘测设计研究院有限公司 基于dem的升压站场坪雨水口设置区域分析方法
CN110968656A (zh) * 2018-09-30 2020-04-07 北京金风科创风电设备有限公司 地形图处理方法及处理装置
CN111951393A (zh) * 2020-07-27 2020-11-17 南京师范大学 水平构造台地的三维模型构造方法及装置
CN112116709A (zh) * 2020-09-17 2020-12-22 滁州学院 一种提高地形表达精度的地形特征线处理方法
CN112861639A (zh) * 2021-01-14 2021-05-28 黄河勘测规划设计研究院有限公司 一种基于无人机遥感的河道变迁监测方法
CN113450456A (zh) * 2021-06-28 2021-09-28 中国水利水电科学研究院 高可识别度的dem制作方法
CN114329575A (zh) * 2022-03-15 2022-04-12 武汉峰岭科技有限公司 基于内插山脊线的等高线加密方法、系统、介质及设备
CN115113228A (zh) * 2022-05-09 2022-09-27 江苏省水利科学研究院 一种基于地理信息技术的退圩还湖工程测验方法
CN111127649B (zh) * 2019-12-30 2023-04-14 重庆市勘测院 构建三维体块模型的方法及装置、服务器

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003150950A (ja) * 2001-11-16 2003-05-23 Asia Air Survey Co Ltd 地形モデルの補間方法、補間プログラム、および補間プログラムを記憶した記憶媒体
JP2003216023A (ja) * 2002-01-24 2003-07-30 Takashi Watanabe 地形図情報を用いた数値標高モデル自動生成システム
JP2005106906A (ja) * 2003-09-29 2005-04-21 Hitachi Software Eng Co Ltd 標高データ生成方法及び標高データ生成プログラム
JP2009163317A (ja) * 2007-12-28 2009-07-23 Asahi Koyo Kk 地形モデル補間装置、方法及びプログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003150950A (ja) * 2001-11-16 2003-05-23 Asia Air Survey Co Ltd 地形モデルの補間方法、補間プログラム、および補間プログラムを記憶した記憶媒体
JP2003216023A (ja) * 2002-01-24 2003-07-30 Takashi Watanabe 地形図情報を用いた数値標高モデル自動生成システム
JP2005106906A (ja) * 2003-09-29 2005-04-21 Hitachi Software Eng Co Ltd 標高データ生成方法及び標高データ生成プログラム
JP2009163317A (ja) * 2007-12-28 2009-07-23 Asahi Koyo Kk 地形モデル補間装置、方法及びプログラム

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104282221A (zh) * 2013-07-01 2015-01-14 北京四维图新科技股份有限公司 一种构建室内地图的道路网络的方法
CN103985249A (zh) * 2014-06-03 2014-08-13 北京航空航天大学 道路的高程信息提取方法及装置
CN107154070A (zh) * 2016-03-04 2017-09-12 高德软件有限公司 矢量要素与数字地面模型叠加方法及装置
CN109523630A (zh) * 2016-06-21 2019-03-26 武汉瑞得信息工程有限责任公司 一种基于张力样条插值的等高线生成方法及装置
CN109523631A (zh) * 2016-06-21 2019-03-26 武汉瑞得信息工程有限责任公司 一种基于b样条拟合的等高线生成方法及装置
CN109523630B (zh) * 2016-06-21 2022-07-29 武汉瑞得信息科技有限责任公司 一种基于张力样条插值的等高线生成方法及装置
US11085776B2 (en) 2016-10-24 2021-08-10 Agco International Gmbh Land mapping and guidance system
WO2018080969A1 (fr) * 2016-10-24 2018-05-03 Agco International Gmbh Système de cartographie et de guidage terrestres
CN110968656B (zh) * 2018-09-30 2023-05-26 北京金风科创风电设备有限公司 地形图处理方法及处理装置
CN110968656A (zh) * 2018-09-30 2020-04-07 北京金风科创风电设备有限公司 地形图处理方法及处理装置
CN109684424A (zh) * 2018-11-09 2019-04-26 广西壮族自治区遥感信息测绘院 基于离散特征线的地貌数据快速生成与优化方法
CN109684424B (zh) * 2018-11-09 2023-03-14 广西壮族自治区遥感信息测绘院 基于离散特征线的地貌数据快速生成与优化方法
CN110955742B (zh) * 2019-11-28 2023-07-25 中国电建集团成都勘测设计研究院有限公司 基于dem的升压站场坪雨水口设置区域分析方法
CN110955742A (zh) * 2019-11-28 2020-04-03 中国电建集团成都勘测设计研究院有限公司 基于dem的升压站场坪雨水口设置区域分析方法
CN111127649B (zh) * 2019-12-30 2023-04-14 重庆市勘测院 构建三维体块模型的方法及装置、服务器
CN111951393B (zh) * 2020-07-27 2024-02-27 南京师范大学 水平构造台地的三维模型构造方法及装置
CN111951393A (zh) * 2020-07-27 2020-11-17 南京师范大学 水平构造台地的三维模型构造方法及装置
CN112116709B (zh) * 2020-09-17 2023-07-07 滁州学院 一种提高地形表达精度的地形特征线处理方法
CN112116709A (zh) * 2020-09-17 2020-12-22 滁州学院 一种提高地形表达精度的地形特征线处理方法
CN112861639A (zh) * 2021-01-14 2021-05-28 黄河勘测规划设计研究院有限公司 一种基于无人机遥感的河道变迁监测方法
CN113450456A (zh) * 2021-06-28 2021-09-28 中国水利水电科学研究院 高可识别度的dem制作方法
CN113450456B (zh) * 2021-06-28 2023-10-03 中国水利水电科学研究院 高可识别度的dem制作方法
CN114329575A (zh) * 2022-03-15 2022-04-12 武汉峰岭科技有限公司 基于内插山脊线的等高线加密方法、系统、介质及设备
CN115113228A (zh) * 2022-05-09 2022-09-27 江苏省水利科学研究院 一种基于地理信息技术的退圩还湖工程测验方法
CN115113228B (zh) * 2022-05-09 2023-10-24 江苏省水利科学研究院 一种基于地理信息技术的退圩还湖工程测验方法

Similar Documents

Publication Publication Date Title
WO2013057760A1 (fr) Système de génération de modèle numérique d'altitude et procédé de génération d'un modèle numérique d'altitude
EP1927954B1 (fr) Modélisation de terrain basée sur une zone de surface courbe
US7932904B2 (en) Generation of constrained voronoi grid in a plane
US8368714B2 (en) Curved surface rendering system and method
US9324186B2 (en) Method for representing terrain, method for creating terrain primitives, and apparatus using the methods
CN111581776B (zh) 一种基于几何重建模型的等几何分析方法
US20140257750A1 (en) System and method for estimating uncertainty for geophysical gridding routines lacking inherent uncertainty estimation
Fisher et al. Semantically correct 2.5 D GIS data—The integration of a DTM and topographic vector data
Pfeifer A subdivision algorithm for smooth 3D terrain models
CN108563915B (zh) 车辆数字化仿真测试模型构建系统及方法、计算机程序
US9188435B2 (en) Method for generating error image and program for generating error image
JP2011033987A (ja) 地図生成装置及び地図生成方法
CN116933607A (zh) 一种矿山资源模型的创建方法、装置、设备及存储介质
KR101814023B1 (ko) 유한차분격자 자료 자동 보정 장치 및 방법
Gold Tessellations in GIS: Part II–making changes
JP4543820B2 (ja) ３次元データ処理装置及びプログラム
JP2007264952A (ja) 地盤解析用メッシュ生成方法及び地盤解析用メッシュ生成プログラム
Lesparre et al. Simplified 3D city models from LiDAR
US20150149127A1 (en) Methods and Systems to Synthesize Road Elevations
JP2020071501A (ja) 地形推定プログラム、地形推定方法、及び、地形推定装置
Petras et al. Mapping gradient fields of landform migration
JP7122901B2 (ja) ブレイクライン抽出プログラム、ブレイクライン抽出システム
CN114329575B (zh) 基于内插山脊线的等高线加密方法、系统、介质及设备
Šljivarić et al. Methods of Modelling the Distortion Caused by Different Amount and Orientation of Coordinate Corrections between Two Coordinate Systems at Various Locations
RU2415381C1 (ru) Способ картографического отображения двумерных распределений, заданных в цифровой форме

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: 11874434

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: 11874434

Country of ref document: EP

Kind code of ref document: A1