CN103971334B

CN103971334B - High-spectrum remote sensing bearing calibration and device

Info

Publication number: CN103971334B
Application number: CN201410171313.3A
Authority: CN
Inventors: 吴远峰; 张�浩; 陈正超; 孙旭; 高连如; 张兵
Original assignee: Institute of Remote Sensing and Digital Earth of CAS
Current assignee: Institute of Remote Sensing and Digital Earth of CAS
Filing date: 2014-04-25
Publication date: 2016-11-30
Anticipated expiration: 2034-04-25

Abstract

The embodiment of the invention discloses a kind of high-spectrum remote sensing bearing calibration and device, from high-spectrum remote sensing to be corrected, determine that a wave band is for reference to wave band, the relative displacement of each row is determined by the minimum error of adjacent lines, and the absolute offset values (the most each row is relative to the side-play amount of target line pixel) of each row is determined by the relative displacement of each row, thus move row corresponding in high-spectrum remote sensing to be corrected, the high-spectrum remote sensing after being corrected according to the absolute offset values of each row by wave band.It can thus be appreciated that, the high-spectrum remote sensing bearing calibration of the embodiment of the present application offer and device, from the feature of high-spectrum remote sensing own, related data between statistics row and row within high-spectrum remote sensing, it is corrected by minimum and position translation line by line the method for row correlated error, do not rely on external attitude data, it is achieved that for the correction of high-spectrum remote sensing geometric distortion.

Description

High-spectrum remote sensing bearing calibration and device

Technical field

The present invention relates to technical field of image processing, more particularly, it relates to a kind of high-spectrum remote sensing bearing calibration And device.

Background technology

Along with developing rapidly of science and technology, remote sensing technology is increasingly paid close attention to by people, classifies according to remote sensing platform, Substantially can be divided into space remote sensing, air remote sensing and ground remote sensing, wherein, air remote sensing refers to and (includes man-machine by aircraft And unmanned plane), dirigible, the aerial platform such as balloon carry the Remote Sensing Technical System of imaging sensor earth observation.

But, the easy wind-engaging of aerial platform, the impact of unstable air-flow and himself mechanical part vibrations etc. and occur Side-sway, the phenomenon rolled or shake, and in high spectrum resolution remote sensing technique, imaging sensor typically uses push-scanning image mode to obtain ground The multispectral data (imaging the most line by line) of thing, the instability of aerial platform will greatly reduce image quality so that the Gao Guang of acquisition There is " burr " or the geometric distortion phenomenon of distortion in spectrum remote-sensing image, as shown in Figure 1.Therefore, it is necessary to geometric distortion occurs High-spectrum remote sensing be corrected.

At present, the method that conventional a kind of high-spectrum remote sensing to there is geometric distortion is corrected is: by sky Middle platform carries attitude recording equipment, utilizes the data of these attitude recording equipment records to enter the high-spectrum remote sensing obtained Row corrects by pixel.But the general small volume of unmanned plane, for weight and the consideration of cost factor, it is high-precision that unmanned plane is difficult to lift-launch The attitude recording equipment of degree (its volume and weight of the most accurate equipment is the biggest), therefore, when aerial platform is unmanned plane, passes through The mode carrying attitude recording equipment is difficult to be corrected the high-spectrum remote sensing of generation geometric distortion.

Therefore, how geometric distortion to the high-spectrum remote sensing obtained by unmanned plane mode is corrected becoming urgently Problem to be solved.

Summary of the invention

It is an object of the invention to provide a kind of high-spectrum remote sensing bearing calibration and device, with to by unmanned plane mode The geometric distortion of the high-spectrum remote sensing obtained is corrected.

For achieving the above object, the technical scheme is that

A kind of high-spectrum remote sensing bearing calibration, including:

Determine all wave bands that reference wave band, described reference wave band are high-spectrum remote sensing to be corrected meet and preset The wave band of condition；

According to default movement rule by the described i+1 row pixel with reference to wave band relative to the i-th of described reference wave band Row pixel moves pixel-by-pixel, often moves a location of pixels, calculates the described i+1 row pixel with reference to wave band according to the first formula With the error amount of described the i-th row pixel with reference to wave band, described first formula is:

E=[∑ (R_p-C_q)²]/m²

Wherein, E is the described i+1 row pixel with reference to wave band and the error amount of described the i-th row pixel with reference to wave band, R_p It is the value of pth pixel, C in the i-th row pixel_qFor the value of q-th pixel in i+1 row pixel, in described i+1 row pixel Q pixel is alignd with pth pixel column in described i-th row pixel；M is in i+1 row pixel to be column alignment with the i-th row pixel The number of pixel；

During by error amount minimum, described i+1 row pixel is defined as described relative to the side-play amount of described i-th row pixel The relative displacement of i+1 row pixel, wherein, described i+1 row pixel is positive offset relative to the side-play amount of described i-th row pixel Amount or negative side-play amount, described i+1 row pixel is preset relative to the positive and negative foundation of the value of the side-play amount of described i-th row pixel Positive offset direction determine；

Calculate the first absolute offset values of described each row pixel with reference to wave band, including: calculate and be positioned at the i-th row pixel with pre- The relative displacement sum of all row between the target line pixel first determined, obtains first and value；The of described i-th row pixel One absolute offset values is the relative displacement sum of described first and value with described i-th row pixel；

According to the first absolute offset values of described each row pixel with reference to wave band, move EO-1 hyperion to be corrected by wave band distant Row corresponding in sense image, it is thus achieved that the high-spectrum remote sensing after the first correction；Wherein, the pixel count that the i-th row pixel moves is The absolute value of the first absolute offset values of described the i-th row pixel with reference to wave band, the side that each wave band the i-th row pixel described moves To identical with described offset direction with reference to corresponding to the first absolute offset values of the i-th row pixel of wave band；Wherein, i is for being more than Or the positive integer equal to 1.

Said method, it is preferred that described according to default movement rule by relative for the described i+1 row pixel with reference to wave band Move pixel-by-pixel in described the i-th row pixel with reference to wave band and include:

Determine starting column location and position, end column that described i+1 row pixel moves relative to described i-th row pixel Put；

By the picture at the pixel maximum for the column position of the described i+1 row pixel starting column location with described i-th row pixel Element column alignment；

Described i+1 row pixel is moved pixel-by-pixel to described end column locality, until described i+1 row pixel The pixel of column position minimum is alignd with the pixel column of the position, end column of described i-th row pixel.

Said method, it is preferred that described determine what described i+1 row pixel moved relative to described i-th row pixel Starting column location and position, end column include:

Determine the minimum pixel column position of the column position of described i-th row pixel be described i+1 row pixel relative to The starting column location that described i-th row pixel moves；

