The computer virtual of remote-sensing distributed image on rolling ground
Technical field
The present invention relates to a kind of remote sensing and geographic information system technology, particularly utilize the method for the computer virtual of computer virtual structure remote-sensing distributed image on rolling ground.
Background technology
Because of landform blocking and shielding to direct sunlight, the direct sunlight Illumination Distribution of comparing on the natural terrain with wide level ground will be changed by the influence of topography, thereby changed atural object sensor information, influenced the quality and the quantitative examination of remote sensing digital image, the precision of application.Landform always is the difficult problem of satellite remote sensing digital image processing and application to the elimination of the influence of satellite remote sensing image and massif shade, clouds layer shadow.Existing satellite remote sensing digital image processing and application technology all fundamentally do not address this problem.This invention is under GIS supports, the method for the computer virtual by computer virtual structure remote-sensing distributed image on rolling ground has solved this problem preferably from theory and practice.
Summary of the invention
The objective of the invention is to provides a kind of computer virtual of remote-sensing distributed image on rolling ground at a difficult problem that exists in the prior art.
In order to reach the foregoing invention purpose, the present invention adopts following technical proposals:
A kind of computer virtual of remote-sensing distributed image on rolling ground is characterized in that:
At first, under the support of Geographic Information System, virtual different position of sun parameter and the geographic coordinate that utilizes each point on digital topography map and the digital topography map, calculate the virtual direct sunlight landform illumination coefficient of each point on the digital topography map by the computing machine pointwise, at last the virtual direct sunlight landform of the each point illumination coefficient that obtains is converted into virtual direct sunlight distribution plan, described position of sun parameter is meant the elevation angle and the position angle of the sun; Described digital topography map is meant the digital terrain model of defending sheet remote measurement digital topography map registration; The geographic coordinate of described each point is meant the corresponding longitude and latitude of each point on the digital topography map.Set up the mathematical model of quantitative relationship between satellite remote sensing information and landform, direct sunlight and the sky scattering light; The parsing satellite remote sensing digital image that continues: generate direct sunlight remote sensing images and sky scattering light remote sensing images, carry out topographic change respectively, the radiation difference on the elimination rugged topography and this species diversity are to the influence of satellite remote sensing digital image; Judge and eliminate the massif and the clouds layer shadow of remote sensing digital image then; Finish the computer virtual of remote-sensing distributed image on rolling ground on this basis at last;
Its concrete operations step is as follows:
(1), input research area digital topography map, the satellite remote sensing digital image in Study of Registration area and digital topography map, and be divided into grid;
(2), sun altitude, the position angle of each point on the digital topography map calculated in pointwise;
(3), the judgement of massif and clouds layer shadow on the calculating analog digital topomap;
(4), the ground direct sunlight landform illumination coefficient of each point on the virtual calculating digital topography map;
(5), the computing machine of direct sunlight distributed image generates with visual on the rolling ground.
Above-described computing machine pointwise is calculated virtual sun altitude, the position angle of each point on the digital topography map and is calculated by following formula:
θ
ij=arcsin(sinφ*sinδ+cosφ*cosδ*cost(i,j)) (1)
A(i,j)=arcsin(sinθ
ij*sinφ-sinδ)/cosθ
ij*cosφ) (2)
δ=arcsin(sinθ*sinφ-cosθ*cosφ*cosA) (3)
t=arcsin(cosθ
ij*sinA/cosδ)+Δλ (4)
θ wherein, A, δ are respectively and defend substar sun altitude in the sheet annotation, position angle and sun Chi Jiao; λ, φ are respectively the geographical longitude and latitude of substar.Δ λ is the longitude increment of pixel point to substar.
The judgement of massif and clouds layer shadow is on the above-mentioned analog digital ground type figure: utilize digital topography map and the location parameter of defending the sun of sheet pixel with the satellite remote sensing digital image registration, be elevation angle, position angle, carry out the judgement of massif and clouds layer shadow, the criterion of its judgement is: the sun altitude that is equal to, or greater than this pixel point in the maximum landform, cloud layer, the elevation angle that shine upon this pixel point of direction promptly: DH (i, j) 〉=θ
Ij, then this pixel point is a shade; Otherwise then not shade.
Direct sunlight landform illumination coefficient in each point ground (domatic) on the above-mentioned virtual calculating digital topography map; Comprise for digital topography map top shadow and unshaded point, calculate its virtual direct sunlight landform illumination coefficient respectively.
(1) calculating of the virtual direct sunlight landform of the some illumination coefficient in the digital topography map top shadow:
Because the virtual direct sunlight illumination of the point in the digital topography map top shadow is 0, so virtual direct sunlight landform illumination coefficient: F
Ij=0;
(2) the virtual direct sunlight landform illumination coefficient of unshaded point is calculated on the digital topography map:
The virtual direct sunlight landform illumination coefficient of non-shadow spots is calculated on the digital topography map, and is as follows:
F
ij=1-tgα
ij·ctgθ
ij·cosω
ij,ω
ij=AL
ij-A
ij (5)
Wherein, pixel ground altitude of the sun θ
Ij, position angle AL
IjProvide pixel ground inclination α by the virtual sheet annotation of defending
Ij, aspect A
IjAnd pixel provides ω through, latitude by the virtual corresponding DTM of sheet (digital terrain model that digital topography map generates) that defends
IjRepresent the angle of solar azimuth and ground aspect.
Above-mentioned again the computing machine of virtual direct sunlight (normalization) distributed image on the rolling ground is generated with visual.Be that the virtual direct sunlight landform illumination coefficient of the each point that will obtain is implemented linear stretch by original ranks series arrangement and pointwise and calculated, realize the one-tenth figure of virtual direct sunlight landform irradiation normalization coefficient and visual.
GN
ij=INT(F
ij×MAX(DN)/F)+0.5);I=1、2、3…M;J=1、2、3…N (6)
F=INT(MAX((F
ij));I=1、2、3…M;J=1、2、3…N) (7)
Wherein DN is the brightness of image value, MAX (DN)=2
K=-1 is maximum brightness value, and the maximum that M, N are respectively image is the pixel number in length and breadth.
Advantage that this invention is compared with prior art had and good effect
Generate ground direct sunlight distribution plan according to the inventive method, expressed the virtual direct sunlight relative exposure of nature ground each point under the influence of topography: virtual direct sunlight distribution (normalization) image that computing machine generates quantitatively, intuitively, has visually reflected the relative distribution of virtual direct sunlight illumination at natural terrain; The inventive method is in the influence of topography correction of the direct sunlight of remote sensing images, the Digital Simulation of radar detection effect, the Computer Simulation of remote sensing digital image and virtual, the dizzy computer manufacture that dyes map, defend the automatic differentiation of topographic shadowing on the sheet and the reallocation research of ground natural light radiation, the aspects such as planning and design of ground radar site all play an important role.
Description of drawings
Fig. 1 is the computer virtual process flow diagram of remote-sensing distributed image on rolling ground of the present invention.
Fig. 2 is 1: 10 ten thousand digital topography map (1024 * 1024) in area, Longyan, Fujian;
Fig. 3 is the linear stretch of true direct sunlight topography profile conversion coefficient.Be expressed as picture moment direct sunlight radiation at the true normalization on the rolling ground (the direct sunlight radiant illumination is 1 on the level ground) distributed image.
Fig. 4 is the linear stretch of virtual direct sunlight topography profile conversion coefficient.Be expressed as picture virtual direct sunlight radiation of moment at the normalization on the rolling ground (the direct sunlight radiant illumination is 1 on the level ground) distributed image.
Embodiment
A preferred embodiment of the present invention is described with reference to the accompanying drawings as follows:
This example is the computer virtual of area, Longyan, Fujian 1: 10 ten thousand digital topography map (1024 * 1024) with remote-sensing distributed image on rolling ground.
Referring to Fig. 2, it has shown this enforcement and has wanted the research area, i.e. 1: 10 ten thousand digital topography map (1024 * 1024) in area, Fujian China Longyan.
Referring to Fig. 3, it has shown area, Fujian China Longyan direct sunlight topography profile image (1024 * 1024, the solar direction southeast).
See also Fig. 1, it is the computer virtual process flow diagram of remote-sensing distributed image on rolling ground of the present invention, and in GIS, operating procedure is as follows:
Step 1. beginning, start-up routine, operation routinely earlier, input research area digital topography map, the satellite remote sensing digital image in Study of Registration area and digital topography map, and be divided into grid.
Step 2. is chosen first point on the digital topography map (i.e. first grid), the sun altitude of virtual this point, position angle on digital topography map: calculate by following formula:
θ
ij=arcsin(sinφ*sinδ+cosφ*cosδ*cost(i,j)) (1)
A(i,j)=arcsin(sinθ
ij*sinφ-sinδ)/cosθ
ij*cosφ) (2)
δ=arcsin(sinθ*sinφ-cosθ*cosφ*cosA) (3)
t=arcsin(cosθ
ij*sinA/cosδ)+Δλ (4)
θ, A, δ: defend substar sun altitude in the sheet annotation, position angle and sun Chi Jiao; λ, φ: be respectively the geographical longitude and latitude of substar.Δ λ is the longitude increment of pixel point to substar.
The judgement of massif and clouds layer shadow on the step 3. computer simulation digital topography map:
Utilize digital topography map and the location parameter of defending the sun of sheet pixel with the satellite remote sensing digital image registration, be elevation angle, position angle, carry out the judgement of massif and clouds layer shadow, the criterion of its judgement is: the sun altitude that is equal to, or greater than this pixel point in the maximum landform, cloud layer, the elevation angle that shine upon this pixel point of direction promptly: DH (i, j) 〉=θ
Ij, then this pixel point is a shade; Otherwise then not shade.
Ground (domatic) the direct sunlight landform illumination coefficient of each point on the virtual calculating digital topography map of step 4.:
For digital topography map top shadow and unshaded point, calculate its virtual direct sunlight landform illumination coefficient respectively.
(1). the calculating of the virtual direct sunlight landform of the some illumination coefficient in the digital topography map top shadow:
Because the virtual direct sunlight illumination of the point in the digital topography map top shadow is 0, so virtual direct sunlight landform illumination coefficient: F
Ij=0
(2). the virtual direct sunlight landform illumination coefficient of unshaded point is calculated on the digital topography map:
The virtual direct sunlight landform illumination coefficient of non-shadow spots is calculated on the digital topography map, and is as follows:
F
ij=1-tgα
ij·ctgθ
ij·cosω
ij ω
ij=AL
ij-A
ij (5)
Pixel ground altitude of the sun θ
Ij, position angle AL
IjProvide pixel ground inclination α by the virtual sheet annotation of defending
Ij, aspect A
IjAnd pixel provides ω through, latitude by the virtual corresponding DTM of sheet (digital terrain model that digital topography map generates) that defends
IjRepresent the angle of solar azimuth and ground aspect.
Whether step 5. is judged to have a few on the digital topography map to calculate and is finished, if do not finish, carries out ' next point ', returns execution in step 2; Otherwise execution in step 6.
The computing machine of direct sunlight on step 6. rolling ground (normalization) distributed image generates with visual.
The virtual direct sunlight landform illumination coefficient of the each point that obtains is implemented linear stretch by original ranks series arrangement and pointwise calculate, realize the one-tenth figure of virtual direct sunlight landform irradiation normalization coefficient and visual.As shown in Figure 4, wherein, GN
Ij=INT (F
Ij* MAX (DN)/F)+0.5); I=1,2,3 ... M; J=1,2,3 ... N (6)
F=INT(MAX((F
ij));I=1、2、3…M;J=1、2、3…N) (7)
Wherein DN is the brightness of image value; MAX (DN)=2
K=-1 maximum that is respectively image for maximum brightness value .M, N is the pixel number in length and breadth.
Step 7. output associated picture, program is finished.