CN106780321A - A kind of overall tight orientation of the satellite HR sensors images of CBERS 02 and correction joining method - Google Patents
A kind of overall tight orientation of the satellite HR sensors images of CBERS 02 and correction joining method Download PDFInfo
- Publication number
- CN106780321A CN106780321A CN201611047929.5A CN201611047929A CN106780321A CN 106780321 A CN106780321 A CN 106780321A CN 201611047929 A CN201611047929 A CN 201611047929A CN 106780321 A CN106780321 A CN 106780321A
- Authority
- CN
- China
- Prior art keywords
- image
- row
- correction
- image blocks
- raw video
- 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
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005304 joining Methods 0.000 title claims abstract description 9
- 238000012937 correction Methods 0.000 title claims description 50
- 238000007670 refining Methods 0.000 claims abstract description 7
- 238000012952 Resampling Methods 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 4
- 238000003491 array Methods 0.000 claims description 3
- 238000007493 shaping process Methods 0.000 claims description 3
- 244000062793 Sorghum vulgare Species 0.000 claims 1
- 235000019713 millet Nutrition 0.000 claims 1
- 238000010276 construction Methods 0.000 abstract description 4
- 238000012544 monitoring process Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 3
- 238000012827 research and development Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000012892 rational function Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- 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
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
The present invention provides a kind of overall tight orientation of each image blocks of satellite HR sensors of CBERS 02 and corrects joining method, is adapted to the satellite image of the star HR sensors of CEERS 02 or similar sensor construction.The method, by matching the tie point between the enough control points of acquisition and image blocks, overall tight orientation is carried out to each raw video block, and raw video is corrected and spliced using the orientation parameter after refining under the support of existing DOM and DEM.0 grade of image blocks without geometric manipulations are processed by tight model The present invention gives a kind of, the seamless spliced method just penetrated and correct image is generated, the quality and efficiency of CBERS 02C geometric manipulations is lifted.
Description
Technical field
The present invention relates to satellite remote sensing image geometry process field, more particularly to CBERS-02 satellites HR sensor images are whole
Body is oriented and corrects joining method.
Background technology
CBERS-02 includes two satellites, i.e. CBERS-02B and CBERS-02C stars.26 points during September in 2007 19 days 11,
By number 02B satellite of resource of mini-bus cooperation research and development, launched in Taiyuan, Shanxi satellite hub.Number 02B star of resource is real
When mode transmission earth observation remote sensing satellite, be China and Brazil cooperation the 3rd landsat.
In November, 2004, after the official signature of government of mini-bus two countries is on number agreement for 02B stars of joint research and development resource, press
According to leaders requirement, our Project R＆D troop is under commander in chief, the leading of chief designer, and Palestine Side's personnel's solidarity and cooperation,
Effort is striven, smoothly completed satellite home with Brazil related work, transmitting the same day given by Long March 4B
Enter space.Number 02B star orbitals road of resource selects subcircular sun synchronization, recurrence, Frozen Orbit, and the average height of orbit is 778km,
Repetition period is 26 days, and the lifetime of satellite is more than 2 years, 1455 kilograms of satellite weight, 1100 watts of end of lifetime power, satellite size
1800 × 2000 × 2000 millimeters.Satellite Payloads include 5 spectral coverage CCD camera, full-color high resolution camera, one
The wide visual field imager of the spectral coverage of platform 2, Data transfer system, solid-state memory, data collecting system and Space environment monitor system.
26 points during 22 days 11 December in 2011, China is delivered in Taiyuan Satellite Launch Center (TSLC) with " No. four second of the Long March "
Rocket, lift-off is succeeded in sending up by " resource one " 02C satellites.The HR cameras of 02C Seeds of First Post-flight two of resource one, spatial resolution
It it is 2.36 meters, the breadth of two splicings reaches 54km；The panchromatic and multispectral camera resolution ratio carried is respectively 5 meters and 10 meters,
Breadth is 60km.Two satellite image qualities of data meet 1:2.5 ten thousand -1:The requirement of 100000 survey of territorial resources monitoring accuracies, most
Small monitoring figure spot area reaches 0.2 mu, meets developed area, pays close attention to region resource present situation high-resolution investigation and monitoring
It is required that, index and conventional use of French SPOT -5, the moral such as attribute accuracy, area precision, minimum monitoring figure spot of fusion evaluation
State RapidEye data are approached, and be can be widely used for landuse dynamic monitoring by remote sensing, used in land use change survey, mineral resources and are opened
Hair is protected and utilization, Geological Hazards Investigation and monitoring, Eco-Geo-Environment are investigated and overseas the territory such as Mineral Resource Survey provides
Source main body business, at the same serve agricultural, forestry, water conservancy, environmental protection, ocean, mapping, traffic, live build, count, the correlation such as earthquake
Sector application.
Image geometry treatment model is that remote sensing image carries out photogrammetric basis, and image geometry model mainly has tight mould
Type (physical model), rational polynominal model and direct linear transformation's model (DLT) etc..The tight geometrical model of image mainly passes through
Aerial triangulation obtains high accuracy orientation parameter, i.e., intersect meter by the corner of the same name of the geometric manipulations to image or imaging point
Calculate to determine the two-dimensionally or three-dimensionally reason coordinate of target subject.The main method of the tight geometry of optical image is still common with basis
Based on line equation, the relevant parameters such as calibration model, model trajectory, the attitude mode of sensor are contacted by collinearity equation
Carry out and realize the solution of these model parameters.With reference to different model trajectorys, attitude mode, self calibration model, so as to occur in that
Several different tight geometrical models.
For space remote sensing image, Ebner is summarized as two kinds of main methods to describe the track of remote sensing aircraft, and one is logical
Discrete orientation point methods are crossed, another is by track restrained method.Orientation point methods provide almanac data by secondary file
To express, i.e., by orientation point information such as regular or irregularly time interval spatial point position, attitudes, to any photography moment
Orientation parameter, by interpolation obtain.Track restrained method assumes that remote sensing satellite is operated on a smooth mathematic curve,
All image exposure moment take the photograph station location all on this curve.Rational function model (RFM) is to apply very wide at present
General non-physical sensor model, is a kind of more broadly expression-form of various sensor geometric models, it be applied to it is various not
Same sensor, correction and the positioning of image are carried out by the rational polynominal coefficient for being given with rational polynominal function.By
Sensor parameters are concealed in it, many satellite image suppliers consider to use rational polynominal coefficient (RPC) several as image
The Transfer Standards of what information.In the case where there is ground control point, RFM models generally obtain a set of refining by block adjustment
Model parameter (such as affine transformation parameter) is realized improving precision to the amendment of geometry location error.The geometric manipulations of remote sensing image
Must with sensor in itself the characteristics of based on, pointedly it is modeled and adaptive processes is to obtain high-precision product
Condition and guarantee.But in the market does not carry out the software of specific aim geometric calibration and geometric manipulations also specifically designed for 02C satellites
Product.
The geometric manipulations of CBERS-02 satellite sensor images are not only the direct requirement of audio and video products production, are also this
The basis of image other application.It is obvious that number CBERS-02 satellite sensor of resource is different from other remote sensing satellite sensors
The characteristics of so that the geometric manipulations of CBERS-02 satellite images are also more different from other satellite remote-sensing images.CBERS-02 satellites
Can image be applied extensively and in depth, there is pass closely with geometric manipulations precision and the efficiency of CBERS-02 satellite images
System.
The present invention can lift CBERS-02 satellite images for the geometric manipulations of CBERS-02 satellite remote sensing satellite images
The quality of geometric product treatment, promotes the application of CBERS-02 satellite images.Currently for CBERS02 star HR sensor images
Geometric manipulations method is more, but mostly spliced on raw video after processed again, cause image joint and splicing after
Treatment mass effect it is not ideal.The present invention carries out entirety using orienting piece model to the HR image blocks of multiple observations simultaneously
Orientation and correction, realize the seamless spliced of image, realize 0 grade to the once sampling for just penetrating correction.
Bibliography：
Shao Jun, Guo Jianning .CBERS-02B satellite HR camera remote sensing image RPC bearing calibrations, space flight return and remote sensing
.2010.31(2):29-37.
The block adjustment remote sensing journals of Yuan Xiuxiao, Wang Tao sun .2012.CBERS-02B satellite remote-sensing images, 16 (2):
310-324.
Song Wei .CBERS-02 star charts are as geometric correction method experimental study land resources remote sensing .2009.79 (1):51-
54.
Cheng Chunquan, Zhang Jixian, Huang Guoman, horse crystalline substance Wuhan University Journal information sciences version .2016.41 (5):605-611.
Kingly way army, Gong Jianhua, it is CBERS-02 CCD non-Ortho images mosaic research ground of the based on geographical feature that horse is friendly
Reason and Geographical Information Sciences .2008.24 (2):1-4.
Xu Zhengnan, Feng Zhong certain herbaceous plants with big flowers are based on the CBERS-02B CCD image geometry error analysis remote sensing of computer simulation technique
Information .2010.6.26-30.
The CBERS-02B satellite HR camera image joining methods that Li Shiwei, Liu Tuanjie, Wang Hong fine jade are based on images match is distant
Sense technology and application .2009.24 (3):374-378.
Hao Yanling, Zhang Jie, the firm oceanographic researches .2010.28 (1) of horse:39-45.
Li Wenhua, Li Shihua, 5 red tinkling of pieces of jade rail-frees parameter high-resolution satellite image high accuracy geometric correction mappings sections
Learn .2009.34. supplementary issues:110-112.
The content of the invention
The invention provides a kind of CBERS-02 satellites HR sensor image global orientations and correction joining method, solve
The problem of image rectification, problem producing cause is：CBERS-02 star HR sensor high score images are detected by a plurality of short CCD linear arrays
Device is spliced and there is larger interval between ccd detector long, but each CCD, causes adjacent C CD image blocks laps same
Atural object exists on the acquisition time and there is obvious parallax between larger difference, image, and direct splicing can have obvious splicing seams.
The present invention does not carry out image processing after splicing to raw video block, but passes through during geometric manipulations
Each image blocks global orientation is directly realized by the correction of image with seamless spliced, the quality and efficiency of lifting geometric manipulations.It is adapted to
The satellite image of CEERS-02 star HR sensors or similar sensor construction.The method is led under the support of existing DOM and DEM
Overmatching obtains the tie point between enough control point and image blocks, and overall tight orientation, and profit are carried out to each raw video block
Image is corrected and spliced with the orientation parameter after refining.The method is realized by following step：
Step 1：Each image blocks are obtained and with reference to the same place between image by Image Matching, is sat according to reference to image picture point
Mark obtains ground point plane coordinates, and according to plane coordinates from DEM interpolation elevation, and then obtain the enough (logical of each image blocks
Often greater than 100) control point；
Step 2：By Image Matching, the tie point of the same name of lap between image blocks is obtained；
Step 3：Based on the collinearity equation of geocentric rectangular coordinate system space remote sensing image, the tight mould of orienting piece is set up
Type, each image blocks set common moment point as orientation point, and time interval is set to 1.0 seconds or so between each orientation point, respectively
Image blocks share a set of orbit refinement parameter, and each image blocks are each refined parameter using a set of independent attitude, by orienting piece
Tight model obtains the Precise Orbit and attitude parameter of each image blocks to each image global orientation；
Step 4：According to the track profile and dem data of each raw video block, the ground of whole scape image after correction splicing is obtained
Manage scope and width, the height of image after correction are calculated according to the Pixel size of image after correction；
Step 5：According to the resolution ratio of image after image capturing range after correction splicing and correction, the width of image after correction is calculated
W and height H, opens up two-dimentional shaping array space I [H] [K], and wherein K is according to the following two kinds situation value 4 or 7；When only HRA shadows
When picture or HRB images participate in splicing, value K=4, the setting image blocks 1 of raw video the 1st are arranged, 4086 row, the 2nd image blocks 4086
Row, the 3rd image blocks 4096 arrange totally 4 images row, to choosing row from left to right numbering 1 to 4, calculate choose all picture point things of row successively
Square coordinate and the ranks coordinate after correcting on image simultaneously round value, if arbitrarily certain is classified as k (1≤k≤4), then k represents HRA or HRB
From left to right k-th image row that image is set, if k arrange coordinate of certain picpointed coordinate on image after correction round value be (r,
C), then to array space I assignment I [r] [k-1]=c；
When HRA and HRB simultaneously participate in correction splicing, K values 7 set the 1st raw video block 1 of raw video HRA
Row, the 1st raw video block 4086 are arranged, and the 2nd raw video block 4086 is arranged, the 3rd raw video block 4000 arranges, belongs to the 1st former of HRB
Beginning image blocks 4086 are arranged, the 2nd raw video block 4086 is arranged, the 3rd raw video block 4096 arranges totally 7 row images, to choosing row a from left side
To right numbering 1 to 7, the ranks coordinate after choosing all picture point object coordinates of row and correcting on image is calculated successively and value is rounded,
Any row are set to k (1≤k≤7), then k represents k-th image row that HRA and HRB images are from left to right set, if k arranges certain picture point
Coordinate of the coordinate on image after correction rounds value for (r, c), then to element assignment I [r] [the k-1]=c of array space I.It is logical
Above-mentioned steps are crossed, each adjacent image blocks lap medium line picpointed coordinate is obtained according to direct remedy, recorded in two-dimemsional number
In group；
Step 6：All image blocks are carried out using indirect remedy unify to correct and splice, be achieved by the steps of：
Step (6.1)：The image walked after will definitely correcting according to the corresponding picture point of storage coordinate in I arrays is divided into K-1 blocks
(from the 1st to K-1 blocks)；
Step (6.2)：Successively since the image the first row after correction, terminate to last column, for each image row i
Perform step (6.3)；
Step (6.3)：Successively since the image first row after correction, terminate to last row, for correction image
Whether the image coordinate (i, j) of picture point p, is belonged to according to condition I [i] [n-1] ＜ j≤I [i] [n] (1≤n≤K-1) decision-point
N blocks are interval, if so, using correspondence the n-th block raw video and its orientation parameter of refining, p points correspondence object space coordinate projection is arrived
The image blocks carry out resampling；
Step (6.4)：4 pixels each to image block borderline pixels after correction and its both sides amount to 9 pixels again
Sampled, if being numbered with 1-9 successively to the 4th, the right pixel from the 1st, left side pixel, 9 pixels belong to weight
Folded area, in may map to two adjacent raw video blocks, if m (1≤m≤9) individual pixel is calculated by indirect method arrives original
Resampling calculated value is set to g in adjacent image blockslAnd gr, then the pixel value g of m-th pixel of image after correctingmBy formula gm
=0.1 [(10-m) gl+mgr] calculate.
The present invention realizes the correction of image with seamless spliced, the quality of lifting geometric manipulations by each image blocks global orientation
With efficiency.
Brief description of the drawings
Fig. 1 is that the CEERS-02 star HR images global orientation that the present invention is given corrects splicing step；
Fig. 2 is image and two-dimensional array I storage content schematic diagrames after HR sensor constructions, raw video block, correction.
Specific embodiment
The present invention gives a kind of CEERS-02 stars HR sensors or the satellite image geometric correction of similar sensor construction
With joining method, can in the case where that need not splice to raw video, directly each image blocks unify orientation,
It is unified to correct and splicing, the efficiency and precision of lifting image processing, it is to avoid multiple resampling brings information loss.As shown in figure 1,
The present invention is realized by following step：
Step 1. obtains each image blocks and with reference to the same place between image by Image Matching, is sat according to reference to image picture point
Mark obtains ground point plane coordinates, and according to plane coordinates from DEM interpolation elevation, and then obtain the enough control of each image blocks
Point, Image Matching lifts Image Matching speed and Image Matching using the matching of pyramid step by step and the matching process of object space constraint
Accuracy rate；
Step 2. passes through the tie point of the same name of lap between Image Matching, acquisition image blocks, the company of superimposed image part
Contact matching is realized using Gray-scale Matching method；
Step 3. sets up the tight mould of orienting piece based on the collinearity equation of geocentric rectangular coordinate system space remote sensing image
Type, each image blocks set common moment point as orientation point, and time interval is arranged on 1.0 seconds or so between each orientation point, respectively
Image blocks share a set of orbit refinement parameter, and each image blocks are each refined parameter using a set of independent attitude, by orienting piece
Tight model obtains the Precise Orbit and attitude parameter of each image blocks to each image global orientation.
Wherein orienting piece model is based on the collinearity equation of geocentric rectangular coordinate system satellite image, if Rc EFor sensor is sat
Mark is tied to the transition matrix of geocentric rectangular coordinate system, and its matrix element is [mij] (i, j=0,1,2), then there is collinearity equation：
Wherein (x, y) is picpointed coordinate, (x0,y0) it is principal point coordinate, (X, Y, Z) is topocentric coordinates, (XS,YS,ZS)
It is the coordinate of sensor projection centre in geocentric rectangular coordinate system, f is sensor focal distance.
Attitude and gps system error are still expressed using lower order polynomial expressions, are had：
Wherein (ai,bi,ci,di,ei,fi,gi) (i=0,1,2) be position and attitude system error polynomial model coefficient,
T is the interpolation moment.
Because Euler's compensation offset angle has initial value and precision information, can still facilitate together with other observations, be joined
Close adjustment.Install inclined with topocentric coordinates increment, GPS location, compensation offset angle and its systematic error multinomial coefficient, sensor
Incrementation parameter, sensor self calibration parameter are put as unknown number undetermined, by (2) respectively as observation error equation and strict side
Cheng Shi, builds two kinds of tight orientation error equation groups of form：
Wherein, Vi,Vp,Vs,VtPicpointed coordinate observation, the shadow according to formula (2) structure for respectively being built according to formula (1)
As photography moment POS measured value, self calibration parameter dummy observation and POS system error polynomial coefficient virtual observation change
Positive number vector；G, p, s, e, r represent the unknown number vector of tie point plane coordinates increment [△ X, △ Y] respectivelyT, photography moment orientation
Unknown number vectorSelf calibration unknown parameters number vector, sensor are installed inclined
Put away from vectorial with biasing angle increment unknown numberAnd the track appearance of formula (2) description
Unknown number vector [a of state systematic error multinomial model coefficienti,bi,ci,ei,fi,gi]T(i=0,1,2)；Bg、Bp、Bs、BeFor not
Know the coefficient matrix of number vector g, p, s, e；Measured value is the design matrix of corresponding unknown parameter coefficient；Li、Lp、Ls、LrRepresent corresponding
Observation error equation constant item vector；Pi、Pp、Ps、PrRepresent the weight matrix of corresponding observation.Above formula passes through a cum rights most young waiter in a wineshop or an inn
Multiply solution each unknown number vector g, p, s, e, r.
Any photography moment position and attitude initial value are formed by the moment neighbouring track and attitude measurement value interpolation, arbitrarily
Photography moment track profile correction value is that the track and attitude increment value interpolation at two neighbouring orienting piece moment of the moment are obtained.
The orientation parameter and dem data of each raw video block that step 4. is obtained according to tight model, obtain and correct splicing
Afterwards the geographic range of whole scape image and according to after correction image Pixel size calculate correct after image width, height.Wherein
Spliced image geographic range is corrected to be realized by step (4.1) and the step of step (4.2) two：
Step (4.1) calculates the corresponding object coordinates scope in each image blocks corner using the orientation parameter after refining；
Step (4.2) chooses the value X of transverse axis minimum in four angular coordinate scope0, the maximum value Y of the longitudinal axis1After splicing is corrected
Image top left co-ordinate (X0,Y1), choose the value X of maximum in transverse axis1, the minimum value Y of the longitudinal axis0As the image lower right corner after splicing
Coordinate (X1,Y0)；
Step 5. calculates the width of image after correction according to the resolution ratio of image after image capturing range after correction splicing and correction
W and height H, opens up two-dimentional shaping array space I [H] [K], and wherein K is according to the following two kinds situation value 4 or 7；When only HRA shadows
When picture or HRB images participate in splicing, value K=4, the setting image blocks 1 of raw video the 1st are arranged, 4086 row, the 2nd image blocks 4086
Row, the 3rd image blocks 4096 arrange totally 4 images row, to choosing row from left to right numbering 1 to 4, calculate choose all picture point things of row successively
Square coordinate and the ranks coordinate after correcting on image simultaneously round value, if being classified as k (1≤k≤4) selected by any, then k represent HRA or
K-th image that HRB images are from left to right set row, round value and are if k arranges coordinate of certain picpointed coordinate on image after correction
(r, c), then to element I [r] [k-1] assignment I [r] [the k-1]=c of array space I；
When HRA and HRB simultaneously participate in correction splicing, K values 7 set the 1st raw video block 1 of raw video HRA
Row, the 1st raw video block 4086 are arranged, and the 2nd raw video block 4086 is arranged, the 3rd raw video block 4000 arranges, belongs to the 1st former of HRB
Beginning image blocks 4086 are arranged, the 2nd raw video block 4086 is arranged, the 3rd raw video block 4096 arranges totally 7 row images, to choosing row a from left side
To right numbering 1 to 7, the ranks coordinate after choosing all picture point object coordinates of row and correcting on image is calculated successively and value is rounded,
If selected is arbitrarily classified as k (1≤k≤7), then k represents k-th image row that HRA and HRB images are from left to right set, if k
Arrange coordinate of certain picpointed coordinate on image after correction and round value for (r, c), then the element I [r] [k-1] to array space I is assigned
Value I [r] [k-1]=c.By above-mentioned steps, each adjacent image blocks lap medium line picture point is obtained according to direct remedy
Coordinate, records in two-dimensional array.
Image block is distinguished and I bis- after image, correction after CBERS-02 star sensors structure, raw video block, correction
Coordinate relation is recorded in dimension group as shown in Figure 2；
Step 6. uses indirect method, each image blocks is carried out with resampling correction, and the image near medium line is become
Weighting resampling.(6.1)-step (6.4) is realized as follows：
Step (6.1)：The image walked after will definitely correcting according to the corresponding picture point of storage coordinate in I arrays is divided into K-1 blocks
(from first to K-1 blocks)；
Step (6.2) since the image the first row after correction, terminates to last column successively, for each image row i
Perform step (6.3)；
Step (6.3) since the image first row after correction, terminates, for the picture point after correction to last row successively
P according to its correction after image coordinate (i, j), according to condition I [i] [n-1] ＜ j≤I [i] [n] (1≤n≤K-1) decision-point whether
Belong to n-th piece of interval, if then using correspondence the n-th block raw video and its orientation parameter of refining, by p points correspondence object coordinates
Projecting to the original image carries out resampling；
Step (6.4) 4 pixels each to image block borderline pixels after correction and its both sides amount to 9 pixels and enter again
Row sampling, if being numbered with 1-9 successively to the 4th, the right pixel from the 1st, left side pixel, 9 pixels belong to overlap
Area, in may map to two adjacent raw video blocks, if m (1≤m≤9) individual pixel is calculated to original phase by indirect method
Resampling calculated value is set to g in adjacent image blockslAnd gr, then the pixel value g of m-th pixel of image after correctingmBy formula gm=
0.1[(10-m)gl+mgr] calculate.By this step, realize weighting resampling to image blocks lap.
Claims (5)
1. a kind of each image blocks global orientation of CBERS-02 satellites HR sensors with correct joining method, it is characterised in that：With
Have with reference to orthography data as control and constrained by the tie point of lap between adjacent image blocks, by tight fixed
To piece model, each image blocks are carried out with global orientation, are corrected and is spliced, realize that the nothing of the whole scape image of 02CHR image blocks is stitched into,
Comprise the following steps：(1) each image blocks are obtained and with reference to the same place between orthography by Image Matching, is just penetrated according to reference
Image picpointed coordinate calculate ground point plane coordinates, and according to plane coordinates from digital elevation dem data interpolation height value, obtain
Obtain the control point of raw video block；(2) by Image Matching, the tie point of the same name of lap between acquisition image blocks；(3) with fixed
Based on to the tight model of piece, the Precise Orbit and attitude parameter of each image blocks are obtained by unified orientation；(4) according to accurate rail
Road attitude parameter and existing dem data, obtain the geographic range of whole scape image after correction splicing, and according to the picture of image after correction
First size calculates the width and height of image；(5) each adjacent image blocks lap medium line picture is obtained according to direct remedy
Point coordinates, and record in array；(6) use indirect method, each image blocks corrected, and to overlapping region medium line near
Image be weighted resampling.
2. CBERS-02 satellites HR sensor image global orientations according to claim 1 with correct joining method, it is special
Levy and be：In the step (3), in the orienting piece model, each image blocks set common time interval point as orientation point,
Each image blocks share a set of orbit refinement parameter, and each image blocks are each refined parameter using a set of independent attitude.
3. CBERS-02 satellites HR sensor image global orientations according to claim 1 with correct joining method, it is special
Levy and be：In the step (3), the three-dimensional coordinate at the control point for obtaining will be matched in the orienting piece model by step (1)
As given value, by step (2) match tie point between the adjacent image blocks for obtaining accordingly millet cake plane coordinates as not
Know number, elevation coordinate interpolation from DEM is obtained as given value.
4. CBERS-02 satellites HR sensor image global orientations according to claim 1 with correct joining method, it is special
Levy and be：In the step (5), according to the resolution ratio of image after image capturing range after correction splicing and correction, shadow after correction is calculated
The width W and height H of picture, open up two-dimentional shaping array space I [H] [K], and wherein K is according to the following two kinds situation value 4 or 7；When
When only HRA images or HRB images participate in correcting splicing, value K=4 sets the row of the 1st image blocks of raw video the 1st, the 4086th
Totally 4 images are arranged for row, the row of the 2nd image blocks the 4086th, the row of the 3rd image blocks the 4096th, to choosing row from left to right numbering 1 to 4, arbitrarily
Row are set to k, wherein 1≤k≤4, then k represents k-th image picture point row that HRA or HRB images are from left to right set, and calculates successively
Choose the ranks coordinate after arranging all picture point object coordinates and correcting on image and round value, if the raw video k for choosing arranges certain
Coordinate of the picpointed coordinate on image after correction rounds value for (r, c), then to element assignment I [r] [the k-1]=c in array I；
When HRA and HRB simultaneously participate in correction splicing, K values 7, the row of the 1st raw video block the 1st of setting raw video HRA,
1st raw video block 4086 is arranged, and the 2nd raw video block 4086 is arranged, the 3rd raw video block 4000 arranges, belongs to the 1st original of HRB
Image blocks 4086 are arranged, the 2nd raw video block 4086 is arranged, the 3rd raw video block 4096 arranges totally 7 row images, to choose row from a left side to
Right numbering 1 to 7, any row are set to k, wherein 1≤k≤7, then k represents k-th image that HRA and HRB images are from left to right set
Picture point is arranged, and the ranks coordinate after choosing all picture point object coordinates of row and correcting on image is calculated successively and value is rounded, if k arranges certain
Coordinate of the picpointed coordinate on image after correction rounds value for (r, c), then to array space I assignment I [r] [k-1]=c.
5. CBERS-02 satellites HR sensor image global orientations according to claim 1 with correct joining method, it is special
Levy and be：In the step (6), the indirect scheme of digital rectification is achieved by the steps of：
Step (6.1)：The image walked after will definitely correcting according to the corresponding picture point of storage coordinate in I arrays is divided into K-1 blocks, then entangle
Each image blocks after just are corresponded with raw video block；
Step (6.2) since the image the first row after correction, terminates to last column successively, is performed for each image row i
Step (6.3)；
Step (6.3)：Successively since the image first row after correction, terminate to last row, for the picture point p shadows after correction
As coordinate (i, j), whether n-th piece of interval is belonged to according to condition I [i] [n-1] ＜ j≤I [i] [n] (1≤n≤K-1) decision-point,
If using correspondence the n-th block raw video and its orientation parameter of refining, p points correspondence object space coordinate projection is former to corresponding n-th
Beginning image blocks image space coordinate carries out resampling；
Step (6.4) 4 pixels each to image block borderline pixels after correction and its both sides amount to 9 pixels and re-start adopts
Sample, if being numbered with 1-9 successively to the 4th, the right pixel from the 1st, left side pixel, 9 pixels belong to overlay region, can
To be mapped in two adjacent raw video blocks, if m (1≤m≤9) individual pixel is calculated to original adjacent left and right by indirect method
Resampling calculated value is set to g in image blockslAnd gr, then the pixel value g of m-th pixel of image after correctingmBy formula gm=
0.1[(10-m)gl+mgr] calculate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611047929.5A CN106780321B (en) | 2016-11-21 | 2016-11-21 | CBERS-02 satellite HR sensor image overall tight orientation and correction splicing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611047929.5A CN106780321B (en) | 2016-11-21 | 2016-11-21 | CBERS-02 satellite HR sensor image overall tight orientation and correction splicing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106780321A true CN106780321A (en) | 2017-05-31 |
| CN106780321B CN106780321B (en) | 2020-08-21 |
Family
ID=58974577
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611047929.5A Expired - Fee Related CN106780321B (en) | 2016-11-21 | 2016-11-21 | CBERS-02 satellite HR sensor image overall tight orientation and correction splicing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106780321B (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107465900A (en) * | 2017-08-22 | 2017-12-12 | 长光卫星技术有限公司 | A kind of preprocess method suitable for large area array sub-meter grade remote sensing video image unit |
| CN108387899A (en) * | 2018-04-17 | 2018-08-10 | 南京师范大学 | Ground control point automatically selecting method in synthetic aperture radar interferometry |
| CN108492334A (en) * | 2018-03-27 | 2018-09-04 | 中国海监南海航空支队 | A method of realizing commercial camera photo geographical calibration based on positioning and directing data |
| CN109063711A (en) * | 2018-07-06 | 2018-12-21 | 航天星图科技（北京）有限公司 | A kind of satellite image based on LLTS frame just penetrates correct algorithm |
| CN109579796A (en) * | 2018-12-24 | 2019-04-05 | 中国科学院遥感与数字地球研究所 | A kind of block adjustment method of image after projection |
| CN112184561A (en) * | 2020-09-25 | 2021-01-05 | 合肥工业大学 | Rocket stage pose deviation rectifying method and device based on multivariate search |
| CN113012083A (en) * | 2021-03-03 | 2021-06-22 | 上海景运信息技术有限公司 | Data processing method for image zero IO in remote sensing ortho image processing process |
| CN113624207A (en) * | 2021-08-06 | 2021-11-09 | 苏州九宇遥感科技有限公司 | Satellite image three-dimensional area network adjustment method based on orbit constraint |
| CN113899386A (en) * | 2021-09-27 | 2022-01-07 | 武汉大学 | Multi-source optical satellite remote sensing image collaborative regional net adjustment method and system based on three-dimensional reference net |
| CN115933652A (en) * | 2022-11-29 | 2023-04-07 | 北京航天飞行控制中心 | Lunar vehicle direct-drive teleoperation driving method based on sequence image splicing and fusion |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101303231A (en) * | 2008-07-08 | 2008-11-12 | 武汉大学 | Automatic selection of orthogonal projecting inlay line and orthogonal projection image seamless inlay method |
| CN101604018A (en) * | 2009-07-24 | 2009-12-16 | 中国测绘科学研究院 | High-definition remote sensing image data disposal route and system thereof |
| CN102202159A (en) * | 2011-03-29 | 2011-09-28 | 段连飞 | Digital splicing method for unmanned aerial photographic photos |
| CN103927731A (en) * | 2014-05-05 | 2014-07-16 | 武汉大学 | Low-altitude remote sensing image rapid and automatic splicing method without POS assisting |
| CN104318540A (en) * | 2014-10-20 | 2015-01-28 | 中南大学 | Aerial image online splicing method capable of using synergy of CPU (central processing unit) and GPU (graphics processing unit) |
| KR20150102525A (en) * | 2014-02-28 | 2015-09-07 | 경북대학교 산학협력단 | Image processing apparatus and method, and system for assisting vehicle driving using the same |
-
2016
- 2016-11-21 CN CN201611047929.5A patent/CN106780321B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101303231A (en) * | 2008-07-08 | 2008-11-12 | 武汉大学 | Automatic selection of orthogonal projecting inlay line and orthogonal projection image seamless inlay method |
| CN101604018A (en) * | 2009-07-24 | 2009-12-16 | 中国测绘科学研究院 | High-definition remote sensing image data disposal route and system thereof |
| CN102202159A (en) * | 2011-03-29 | 2011-09-28 | 段连飞 | Digital splicing method for unmanned aerial photographic photos |
| KR20150102525A (en) * | 2014-02-28 | 2015-09-07 | 경북대학교 산학협력단 | Image processing apparatus and method, and system for assisting vehicle driving using the same |
| CN103927731A (en) * | 2014-05-05 | 2014-07-16 | 武汉大学 | Low-altitude remote sensing image rapid and automatic splicing method without POS assisting |
| CN104318540A (en) * | 2014-10-20 | 2015-01-28 | 中南大学 | Aerial image online splicing method capable of using synergy of CPU (central processing unit) and GPU (graphics processing unit) |
Non-Patent Citations (3)
| Title |
|---|
| J. BERMUDEZ-CAMEO ET AL: ""Hypercatadioptric line images for 3D orientation and image rectification"", 《ROBOTICS AND AUTONOMOUS SYSTEMS》 * |
| 常军: ""基于严格物理成像模型的卫星影像纠正"", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 林铁 等: ""CBERS-02C卫星HR正射影像的快速拼接方法"", 《测绘地理信息》 * |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107465900A (en) * | 2017-08-22 | 2017-12-12 | 长光卫星技术有限公司 | A kind of preprocess method suitable for large area array sub-meter grade remote sensing video image unit |
| CN108492334A (en) * | 2018-03-27 | 2018-09-04 | 中国海监南海航空支队 | A method of realizing commercial camera photo geographical calibration based on positioning and directing data |
| CN108387899B (en) * | 2018-04-17 | 2020-07-31 | 南京师范大学 | Ground control point automatic selection method in synthetic aperture radar interferometry |
| CN108387899A (en) * | 2018-04-17 | 2018-08-10 | 南京师范大学 | Ground control point automatically selecting method in synthetic aperture radar interferometry |
| CN109063711B (en) * | 2018-07-06 | 2021-10-29 | 中科星图股份有限公司 | Satellite image orthorectification algorithm based on LLTS framework |
| CN109063711A (en) * | 2018-07-06 | 2018-12-21 | 航天星图科技（北京）有限公司 | A kind of satellite image based on LLTS frame just penetrates correct algorithm |
| CN109579796A (en) * | 2018-12-24 | 2019-04-05 | 中国科学院遥感与数字地球研究所 | A kind of block adjustment method of image after projection |
| CN109579796B (en) * | 2018-12-24 | 2020-11-10 | 中国科学院遥感与数字地球研究所 | Area network adjustment method for projected image |
| CN112184561A (en) * | 2020-09-25 | 2021-01-05 | 合肥工业大学 | Rocket stage pose deviation rectifying method and device based on multivariate search |
| CN112184561B (en) * | 2020-09-25 | 2024-02-20 | 合肥工业大学 | Rocket stage pose correction method and device based on multiple search |
| CN113012083A (en) * | 2021-03-03 | 2021-06-22 | 上海景运信息技术有限公司 | Data processing method for image zero IO in remote sensing ortho image processing process |
| CN113624207A (en) * | 2021-08-06 | 2021-11-09 | 苏州九宇遥感科技有限公司 | Satellite image three-dimensional area network adjustment method based on orbit constraint |
| CN113899386A (en) * | 2021-09-27 | 2022-01-07 | 武汉大学 | Multi-source optical satellite remote sensing image collaborative regional net adjustment method and system based on three-dimensional reference net |
| CN113899386B (en) * | 2021-09-27 | 2023-11-21 | 武汉大学 | Multi-source optical satellite remote sensing image collaborative regional network adjustment method and system based on three-dimensional reference network |
| CN115933652A (en) * | 2022-11-29 | 2023-04-07 | 北京航天飞行控制中心 | Lunar vehicle direct-drive teleoperation driving method based on sequence image splicing and fusion |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106780321B (en) | 2020-08-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106780321A (en) | A kind of overall tight orientation of the satellite HR sensors images of CBERS 02 and correction joining method | |
| CN105627991B (en) | A kind of unmanned plane image real time panoramic joining method and system | |
| CN110310248B (en) | A kind of real-time joining method of unmanned aerial vehicle remote sensing images and system | |
| CN103914808B (en) | Method for splicing ZY3 satellite three-line-scanner image and multispectral image | |
| Zongjian | UAV for mapping—low altitude photogrammetric survey | |
| CN103674063B (en) | A kind of optical remote sensing camera geometric calibration method in-orbit | |
| CN102506824B (en) | Method for generating digital orthophoto map (DOM) by urban low altitude unmanned aerial vehicle | |
| CN106767714B (en) | Improve the equivalent mismatch model multistage Calibration Method of satellite image positioning accuracy | |
| CN104748728B (en) | Intelligent machine attitude matrix calculation method and its applied to photogrammetric method | |
| CN103115614B (en) | Associated parallel matching method for multi-source multi-track long-strip satellite remote sensing images | |
| CN106895851B (en) | A kind of sensor calibration method that the more CCD polyphasers of Optical remote satellite are uniformly processed | |
| CN103900539B (en) | A kind of aerial cube panoramic imagery object localization method | |
| CN110111274B (en) | Method for calibrating exterior orientation elements of satellite-borne push-broom optical sensor | |
| CN103047985A (en) | Rapid positioning method for space target | |
| CN102968631A (en) | Automatic geometric correction and orthorectification method for multispectral remote sensing satellite images of mountainous area | |
| CN108305237A (en) | Consider more stereopsis fusion drafting method of different illumination image-forming conditions | |
| CN101551275A (en) | Technical method of vehicular multispectral scanner for monitoring industrial warm discharge water | |
| KR102075028B1 (en) | Unmanned High-speed Flying Precision Position Image Acquisition Device and Accurate Position Acquisition Method Using the same | |
| CN100443919C (en) | Automatic assimilation method for multi-source thermal infrared wave band data of polar-orbit meteorological satellite | |
| CN105444778B (en) | A kind of star sensor based on imaging geometry inverting is in-orbit to determine appearance error acquisition methods | |
| CN114972545B (en) | On-orbit data rapid preprocessing method for hyperspectral satellite | |
| CN106873004A (en) | The in-orbit geometry calibration method of rail level array camera high based on sun altitude self adaptation | |
| Tang et al. | Geometric accuracy analysis model of the Ziyuan-3 satellite without GCPs | |
| CN107316280A (en) | Li Island satellite image RPC models high accuracy geometry location method | |
| CN108253942B (en) | Method for improving oblique photography measurement space-three quality |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200821 Termination date: 20211121 |