Determine the maximum pixel column position of the column position of described i-th row pixel be described i+1 row pixel relative to The position, end column that described i-th row pixel moves.

Determine that the first row position that user inputs is that described i+1 row pixel moves relative to described i-th row pixel Starting column location；

Determine that the secondary series position that user inputs is that described i+1 row pixel moves relative to described i-th row pixel Position, end column；

Wherein, described first row position is different from the pixel column position of the column position of described i-th row pixel minimum, institute The pixel column position stating secondary series position and the column position of described i-th row pixel maximum is different, described first row position with Described secondary series position is different.

Determine that described i+1 row pixel is entered relative to described i-th row pixel according to the degree of overlapping of predetermined adjacent lines The starting column location of row movement and position, end column, wherein,

The starting column location that described i+1 row pixel moves relative to described i-th row pixel is:

The position, end column that described i+1 row pixel moves relative to described i-th row pixel is:

Wherein, n is the number of pixels of described i-th row pixel；R is the degree of overlapping of described predetermined adjacent lines.

Said method, it is preferred that after the high-spectrum remote sensing after obtaining the first correction, also include:

Reference wave band according to first user operation high-spectrum remote sensing after described first correction draws target The distortion contour curve of atural object；

Reference wave band according to second user operation high-spectrum remote sensing after described first correction is drawn described The correct contour curve of Target scalar；

Determine the 3rd column position of described the i-th row pixel with reference to wave band and the intersection point of described distortion contour curve, and really 4th column position of the intersection point of fixed described the i-th row pixel with reference to wave band and described correct contour curve；

Obtain the second absolute offset values of described the i-th row pixel with reference to wave band, described the i-th row pixel with reference to wave band Second absolute offset values is the difference of described 3rd column position and described 4th column position；

According to the second absolute offset values of described i-th row pixel, by wave band by the high-spectrum remote-sensing after described first correction I-th row pixel of image moves to the direction of described correct contour curve, it is thus achieved that the high-spectrum remote sensing after the second correction；Institute State the i-th row pixel and move the absolute value that amount of movement is described second absolute offset values to the direction of described correct contour curve.

A kind of high-spectrum remote sensing correcting unit, including:

First determines module, is used for determining with reference to wave band, and described is high-spectrum remote sensing to be corrected with reference to wave band All wave bands meet pre-conditioned wave band；

First computing module, for according to default movement rule by the described i+1 row pixel with reference to wave band relative to Described the i-th row pixel with reference to wave band moves pixel-by-pixel, often moves a location of pixels, calculates described reference according to the first formula The i+1 row pixel of wave band and the error amount of described the i-th row pixel with reference to wave band, described first formula is:

E=[∑ (R_p-C_q)²]/m²

Second determines module, and when by error amount minimum, described i+1 row pixel is relative to described i-th row pixel Side-play amount is defined as the relative displacement of described i+1 row pixel, and wherein, described i+1 row pixel is relative to described i-th row picture The side-play amount of element is positive offset amount or negative side-play amount, and described i+1 row pixel is relative to the taking of side-play amount of described i-th row pixel The positive offset direction that the positive and negative foundation of value is preset determines；

Second computing module, for calculating the first absolute offset values of described each row pixel with reference to wave band, including: calculate The relative displacement sum of all row between the i-th row pixel and predetermined target line pixel, obtains first and value；Described First absolute offset values of the i-th row pixel is the relative displacement sum of described first and value with described i-th row pixel；

First correction module, for the first absolute offset values according to described each row pixel with reference to wave band, moves by wave band Row corresponding in dynamic high-spectrum remote sensing to be corrected, it is thus achieved that the high-spectrum remote sensing after the first correction；Wherein, the i-th row The pixel count that pixel moves is the absolute value of the first absolute offset values of described the i-th row pixel with reference to wave band, each ripple described The skew side corresponding to first absolute offset values of the direction moved of section the i-th row pixel and described the i-th row pixel with reference to wave band To identical；Wherein, i is the positive integer more than or equal to 1.

Said apparatus, it is preferred that described first computing module includes: mover module and calculating sub module；Wherein,

Described mover module for according to default movement rule by the described i+1 row pixel with reference to wave band relative to Described the i-th row pixel with reference to wave band moves pixel-by-pixel；Including:

First determines unit, for determining that described i+1 row pixel is initial relative to what described i-th row pixel moved Column position and position, end column；

Mobile unit, for by the pixel of column position maximum of described i+1 row pixel and initiateing of described i-th row pixel Pixel column alignment at column position, moves to described end column locality pixel-by-pixel by described i+1 row pixel, until described The pixel of the column position minimum of i+1 row pixel is alignd with the pixel column of the position, end column of described i-th row pixel；

Described calculating sub module is used for often moving a location of pixels in the described i+1 row pixel with reference to wave band, according to First formula calculates the error amount of the described i+1 row pixel with reference to wave band and described the i-th row pixel with reference to wave band, and described the One formula is:

E=[∑ (R_p-C_q)²]/m²

Wherein, E is the described i+1 row pixel with reference to wave band and the error amount of described the i-th row pixel with reference to wave band, R_p It is the value of pth pixel, C in the i-th row pixel_qFor the value of q-th pixel in i+1 row pixel, in described i+1 row pixel Q pixel is alignd with pth pixel column in described i-th row pixel；M is in i+1 row pixel to be column alignment with the i-th row pixel The number of pixel.

Said apparatus, it is preferred that described first determines that unit includes:

First determines subelement, is institute for determining the pixel column position of the column position minimum of described i-th row pixel State the starting column location that i+1 row pixel moves relative to described i-th row pixel；

Second determines subelement, is institute for determining the pixel column position of the column position maximum of described i-th row pixel State the position, end column that i+1 row pixel moves relative to described i-th row pixel.

3rd determines subelement, for determining that the first row position that user inputs is that described i+1 row pixel is relative to institute State the starting column location that the i-th row pixel moves；

4th determines subelement, for determining that the secondary series position that user inputs is that described i+1 row pixel is relative to institute State the position, end column that the i-th row pixel moves；

5th determines subelement, for determining described i+1 row pixel phase according to the degree of overlapping of predetermined adjacent lines The starting column location that described i-th row pixel is moved and position, end column, wherein,

Said apparatus, it is preferred that also include:

First drafting module, the ginseng of the high-spectrum remote sensing after correcting described first according to first user operation Examine the distortion contour curve drawing Target scalar in wave band；

Second drafting module, the ginseng of the high-spectrum remote sensing after correcting described first according to the second user operation Examine the correct contour curve drawing described Target scalar in wave band；

