CN104766328A - Remote-sensing image Bowtie effect correction method based on path tracing - Google Patents
Remote-sensing image Bowtie effect correction method based on path tracing Download PDFInfo
- Publication number
- CN104766328A CN104766328A CN201510181391.6A CN201510181391A CN104766328A CN 104766328 A CN104766328 A CN 104766328A CN 201510181391 A CN201510181391 A CN 201510181391A CN 104766328 A CN104766328 A CN 104766328A
- Authority
- CN
- China
- Prior art keywords
- image
- longitude
- latitude
- scanning strip
- 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
Abstract
The invention discloses a remote-sensing image Bowtie effect correction method based on path tracing. The method includes the steps that interpolation is conducted, and longitude and latitude data have the same distinguishability as a raw image; a search starting point of a first correction image element is determined according to the lifting and falling track direction and the scanning direction of a satellite, a resampling sequence conducted on a correction image is set, locations of the image elements of the correction image in the raw image are positioned in sequence by using a path tracing mode, namely the locations of the image elements in a longitude and latitude grid are searched after the interpolation is conducted, and a plurality of sampling points, most close to the image elements of the correction image, are acquired to be used for gray level resampling to remove the image overlapping phenomenon. By means of the method, the calculating efficiency is good, the raw image information is taken full advantage, and the method is a strict high-precision low-medium distinguishability satellite image geometric correction method.
Description
Technical field
The present invention relates to remote sensing image technical field, especially relate to a kind of remote sensing image Bowtie effect correcting method based on path tracing.
Background technology
Bow-tie effect is again Bowtie phenomenon, is the data overlap phenomenon occurred in a kind of remote sensing images.Bowtie effect is prevalent in MODIS, AVHRR, FY-3 etc. and adopts and morely visit first cross rail and to sweep in the polar orbit of mode imaging in low resolution satellite image, it is when detector transversal scanning exceedes certain angle, because detector is to the joint effect of the factors such as the shake in the curvature of the visual field geometrical property of earth observation, earth surface, topographic relief and detector motion, between adjacent scan bands ground coverage, there is overlapping phenomenon.
At present, Bowtie effect minimizing technology mainly contains ephemeris method and non-ephemeris method.Ephemeris method mainly generates study plot reason grid according to the ephemeris of satellite, then data is projected on this grid according to its geographic coordinate and mates, while geography calibration, eliminate Bowtie effect.This method needs the ephemeris of satellite, and Measures compare is complicated, and the data product that user takes usually does not comprise ephemeris, therefore less employing.
Under the prerequisite not having satellite ephemeris, scholar both domestic and external has done a lot of relevant research to the Bowtie phenomenon removing satellite image.Guo Guangmeng (2003), Yu Junhui etc. (2004), Zhang Juan etc. (2009), Li Zhen etc. (2010), Wang Hanyu etc. (2014) propose to utilize the method for the correction of image to determine the repetition line number of MODIS image adjacent scan bands, then remove the Bowtie effect of image by the method for resampling; Liu Liangming etc. (2007) propose the method for the quick removal Bowtie effect of Corpus--based Method, utilize detector from left to right in scanning process the repetition line number on both sides about the hypothesis of substar almost symmetry, adopt a symmetric function to carry out matching, then remove the Bowtie effect of other image according to this rule.Not strict on these theoretical methods, and the problem of repetition line number instability cannot be avoided.Some researcher proposes to utilize longitude and latitude data subsidiary in satellite data file to assist and removes Bowtie effect.Cheng Liang etc. (2005,2007) determine the repetition line number of MODIS image according to MOD03 latitude data, then adopt fractal interpolation, Kriging interpolation method carries out resampling to image data.Song Shasha etc. (2010) then adopt net function interpolation to carry out interpolation to MODIS latitude data, calculate the repetition line number of image according to the latitude data after interpolation.These methods all only remove the pixel of repetition on raw video, overlay region is filled up again by the method for gray scale interpolation, there is no to solve the problem from image central authorities to the non-linear increase of image both sides pixel ground coverage, and subsequent applications still needs to carry out geometry correction process.The Modistools module (2004) of R & D ScanEx company of Russia exploitation utilizes geometry correction to remove Bowtie effect, but must use under ENVI environment, and travelling speed is very slow; Jiang Gengming etc. (2004), Xu Meng etc. (2005), after Liang Zhihua (2012) all adopts pair warp and weft degrees of data to carry out interpolation, the method that recycling forward projection and rear orientation projection combine, Bowtie effect is removed while carrying out geometry correction, but calculation of complex, needs to consume a large amount of time.
Summary of the invention
For not good, the slow-footed shortcoming of existing method Bowtie effect removal effect, the invention provides a kind of remote sensing image Bowtie effect correcting method based on path tracing.
Technical solution of the present invention provides a kind of remote sensing image Bowtie effect correcting method based on path tracing, raw video for having Bowtie effect to arbitrary width is corrected, generate the correction image eliminating Bowtie effect, described raw video adopts many spy units and sweeps mode imaging, be spliced up and down by several scanning strips, the line number of each scanning strip equals the number visiting unit, comprises the following steps
Step 1, to the situation of longitude and latitude data resolution subsidiary in satellite data file lower than the resolution of satellite image, carries out interpolation, makes longitude and latitude data and raw video have identical resolution;
Step 2, determines first search starting point S0 correcting pixel according to the lift rail direction of satellite and direction of scanning, comprises the center of the longitude and latitude graticule mesh corresponding to maximum for minimum for longitude in raw video, latitude a jiao as S0; The current location S=S0 of order search, the longitude and latitude graticule mesh claiming S place is current mesh, and the scanning strip at S place is Strip1;
Step 3, set and adopt the order of serpentine to carry out resampling by column to correcting image, getting according to order the pixel correcting image is current correction pixel, with the impact point T of the latitude and longitude coordinates of this pixel for search;
Step 4, from S in the scanning strip Strip1 of raw video, carries out first time path tracing to T; If follow the trail of result be S and T in same longitude and latitude graticule mesh, then find the target mesh at T place, forward Step.5 to, if S follows the trail of the border of arrival Strip1 longitude and latitude graticule mesh and T outside border, forward Step.9 to;
Step 5, according to lift rail direction and the direction of scanning of satellite, and current resampling pixel is in the parity of correcting the row in image, determines new scanning strip Strip2 T being carried out to second time path tracing, and the current location S of more new search;
Step 6, from S in the scanning strip Strip2 of raw video, carries out second time path tracing to T; If follow the trail of result be S and T in same longitude and latitude graticule mesh, then find the target mesh at T place, forward Step.7 to, if the border of S arrival Strip2 longitude and latitude graticule mesh and T, outside border, forward Step.8 to;
Step 7, utilizes 8 pixels the most contiguous with T, and after carrying out gray resample by inverse distance weighted method, assignment is to correcting image, and upgrades the center that S is the target mesh at T place, forwards Step.10 to;
Described 8 pixels the most contiguous with T, be T in the overlay region of raw video, the target mesh respectively searched in Strip1 and Strip2, gets summit totally 8 points of these two target mesh;
Step 8, utilizes 4 pixels the most contiguous with T, and after carrying out gray resample by inverse distance weighted method, assignment is to correcting image, and upgrades the center that S is the target mesh at T place, forwards Step.10 to;
Described 4 pixels the most contiguous with T, be T in the non-overlapped district of raw video, only in Strip1, search a target mesh, get summit totally 4 points of this target mesh;
Step 9, if T is outside raw video border, gives and corrects image tax invalid value, and retain the position of S in scanning strip Strip1, forward Step.10 to; If T is on other scanning strip, then Strip1 is updated to the scanning strip that its next one is new, redefines the position of S, return step 4;
Step 10, judge whether the correction of view picture image, if do not complete, the current location S searched at the end of retaining current correction pixel resampling is constant, return step 3 and get that to correct the next pixel of image be current correction pixel in order, with the impact point T of the latitude and longitude coordinates of this pixel for search, carry out the resampling of next pixel, if complete, end loop.
And when carrying out interpolation in step 1, utilize 4 adjacent points to carry out the encryption of cubic curve interpolation in X-direction, before and after utilizing in the Y direction, 2 points carry out linear interpolation encryption.
And, in step 4 and step 6, path trace mode is, judge the line segment ST working direction of determining searching route whether crossing with four limits of current mesh, if with wherein crossing on one side, searching route enters the adjacent graticule mesh on this limit, and upgrade the center that S is adjacent graticule mesh, repeatedly carry out judgement search operation.
And path trace mode is in step 5 and step 9,
New scanning strip determination mode is as follows,
1. judge the characteristic of raw video, comprising the some life in the upper left corner in raw video is No. 0 point, and the point in the upper right corner, the lower right corner, the lower left corner is followed successively by 1,2, No. 3 point in the direction of the clock; If the maximum point of raw video angle point middle latitude is No. 0 point or No. 1 point, make variable Nextscan=1; If the maximum point of raw video angle point middle latitude is No. 2 points or No. 3 points, make variable Nextscan=-1;
2. judge that current correction pixel is in the parity of correcting column in image, if at odd column, makes variable order=1, if at even column, makes variable order=-1;
3. determine that new scanning strip is Strip1+order × Nextscan;
Current location S update mode is as follows,
If search in the above-mentioned scanning strip determined for the first time, then in two kinds of situation,
If 1. order × Nextscan=1, then upgrading S is the graticule mesh center that in the first row of this scanning strip, longitude is minimum;
If 2. order × Nextscan=-1, then upgrading S is the graticule mesh center that in last column of this scanning strip, longitude is minimum;
If carried out search in the above-mentioned scanning strip determined, then upgrading S was the last position that S finally reaches when searching in this scanning strip.
The present invention proposes a kind of quick minimizing technology of Bowtie effect based on path tracing, comprise the interpolation of coordinate to original longitude and latitude data, longitude and latitude data are made to have identical resolution with raw video data, the mode of path tracing is adopted to locate the position of each pixel in raw video of correcting image, namely search for the position in its longitude and latitude graticule mesh after interpolation, obtain several sampled points the most contiguous with correcting image picture element and be used for the overlapping phenomenon that gray resample removes image.This method counting yield is good, makes full use of raw video information, is low resolution satellite image geometric correction method in a kind of strict high precision.
Accompanying drawing explanation
Fig. 1 schematic diagram that to be the embodiment of the present invention search for based on the some position of path tracing.
Fig. 2 is embodiment of the present invention process flow diagram.
Specific implementation method
Understand for the ease of those of ordinary skill in the art and implement the present invention, below in conjunction with drawings and Examples, the present invention is described in further detail, should be appreciated that exemplifying embodiment described herein is only for instruction and explanation of the present invention, is not intended to limit the present invention.
Embodiment is corrected the raw video that arbitrary width has Bowtie effect, generates the correction image eliminating Bowtie effect.Raw video adopts many spy units and sweeps mode imaging, is spliced up and down by several scanning strips, and the line number of each scanning strip equals the number visiting unit.Such as, the MODIS image of 1km resolution adopts 10 visit unit and sweep imaging, and the size of a scape raw video is 2030 × 1354, is spliced up and down by 203 scanning strips, each scanning with 10 row, often row 1354 pixels.Raw video the 1st row is to the 10th behavior the 1st scanning strip, and the 11st row is to the 20th behavior the 2nd scanning strip, and the rest may be inferred.
See Fig. 2, the method that embodiment provides comprises the following steps, and during concrete enforcement, those skilled in the art can adopt software mode to realize automatically running:
Step.1, to the situation of longitude and latitude data resolution subsidiary in satellite data file lower than the resolution of satellite image, carries out interpolation, makes longitude and latitude data and raw video have identical resolution;
At present, low resolution satellite provides 250m, 500m and 1000m image of tri-kinds of resolution usually, because the longitude and latitude graticule mesh that satellite image provides is 1000m resolution, therefore first can judge before search that longitude and latitude data resolution subsidiary in satellite data file is whether lower than the situation of the resolution of satellite image, that longitude and latitude data interpolating be become have identical resolution with raw video, such as, be 250m or 500m.Consider the imaging characteristic of middle low resolution satellite, satellite flies at a constant speed in orbital direction, and the longitude and latitude of Y-direction is changed to linear change; And the impact in vertical track direction due to earth curvature, from image central authorities to image both sides, the rate of change of longitude and latitude becomes large gradually, and namely the longitude and latitude change of X-direction is nonlinearities change.The present invention proposes further, during interpolation, utilizes 4 adjacent points to carry out the encryption of cubic curve interpolation in X-direction, and before and after utilizing in the Y direction, 2 points carry out linear interpolation encryption, and concrete interpolation is embodied as prior art, and it will not go into details in the present invention.After X, Y-direction are encrypted to graticule mesh, the point of graticule mesh inside is adopted to the method interpolation of net function (Song Shasha, 2010).Consider the phenomenon that there is picture-dot interlacing between scanning strip, when carrying out row encryption to each scanning strip last column, the data of next scanning strip can not be utilized to carry out interpolation, and the data of Current Scan band should be utilized to carry out linear extrapolation.
Step.2 determines first search starting point S0 correcting pixel according to the lift rail direction of satellite and direction of scanning, the current location S=S0 of order search, and the longitude and latitude graticule mesh claiming S place is current mesh, and the scanning strip at S place is Strip1;
Due to the difference of satellite lift rail direction and direction of scanning, may cause between the impact point T of S0 and search as S0 distant using the center of the longitude and latitude graticule mesh corresponding to a jiao (upper left corner as image) that raw video is fixing, affect the counting yield of method.Therefore the present invention proposes further, using the center of the longitude and latitude graticule mesh corresponding to maximum for minimum for longitude in raw video, latitude a jiao as S0, specifically, be the upper right corner of image for falling the satellite that rail scans from left to right, be the upper left corner of image for falling the satellite that rail scans from right to left, the satellite scanned from left to right for rail lift is the lower left corner of image, and the satellite scanned from right to left for rail lift is the lower right corner of image.Such selection S0, then the correction image picture element correction order coordinating serpentine in Step.3, just can make the shortest path searched for.
Step.3 setting adopts the order of serpentine to carry out resampling by column to correcting image, getting according to order the pixel correcting image is current resampling handling object (i.e. current correction pixel), with the impact point T of the latitude and longitude coordinates of this pixel for search;
When the row correcting image carry out resampling from top to bottom, the latitude of T reduces gradually, and the latitude of S also reduces gradually in the process of path tracing, so when entering next column resampling, in order to shorten ST, the correction of pixel just should adopt mode from bottom to top to carry out.Therefore, the present invention adopts the mode of serpentine to correct correction image picture element by column, and first row is corrected from top to bottom, and secondary series is corrected from bottom to top, and the 3rd row are corrected from top to bottom again, and the rest may be inferred.
Step.4 from S, carries out first time path tracing to T in the scanning strip Strip1 of raw video.Judge the line segment ST working direction of determining searching route whether crossing with four limits of current mesh, if with wherein crossing on one side, searching route enters the adjacent graticule mesh on this limit, and upgrade the center that S is adjacent graticule mesh, repeatable operation is until S and T (namely finds the target mesh at T place) in same longitude and latitude graticule mesh, forward Step.5 to, if S arrives the border of Strip1 longitude and latitude graticule mesh and T, outside border, (namely T is not in the areas imaging of scanning strip Strip1, S and T cannot be made in same longitude and latitude graticule mesh), forward Step.9 to;
The principle of path tracing is that if with wherein crossing, searching route enters the adjacent graticule mesh on this limit, and upgrades the center that S is adjacent graticule mesh, as shown in Figure 1 by judging the line segment ST working direction determining searching route whether crossing with four limits of current mesh.
Wherein, line segment ST and current mesh P is judged
1p
2p
3p
4certain graticule mesh limit P
1p
2whether there is the condition of intersection point to be formula (1) and formula (2), if meet the requirement of formula (1) and formula (2) simultaneously, then line segment ST and line segment P is described
1p
2there is intersection point, in like manner calculate the relation on ST and other graticule mesh limit.
Search operation is judged more than repeating, constantly update the current location S of search, S levels off to T gradually along searching route, until S and T (namely finds the target mesh at T place) in same longitude and latitude graticule mesh or S arrives the border of Strip1 longitude and latitude graticule mesh and T (namely T is not in areas imaging of scanning strip Strip1) outside border.Finally, the path of S process is the searching route of T, path as illustrated by the arrows in fig. 1.
Step.5 is according to the lift rail direction of satellite and direction of scanning, and current resampling pixel is in the parity of correcting the row in image, determine the new scanning strip (being designated as Strip2) T being carried out to second time path tracing, and the current location S of more new search;
(1) determination of second time path tracing scanning strip
Due to satellite lift rail direction and the difference of correcting image picture element correction order, the scanning strip position of second time path tracing is not fixed, but can be judged by the characteristic of raw video and the parity of correcting image picture element column, and concrete determination methods is as follows:
1. the characteristic of raw video is judged;
Be No. 0 point by the some life in the upper left corner in raw video, the point in the upper right corner, the lower right corner, the lower left corner is followed successively by 1,2, No. 3 point in the direction of the clock.If the maximum point of raw video angle point middle latitude is No. 0 point or No. 1 point, illustrate raw video from top to bottom latitude reduce gradually, make variable Nextscan=1; If the maximum point of raw video angle point middle latitude is No. 2 points or No. 3 points, illustrate raw video from top to bottom latitude increase gradually, make variable Nextscan=-1.
2. judge that current correction pixel is in the parity of correcting column in image;
Correct because correction image have employed serpentine order, the parity of scanning strip position also with current correction pixel column of second time path tracing is relevant.If current correction pixel, correcting the odd column of image, makes variable order=1, from top to bottom, the latitude of T reduces correction order gradually; If current correction pixel, correcting the even column of image, makes variable order=-1, from bottom to top, the latitude of T increases correction order gradually.
3. determine the scanning strip of second time search, general scanning strip adopts reel number mark, i.e. Strip2=Strip1+order × Nextscan.
(2) determination of the search current location S of second time path tracing
If search in the above-mentioned scanning strip determined for the first time, then in two kinds of situation:
If 1. order × Nextscan=1, then upgrade the graticule mesh center that S is longitude minimum (judging according to the angle point that longitude is minimum) in the first row of scanning strip Strip2;
If 2. order × Nextscan=-1, then upgrading S is the graticule mesh center that in last column of scanning strip Strip2, longitude is minimum.
If carried out search in scanning strip Strip2, then upgrading S was the last position that S finally reaches when searching in Strip2.
Step.6 from S, carries out second time path tracing to T in the scanning strip Strip2 of raw video.The method of path tracing and identical in Step.4, until S and T (namely finds the target mesh at T place) in same longitude and latitude graticule mesh, forward Step.7 to.If S arrives the border of Strip2 longitude and latitude graticule mesh and T (namely T is not in areas imaging of scanning strip Strip2) outside border, forward Step.8 to.
Step.7 utilizes 8 pixels the most contiguous with T, carries out gray resample by inverse distance weighted method, its assignment is given and corrects image, and upgrade the center that S is the target mesh at T place, forward Step.10 to;
Pixel is corrected and is performed Step.7, illustrate that T is in the overlay region of raw video, a target mesh can be respectively searched in Strip1 and Strip2, get the summit of these two target mesh, press inverse distance weighted method and carry out gray resample for totally 8, inverse distance weighted method is prior art, and it will not go into details in the present invention.
Step.8 utilizes 4 pixels the most contiguous with T, carries out gray resample by inverse distance weighted method, its assignment is given and corrects image, and upgrade the center that S is the target mesh at T place, forward Step.10 to;
Pixel is corrected and is performed to Step.8, illustrates that T is in the non-overlapped district of raw video, only can search a target mesh in Strip1, get the summit of this target mesh, presses and carries out gray resample against distance weighted method for totally 4.
If Step.9 T is outside raw video border, gives and correct image tax invalid value, and retain the position of S in scanning strip Strip1, forward Step.10 to.If T is on other scanning strip, then Strip1 is updated to the scanning strip that its next one is new, redefines the position of S, return step 4;
If T is on other scanning strip, the defining method of the scanning strip that the next one of Strip1 is new is identical with the defining method of Strip2 in Step.5, i.e. identical also with Step.5 of the method that redefines of Strip1=Strip1+order × Nextscan, S.
If 1. order × Nextscan=1, then upgrade the graticule mesh center that S is longitude minimum (judging according to the angle point that longitude is minimum) in the first row of new scanning strip Strip1;
If 2. order × Nextscan=-1, then upgrading S is the graticule mesh center that in last column of new scanning strip Strip1, longitude is minimum.
If carried out search in new scanning strip Strip1, then upgrading S was the last position that S finally reaches when searching in Strip1.
Step.10 has judged whether the correction of view picture image, if do not complete, the current location S searched at the end of retaining current correction pixel resampling is constant, return Step.3 and get that to correct the next pixel of image be current correction pixel in order, with the impact point T of its latitude and longitude coordinates for search, carry out the resampling of next pixel, if complete, end loop.
Specific embodiment described herein is only to the explanation for example of the present invention's spirit.Those skilled in the art can make various amendment or supplement or adopt similar mode to substitute to described specific embodiment, but can't depart from spirit of the present invention or surmount the scope that appended claims defines.
Claims (5)
1. the remote sensing image Bowtie effect correcting method based on path tracing, raw video for having Bowtie effect to arbitrary width is corrected, generate the correction image eliminating Bowtie effect, described raw video adopts many spy units and sweeps mode imaging, be spliced up and down by several scanning strips, the line number of each scanning strip equals the number visiting unit, it is characterized in that: comprise the following steps
Step 1, to the situation of longitude and latitude data resolution subsidiary in satellite data file lower than the resolution of satellite image, carries out interpolation, makes longitude and latitude data and raw video have identical resolution;
Step 2, determines first search starting point S0 correcting pixel according to the lift rail direction of satellite and direction of scanning, comprises the center of the longitude and latitude graticule mesh corresponding to maximum for minimum for longitude in raw video, latitude a jiao as S0; The current location S=S0 of order search, the longitude and latitude graticule mesh claiming S place is current mesh, and the scanning strip at S place is Strip1;
Step 3, set and adopt the order of serpentine to carry out resampling by column to correcting image, getting according to order the pixel correcting image is current correction pixel, with the impact point T of the latitude and longitude coordinates of this pixel for search;
Step 4, from S in the scanning strip Strip1 of raw video, carries out first time path tracing to T; If follow the trail of result be S and T in same longitude and latitude graticule mesh, then find the target mesh at T place, forward 0 to, if S follows the trail of the border of arrival Strip1 longitude and latitude graticule mesh and T outside border, forward 0 to;
Step 5, according to lift rail direction and the direction of scanning of satellite, and current resampling pixel is in the parity of correcting the row in image, determines new scanning strip Strip2 T being carried out to second time path tracing, and the current location S of more new search;
Step 6, from S in the scanning strip Strip2 of raw video, carries out second time path tracing to T; If follow the trail of result be S and T in same longitude and latitude graticule mesh, then find the target mesh at T place, forward 0 to, if the border of S arrival Strip2 longitude and latitude graticule mesh and T, outside border, forward 0 to;
Step 7, utilizes 8 pixels the most contiguous with T, and after carrying out gray resample by inverse distance weighted method, assignment is to correcting image, and upgrades the center that S is the target mesh at T place, forwards 0 to;
Described 8 pixels the most contiguous with T, be T in the overlay region of raw video, the target mesh respectively searched in Strip1 and Strip2, gets summit totally 8 points of these two target mesh;
Step 8, utilizes 4 pixels the most contiguous with T, and after carrying out gray resample by inverse distance weighted method, assignment is to correcting image, and upgrades the center that S is the target mesh at T place, forwards 0 to;
Described 4 pixels the most contiguous with T, be T in the non-overlapped district of raw video, only in Strip1, search a target mesh, get summit totally 4 points of this target mesh;
Step 9, if T is outside raw video border, gives and corrects image tax invalid value, and retain the position of S in scanning strip Strip1, forward 0 to; If T is on other scanning strip, then Strip1 is updated to the scanning strip that its next one is new, redefines the position of S, return step 4;
Step 10, judge whether the correction of view picture image, if do not complete, the current location S searched at the end of retaining current correction pixel resampling is constant, return step 3 and get that to correct the next pixel of image be current correction pixel in order, with the impact point T of the latitude and longitude coordinates of this pixel for search, carry out the resampling of next pixel, if complete, end loop.
2. according to claim 1 based on the remote sensing image Bowtie effect correcting method of path tracing, it is characterized in that: when carrying out interpolation in step 1, utilize 4 adjacent points to carry out the encryption of cubic curve interpolation in X-direction, before and after utilizing in the Y direction, 2 points carry out linear interpolation encryption.
3. according to claim 1 or 2 based on the remote sensing image Bowtie effect correcting method of path tracing, it is characterized in that: in step 4 and step 6, path trace mode is, judge the line segment ST working direction of determining searching route whether crossing with four limits of current mesh, if with wherein crossing on one side, searching route enters the adjacent graticule mesh on this limit, and upgrade the center that S is adjacent graticule mesh, repeatedly carry out judgement search operation.
4. according to claim 1 or 2 based on the remote sensing image Bowtie effect correcting method of path tracing, it is characterized in that: in step 5 and step 9, path trace mode is,
New scanning strip determination mode is as follows,
1. judge the characteristic of raw video, comprising the some life in the upper left corner in raw video is No. 0 point, and the point in the upper right corner, the lower right corner, the lower left corner is followed successively by 1,2, No. 3 point in the direction of the clock; If the maximum point of raw video angle point middle latitude is No. 0 point or No. 1 point, make variable Nextscan=1; If the maximum point of raw video angle point middle latitude is No. 2 points or No. 3 points, make variable Nextscan=-1;
2. judge that current correction pixel is in the parity of correcting column in image, if at odd column, makes variable order=1, if at even column, makes variable order=-1;
3. determine that new scanning strip is Strip1+order × Nextscan;
Current location S update mode is as follows,
If search in the above-mentioned scanning strip determined for the first time, then in two kinds of situation,
If 1. order × Nextscan=1, then upgrading S is the graticule mesh center that in the first row of this scanning strip, longitude is minimum;
If 2. order × Nextscan=-1, then upgrading S is the graticule mesh center that in last column of this scanning strip, longitude is minimum;
If carried out search in the above-mentioned scanning strip determined, then upgrading S was the last position that S finally reaches when searching in this scanning strip.
5. according to claim 3 based on the remote sensing image Bowtie effect correcting method of path tracing, it is characterized in that: in step 5 and step 9, path trace mode is,
New scanning strip determination mode is as follows,
1. judge the characteristic of raw video, comprising the some life in the upper left corner in raw video is No. 0 point, and the point in the upper right corner, the lower right corner, the lower left corner is followed successively by 1,2, No. 3 point in the direction of the clock; If the maximum point of raw video angle point middle latitude is No. 0 point or No. 1 point, make variable Nextscan=1; If the maximum point of raw video angle point middle latitude is No. 2 points or No. 3 points, make variable Nextscan=-1;
2. judge that current correction pixel is in the parity of correcting column in image, if at odd column, makes variable order=1, if at even column, makes variable order=-1;
3. determine that new scanning strip is Strip1+order × Nextscan;
Current location S update mode is as follows,
If search in the above-mentioned scanning strip determined for the first time, then in two kinds of situation,
If 1. order × Nextscan=1, then upgrading S is the graticule mesh center that in the first row of this scanning strip, longitude is minimum;
If 2. order × Nextscan=-1, then upgrading S is the graticule mesh center that in last column of this scanning strip, longitude is minimum;
If carried out search in the above-mentioned scanning strip determined, then upgrading S was the last position that S finally reaches when searching in this scanning strip.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510181391.6A CN104766328B (en) | 2015-04-16 | 2015-04-16 | Remote sensing image Bowtie effect correcting methods based on path tracing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510181391.6A CN104766328B (en) | 2015-04-16 | 2015-04-16 | Remote sensing image Bowtie effect correcting methods based on path tracing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104766328A true CN104766328A (en) | 2015-07-08 |
| CN104766328B CN104766328B (en) | 2017-07-25 |
Family
ID=53648137
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510181391.6A Expired - Fee Related CN104766328B (en) | 2015-04-16 | 2015-04-16 | Remote sensing image Bowtie effect correcting methods based on path tracing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104766328B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108154471A (en) * | 2017-11-15 | 2018-06-12 | 北京海致网聚信息技术有限公司 | A kind of method and apparatus for remote sensing image splicing |
| CN110956588A (en) * | 2019-09-23 | 2020-04-03 | 武汉大学 | Image high-precision geometric correction method based on shortest distance of encrypted points |
| CN111899183A (en) * | 2019-05-06 | 2020-11-06 | 中国海洋大学 | Remote sensing image geometric fine correction method based on geographic positioning data |
| CN113469899A (en) * | 2021-06-04 | 2021-10-01 | 中国资源卫星应用中心 | Optical remote sensing satellite relative radiation correction method based on radiant energy reconstruction |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20110094902A (en) * | 2010-02-18 | 2011-08-24 | 공주대학교 산학협력단 | Normalized difference vegetation index correction method using spatio-temporal continuity |
| CN103218787A (en) * | 2013-04-23 | 2013-07-24 | 国家测绘地理信息局卫星测绘应用中心 | Multi-source heterogeneous remote-sensing image control point automatic collecting method |
| KR20130100851A (en) * | 2012-02-08 | 2013-09-12 | 가이아쓰리디 주식회사 | Method for processing satellite image and system for processing the same |
| CN103530499A (en) * | 2013-08-29 | 2014-01-22 | 西南林业大学 | Method for building mountainous area surface temperature base line and application |
| CN104021556A (en) * | 2014-06-13 | 2014-09-03 | 西南交通大学 | Heterological remote-sensing image registration method based on geometric structure similarity |
| US20150049855A1 (en) * | 2010-03-19 | 2015-02-19 | Triple Ring Technologies, Inc. | Method and apparatus for tomographic x-ray imaging and source configuration |
-
2015
- 2015-04-16 CN CN201510181391.6A patent/CN104766328B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20110094902A (en) * | 2010-02-18 | 2011-08-24 | 공주대학교 산학협력단 | Normalized difference vegetation index correction method using spatio-temporal continuity |
| US20150049855A1 (en) * | 2010-03-19 | 2015-02-19 | Triple Ring Technologies, Inc. | Method and apparatus for tomographic x-ray imaging and source configuration |
| KR20130100851A (en) * | 2012-02-08 | 2013-09-12 | 가이아쓰리디 주식회사 | Method for processing satellite image and system for processing the same |
| CN103218787A (en) * | 2013-04-23 | 2013-07-24 | 国家测绘地理信息局卫星测绘应用中心 | Multi-source heterogeneous remote-sensing image control point automatic collecting method |
| CN103530499A (en) * | 2013-08-29 | 2014-01-22 | 西南林业大学 | Method for building mountainous area surface temperature base line and application |
| CN104021556A (en) * | 2014-06-13 | 2014-09-03 | 西南交通大学 | Heterological remote-sensing image registration method based on geometric structure similarity |
Non-Patent Citations (2)
| Title |
|---|
| 刘良明 等: "MODIS数据Bowtie效应快速消除算法研究", 《国土资源遥感》 * |
| 王汉禹 等: "MODIS数据Bowtie效应消除算法", 《计算机工程》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108154471A (en) * | 2017-11-15 | 2018-06-12 | 北京海致网聚信息技术有限公司 | A kind of method and apparatus for remote sensing image splicing |
| CN111899183A (en) * | 2019-05-06 | 2020-11-06 | 中国海洋大学 | Remote sensing image geometric fine correction method based on geographic positioning data |
| CN110956588A (en) * | 2019-09-23 | 2020-04-03 | 武汉大学 | Image high-precision geometric correction method based on shortest distance of encrypted points |
| CN110956588B (en) * | 2019-09-23 | 2022-08-05 | 武汉大学 | Image high-precision geometric correction method based on shortest distance of encrypted points |
| CN113469899A (en) * | 2021-06-04 | 2021-10-01 | 中国资源卫星应用中心 | Optical remote sensing satellite relative radiation correction method based on radiant energy reconstruction |
| CN113469899B (en) * | 2021-06-04 | 2023-12-29 | 中国资源卫星应用中心 | Optical remote sensing satellite relative radiation correction method based on radiation energy reconstruction |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104766328B (en) | 2017-07-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107504981B (en) | Satellite attitude error correction method and device based on laser height measurement data | |
| US8315477B2 (en) | Method and apparatus of taking aerial surveys | |
| JP3869814B2 (en) | Orthorectification processing method for satellite images | |
| CN104766328A (en) | Remote-sensing image Bowtie effect correction method based on path tracing | |
| JP2013080467A (en) | Enhancing video using super-resolution | |
| CN102735216B (en) | CCD stereoscopic camera three-line imagery data adjustment processing method | |
| CN102968788A (en) | Wave band registering method based on regular grid surface element | |
| CN106709944B (en) | Satellite remote sensing image registration method | |
| CN103822615A (en) | Unmanned aerial vehicle ground target real-time positioning method with automatic extraction and gathering of multiple control points | |
| JP2009223220A (en) | Road surface marking map creating method | |
| CN101609149A (en) | A kind of method that improves attitude determination precision of airborne laser radar | |
| CN110246082A (en) | A kind of remote sensing Panorama Mosaic method | |
| CN113706702A (en) | Mining area three-dimensional map construction system and method | |
| CN111709876B (en) | Image splicing method, device, equipment and storage medium | |
| CN102147249B (en) | Method for precisely correcting satellite-borne optical linear array image based on linear characteristic | |
| CN103398701B (en) | Satellite-borne non-colinear TDI (time delay integral) CCD (charge coupled device) image splicing method based on object space projection plane | |
| CN103778610A (en) | Geometric pretreatment method for vertical rail swing images of satellite-borne linear array sensor | |
| US11361502B2 (en) | Methods and systems for obtaining aerial imagery for use in geospatial surveying | |
| CN105335944A (en) | Onboard staggered TDI (Transport Driver Interface) infrared image recovering method and system based on geometric correction | |
| Li | Coastline mapping and change detection using one-meter resolution satellite imagery | |
| CN112863183B (en) | Traffic flow data fusion method and system | |
| Kombiadou et al. | EVREST Project Report: Remote Sensing Database Report | |
| CN116681611A (en) | Linear array swing scanning type optical satellite image geometric correction method and system | |
| CN117170401A (en) | Aerial photography ground control point position planning method and related equipment | |
| CN105279341A (en) | Method for intercepting google earth image and directly applying in engineering design |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate 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: 20170725 Termination date: 20190416 |