CN101923729A - Reconstruction method of three-dimensional shape of lunar surface based on single gray level image - Google Patents
Reconstruction method of three-dimensional shape of lunar surface based on single gray level image Download PDFInfo
- Publication number
- CN101923729A CN101923729A CN 201010262508 CN201010262508A CN101923729A CN 101923729 A CN101923729 A CN 101923729A CN 201010262508 CN201010262508 CN 201010262508 CN 201010262508 A CN201010262508 A CN 201010262508A CN 101923729 A CN101923729 A CN 101923729A
- Authority
- CN
- China
- Prior art keywords
- menology
- pixel
- value
- dimensional shape
- reconstruction method
- 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.)
- Granted
Links
Images
Landscapes
- Processing Or Creating Images (AREA)
- Image Generation (AREA)
Abstract
The invention relates to a reconstruction method of a three-dimensional shape of a lunar surface based on a single gray level image, belonging to the technical field of computer simulation. The method comprises the following steps: firstly, extracting the gray level value information of each pixel point in the single gray level image of the lunar surface, and finding out the pixel point with the maximum gray level value; dividing the maximum gray level value by the gray level value of each pixel point in the image to acquire the cosine of the inclination angle of each pixel point; carrying out an approximate spherical hypothesis on the lunar surface to acquire the cosine of the deflection angle of each pixel point; computing the inclination angle and the deflection angle of each pixel point, then computing the surface normal vector of each pixel point, and converting the surface normal vector into the height coordinate value of each pixel point; carrying out filtering treatment; and reconstructing the three-dimensional shape of the lunar surface by utilizing the OpenGL technology according to each acquired surface normal vector and each acquired height coordinate value. The method of the invention can reduce the requirements for implementation of three-dimensional surface restoration, improve the speed and the efficiency of restoration and simplify the steps of restoration; and by using the method for carrying out treatment, the three-dimensional shape of the lunar surface can be more accurately reconstructed.
Description
Technical field
The present invention relates to menology reconstruction method of three-dimensional shape, belong to the computer simulation technique field based on individual gray level image.
Background technology
By menology 3D shape that the lunar map picture of taking is analyzed and researched is one of fundamental research of moon exploration.In the various countries' lunar exploration engineering that comprises China's " Chang'e I " and in the works, the research of menology dimensional topography is in critical positions all the time.Remain a kind of important menology dimensional topography analytical approach at present by image reconstruction menology 3D shape.Six kinds of contained measurements of the LRO of the U.S.'s in October, 2008 emission have two kinds in order to measure moon landform in equipping, and one of them takes the lunar map picture exactly.As far back as nineteen fifty-one, just find that menology has the characteristic of approximate lambert's body (Lambertian) reflection model before the human moonfall.The domestic nearest moon 3D shape restoration methods that has proposed based on binocular stereo vision, but there are two width of cloth image registration problems in the binocular stereo vision method.The single width gray level image light and shade that the seventies in 20th century, Horn proposed is recovered 3D shape, and (Shape-from-Shading, SFS) method obtains many scholars' attention.The principle of SFS algorithm is to calculate the intensity of reflected light of respective point on the three-dimensional object surface according to each point gray-scale value on the image, and by the analysis of physics and geometrical optics as can be known, the geometric configuration of catoptrical intensity and body surface character and body surface has relation.In order to eliminate the pathosis of problem, typical SFS algorithm has minimization algorithm, finite element algorithm, neural network algorithm and viscosity to separate theory etc., yet the ubiquitous problem of algorithm has following two at present: the one, and the reflection model of selecting does not meet the reflection characteristic of body surface; The 2nd, constraint condition of introducing and solution procedure are too complicated, and it is slow to find the solution speed, and efficient is low.
Summary of the invention
The purpose of this invention is to provide a kind of menology reconstruction method of three-dimensional shape, find the solution slow, the inefficient problem of speed to solve existing reconstructing method based on individual gray level image.
For achieving the above object, the menology reconstruction method of three-dimensional shape step based on individual gray level image of the present invention is as follows:
(1) extracts each gray values of pixel points information in individual menology gray level image, find out the wherein pixel of maximum gradation value;
(2) each gray values of pixel points and maximum gradation value in the image are divided by, obtain the cosine at each pixel inclination angle; According to the menology reflection characteristic, menology is done the almost spherical hypothesis, try to achieve the cosine of each pixel drift angle;
(3) according to the inclination angle cosine value of each pixel, the surface normal vector that the drift angle cosine value is tried to achieve each pixel;
(4) the surface normal vector of each pixel that will try to achieve is converted to the height coordinate value of each pixel;
(5) utilize medium filtering that each the height coordinate value that calculates is carried out Filtering Processing;
(6) according to the surface normal vector and the height coordinate value that obtain, utilize OpenGL technology reconstruct menology 3D shape.
Further, maximum gradation value is in the described step (1)
, the pairing pixel of this maximum gradation value is designated as
Further, pixel in described step (2) image
Gray-scale value
With maximum gradation value
Be divided by, obtain pixel
The inclination angle
Cosine
According to the menology reflection characteristic, menology is done the almost spherical hypothesis, try to achieve pixel
The drift angle
Cosine
, wherein
,
Be x for the gradation of image value, the calculus of differences of y axle obtains.
Further, described step (3) is according to pixel
Inclination angle cosine value, the surface normal vector that the drift angle cosine value is tried to achieve this pixel,
, wherein
Mould for vector.
Further, described step (4) is with the pixel that obtains
The surface normal vector be converted to the height coordinate value, formula is
Further, the medium filtering in the described step (5) is to establish the height sequence of sets
, with length be
Moving window it is carried out medium filtering, obtain the window correspondence
Number
, wherein,
Be the value of window center point,
With this
Number is pressed the size ordering of numerical value, and getting its sequence number is several conducts of middle
Filtering output, establish
Be filtering output, then
Further, also add illumination model after the reconstruct menology 3D shape in the described step (6).
Further, described illumination model is the Hapke illumination model, for
Further, the menology 3D shape stack texture to reconstruct shows.
The present invention is for emulation menology 3D shape, based on traditional SFS method, adopted the lambert's body illumination reflection model that meets the menology reflection characteristic, menology is done spherical hypothesis, then image is done the approximate differential computing and obtain height function, again the result is done medium filtering, the correctness that assurance is found the solution, thus realized utilizing the single width gray level image to recover the menology 3D shape.Adopt method of the present invention can reduce the requirement that three-dimensional surface recover to be implemented, improve speed and the efficient recovered, simplify the step of recovering, after handling through this method, the 3D shape of reconstruct menology more accurately.
Further, add the Hapke illumination model, because the Hapke model mainly is at the moon reflection case that is caused by single source, according to the geometrical property of light in the menology reflection, with the reflection coefficient of menology to light, calculate the radiation intensity that arrives the observation place by the light of menology reflection, better actual conditions of pressing close to menology.
Description of drawings
Fig. 1 is a process flow diagram of the present invention;
Fig. 2 is the moonscape emulation gray-scale map of embodiment;
Fig. 3 is to the cloud data display effect behind Fig. 1 reconstruct three-dimensional surface among the embodiment;
Fig. 4 is the displayed map of embodiment after to Fig. 1 reconstruct three-dimensional surface;
Fig. 5 is embodiment to Fig. 1 reconstruct 3D shape and adds the design sketch of Hapke illumination;
Fig. 6 be embodiment to Fig. 1 reconstruct 3D shape after the design sketch of stack Hapke illumination and texture;
Fig. 7 is the menology image that NASA takes;
Fig. 8 be embodiment to Fig. 7 reconstruct shape after the design sketch of stack texture and illumination.
Embodiment
The SFS algorithm is an intensity of reflected light of calculating respective point on the three-dimensional object surface according to each point gray-scale value on the image, and by the analysis of physics and geometrical optics as can be known, the geometric configuration of catoptrical intensity and body surface character and body surface has relation.
If will show the menology shape, the half-tone information in the image need be converted to the amount relevant with moon geological information.According to lambert's body reflection law, the diffuse reflection light intensity that desirable diffuse reflection body surface reflects is directly proportional with the cosine of the angle between incident light and the body surface normal vector, that is:
,
In the formula,
Be the intensity that diffuses of body surface,
Be the intensity of light source,
Be the body surface reflection coefficient,
Be the angle between incident light and the surface normal.
If the unit normal vector of body surface at illuminated some place is
, reflection spot to the vector of unit length of pointolite is
, the vector form that can be expressed as then:
,
When pointolite is far from illuminated surface,
Change very for a short time, thereby it can be seen as a constant.The light that this moment, pointolite sent is directional light, by vector
Definite fully.Under gradient space, suppose
Be light source incident vector,
Method vector for the body surface each point.For natural light, the incident vector is a definite value, and for the body surface normal vector, according to higher mathematics knowledge, provides functional form
, that is:
Hence one can see that
Be the object surfaces gradient.On mathematics, the dot product of two vectors equals the product of their moulds and included angle cosine, that is:
Can set up with any point on the gray level image by following formula
The intensity that diffuses at place is:
The intensity that diffuses of body surface
Being reflected on the image is exactly this dot image brightness.By following formula as can be known, when
Value is 1, when promptly light source direction is identical with surperficial direction of normal,
Be maximal value.Finding the solution the SFS problem is exactly according to following formula, utilizes the brightness of known image and light source direction to determine the body surface gradient
, and by
With
Relation can further obtain the object surfaces height
According to above-mentioned formula, each point has two surface graded components to be asked in the scene, and every had only a known gray scale, so it is the inverse process of a morbid state.
The process flow diagram of the menology reconstruction method of three-dimensional shape based on individual gray level image of the present invention as shown in Figure 1, specific embodiment is an example with Fig. 2 moonscape emulation gray level image, table 1 is part pixel gray-scale value tabulation in the image of Fig. 2.
1, establishes moon view data point
Gray-scale value be
, at first obtain the point of gray-scale value maximum in the view data
, its gray-scale value
2, will
The gray-scale value of point
With
Do division, its result is as the inclination angle cosine of this pixel
。As shown in table 2 below is the inclination angle cosine value data list that table 1 image slices vegetarian refreshments gray-scale value calculates.
According to the menology reflection characteristic, can do almost spherical hypothesis, the then drift angle of pixel to menology
Cosine can calculate by following formula,
, wherein
,
The representative image gray-scale value is x, and the calculus of differences of y axle is promptly differentiated to x, y.Inclination angle cosine value data list as shown in table 3 below, that table 1 image slices of serving as reasons vegetarian refreshments gray-scale value calculates.
3, obtain just can calculating the surface normal vector at this some place after the inclination angle and drift angle of pixel,
4, can draw the three-dimensional appearance of menology according to the surface normal vector, be to the cloud data display effect figure behind Fig. 1 reconstruct three-dimensional surface as shown in Figure 3; Fig. 4 is to the displayed map behind Fig. 1 reconstruct three-dimensional surface.The surface normal vector can only be described the features of shape of menology, can not obtain the relative height value of surface point, so, need carry out of the conversion of surface normal vector to the height coordinate value,
Height value data list as shown in table 4 below, that table 1 image slices of serving as reasons vegetarian refreshments gray-scale value recovers.
5, obtain altitude information after owing to there is the influence of picture noise, if especially imaging surface exists under the situation of spike sudden change, the height value that calculates has than mistake, so also need utilize medium filtering to carry out the processing of following filtering:
If sequence of sets
, with length be
Moving window it is carried out medium filtering, obtain the window correspondence
Number
, wherein,
Be the value of window center point,
With this
Number is pressed the size ordering of numerical value, and getting its sequence number is several conducts of middle
Filtering output, establish
Be filtering output, then
As shown in table 5 below, for table 1 picture point is recovered and filtered height value data list.
Table 5 table 1 picture point is recovered and filtered height value
6, pass through said method, utilize the OpenGL technology menology dimensional topography of emulation can be shown, but for effect is close to the real terrain of menology, also need to add illumination model, be to Fig. 1 reconstruct 3D shape as shown in Figure 5 and add the design sketch of Hapke illumination.Method has adopted the Hapke illumination model.The Hapke model mainly is at the moon reflection case that is caused by single source, according to the geometrical property of light in the menology reflection, with the reflection coefficient of menology, calculate the radiation intensity that arrives the observation place by the light of menology reflection, better actual conditions of pressing close to menology to light.The Hapke model can be expressed as:
Wherein
The albedo of expression celestial body surface atural object correspondence,
The condition that the expression viewpoint direction overlaps with light source direction,
Representation unit angle scattering equation,
The bidirectional reflectance coefficient of expression scattered beam in infinitely-great medium.
With
Can represent by the product of celestial body surface gradient vector respectively, promptly
Wherein
The unit vector of expression light source direction,
The unit vector of expression viewpoint direction.
Because when the reflection angle of light and incident angle were approaching, the direction vector of light source and the direction vector of viewpoint were very approaching, singularity can appear in the Hapke model at this moment, and bigger error can appear in result calculated.Therefore, adopt a kind of improved model:
, wherein
Be the normal vector of light source direction,
Be light source direction vector and
Be the vector of viewpoint direction,
,
,
,
,
,
,
The albedo of expression celestial body surface atural object correspondence.
6.1, improved model simplified the expression formula of Hapke model, make it singularity can not occur, and kept the fundamental property of Hapke model.
7, under this illumination model, scene stack texture is shown that effect is better true to nature.The stack texture is to utilize texture technology among the OpenGL, the general step that texture is drawn is: definition texture mapping, control texture, explanation texture mapping mode and definition texture coordinate, it realizes that function is glTexImage2D (), is the design sketch to stack Hapke illumination and texture after Fig. 1 reconstruct 3D shape as shown in Figure 6.
Shown in Figure 7 is the menology image that NASA takes; Fig. 8 be reconstructing method according to the present invention to Fig. 7 reconstruct shape after the design sketch of stack texture and illumination.
Claims (10)
1. menology reconstruction method of three-dimensional shape based on individual gray level image is characterized in that this method step is as follows:
(1) extracts each gray values of pixel points information in individual menology gray level image, find out the wherein pixel of maximum gradation value;
(2) each gray values of pixel points and maximum gradation value in the image are divided by, obtain the cosine at each pixel inclination angle; According to the menology reflection characteristic, menology is done the almost spherical hypothesis, try to achieve the cosine of each pixel drift angle;
(3) according to the inclination angle cosine value of each pixel, the surface normal vector that the drift angle cosine value is tried to achieve each pixel;
(4) the surface normal vector of each pixel that will try to achieve is converted to the height coordinate value of each pixel;
(5) utilize medium filtering that each the height coordinate value that calculates is carried out Filtering Processing;
(6) according to the surface normal vector and the height coordinate value that obtain, utilize OpenGL technology reconstruct menology 3D shape.
3. menology reconstruction method of three-dimensional shape according to claim 2 is characterized in that: pixel in described step (2) image
Gray-scale value
With maximum gradation value
Be divided by, obtain pixel
The inclination angle
Cosine
According to the menology reflection characteristic, menology is done the almost spherical hypothesis, try to achieve pixel
The drift angle
Cosine
, wherein
,
Be x for the gradation of image value, the calculus of differences of y axle obtains.
4. menology reconstruction method of three-dimensional shape according to claim 3 is characterized in that: described step (3) is according to pixel
Inclination angle cosine value, the surface normal vector that the drift angle cosine value is tried to achieve this pixel,
, wherein
Mould for vector.
6. menology reconstruction method of three-dimensional shape according to claim 5 is characterized in that: the medium filtering in the described step (5) is to establish the height sequence of sets
, with length be
Moving window it is carried out medium filtering, obtain the window correspondence
Number
, wherein,
Be the value of window center point,
With this
Number is pressed the size ordering of numerical value, and getting its sequence number is several conducts of middle
Filtering output, establish
Be filtering output, then
7. menology reconstruction method of three-dimensional shape according to claim 6 is characterized in that: also add illumination model after the reconstruct menology 3D shape in the described step (6).
10. menology reconstruction method of three-dimensional shape according to claim 9 is characterized in that: the menology 3D shape stack texture to reconstruct shows.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102625080A CN101923729B (en) | 2010-08-25 | 2010-08-25 | Reconstruction method of three-dimensional shape of lunar surface based on single gray level image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102625080A CN101923729B (en) | 2010-08-25 | 2010-08-25 | Reconstruction method of three-dimensional shape of lunar surface based on single gray level image |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101923729A true CN101923729A (en) | 2010-12-22 |
CN101923729B CN101923729B (en) | 2012-01-25 |
Family
ID=43338636
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102625080A Expired - Fee Related CN101923729B (en) | 2010-08-25 | 2010-08-25 | Reconstruction method of three-dimensional shape of lunar surface based on single gray level image |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101923729B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102331253A (en) * | 2011-08-09 | 2012-01-25 | 中国科学院西安光学精密机械研究所 | Moon-oriented high-resolution common-rail three-dimensional imaging method and device |
CN102928201A (en) * | 2012-10-24 | 2013-02-13 | 北京控制工程研究所 | Target simulating system of dynamic selenographic imaging sensor |
CN108352083A (en) * | 2015-11-06 | 2018-07-31 | 微软技术许可有限责任公司 | 2D image procossings for being drawn into 3D objects |
CN105869206B (en) * | 2016-04-12 | 2018-09-18 | 广州华欣电子科技有限公司 | three-dimensional rebuilding method and device |
CN108961391A (en) * | 2018-06-12 | 2018-12-07 | 温州大学激光与光电智能制造研究院 | A kind of surface reconstruction method based on curvature filtering |
CN109003332A (en) * | 2018-06-25 | 2018-12-14 | 重庆交通大学 | Bituminous pavement surface texture analogue system and its emulation mode |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6229913B1 (en) * | 1995-06-07 | 2001-05-08 | The Trustees Of Columbia University In The City Of New York | Apparatus and methods for determining the three-dimensional shape of an object using active illumination and relative blurring in two-images due to defocus |
CN1747559A (en) * | 2005-07-29 | 2006-03-15 | 北京大学 | Three-dimensional geometric mode building system and method |
CN101615072A (en) * | 2009-06-18 | 2009-12-30 | 东南大学 | Based on method for reproducing texture force touch from the shading shape technology |
CN101719287A (en) * | 2009-11-26 | 2010-06-02 | 上海大学 | Method for rebuilding shape of hemisphere three-dimensional surface with control point information |
-
2010
- 2010-08-25 CN CN2010102625080A patent/CN101923729B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6229913B1 (en) * | 1995-06-07 | 2001-05-08 | The Trustees Of Columbia University In The City Of New York | Apparatus and methods for determining the three-dimensional shape of an object using active illumination and relative blurring in two-images due to defocus |
CN1747559A (en) * | 2005-07-29 | 2006-03-15 | 北京大学 | Three-dimensional geometric mode building system and method |
CN101615072A (en) * | 2009-06-18 | 2009-12-30 | 东南大学 | Based on method for reproducing texture force touch from the shading shape technology |
CN101719287A (en) * | 2009-11-26 | 2010-06-02 | 上海大学 | Method for rebuilding shape of hemisphere three-dimensional surface with control point information |
Non-Patent Citations (1)
Title |
---|
《计算机与数字工程》 20080229 徐东平, 张忭 基于SFS方法的三维表面重建算法研究 7-9页 , 2 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102331253A (en) * | 2011-08-09 | 2012-01-25 | 中国科学院西安光学精密机械研究所 | Moon-oriented high-resolution common-rail three-dimensional imaging method and device |
CN102331253B (en) * | 2011-08-09 | 2013-03-27 | 中国科学院西安光学精密机械研究所 | Moon-oriented high-resolution common-rail three-dimensional imaging method and device |
CN102928201A (en) * | 2012-10-24 | 2013-02-13 | 北京控制工程研究所 | Target simulating system of dynamic selenographic imaging sensor |
CN102928201B (en) * | 2012-10-24 | 2014-12-17 | 北京控制工程研究所 | Target simulating system of dynamic selenographic imaging sensor |
CN108352083A (en) * | 2015-11-06 | 2018-07-31 | 微软技术许可有限责任公司 | 2D image procossings for being drawn into 3D objects |
CN108352083B (en) * | 2015-11-06 | 2022-06-07 | 微软技术许可有限责任公司 | 2D image processing for stretching into 3D objects |
CN105869206B (en) * | 2016-04-12 | 2018-09-18 | 广州华欣电子科技有限公司 | three-dimensional rebuilding method and device |
CN108961391A (en) * | 2018-06-12 | 2018-12-07 | 温州大学激光与光电智能制造研究院 | A kind of surface reconstruction method based on curvature filtering |
CN109003332A (en) * | 2018-06-25 | 2018-12-14 | 重庆交通大学 | Bituminous pavement surface texture analogue system and its emulation mode |
CN109003332B (en) * | 2018-06-25 | 2022-12-06 | 重庆交通大学 | Asphalt pavement surface texture simulation system and simulation method thereof |
Also Published As
Publication number | Publication date |
---|---|
CN101923729B (en) | 2012-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101923729B (en) | Reconstruction method of three-dimensional shape of lunar surface based on single gray level image | |
CN109035388B (en) | Three-dimensional face model reconstruction method and device | |
CN100594519C (en) | Method for real-time generating reinforced reality surroundings by spherical surface panoramic camera | |
CN107341853B (en) | Virtual-real fusion method and system for super-large virtual scene and dynamic screen shooting | |
Gallego et al. | A variational stereo method for the three-dimensional reconstruction of ocean waves | |
CN106169082A (en) | Training grader is with the method and system of the object in detection target environment image | |
CN108805946B (en) | Method and system for shading two-dimensional ultrasound images | |
CN103530907B (en) | Complicated three-dimensional model drawing method based on images | |
CN106846392A (en) | The method and apparatus of three-dimensional modeling | |
CN103268629B (en) | Unmarked some real time restoration method of 3 D human body form and attitude | |
CN111043988B (en) | Single stripe projection measurement method based on graphics and deep learning | |
CN106447763A (en) | Face image three-dimensional reconstruction method for fusion of sparse deformation model and principal component regression algorithm | |
CN104463968B (en) | The matching of remote sensing image binocular stereo vision and three-dimensional rebuilding method based on power grid GIS three-dimensional platform | |
CN104167011A (en) | Micro-structure surface global lighting drawing method based on direction light radiation intensity | |
Lyu et al. | Structured light-based underwater 3-D reconstruction techniques: A comparative study | |
CN108734772A (en) | High accuracy depth image acquisition methods based on Kinect fusion | |
Wu et al. | 3D film animation image acquisition and feature processing based on the latest virtual reconstruction technology | |
CN102722906B (en) | Feature-based top-down image modeling method | |
CN112150609A (en) | VR system based on indoor real-time dense three-dimensional reconstruction technology | |
Bonaventura et al. | Information measures for terrain visualization | |
KR100942026B1 (en) | Makeup system and method for virtual 3D face based on multiple sensation interface | |
CN104915992A (en) | A real-time shadow volume rendering method based on femur CT images | |
Ng et al. | Surface-from-gradients with incomplete data for single view modeling | |
Ryoo et al. | Simulation of 3D Virtual Space Model Based on Texture Image & Structure Analysis | |
van de Grind | Vision and the graphical simulation of spatial structure |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120125 Termination date: 20180825 |
|
CF01 | Termination of patent right due to non-payment of annual fee |