3rd determines module, for determining described the i-th row pixel with reference to wave band and the intersection point of described distortion contour curve The 3rd column position, and determine the 4th row position of described the i-th row pixel with reference to the wave band intersection point with described correct contour curve Put；

Acquisition module, for obtaining the second absolute offset values of described the i-th row pixel with reference to wave band, described with reference to wave band The difference that the second absolute offset values is described 3rd column position and described 4th column position of the i-th row pixel；

Second correction module, for the second absolute offset values according to described i-th row pixel, by wave band by described first school I-th row pixel of the high-spectrum remote sensing after just moves to the direction of described correct contour curve, it is thus achieved that after the second correction High-spectrum remote sensing；It is described second absolute that described i-th row pixel moves amount of movement to the direction of described correct contour curve The absolute value of side-play amount.

By above scheme, a kind of high-spectrum remote sensing bearing calibration of the application offer and device, from treating school Positive high-spectrum remote sensing determining, a wave band is for reference to wave band, determining the relative of each row by the minimum error of adjacent lines Side-play amount, and (the most each row is inclined relative to target line pixel to determine the absolute offset values of each row by the relative displacement of each row Shifting amount), thus move row corresponding in high-spectrum remote sensing to be corrected according to the absolute offset values of each row by wave band, obtain High-spectrum remote sensing after correction.It follows that the high-spectrum remote sensing bearing calibration of the embodiment of the present application offer and dress Put, from the feature of high-spectrum remote sensing own, the related data between statistics row and row within high-spectrum remote sensing, It is corrected by minimum and position translation line by line the method for row correlated error, does not relies on external attitude data, it is achieved that pin Correction to high-spectrum remote sensing geometric distortion.

Accompanying drawing explanation

In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing In having technology to describe, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only this Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to Other accompanying drawing is obtained according to these accompanying drawings.

Fig. 1 is the exemplary plot of the unmanned plane high-spectrum remote sensing that geometric distortion occurs；

A kind of flowchart of the high-spectrum remote sensing bearing calibration that Fig. 2 provides for the embodiment of the present application；

The exemplary plot of the i+1 row pixel that Fig. 3 a provides for the embodiment of the present application and a kind of relative position of the i-th row pixel；

The example of the i+1 row pixel that Fig. 3 b provides for the embodiment of the present application and the another kind of relative position of the i-th row pixel Figure；

The schematic diagram that the i+1 row pixel that Fig. 4 provides for the embodiment of the present application moves pixel-by-pixel relative to the i-th row pixel；

The another kind of flowchart of the high-spectrum remote sensing bearing calibration that Fig. 5 provides for the embodiment of the present application；

A kind of structural representation of the high-spectrum remote sensing correcting unit that Fig. 6 provides for the embodiment of the present application；

A kind of structural representation of the first computing module that Fig. 7 provides for the embodiment of the present application；

The first a kind of structural representation determining unit that Fig. 8 provides for the embodiment of the present application；

The first another kind of structural representation determining unit that Fig. 9 provides for the embodiment of the present application；

First another structural representation determining unit that Figure 10 provides for the embodiment of the present application；

The another kind of structural representation of the high-spectrum remote sensing correcting unit that Figure 11 provides for the embodiment of the present application；

Figure 12 carries out school by the embodiment of the present application to high-spectrum remote sensing shown in Fig. 1 for what the embodiment of the present application provided Design sketch after the first correction obtained the most afterwards；

The distortion profile that the user that Figure 13 provides for the embodiment of the present application sketches in high-spectrum remote sensing shown in Figure 12 Curve and correct contour curve；

Figure 14 for the embodiment of the present application provide the high-spectrum remote sensing shown in Figure 12 is corrected after obtain Design sketch after two corrections.

Term " first " in specification and claims and above-mentioned accompanying drawing, " second ", " the 3rd " " the 4th " etc. (if Exist) it is for distinguishing similar part, without being used for describing specific order or precedence.Should be appreciated that so use Data can exchange in the appropriate case, in order to embodiments herein described herein can be with except here illustrating Order in addition is implemented.

Detailed description of the invention

Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete Describe, it is clear that described embodiment is only a part of embodiment of the present invention rather than whole embodiments wholely.Based on Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under not making creative work premise Embodiment, broadly falls into the scope of protection of the invention.

The one referring to the high-spectrum remote sensing bearing calibration that Fig. 2, Fig. 2 provide for the embodiment of the present application realizes flow process Figure, may include that

Step S21: determine that, with reference to wave band, described reference wave band is in all wave bands of high-spectrum remote sensing to be corrected Meet pre-conditioned wave band；

High-spectrum remote sensing generally comprises tens, the data of up to a hundred wave bands.In the embodiment of the present application, distant from EO-1 hyperion Selecting a wave band as with reference to wave band in all wave bands of sense image, wherein, described can be that signal to noise ratio meets with reference to wave band Preset the wave band of signal to noise ratio condition；

Preferably, the image of each wave band can be shown, artificial visual determine with reference to wave band, i.e. by manually selecting Select a noise fewer, and compare clearly that wave band is as with reference to wave band, empirically, corresponding high-spectrum remote sensing, front The signal to noise ratio of the several wave band in face and the most several wave band is the most poor, so, it will usually in the middle of selecting, some wave band is as ginseng Examine wave band；Can also determine with reference to wave band by calculating, concrete, high-spectrum remote sensing to be corrected can be calculated by wave band The signal to noise ratio of each wave band, select the maximum wave band of signal to noise ratio as with reference to wave band, the wave band maximum when signal to noise ratio exists many Time individual, a wave band can be randomly choosed from the wave band of signal to noise ratio maximum as with reference to wave band.

Step S22: according to default movement rule by the described i+1 row pixel with reference to wave band relative to described reference wave Section the i-th row pixel move pixel-by-pixel, often move a location of pixels, according to first formula calculate described with reference to wave band i-th+ 1 row pixel and the error amount of described the i-th row pixel with reference to wave band, described first formula is:

E=[∑ (R_p-C_q)²]/m²

When calculating, the value of pixel can be original gray value, it is also possible to obtain after radiant correction Spoke brightness value, it is also possible to be the reflectance value obtained through atmospheric correction.

In the embodiment of the present application, by the described i+1 row pixel with reference to wave band relative to described the i-th row picture with reference to wave band Element moves pixel-by-pixel and refers to that the position of the i-th row pixel keeps constant, and mobile i+1 row pixel, i+1 row pixel is moving Time, mobile step-length once is a pixel.

It should be noted that work as described the i-th+2 row pixel with reference to wave band relative to the described i+1 row with reference to wave band When pixel moves pixel-by-pixel, the position of i+1 row pixel is described i+1 row pixel raw bits in described reference wave band Put rather than i+1 row pixel move after position.

