CN116045920A - Method for extracting DEM based on GF7 stereopair - Google Patents
Method for extracting DEM based on GF7 stereopair Download PDFInfo
- Publication number
- CN116045920A CN116045920A CN202211604277.6A CN202211604277A CN116045920A CN 116045920 A CN116045920 A CN 116045920A CN 202211604277 A CN202211604277 A CN 202211604277A CN 116045920 A CN116045920 A CN 116045920A
- Authority
- CN
- China
- Prior art keywords
- image
- dem
- double
- images
- stereoscopic
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
- G01C11/04—Interpretation of pictures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
- G01C11/02—Picture taking arrangements specially adapted for photogrammetry or photographic surveying, e.g. controlling overlapping of pictures
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/97—Determining parameters from multiple pictures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Software Systems (AREA)
- Radar, Positioning & Navigation (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Processing Or Creating Images (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a method for extracting DEM based on GF7 stereopair. Step 1: obtaining GF7 stereopair remote sensing satellite images; step 2: constructing a adjustment regional network and completing regional network adjustment to realize high-precision orientation of all images; step 3: constructing a plurality of pairs of double-image stereoscopic images by utilizing index values such as the overlapping rate between the front view and the rear view of the GF7 stereoscopic image and the intersection angle of homonymous rays; step 4: extracting a Digital Elevation Model (DEM) by adopting an image dense matching algorithm; step 5: calculating the weight value of each double-image stereoscopic image through the intersection angle, the relative precision and the time phase difference index of the constituent images of the double-image stereoscopic images; step 6: and carrying out weighted fusion on the basis of the weight value to generate a fused DEM product, and editing and correcting the product to generate a high-quality and high-precision double-image stereoscopic relative DEM product. The invention greatly improves the resolution of the DEM and the DEM matching quality of areas with complex topography, detail and the like.
Description
Technical Field
The invention relates to the technical field of topographic and geomorphic production in mapping geographic information products, in particular to a method for extracting DEM (digital elevation model) based on GF7 stereopair.
Background
With the continuous development of national economy, mapping geographic information resources plays an important role in government management decisions, industrial development, people living standard and the like. The DEM is a physical ground model for realizing digital simulation of ground topography through limited topography elevation data and representing ground elevation in the form of a group of ordered value arrays. The DEM is used as an important geographical information mapping result form and is widely applied to the fields of national economy and national defense construction, such as mapping, hydrology, weather, landform, geology, soil, engineering construction, communication, military and the like, and humanity and natural science.
In recent years, obtaining DEM by using remote sensing satellite optical stereo image and space photography measurement means has become an important production method of high-quality DEM products in a large range. The method comprises the following main steps of: and (3) three-dimensional image area network adjustment, three-dimensional image three-dimensional matching extraction DEM, DEM editing and the like. The satellite stereo image is constructed by two or more remote sensing satellite plane images of the same region, which are taken from different imaging angles. Currently, theories, techniques and methods for developing DEM production using two-image stereoscopic images constructed from two images have gradually tended to be mature. In recent years, however, mapping remote sensing satellites capable of acquiring stereoscopic images are sequentially transmitted (such as an ALOS satellite in japan, a satellite series No. three, a satellite No. GF7, a satellite No. gilin, etc.), stereoscopic relative satellite images have become important stereoscopic mapping application image data sources, and provide high-precision and high-quality DEM products for topographic and geomorphic products in the geographic mapping industry.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention aims to provide a method for extracting a DEM (digital elevation model) based on a GF7 stereopair, which utilizes indexes such as the overlapping rate between a front view and a rear view, the intersection angle of homonymous rays and the like, and provides a method for weighting and fusing a DEM product by utilizing a high-resolution remote sensing satellite stereopsis, so that the situation that the DEM on a local image cannot be successfully matched due to terrain shielding or shadow and the like can be effectively reduced, and the resolution of the DEM, the DEM matching quality of areas with complex terrain details and the like can be greatly improved.
In order to achieve the above purpose, the technical scheme adopted by the invention is as follows: a method for extracting DEM based on GF7 stereopair, comprising the steps of:
step 1: acquiring front-rear vision full-color image data of two GF7 stereopair remote sensing satellite images in the same orbit;
step 2: constructing a leveling regional network for the GF7 stereopair remote sensing satellite images and finishing regional network leveling, so as to realize the relative and absolute orientation of high-precision quality of all images;
step 3: constructing a plurality of pairs of double-image stereoscopic images from the GF7 stereoscopic image by utilizing the overlapping rate between the front view and the rear view of the GF7 stereoscopic image and index values such as the intersection angle of homonymous rays;
step 4: generating an approximate epipolar line stereoscopic image by a double-image stereoscopic image firstly by adopting a projection track method, then generating a parallax image by adopting a dense matching algorithm, secondly acquiring the coordinates of each pixel ground object by utilizing an approximate epipolar line stereoscopic image imaging geometric model, generating three-dimensional point cloud data of the ground object in a task area, and finally generating a regular grid DEM product by point cloud rasterization;
step 5: calculating the weight value of each double-image stereoscopic image through indexes such as the intersection angle, the relative precision, the time phase difference forming the image and the like of the double-image stereoscopic image;
step 6: and based on the weight value, carrying out weighted fusion on the DEM extracted from all the double-image stereoscopic images in the region to generate a fused DEM product of the whole task region, and editing and correcting the fused DEM product to generate a high-quality and high-precision double-image stereoscopic relative DEM product. According to the invention, the resolution of the DEM and the DEM matching quality of areas with complex topography and detail and the like are effectively and greatly improved through the high-resolution stereopair remote sensing satellite images.
And in the step 1, front and rear vision full-color image data of two same-orbit GF7 stereopair remote sensing satellite images are obtained.
In the step 2, a leveling regional network is constructed for the remote sensing satellite images of the GF7 stereopair, the dense connection points in the leveling regional network are measured firstly, and a proper amount of image control points are measured according to the requirement; secondly, selecting a proper regional network adjustment model, uniformly carrying out adjustment iterative computation in the whole region until adjustment precision requirements are met, and acquiring error compensation parameters of all image imaging geometric models; and finally, outputting a new image imaging geometric model to finish high-precision relative and absolute orientation of the satellite images in the region.
In the step 3, the GF7 stereoscopic image is used to construct a plurality of pairs of double-image stereoscopic images by using the index values such as the overlapping rate between the front view and the rear view of the GF7 stereoscopic image and the intersection angle of the same name light.
In the step 4, through the double-image stereoscopic image, firstly, a projection track method is adopted to generate an approximate epipolar stereoscopic image, then, a dense matching algorithm is adopted to generate a parallax image, secondly, an approximate epipolar stereoscopic image imaging geometric model is utilized to perform intersection to acquire the coordinates of each pixel ground object, three-dimensional point cloud data of the ground object in a task area are generated, and finally, a regular grid DEM product is generated through point cloud rasterization.
In the step 5, the weight value of each double-image stereoscopic image is calculated by the indexes such as the intersection angle, the relative precision, the phase difference of the constituent images and the like of the double-image stereoscopic image.
And step 6, based on the weight value, carrying out weighted fusion on the DEM extracted from all the double-image stereoscopic images in the region to generate a fused DEM product of the whole task region, and editing and correcting the fused DEM product to generate a high-quality and high-precision double-image stereoscopic relative DEM product.
Compared with the prior art, the invention has the following beneficial effects:
more landform information in the region can be obtained through the GF7 stereoscopic image, and the method adopts the GF7 stereoscopic image to develop the DEM generation, so that the risk that the local region cannot be successfully matched with the DEM due to relief shielding and the like can be effectively reduced, and the DEM matching success rate of the whole region is improved; the information redundancy advantage of the multi-degree overlapped coverage stereoscopic image is fully utilized, the same-track stereoscopic image and different-track stereoscopic image in the task area are split into a plurality of double-image stereoscopic images, then the two images are respectively matched with the DEM, and then the weighted fusion is carried out to generate a total DEM, so that the resolution of a DEM product can be effectively improved, and the DEM matching quality of the area with complex landform detail can also be improved.
Drawings
Fig. 1 is a flowchart provided in an embodiment of the present invention.
Fig. 2 is a high quality, high precision biparatopic stereoscopic DEM outcome map.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the following embodiments are used to further describe a method for extracting DEM based on GF7 stereopair according to the present invention in detail with reference to the accompanying drawings.
Fig. 1 is a flowchart of a method for extracting DEM based on GF7 stereopair, mainly comprising the steps of:
1) Step 1: acquiring front-rear vision full-color image data of two GF7 stereopair remote sensing satellite images in the same orbit;
2) Step 2: constructing a leveling regional network for the GF7 stereopair remote sensing satellite images and finishing regional network leveling, so as to realize the relative and absolute orientation of high-precision quality of all images;
3) Step 3: constructing a plurality of pairs of double-image stereoscopic images from the GF7 stereoscopic image by utilizing the overlapping rate between the front view and the rear view of the GF7 stereoscopic image and index values such as the intersection angle of homonymous rays;
4) Step 4: generating an approximate epipolar line stereoscopic image by a double-image stereoscopic image firstly by adopting a projection track method, then generating a parallax image by adopting a dense matching algorithm, secondly acquiring the coordinates of each pixel ground object by utilizing an approximate epipolar line stereoscopic image imaging geometric model, generating three-dimensional point cloud data of the ground object in a task area, and finally generating a regular grid DEM product by point cloud rasterization;
5) Step 5: calculating the weight value of each double-image stereoscopic image through indexes such as the intersection angle, the relative precision, the time phase difference forming the image and the like of the double-image stereoscopic image;
6) Step 6: and based on the weight value, carrying out weighted fusion on the DEM extracted from all the double-image stereoscopic images in the region to generate a fused DEM product of the whole task region, and editing and correcting the fused DEM product to generate a high-quality and high-precision double-image stereoscopic relative DEM product. According to the invention, the resolution of the DEM and the DEM matching quality of areas with complex topography and detail and the like are effectively and greatly improved through the high-resolution stereopair remote sensing satellite images.
Further elaboration is made on steps 1-6:
and further, in the step 1, front and rear vision full-color image data of two same-orbit GF7 stereopair remote sensing satellite images are obtained.
In step 2, a adjustment regional network is constructed for the remote sensing satellite images of the GF7 stereopair, the dense connection points in the adjustment regional network are measured firstly, and a proper amount of image control points are measured according to the requirement; secondly, selecting a proper regional network adjustment model, uniformly carrying out adjustment iterative computation in the whole region until adjustment precision requirements are met, and acquiring error compensation parameters of all image imaging geometric models; and finally, outputting a new image imaging geometric model to finish high-precision relative and absolute orientation of the satellite images in the region.
In the step 3, the GF7 stereoscopic image is constructed into a plurality of pairs of double-image stereoscopic images by using the index values such as the overlapping rate between the front view and the rear view of the GF7 stereoscopic image and the intersection angle of the same name rays.
In the step 4, a projection trajectory method is adopted to generate an approximate epipolar line stereoscopic image, then a dense matching algorithm is adopted to generate a parallax image, then an approximate epipolar line stereoscopic image imaging geometric model is utilized to perform intersection to obtain coordinates of each pixel ground object, three-dimensional point cloud data of the ground object in a task area are generated, and finally a regular grid DEM product is generated through point cloud rasterization.
In step 5, the weight value of each double-image stereoscopic image is calculated by the indexes such as the intersection angle, the relative precision, the phase difference of the constituent images, and the like of the double-image stereoscopic image.
Further, in the step 6, based on the weight value, weighting and fusing all DEMs extracted from the double-image stereoscopic images in the region to generate a fused DEM product of the whole task region, and performing manual intervention, editing and correction on the fused DEM product; and finally, a high-quality and high-precision double-image stereoscopic relative DEM result diagram shown in fig. 2 is generated.
In summary, the method for weighting and fusing the DEM product by utilizing the high-resolution remote sensing satellite stereoscopic images is provided by utilizing indexes such as the overlapping rate between the front view and the rear view and the intersection angle of the same name rays, so that the situation that the DEM on the local image cannot be successfully matched due to terrain shielding or shadow and the like can be effectively reduced, and the resolution of the DEM and the DEM matching quality of areas with complex topography and detail and the like are greatly improved.
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential characteristics thereof. The present embodiments are, therefore, to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein.
Claims (7)
1. A method for extracting DEM based on GF7 stereopair, which is characterized in that: the method comprises the following steps:
step 1: acquiring front-rear vision full-color image data of two GF7 stereopair remote sensing satellite images in the same orbit;
step 2: constructing a leveling regional network for the GF7 stereopair remote sensing satellite images and finishing regional network leveling, so as to realize the relative and absolute orientation of high-precision quality of all images;
step 3: constructing a plurality of pairs of double-image stereoscopic images from the GF7 stereoscopic image by utilizing the overlapping rate between the front view and the rear view of the GF7 stereoscopic image and the index value of the homonymous ray intersection angle;
step 4: generating an approximate epipolar line stereoscopic image by a double-image stereoscopic image firstly by adopting a projection track method, then generating a parallax image by adopting a dense matching algorithm, secondly acquiring the coordinates of each pixel ground object by utilizing an approximate epipolar line stereoscopic image imaging geometric model, generating three-dimensional point cloud data of the ground object in a task area, and finally generating a regular grid DEM product by point cloud rasterization;
step 5: calculating the weight value of each double-image stereoscopic image through the intersection angle, the relative precision and the time phase difference index of the constituent images of the double-image stereoscopic images;
step 6: and based on the weight value, carrying out weighted fusion on the DEM extracted from all the double-image stereoscopic images in the region to generate a fused DEM product of the whole task region, and editing and correcting the fused DEM product to generate a high-quality and high-precision double-image stereoscopic relative DEM product.
2. The method for extracting DEM based on GF7 stereopair according to claim 1, wherein:
and in the step 1, front and rear vision full-color image data of two same-orbit GF7 stereopair remote sensing satellite images are obtained.
3. The method for extracting DEM based on GF7 stereopair according to claim 1, wherein:
in the step 2, a leveling regional network is constructed for the remote sensing satellite images of the GF7 stereopair, the dense connection points in the leveling regional network are measured firstly, and a proper amount of image control points are measured according to the requirement; secondly, selecting a proper regional network adjustment model, uniformly carrying out adjustment iterative computation in the whole region until adjustment precision requirements are met, and acquiring error compensation parameters of all image imaging geometric models; and finally, outputting a new image imaging geometric model to finish high-precision relative and absolute orientation of the satellite images in the region.
4. The method for extracting DEM based on GF7 stereopair according to claim 1, wherein:
in the step 3, the GF7 stereoscopic image is used to construct a plurality of pairs of double-image stereoscopic images by using the index values such as the overlapping rate between the front view and the rear view of the GF7 stereoscopic image and the intersection angle of the same name light.
5. The method for extracting DEM based on GF7 stereopair according to claim 1, wherein:
in the step 4, through the double-image stereoscopic image, firstly, a projection track method is adopted to generate an approximate epipolar stereoscopic image, then, a dense matching algorithm is adopted to generate a parallax image, secondly, an approximate epipolar stereoscopic image imaging geometric model is utilized to perform intersection to acquire the coordinates of each pixel ground object, three-dimensional point cloud data of the ground object in a task area are generated, and finally, a regular grid DEM product is generated through point cloud rasterization.
6. The method for extracting DEM based on GF7 stereopair according to claim 1, wherein:
in the step 5, the weight value of each double-image stereoscopic image is calculated by the indexes such as the intersection angle, the relative precision, the phase difference of the constituent images and the like of the double-image stereoscopic image.
7. The method for extracting DEM based on GF7 stereopair according to claim 1, wherein:
and step 6, based on the weight value, carrying out weighted fusion on the DEM extracted from all the double-image stereoscopic images in the region to generate a fused DEM product of the whole task region, and editing and correcting the fused DEM product to generate a final DEM product.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211604277.6A CN116045920A (en) | 2022-12-14 | 2022-12-14 | Method for extracting DEM based on GF7 stereopair |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211604277.6A CN116045920A (en) | 2022-12-14 | 2022-12-14 | Method for extracting DEM based on GF7 stereopair |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116045920A true CN116045920A (en) | 2023-05-02 |
Family
ID=86115394
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211604277.6A Pending CN116045920A (en) | 2022-12-14 | 2022-12-14 | Method for extracting DEM based on GF7 stereopair |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116045920A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117649611A (en) * | 2024-01-30 | 2024-03-05 | 西安宇速防务集团有限公司 | DEM data production processing method based on twice orientation |
-
2022
- 2022-12-14 CN CN202211604277.6A patent/CN116045920A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117649611A (en) * | 2024-01-30 | 2024-03-05 | 西安宇速防务集团有限公司 | DEM data production processing method based on twice orientation |
| CN117649611B (en) * | 2024-01-30 | 2024-04-30 | 西安宇速防务集团有限公司 | DEM data production processing method based on twice orientation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111126148B (en) | DSM (digital communication system) generation method based on video satellite images | |
| CN102506824B (en) | Method for generating digital orthophoto map (DOM) by urban low altitude unmanned aerial vehicle | |
| CN113358091B (en) | Method for producing digital elevation model DEM (digital elevation model) by using three-linear array three-dimensional satellite image | |
| CN103198524B (en) | A kind of three-dimensional reconstruction method for large-scale outdoor scene | |
| CN108168521A (en) | One kind realizes landscape three-dimensional visualization method based on unmanned plane | |
| CN104931022A (en) | Satellite image three-dimensional area network adjustment method based on satellite-borne laser height measurement data | |
| CN105160702A (en) | Stereoscopic image dense matching method and system based on LiDAR point cloud assistance | |
| Toutin | Error tracking of radargrammetric DEM from RADARSAT images | |
| CN110889899B (en) | Digital earth surface model generation method and device | |
| CN114241125B (en) | Multi-view satellite image-based fine three-dimensional modeling method and system | |
| US20230186561A1 (en) | Method for 3d reconstruction from satellite imagery | |
| CN114998399B (en) | Heterogeneous optical remote sensing satellite image stereopair preprocessing method | |
| CN116045920A (en) | Method for extracting DEM based on GF7 stereopair | |
| CN117576343B (en) | Three-dimensional MESH model manufacturing method based on high-resolution satellite stereoscopic image | |
| KR100732915B1 (en) | Method for three-dimensional determining of basic design road route using digital photommetry and satellite image | |
| CN111986320B (en) | Smart city application-oriented DEM and oblique photography model space fitting optimization method | |
| CN111611525B (en) | Remote sensing data elevation calculation method based on object space matching elevation deviation iterative correction | |
| CN112785686A (en) | Forest map construction method based on big data and readable storage medium | |
| CN113324527B (en) | Co-rail laser height measurement point and three-linear array three-dimensional image combined surveying and mapping processing method | |
| CN107784666A (en) | The detection of terrain and its features three dimensional change and update method based on stereopsis | |
| CN114252058A (en) | Aviation fine large scene manufacturing and measuring method | |
| Previtali et al. | An automatic multi-image procedure for accurate 3D object reconstruction | |
| Rosado | Introduction to applied photogrammetry and cartography for civil engineering | |
| Picu | Updating geospatial data by creating a high resolution digital surface model | |
| Wang et al. | Digital Surface Model Refinement Based on Projected Images |
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 |