In order to better illustrate the computational methods of i+1 row pixel and the error amount of the i-th row pixel, refer to Fig. 3 a, figure The exemplary plot of the i+1 row pixel that 3a provides for the embodiment of the present application and a kind of relative position of the i-th row pixel；

In Fig. 3 a, having 3 pixels in i+1 row pixel with the i-th row pixel is column alignment, wherein, in i+1 row pixel R1 column alignment in C6 and the i-th row pixel, R2 column alignment, C8 in i+1 row pixel in C7 and the i-th row pixel in i+1 row pixel With R3 column alignment in the i-th row pixel；Then, the i+1 row pixel shown in Fig. 3 a and the error value E of the i-th row pixel_example1Calculating Method is:

E_example1=[(R1-C6)²+(R2-C7)²+(R3-C8)²]/3²。

On the basis of Fig. 3 a, after i+1 row pixel moves right a location of pixels, i+1 row pixel and the i-th row The exemplary plot of the another kind of position relatively of pixel is as shown in Figure 3 b；

In Fig. 3 b, having 4 pixels in i+1 row pixel with the i-th row pixel is column alignment, wherein, in i+1 row pixel R1 column alignment in C5 and the i-th row pixel, R2 column alignment, C7 in i+1 row pixel in C6 and the i-th row pixel in i+1 row pixel With R3 column alignment in the i-th row pixel, R4 column alignment in C8 and the i-th row pixel in i+1 row pixel；The then i+1 shown in Fig. 3 b Row pixel and the error value E of the i-th row pixel_example2Computational methods be:

E_example2=[(R1-C5)²+(R2-C6)²+(R3-C7)²+(R4-C8)²]/4²。

Step S23: during by error amount minimum, described i+1 row pixel determines relative to the side-play amount of described i-th row pixel For the relative displacement of described i+1 row pixel, wherein, described i+1 row pixel is relative to the side-play amount of described i-th row pixel For positive offset amount or negative side-play amount, described i+1 row pixel is positive and negative relative to the value of the side-play amount of described i-th row pixel Determine according to the positive offset direction preset；

Owing to often moving an i+1 row pixel, calculate the error amount of an i+1 row pixel and the i-th row pixel, institute With, i+1 row pixel has multiple with the error amount of the i-th row pixel, in the embodiment of the present application, from the i+1 row pixel obtained with In multiple error amounts of the i-th row pixel, determine the error amount that value is minimum, and by corresponding for error amount that this value is minimum i-th+ 1 row pixel is defined as the relative displacement of described i+1 row pixel relative to the side-play amount of described i-th row pixel.

In the embodiment of the present application, i+1 row pixel includes i+1 row pixel relative to the side-play amount of described i-th row pixel L row deviate described i-th row pixel L row pixel count and offset direction.Therefore, in the embodiment of the present application, described i-th + 1 row pixel is positive offset amount or negative side-play amount relative to the side-play amount of described i-th row pixel, described i+1 row pixel relative to The positive offset direction that the positive and negative foundation of the value of the side-play amount of described i-th row pixel is preset determines；Concrete, when i+1 row picture When the offset direction of the L row that the L row of element deviate described i-th row pixel is positive offset direction, i+1 row pixel is relative to institute The side-play amount stating the i-th row pixel is positive offset amount, and otherwise, i+1 row pixel is negative relative to the side-play amount of described i-th row pixel Side-play amount.

Concrete, when the L of i+1 row pixel arranges the first side of the L row being positioned at described i-th row pixel, Ke Yiding Justice is positive offset direction；When the L of i+1 row pixel arranges the second side of the L row being positioned at described i-th row pixel, it is defined as Negative offset direction.Such as, as shown in Figure 3 a, can define when the 1st row (i.e. the row at C1 place) of i+1 row pixel are positioned at described It is positive offset direction time on the right side of the of 1st row (i.e. the row at R1 place) of the i-th row pixel, then, now, i+1 row pixel phase Relative displacement for described i-th row pixel is-5.

It should be noted that in the embodiment of the present application, the i-th row pixel is 0 relative to the relative displacement of the i-th row pixel.

Step S24: calculate the first absolute offset values of described each row pixel with reference to wave band, including: calculate and be positioned at the i-th row The relative displacement sum of all row between pixel and predetermined target line pixel, obtains first and value；Described i-th row First absolute offset values of pixel is the relative displacement sum of described first and value with described i-th row pixel；

Described predetermined target line pixel can be described with reference to one-row pixels any in wave band；The embodiment of the present application In, the first absolute offset values of the i-th row pixel is substantially the i-th row pixel phase relative to described predetermined target line pixel To side-play amount.

In the embodiment of the present application, the absolute offset values of predetermined target line pixel is 0.

For example, it is assumed that predetermined target line pixel is the 6th row pixel with reference to wave band, then, the 6th row pixel exhausted It is 0 to side-play amount；The relative displacement that absolute offset values is the 5th row pixel of the 5th row pixel, the absolute drift of the 4th row pixel Amount is then relative displacement and the relative displacement sum of the 5th row pixel of the 4th row pixel, the absolute offset values of the 3rd row pixel It is then the relative displacement of the 3rd row pixel, the relative displacement of the 4th row pixel and the relative displacement sum of the 5th row pixel；With Reason, the relative displacement that absolute offset values is the 7th row pixel of the 7th row pixel, the absolute offset values of eighth row pixel is then the 8th The relative displacement of row pixel and the relative displacement sum of the 7th row pixel, the absolute offset values of the 9th row pixel is then the 9th row The relative displacement of pixel, the relative displacement of eighth row pixel and the relative displacement sum of the 7th row pixel.

Preferably, for ease of calculating, described target line pixel is the described the first row pixel with reference to wave band, then the i-th row picture Element the relative displacement that the first absolute offset values is the i-th row pixel and the i-th row pixel before all row relative displacement it With.

Step S25: according to the first absolute offset values of described each row pixel with reference to wave band, move to be corrected by wave band Row corresponding in high-spectrum remote sensing, it is thus achieved that the high-spectrum remote sensing after the first correction；Wherein, the i-th row pixel moves Pixel count is the absolute value of the first absolute offset values of described the i-th row pixel with reference to wave band, each wave band the i-th row pixel described The direction of movement is identical with described offset direction with reference to corresponding to the first absolute offset values of the i-th row pixel of wave band；Wherein, I is the positive integer more than or equal to 1, i.e. i=1,2,3 ..., M-1；M is the line number of the pixel of high-spectrum remote sensing.

After mobile i-th row pixel, can be 0 by the pixel assignment removing position of the i-th row pixel.

In the embodiment of the present application, in each wave band of high-spectrum remote sensing to be corrected, the move mode of the i-th row pixel All identical with the move mode with reference to the i row pixel in wave band.For example, it is assumed that the first absolute offset values of the i-th row pixel is-3, Then the move mode of the i-th row pixel with reference to wave band is: move 3 location of pixels to negative offset direction set in advance, in like manner, In other wave band in high-spectrum remote sensing to be corrected, the i-th row pixel all moves 3 to described negative offset direction set in advance Individual location of pixels.

It is to say, in the embodiment of the present application, first determine the absolute offset values with reference to wave band each row pixel, then according to ginseng The side-play amount corresponding row of wave band each to high-spectrum remote sensing to be corrected of each row pixel examining wave band moves, thus realizes Correction to high-spectrum remote sensing geometric distortion.

A kind of high-spectrum remote sensing bearing calibration that the embodiment of the present application provides, from high-spectrum remote sensing to be corrected In determine that a wave band, for reference to wave band, determines the relative displacement of each row by the minimum error of adjacent lines, and by each row Relative displacement determine the absolute offset values (the most each row is relative to the side-play amount of target line pixel) of each row, thus according to each row Absolute offset values move in high-spectrum remote sensing to be corrected corresponding row, the high-spectrum remote-sensing after being corrected by wave band Image.It follows that the high-spectrum remote sensing bearing calibration that the embodiment of the present application provides, special from high-spectrum remote sensing itself Levy and set out, the related data between statistics row high-spectrum remote sensing within and row, by row correlated error minimum with line by line The method of position translation is corrected, and does not relies on external attitude data, it is achieved that for high-spectrum remote sensing geometric distortion Correction.

In embodiment illustrated in fig. 2, it is preferred that described according to default movement rule by the described i+1 row with reference to wave band Pixel is mobile pixel-by-pixel relative to described the i-th row pixel with reference to wave band be may include that

Wherein it is possible to determine that the pixel column position of the column position minimum of described i-th row pixel is described i+1 row picture The starting column location that element moves relative to described i-th row pixel；

That is, it may be determined that the column position at first pixel place of the i-th row pixel be i+1 row pixel relative to The starting column location that described i-th row pixel moves；The column position determining last pixel place of the i-th row pixel is The position, end column that i+1 row pixel moves relative to described i-th row pixel.

In the embodiment of the present application, by the pixel of column position maximum of described i+1 row pixel and rising of described i-th row pixel After pixel column alignment at beginning column position, move i+1 row pixel pixel-by-pixel.

In like manner, it is also possible to by the end column of pixel minimum for the column position of described i+1 row pixel with described i-th row pixel The pixel column alignment of position；

Described i+1 row pixel is moved pixel-by-pixel to described starting column location direction, until described i+1 row pixel The pixel of column position maximum is alignd with the pixel column at the starting column location of described i-th row pixel.

In above-described embodiment, it is preferred that in order to reduce amount of calculation, described determine that described i+1 row pixel is relative to described Starting column location and the another kind of implementation of position, end column that i-th row pixel moves can be:

Determine that the secondary series position that user inputs is that described i+1 row pixel moves relative to described i-th row pixel Position, end column.

Wherein, described first row position and secondary series position are empirically determined by expert of the art, described first row position Put the pixel column position minimum from the column position of described i-th row pixel different, described secondary series position and described i-th row picture The pixel column position difference that the column position of element is maximum, described first row position is different from described secondary series position.

In above-described embodiment, it is preferred that in order to reduce amount of calculation, described determine that described i+1 row pixel is relative to described Starting column location and another implementation of position, end column that i-th row pixel moves be:

Wherein, n is the number of pixels of the i-th row pixel；R is the degree of overlapping of described predetermined adjacent lines, 0 < r < 1.

Represent downward rounding operation；

The degree of overlapping of adjacent lines is in the case of high-spectrum remote sensing does not occurs geometric distortion, can ground between adjacent lines The number of pixels of thing correspondence continuously is divided by the width (i.e. the number of pixel in one-row pixels) of image.In the embodiment of the present application, institute State the degree of overlapping of predetermined adjacent lines can the corresponding continuously number of pixels of atural object to be removed by between the adjacent lines pre-estimated Obtain with the width of high-spectrum remote sensing.

It is exemplified below and i+1 row pixel is implemented process relative to what the i-th row pixel moved pixel-by-pixel, please join See the schematic diagram that Fig. 4, Fig. 4 move pixel-by-pixel relative to the i-th row pixel for the i+1 row pixel that the embodiment of the present application provides；

The number of pixels assuming every one-row pixels is 8, and the degree of overlapping of adjacent lines is 0.4, then, it may be determined that initial row position It is set toPosition, end column is 8-3+1=6, and therefore, in embodiment illustrated in fig. 4, starting column location is the i-th row picture The column position at R3 place in element, position, end column is the column position at R6 place in the i-th row pixel；

When mobile i+1 row pixel, the original position of i+1 row pixel is: the C8 pixel in i+1 row pixel and the R3 column alignment in i row pixel, then, i+1 row pixel pixel-by-pixel in the i-th row pixel the column position at R6 place move, directly R6 column alignment in R1 in i+1 row pixel and the i-th row pixel, the most no longer moves i+1 row pixel.

As seen from Figure 4, the number of times that i+1 row pixel moves is 11 times, and the corresponding number of times calculating error is also 11 Secondary；

And if using R1 column position in the i-th row pixel as starting column location, using R8 column position as termination Column position, then i+1 row pixel needs the number of times of movement to be 8*2=16 time, and the number of times of corresponding Error Calculation is also 16 times.

It should be noted that in Fig. 4 example shown, i+1 row only has a line, after this example is in order to illustrate to move, Position relationship between i+1 row pixel and the i-th row pixel, it is shown that multiple i+1 row pixels.

In order to optimize above-described embodiment further, high-spectrum remote sensing bearing calibration another that the embodiment of the present application provides A kind of flowchart as it is shown in figure 5, obtain first correction after high-spectrum remote sensing after, it is also possible to including:

Step S51: in the reference wave band of the high-spectrum remote sensing after correcting described first according to first user operation Draw the distortion contour curve of Target scalar；

Step S52: in the reference wave band of the high-spectrum remote sensing after correcting described first according to the second user operation Draw the correct contour curve of described Target scalar；

In the embodiment of the present application, by user with reference to wave band sketches out the atural object profile after being corrected, and correctly Atural object profile.Wherein, correct atural object profile can be checked according to high-resolution satellite image comparison by domain expert and know, Or can be determined by the photo gathered on the spot.And when the atural object profile of reality is more complicated, one can be laid in advance directly Wire tag band, so can delineate the profile being distorted of linear mark band and correct profile in remote sensing images.

Step S53: determine the 3rd row position of the intersection point of described the i-th row pixel with reference to wave band and described distortion contour curve Put, and determine the 4th column position of the intersection point of described the i-th row pixel with reference to wave band and described correct contour curve；

Step S54: obtain the second absolute offset values of described the i-th row pixel with reference to wave band, described with reference to the i-th of wave band Second absolute offset values of row pixel is the difference of described 3rd column position and described 4th column position；

The absolute value of the second absolute offset values of the i-th row pixel is i.e. the pixel count of the pixel deviation tram of the i-th row.

Step S55: according to the second absolute offset values of described i-th row pixel, by wave band by the height after described first correction I-th row pixel of spectral remote sensing image moves to the direction of described correct contour curve, it is thus achieved that the EO-1 hyperion after the second correction is distant Sense image；It is described second absolute offset values that described i-th row pixel moves amount of movement to the direction of described correct contour curve Absolute value.

Corresponding with embodiment of the method, the embodiment of the present application also provides for a kind of high-spectrum remote sensing correcting unit, this Shen Please embodiment provide high-spectrum remote sensing correcting unit a kind of structural representation as shown in Figure 6, may include that

First determines module 61, the first computing module 62, and second determines module 63, the second computing module 64 and the first correction Module 65；Wherein,

First determines that module 61 is for determining that described reference wave band is high-spectrum remote sensing to be corrected with reference to wave band All wave bands meet pre-conditioned wave band；

First computing module 62 for according to default movement rule by the described i+1 row pixel with reference to wave band relative to Described the i-th row pixel with reference to wave band moves pixel-by-pixel, often moves a location of pixels, calculates described reference according to the first formula The i+1 row pixel of wave band and the error amount of described the i-th row pixel with reference to wave band, described first formula is:

E=[∑ (R_p-C_q)²]/m²

Second determine module 63 for during by error amount minimum, described i+1 row pixel is relative to described i-th row pixel Side-play amount is defined as the relative displacement of described i+1 row pixel, and wherein, described i+1 row pixel is relative to described i-th row picture The side-play amount of element is positive offset amount or negative side-play amount, and described i+1 row pixel is relative to the taking of side-play amount of described i-th row pixel The positive offset direction that the positive and negative foundation of value is preset determines；

Second computing module 64 is used for calculating the first absolute offset values of described each row pixel with reference to wave band, including: meter The relative displacement sum of the calculation all row between the i-th row pixel and predetermined target line pixel, obtains the first He Value；First absolute offset values of described i-th row pixel be described first and the relative displacement of value and described i-th row pixel it With；

First correction module 65, for the first absolute offset values according to described each row pixel with reference to wave band, moves by wave band Row corresponding in dynamic high-spectrum remote sensing to be corrected, it is thus achieved that the high-spectrum remote sensing after the first correction；Wherein, the i-th row The pixel count that pixel moves is the absolute value of the first absolute offset values of described the i-th row pixel with reference to wave band, each ripple described The skew side corresponding to first absolute offset values of the direction moved of section the i-th row pixel and described the i-th row pixel with reference to wave band To identical；Wherein, i is the positive integer more than or equal to 1.

A kind of high-spectrum remote sensing correcting unit that the embodiment of the present application provides, from high-spectrum remote sensing to be corrected In determine that a wave band, for reference to wave band, determines the relative displacement of each row by the minimum error of adjacent lines, and by each row Relative displacement determine the absolute offset values (the most each row is relative to the side-play amount of the first row) of each row, thus exhausted according to each row Side-play amount is moved row corresponding in high-spectrum remote sensing to be corrected, the high-spectrum remote-sensing figure after being corrected by wave band Picture.It follows that the high-spectrum remote sensing correcting unit that the embodiment of the present application provides, from the feature of high-spectrum remote sensing own Setting out, the related data between statistics row and row within high-spectrum remote sensing, by row correlated error minimum and position line by line The method that horizontalization moves is corrected, and does not relies on external attitude data, it is achieved that for high-spectrum remote sensing geometric distortion Correction.

In above-described embodiment, it is preferred that a kind of structural representation of described first computing module 62 is as it is shown in fig. 7, permissible Including:

Mover module 71 and calculating sub module 72；Wherein,

Described mover module 71 is used for relative for the described i+1 row pixel with reference to wave band according to default movement rule Move pixel-by-pixel in described the i-th row pixel with reference to wave band；Specifically may include that

First determines unit 711, for determining what described i+1 row pixel moved relative to described i-th row pixel Starting column location and position, end column；

Mobile unit 712, for the pixel that the column position of described i+1 row pixel is maximum and described i-th row pixel Pixel column alignment at starting column location, moves to described end column locality pixel-by-pixel by described i+1 row pixel, until The pixel of the column position minimum of described i+1 row pixel is alignd with the pixel column of the position, end column of described i-th row pixel；

Described calculating sub module 72, for often moving a location of pixels in the described i+1 row pixel with reference to wave band, is pressed The described i+1 row pixel with reference to wave band and the error amount of described the i-th row pixel with reference to wave band is calculated according to the first formula, described First formula is:

E=[∑ (R_p-C_q)²]/m²

Above-described embodiment, it is preferred that the first of the embodiment of the present application offer determines a kind of structural representation of unit 711 such as Shown in Fig. 8, may include that

First determines subelement 81, for determining that the pixel column position of the column position minimum of described i-th row pixel is The starting column location that described i+1 row pixel moves relative to described i-th row pixel；

Second determines subelement 82, for determining that the pixel column position of the column position maximum of described i-th row pixel is The position, end column that described i+1 row pixel moves relative to described i-th row pixel.

Above-described embodiment, it is preferred that the first of the embodiment of the present application offer determines the another kind of structural representation of unit 711 As shown in Figure 9, it is also possible to including:

3rd determines subelement 91, for determine the first row position that user inputs be described i+1 row pixel relative to The starting column location that described i-th row pixel moves；

4th determines subelement 92, for determine the secondary series position that user inputs be described i+1 row pixel relative to The column position minimum of the position, end column that described i-th row pixel moves, described first row position and described i-th row pixel Pixel column position is different, the pixel column position that described secondary series position is maximum with the column position of described i-th row pixel Difference, described first row position is different from described secondary series position.

Above-described embodiment, it is preferred that the first of the embodiment of the present application offer determines another structural representation of unit 711 As shown in Figure 10, may include that

5th determines subelement 101, for determining described i+1 row pixel according to the degree of overlapping of predetermined adjacent lines The starting column location moved relative to described i-th row pixel and position, end column, wherein,

Wherein, n is the number of pixels of the i-th row pixel；R is the degree of overlapping of described adjacent lines, 0 < r < 1.

Above-described embodiment, it is preferred that the another kind of knot of the high-spectrum remote sensing correcting unit that the embodiment of the present application provides Structure schematic diagram is logical as shown in figure 11, it is also possible to including:

First drafting module 111, the second drafting module 112, the 3rd determines module 113, acquisition module 114 and the second correction Module 115；Wherein,

First drafting module 111 high-spectrum remote sensing after correcting described first according to first user operation With reference to the distortion contour curve drawing Target scalar in wave band；

Second drafting module 112 high-spectrum remote sensing after correcting described first according to the second user operation With reference to the correct contour curve drawing described Target scalar in wave band；

3rd determines that module 113 is for determining described the i-th row pixel with reference to wave band and the friendship of described distortion contour curve 3rd column position of point, and determine the 4th row of the intersection point of described the i-th row pixel with reference to wave band and described correct contour curve Position；

Acquisition module 114 is for obtaining the second absolute offset values of described the i-th row pixel with reference to wave band, described reference wave The difference that second absolute offset values is described 3rd column position and described 4th column position of the i-th row pixel of section；

Second correction module 115 is for the second absolute offset values according to described i-th row pixel, by wave band by described first I-th row pixel of the high-spectrum remote sensing after correction moves to the direction of described correct contour curve, it is thus achieved that after the second correction High-spectrum remote sensing；It is described second exhausted that described i-th row pixel moves amount of movement to the direction of described correct contour curve Absolute value to side-play amount.

It should be noted that the high-spectrum remote sensing bearing calibration of the embodiment of the present application offer and device, except permissible Carry out on application unmanned plane outside the geometric correction of high-spectrum remote sensing, it is also possible to be applied to other aerial platform (if any man-machine, Dirigible, balloon etc.) geometric correction of upper high-spectrum remote sensing.

It addition, the high-spectrum remote sensing after being corrected by the embodiment of the present application can also utilize attitude recording equipment record Data further corrected, to improve calibration result further.

Below as a example by high-spectrum remote sensing shown in Fig. 1, the calibration result of the embodiment of the present application is described, such as Figure 12 institute Show, Figure 12 be by the embodiment of the present application to obtain after shown in Fig. 1, high-spectrum remote sensing is corrected first correction after Design sketch.

Figure 13 is that the distortion contour curve sketched in high-spectrum remote sensing shown in Figure 12 of user is (such as song thicker in figure Line, i.e. L1) and correct contour curve (such as straight line thinner in figure, i.e. L2)；

Described above to the disclosed embodiments, makes professional and technical personnel in the field be capable of or uses the present invention. Multiple amendment to these embodiments will be apparent from for those skilled in the art, as defined herein General Principle can realize without departing from the spirit or scope of the present invention in other embodiments.Therefore, the present invention It is not intended to be limited to the embodiments shown herein, and is to fit to and principles disclosed herein and features of novelty phase one The widest scope caused.

Claims

1. a high-spectrum remote sensing bearing calibration, it is characterised in that including:

Determine with reference to wave band, described with reference to wave band be high-spectrum remote sensing to be corrected all wave bands in meet pre-conditioned Wave band；

According to default movement rule by the described i+1 row pixel with reference to wave band relative to described the i-th row picture with reference to wave band Element moves pixel-by-pixel, often moves a location of pixels, calculates the described i+1 row pixel with reference to wave band and institute according to the first formula Stating the error amount of the i-th row pixel with reference to wave band, described first formula is:

E=[∑ (R_p-C_q)²]/m²

Wherein, E is the described i+1 row pixel with reference to wave band and the error amount of described the i-th row pixel with reference to wave band, R_pIt is i-th The value of pth pixel, C in row pixel_qFor the value of q-th pixel in i+1 row pixel, q-th picture in described i+1 row pixel In plain and described i-th row pixel, pth pixel column aligns；M is in i+1 row pixel to be the pixel of column alignment with the i-th row pixel Number；

During by error amount minimum, described i+1 row pixel is defined as described i+1 relative to the side-play amount of described i-th row pixel The relative displacement of row pixel, wherein, described i+1 row pixel is positive offset amount relative to the side-play amount of described i-th row pixel Or negative side-play amount, described i+1 row pixel is preset relative to the positive and negative foundation of the value of the side-play amount of described i-th row pixel Positive offset direction determines；

Calculate the first absolute offset values of described each row pixel with reference to wave band, including: calculate and be positioned at the i-th row pixel with the most true The relative displacement sum of the fixed all row between target line pixel, obtains first and value；The first of described i-th row pixel is exhausted It it is the relative displacement sum of described first and value with described i-th row pixel to side-play amount；

According to the first absolute offset values of described each row pixel with reference to wave band, move high-spectrum remote-sensing figure to be corrected by wave band Row corresponding in Xiang, it is thus achieved that the high-spectrum remote sensing after the first correction；Wherein, the pixel count that the i-th row pixel moves is described The absolute value of the first absolute offset values of the i-th row pixel of reference wave band, the direction that each wave band the i-th row pixel moves is with described Identical with reference to the offset direction corresponding to the first absolute offset values of the i-th row pixel of wave band；Wherein, i is more than or equal to 1 Positive integer.

Method the most according to claim 1, it is characterised in that described according to default movement rule by described with reference to wave band I+1 row pixel mobile pixel-by-pixel relative to described the i-th row pixel with reference to wave band include:

Determine starting column location and position, end column that described i+1 row pixel moves relative to described i-th row pixel；

By the pixel column at the pixel maximum for the column position of the described i+1 row pixel starting column location with described i-th row pixel Alignment；

Described i+1 row pixel is moved pixel-by-pixel to described end column locality, until the row position of described i+1 row pixel The pixel putting minimum is alignd with the pixel column of the position, end column of described i-th row pixel.

Method the most according to claim 2, it is characterised in that described determine that described i+1 row pixel is relative to described i-th Starting column location and position, end column that row pixel moves include:

Determine that the minimum pixel column position of the column position of described i-th row pixel is that described i+1 row pixel is relative to described The starting column location that i-th row pixel moves；

Determine that the maximum pixel column position of the column position of described i-th row pixel is that described i+1 row pixel is relative to described The position, end column that i-th row pixel moves.

Determine that the first row position that user inputs is that described i+1 row pixel is relative to rising that described i-th row pixel moves Beginning column position；

Determine that the secondary series position that user inputs is the end that described i+1 row pixel moves relative to described i-th row pixel Only column position；

Wherein, the minimum pixel column position of described first row position and the column position of described i-th row pixel is different, and described the Two column positions are different from the pixel column position of the column position of described i-th row pixel maximum, and described first row position is with described Secondary series position is different.

Determine that described i+1 row pixel is moved relative to described i-th row pixel according to the degree of overlapping of predetermined adjacent lines Dynamic starting column location and position, end column, wherein,

6. according to the method described in claim 1-5 any one, it is characterised in that the EO-1 hyperion after obtaining the first correction is distant After sense image, also include:

Reference wave band according to first user operation high-spectrum remote sensing after described first correction draws Target scalar Distortion contour curve；

Reference wave band according to second user operation high-spectrum remote sensing after described first correction draws described target The correct contour curve of atural object；

Determine the 3rd column position of the intersection point of described the i-th row pixel with reference to wave band and described distortion contour curve, and determine institute State the 4th column position of the i-th row pixel with reference to wave band and the intersection point of described correct contour curve；

Obtain the second absolute offset values of described the i-th row pixel with reference to wave band, the second of described the i-th row pixel with reference to wave band Absolute offset values is the difference of described 3rd column position and described 4th column position；

According to the second absolute offset values of described i-th row pixel, by wave band by the high-spectrum remote sensing after described first correction The i-th row pixel move to the direction of described correct contour curve, it is thus achieved that second correction after high-spectrum remote sensing；Described I row pixel moves the absolute value that amount of movement is described second absolute offset values to the direction of described correct contour curve.

7. a high-spectrum remote sensing correcting unit, it is characterised in that including:

First determines module, is used for determining with reference to wave band, and described is all of high-spectrum remote sensing to be corrected with reference to wave band Wave band meets pre-conditioned wave band；

First computing module, for according to default movement rule by the described i+1 row pixel with reference to wave band relative to described I-th row pixel of reference wave band moves pixel-by-pixel, often moves a location of pixels, calculates described with reference to wave band according to the first formula The error amount of i+1 row pixel and described the i-th row pixel with reference to wave band, described first formula is:

E=[∑ (R_p-C_q)²]/m²

Second determines module, and when by error amount minimum, described i+1 row pixel is relative to the skew of described i-th row pixel Amount is defined as the relative displacement of described i+1 row pixel, and wherein, described i+1 row pixel is relative to described i-th row pixel Side-play amount is positive offset amount or negative side-play amount, and described i+1 row pixel is relative to the value of the side-play amount of described i-th row pixel The positive offset direction that positive and negative foundation is preset determines；

Second computing module, for calculating the first absolute offset values of described each row pixel with reference to wave band, including: calculate i-th The relative displacement sum of all row between row pixel and predetermined target line pixel, obtains first and value；Described i-th First absolute offset values of row pixel is the relative displacement sum of described first and value with described i-th row pixel；

First correction module, for the first absolute offset values according to described each row pixel with reference to wave band, moves by wave band and treats Row corresponding in the high-spectrum remote sensing of correction, it is thus achieved that the high-spectrum remote sensing after the first correction；Wherein, the i-th row pixel The pixel count of movement is the absolute value of the first absolute offset values of described the i-th row pixel with reference to wave band, each wave band the i-th row picture The direction of element movement is identical with described offset direction with reference to corresponding to the first absolute offset values of the i-th row pixel of wave band；Its In, i is the positive integer more than or equal to 1.

Device the most according to claim 7, it is characterised in that described first computing module includes: mover module and meter Operator module；Wherein,

Described mover module for according to default movement rule by the described i+1 row pixel with reference to wave band relative to described I-th row pixel of reference wave band moves pixel-by-pixel；Including:

First determines unit, for determining the initial row position that described i+1 row pixel moves relative to described i-th row pixel Put and position, end column；

Mobile unit, for the initial row position of the pixel that the column position of described i+1 row pixel is maximum with described i-th row pixel Put the pixel column alignment at place, described i+1 row pixel is moved pixel-by-pixel to described end column locality, until described i-th+ The pixel of the column position minimum of 1 row pixel is alignd with the pixel column of the position, end column of described i-th row pixel；

Described calculating sub module is for often moving a location of pixels, according to first in the described i+1 row pixel with reference to wave band Formula calculates the described i+1 row pixel with reference to wave band and the error amount of described the i-th row pixel with reference to wave band, and described first is public Formula is:

E=[∑ (R_p-C_q)²]/m²

Wherein, E is the described i+1 row pixel with reference to wave band and the error amount of described the i-th row pixel with reference to wave band, R_pIt is i-th The value of pth pixel, C in row pixel_qFor the value of q-th pixel in i+1 row pixel, q-th picture in described i+1 row pixel In plain and described i-th row pixel, pth pixel column aligns；M is in i+1 row pixel to be the pixel of column alignment with the i-th row pixel Number.

Device the most according to claim 8, it is characterised in that described first determines that unit includes:

First determines subelement, is described i-th for determining the pixel column position of the column position minimum of described i-th row pixel The starting column location that+1 row pixel moves relative to described i-th row pixel；

Second determines subelement, is described i-th for determining the pixel column position of the column position maximum of described i-th row pixel The position, end column that+1 row pixel moves relative to described i-th row pixel.

3rd determines subelement, for determining that the first row position that user inputs is that described i+1 row pixel is relative to described i-th The starting column location that row pixel moves；

4th determines subelement, for determining that the secondary series position that user inputs is that described i+1 row pixel is relative to described i-th The position, end column that row pixel moves；

11. devices according to claim 8, it is characterised in that described first determines that unit includes:

5th determines subelement, for according to the degree of overlapping of predetermined adjacent lines determine described i+1 row pixel relative to Starting column location that described i-th row pixel moves and position, end column, wherein,

12. according to the device described in claim 7-11 any one, it is characterised in that also include:

First drafting module, the reference wave of the high-spectrum remote sensing after correcting described first according to first user operation The distortion contour curve of Target scalar is drawn in Duan；

Second drafting module, the reference wave of the high-spectrum remote sensing after correcting described first according to the second user operation The correct contour curve of described Target scalar is drawn in Duan；

3rd determines module, for determining the of described the i-th row pixel with reference to the wave band intersection point with described distortion contour curve Three column positions, and determine the 4th column position of the intersection point of described the i-th row pixel with reference to wave band and described correct contour curve；

Acquisition module, for obtaining the second absolute offset values of described the i-th row pixel with reference to wave band, described with reference to the of wave band Second absolute offset values of i row pixel is the difference of described 3rd column position and described 4th column position；

Second correction module, for the second absolute offset values according to described i-th row pixel, by wave band by after described first correction The i-th row pixel of high-spectrum remote sensing move to the direction of described correct contour curve, it is thus achieved that the Gao Guang after the second correction Spectrum remote-sensing image；It is described second absolute drift that described i-th row pixel moves amount of movement to the direction of described correct contour curve The absolute value of amount.