CN102682441A - Hyperspectral image super-resolution reconstruction method based on subpixel mapping - Google Patents
Hyperspectral image super-resolution reconstruction method based on subpixel mapping Download PDFInfo
- Publication number
- CN102682441A CN102682441A CN2012100513658A CN201210051365A CN102682441A CN 102682441 A CN102682441 A CN 102682441A CN 2012100513658 A CN2012100513658 A CN 2012100513658A CN 201210051365 A CN201210051365 A CN 201210051365A CN 102682441 A CN102682441 A CN 102682441A
- Authority
- CN
- China
- Prior art keywords
- pixel
- atural object
- image
- resolution
- inferior
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
The invention discloses a hyperspectral image super to resolution reconstruction method based on the subpixel mapping. The method comprises the following steps of: carrying out the subpixel mapping to acquire proportional values of various ground objects in each pixel element in the original hyperspectral remote sensing image; dividing each pixel element, namely a rough pixel element, in the original hyperspectral remote sensing image into nPoundn subpixels, wherein n is a predefined resolution amplifying factor; assigning each subpixel element to some ground object at random according to the abundance value of each known ground object in each pixel element, so as to obtain the initialized ground object distribution image; then obtaining a ground object distribution image of which resolution is n times of that of the original image by using the simulated annealing algorithm according to the iterative optimization rule; and finally, performing the super to resolution reconstruction based on the subpixel mapping. By considering the partial and overall distribution properties of the ground object at the same time, the hyperspectral image super-resolution reconstruction method based on the subpixel mapping avoids the inauthenticity caused by solving the partial property for the ground object distribution, and thus the hyperspectral remote sensing image can be obtained from the ground object distribution image.
Description
Technical field
The invention belongs to technical field of image processing, be applicable to the high-spectrum remote sensing reconstruction, be specifically related to a kind of high spectrum image super-resolution method for reconstructing based on inferior pixel mapping.
Background technology
In recent years; Satellite remote sensing technology has obtained development at full speed; The new high-spectrum remote-sensing imaging technique that rises can be a lot, obtain face of land observation data on the very narrow and quasi-continuous spectral band; For people provide abundant more atural object observation information, strengthened directly identification and the ability of analyzing the atural object situation from remote sensing images greatly.But high light spectrum image-forming technology spatial resolution when obtaining very high spectral resolution is restricted, and the spatial resolution that therefore how to promote high-spectrum remote sensing becomes problem demanding prompt solution.
For with treat the imaging space plane at a distance of the space remote sensing platform of R, ignoring under the condition that optical system nonlinear distortion and noise influence, its spatial resolution Δ L can be expressed as
Wherein, W is a CCD device array element width, and f is the focal length of optical system.So the spatial resolution that wants to improve remote sensing images can be passed through three kinds of modes, a kind of is to reduce remote-sensing flatform flight track height, but can make the lost of life of this platform like this; Keeping under the constant situation of R, other two kinds of methods that improve resolution are to increase focal length and dwindle array element width, yet; Increasing focal length can make the difficulty of processing of optical element increase; Expense increases, and causes the big Heavy Weight of volume of remote sensor, for the application of reality brings difficulty; The size of CCD array element receives the restriction of technology, because external blockade on new techniques, it is enough little that the CCD aperture that we can access can not be done.To sum up, how can be in distance R, making the remote sensing images spatial resolution get a promotion under the big or small constant prerequisite of focal distance f and CCD aperture W is problem demanding prompt solution.The super-resolution image reconstruction technology is the effective way that addresses this problem, and these techniques make use one or more from different perspectives, the low-resolution image that obtains such as diverse location, different sensors reconstructs the one or more high-definition picture.
The high light spectrum image-forming mode receive the sensor spatial resolution restriction and atural object complex distribution property influence and make that the spatial resolution of gained image is not high; This just causes single pixel institute numerical values recorded information to be mostly the amalgam that multiple type of ground objects spectral signature is constituted, and is called as mixed pixel.The existence of mixed pixel phenomenon has brought great difficulty to the quantitative decipher of image, also brings difficulty for the spatial resolution that improves high spectrum image.Mixed pixel decomposition method is to solve one of this difficult effective way; Its basic meaning is for to extract typical case " pure " object spectrum that is called as end member through mixed pixel spectrum; Estimate each end member shared ratio in mixed pixel simultaneously, this ratio is called the abundance value.For every type of atural object, its abundance value in each pixel forms a width of cloth abundance figure.But this type mixed pixel decomposition method only provides the ratio of every kind of atural object in each pixel, does not provide them and in this pixel, how to distribute.The abundance figure that inferior pixel mapping method utilizes mixed pixel decomposition method to obtain is estimating the distribution of every kind of atural object in each pixel, thereby can access the terrain classification image higher than original image resolution under some restrictive condition.Existing inferior pixel mapping method mainly is divided into two types; Based on study and based on space continuity; The former collects some known high resolving power atural object distributed images as training set; Through being reduced resolution, they obtain abundance figure; Adopt the means training of machine learning to obtain high resolving power atural object distributed image and the relation between the abundance figure accordingly then, and with this relational application on pending high spectrum image abundance figure, thereby obtain needed high resolving power atural object distributed image; Mostly the latter is that based on the algorithm of space continuity distribution has this characteristics proposition different model of local continuity to atural object, distributes through optimizing the optimum atural object that is met this model then.
Summary of the invention
In order to overcome the deficiency of above-mentioned prior art, the object of the present invention is to provide a kind of high spectrum image super-resolution method for reconstructing based on inferior pixel mapping, can access the high resolving power high-spectrum remote sensing higher than original image resolution.
To achieve these goals, the technical scheme of the present invention's employing is:
High spectrum image super-resolution method for reconstructing based on inferior pixel mapping may further comprise the steps:
Step 1, inferior pixel mapping;
Step 1.1 adopts mixed pixel decomposition method to obtain in the original high-spectrum remote sensing ratio value of various atural objects in each pixel, is called the abundance value;
Step 1.2 is called thick pixel with each pixel in the original high-spectrum remote sensing, and each thick pixel is divided into n £ n inferior pixel, and wherein n is predefined resolution enlargement factor;
Step 1.3 according to the abundance value of every kind of known atural object in each pixel, is a certain atural object with each inferior pixel Random assignment, promptly obtains initialized atural object distributed image;
Step 1.4 according to the iteration optimization criterion, utilizes simulated annealing that initialized atural object distributed image is carried out iteration optimization; Until satisfying stopping criterion for iteration; Described iteration optimization criterion is based on the local continuity of atural object distribution and the overall similarity of atural object distribution, wherein
The local continuity that atural object distributes, the atural object that is meant same classification in subrange than other atural object of variety classes distribute near;
The overall similarity that atural object distributes is meant that the atural object distribution in the thick pixel that in global scope, has similar proportion atural object composition has similarity, and arbitrary thus pixel can be by representing with the similar pixel of its atural object distribution;
Following objective function is as described iteration optimization criterion:
Wherein, t is used for belonging in pending pixel of index and the neighborhood thereof the inferior pixel of atural object i, and c refers to atural object classification number,
Finger belongs to the horizontal ordinate and the ordinate of the inferior pixel of atural object classification i, μ
Ix, μ
IyRefer to respectively
Mean value, η is the regularization factor, q
jBe the pixel that has similarity with pixel p, it is selected according to being
The proper vector that f (p) is made up of the abundance value of pixel p and neighborhood thereof, th is a preset threshold value, if a certain pixel q
jProper vector belong to preceding J (j≤J) individual with the similar vector of f (p), this pixel q so
jThe similar pixel that just is considered to pixel p, therefore the inferior pixel distribution in pixel p can be by pixel q
jIn the linear combination that distributes of inferior pixel represent, promptly
α
jBe er
jNormalization reciprocal.The local continuity that on behalf of atural object, this objective function first
distribute, second
represents the overall similarity of atural object distribution;
Obtaining a width of cloth resolution thus is original image n atural object distributed image doubly;
Step 2 is based on the super-resolution method for reconstructing of inferior pixel mapping;
It is m £ m that the atural object distributed image that step 1 is obtained is divided into a series of sizes; The subimage block of 1m<n; Recomputate the ratio that wherein comprises various atural objects for each subimage block; Multiply by the pure end member object spectrum that obtains in the mixed pixel decomposition method with this ratio value then; And linear, additive obtains the curve of spectrum of this subimage block, thereby obtains the n=m high-spectrum remote sensing doubly that resolution is original image.
Mixed pixel decomposition method in the said step 1.1 is meant the ratio of the every kind of typical feature that obtains forming mixed pixel, and wherein mixed pixel is meant the pixel that includes dissimilar atural objects in the high spectrum image of acquisition.
The iteration optimization process is the distribution that changes atural object through continuous iteration in the said step 1.4; Wherein each iterative process all can make objective function
value reduce, and said end condition is that target function value no longer descends or reached the iterations of being scheduled to.
Because the present invention will improve the resolution of high spectrum image, so can claim that this image of waiting to improve resolution is a low-resolution image.
Compared with prior art, advantage of the present invention is:
1) the present invention sets up in the local continuity that atural object is distributed and has adopted the minimum criterion of dispersion in the class that makes between identical atural object on the model;
2) the present invention has considered the overall similarity that atural object distributes when considering the local continuity that atural object distributes, and has avoided the untrue property of bringing when only finding the solution to the local characteristics of atural object distribution;
3) the present invention proposes the super resolution ratio reconstruction method that shines upon based on inferior pixel, can obtain high-spectrum remote sensing from the atural object distributed image.
Description of drawings
Fig. 1 is the dividing mode synoptic diagram of thick pixel among the present invention.
Fig. 2 is the synoptic diagram as a result that obtains of inferior pixel mapping method.
Embodiment
Below in conjunction with embodiment the present invention is explained further details.
High spectrum image super-resolution method for reconstructing based on inferior pixel mapping comprises the steps:
Step 1, inferior pixel mapping;
Step 1.1 adopts mixed pixel decomposition method to obtain in the original high-spectrum remote sensing ratio value of various atural objects in each pixel, is called the abundance value;
Mixed pixel decomposition method is a kind of treatment technology commonly used in this area, provides one here than mixed pixel decomposition method commonly used.This method is at first chosen the pure atural object of suitable composition high spectrum image, is referred to as end member, secondly extracts the spectrum of end member from image or in the wave spectrum storehouse, adopts linear mixed model to be used as the spectrum mixture model then, i.e. the r=∑
ia
is
i+ w 0≤a
i≤1, ∑
ia
i=1, wherein, r is the spectrum of each pixel in the high spectrum image, s
iBe the spectrum of each end member, a
iBe the ratio value of each end member in this pixel, i.e. abundance value, 0≤a
i≤1 is meant that the ratio value of every type of atural object should be between 0 and 1, ∑
ia
i=1 is meant that the ratio value sum of all atural objects in this pixel should be 1, and w is an error.Thus, according to making the minimum principle of error w, the abundance value of each atural object in each pixel can obtain through finding the solution aforesaid equation.Present embodiment is that example specifies embodiment of the present invention with a simple analog image, chooses two types of atural objects in this example, i.e. target " 0 " and background " 1 ".Pending image size is 9 pixels * 9 pixels, is not that real high-spectrum remote sensing and mixed pixel decomposition method do not belong to content of the present invention owing to what adopt, so set the abundance value of two types of atural objects in each pixel here shown in table 1 and table 2.
Table 1, the abundance value of atural object " 0 " in each pixel
? | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
2 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
3 | 0 | 0 | 0 | 0.0625 | 0.2500 | 0.0625 | 0 | 0 | 0 |
[0039]?
4 | 0 | 0 | 0.1250 | 0.9375 | 1 | 0.9375 | 0.1250 | 0 | 0 |
5 | 0 | 0 | 0.5000 | 1 | 1 | 1 | 0.5000 | 0 | 0 |
6 | 0 | 0 | 0.3750 | 1 | 1 | 1 | 0.3750 | 0 | 0 |
7 | 0 | 0 | 0 | 0.5000 | 0.7500 | 0.5000 | 0 | 0 | 0 |
8 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
9 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
Table 2, the abundance value of atural object " 1 " in each pixel
? | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
2 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
3 | 1 | 1 | 1 | 0.9375 | 0.7500 | 0.9375 | 1 | 1 | 1 |
4 | 1 | 1 | 0.8750 | 0.0625 | 0 | 0.0625 | 0.8750 | 1 | 1 |
5 | 1 | 1 | 0.5000 | 0 | 0 | 0 | 0.5000 | 1 | 1 |
6 | 1 | 1 | 0.6250 | 0 | 0 | 0 | 0.6250 | 1 | 1 |
7 | 1 | 1 | 1 | 0.5000 | 0.2500 | 0.5000 | 1 | 1 | 1 |
8 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
9 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
Step 1.2 is called thick pixel with each pixel in the original high-spectrum remote sensing, and each thick pixel is divided into n £ n inferior pixel, and wherein n is predefined resolution enlargement factor;
Among the present invention each pixel in the pending high spectrum image is called thick pixel, promptly 9 * 9 pixels all are called thick pixel in the table 1.Setting the resolution enlargement factor in the present embodiment is 4, and promptly n=4 is equally divided into 4 * 4 inferior pixels with each thick pixel then, and is as shown in Figure 1.
Step 1.3 according to the abundance value of every kind of known atural object in each pixel, is a certain atural object with each inferior pixel Random assignment, promptly obtains initialized atural object distributed image;
With the abundance value of the 5th row target " 0 " of the third line in the table 1 is that 0.25 pixel is an example.After thick pixel is divided into 4 * 4 inferior pixels with each, should there be 4 * 4 * 0.25=4 inferior pixel to belong to target " 0 ", in like manner, there are 12 inferior pixels to belong to background " 1 ".Then in these 4 * 4 inferior pixels at random 4 of pickings be assigned as target " 0 ", remaining 12 inferior pixel is assigned as background " 1 ".
Step 1.4 according to the iteration optimization criterion, utilizes simulated annealing that initialized atural object distributed image is carried out iteration optimization; Until satisfying stopping criterion for iteration; Described iteration optimization criterion is based on the local continuity of atural object distribution and the overall similarity of atural object distribution, wherein
The local continuity that atural object distributes, the atural object that is meant same classification in subrange than other atural object of variety classes distribute near;
The overall similarity that atural object distributes is meant that the atural object distribution in the thick pixel that in global scope, has similar proportion atural object composition has similarity, and arbitrary thus pixel can be by representing with the similar pixel of its atural object distribution;
Following objective function is as described iteration optimization criterion:
Wherein, t is used for belonging in pending pixel of index and the neighborhood thereof the inferior pixel of atural object i, and c refers to atural object classification number,
Finger belongs to the horizontal ordinate and the ordinate of the inferior pixel of atural object classification i, μ
Ix, μ
IyRefer to respectively
Mean value, η is the regularization factor, q
jBe the pixel that has similarity with pixel p, it is selected according to being
The proper vector that f (p) is made up of the abundance value of pixel p and neighborhood thereof, th is a preset threshold value, if a certain pixel q
jProper vector belong to preceding J (j≤J) individual with the similar vector of f (p), this pixel q so
jThe similar pixel that just is considered to pixel p, therefore the inferior pixel distribution in pixel p can be by pixel q
jIn the linear combination that distributes of inferior pixel represent, promptly
α
jBe er
jNormalization reciprocal.The local continuity that on behalf of atural object, this objective function first
distribute, second
represents the overall similarity of atural object distribution;
Simulated annealing is general optimized Algorithm, is used to seek the optimum solution of objective function.In the present embodiment; One by one thick pixel is optimized calculating; Promptly at first calculate the value of objective function under the initial situation; Change the atural object classification ownership of this thick pixel Nei Geya pixel then according to simulated annealing, so constantly repeat this iterative step, no longer descend or reached predetermined iterations up to target function value.
Obtain the atural object distributed image that a width of cloth resolution is 4 times of original images thus, as shown in Figure 2, wherein target " 0 " is represented with black, and background 1 usefulness white is represented.
Step 2,
In the present embodiment; Per 2 * 2 pixels in the high-resolution atural object distributed image that step 1 is obtained are looked as a whole, i.e. m=2, thus formed a series of subimage block; The ratio of two types of atural objects that can wherein be comprised again through simple calculating each subimage block; To multiply by the pure atural object that obtains in the mixed pixel decomposition method be the spectrum of end member with this ratio value then, and linear, additive obtains the curve of spectrum of this subimage block, promptly
B wherein
iBe ratio value, s
iBe the end member spectrum that obtains in the mixed pixel decomposition method.Because present embodiment is example with the analog image, so do not provide the high-resolution remote sensing images that obtain at last.
Claims (3)
1. based on the high spectrum image super-resolution method for reconstructing of inferior pixel mapping, it is characterized in that, may further comprise the steps:
Step 1, inferior pixel mapping;
Step 1.1 adopts mixed pixel decomposition method to obtain in the original high-spectrum remote sensing ratio value of various atural objects in each pixel, is called the abundance value;
Step 1.2 is called thick pixel with each pixel in the original high-spectrum remote sensing, and each thick pixel is divided into n £ n inferior pixel, and wherein n is predefined resolution enlargement factor;
Step 1.3 according to the abundance value of every kind of known atural object in each pixel, is a certain atural object with each inferior pixel Random assignment, promptly obtains initialized atural object distributed image;
Step 1.4 according to the iteration optimization criterion, utilizes simulated annealing that initialized atural object distributed image is carried out iteration optimization; Until satisfying stopping criterion for iteration; Described iteration optimization criterion is based on the local continuity of atural object distribution and the overall similarity of atural object distribution, wherein
The local continuity that atural object distributes, the atural object that is meant same classification in subrange than other atural object of variety classes distribute near;
The overall similarity that atural object distributes is meant that the atural object distribution in the thick pixel that in global scope, has similar proportion atural object composition has similarity, and arbitrary thus pixel can be by representing with the similar pixel of its atural object distribution;
Following objective function is as described iteration optimization criterion:
Wherein, t is used for belonging in pending pixel of index and the neighborhood thereof the inferior pixel of atural object i, and c refers to atural object classification number,
Finger belongs to the horizontal ordinate and the ordinate of the inferior pixel of atural object classification i, μ
Ix, μ
IyRefer to respectively
Mean value, η is the regularization factor, q
jBe the pixel that has similarity with pixel p, it is selected according to being
The proper vector that f (p) is made up of the abundance value of pixel p and neighborhood thereof, th is a preset threshold value, if a certain pixel q
jProper vector belong to preceding J (j≤J) individual with the similar vector of f (p), this pixel q so
jThe similar pixel that just is considered to pixel p, therefore the inferior pixel distribution in pixel p can be by pixel q
jIn the linear combination that distributes of inferior pixel represent, promptly
α
jBe er
jNormalization reciprocal.The local continuity that on behalf of atural object, this objective function first
distribute, second
represents the overall similarity of atural object distribution;
Obtaining a width of cloth resolution thus is original image n atural object distributed image doubly;
Step 2 is based on the super-resolution method for reconstructing of inferior pixel mapping;
It is m £ m that the atural object distributed image that step 1 is obtained is divided into a series of sizes; The subimage block of 1m<n; Recomputate the ratio that wherein comprises various atural objects for each subimage block; Multiply by the pure end member object spectrum that obtains in the mixed pixel decomposition method with this ratio value then; And linear, additive obtains the curve of spectrum of this subimage block, thereby obtains the n=m high-spectrum remote sensing doubly that resolution is original image.
2. according to the said high spectrum image super-resolution method for reconstructing of claim 1 based on inferior pixel mapping; It is characterized in that; Mixed pixel decomposition method in the said step 1.1 is meant the ratio of the every kind of typical feature that obtains forming mixed pixel, and wherein mixed pixel is meant the pixel that includes dissimilar atural objects in the high spectrum image of acquisition.
3. according to the said high spectrum image super-resolution method for reconstructing of claim 1 based on inferior pixel mapping; It is characterized in that; The iteration optimization process is the distribution that changes atural object through continuous iteration in the said step 1.4; Wherein each iterative process all can make objective function
value reduce, and said end condition is that target function value no longer descends or reached the iterations of being scheduled to.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210051365.8A CN102682441B (en) | 2012-03-01 | 2012-03-01 | Hyperspectral image super-resolution reconstruction method based on subpixel mapping |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210051365.8A CN102682441B (en) | 2012-03-01 | 2012-03-01 | Hyperspectral image super-resolution reconstruction method based on subpixel mapping |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102682441A true CN102682441A (en) | 2012-09-19 |
CN102682441B CN102682441B (en) | 2014-11-26 |
Family
ID=46814308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210051365.8A Expired - Fee Related CN102682441B (en) | 2012-03-01 | 2012-03-01 | Hyperspectral image super-resolution reconstruction method based on subpixel mapping |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102682441B (en) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103559506A (en) * | 2013-11-19 | 2014-02-05 | 中国科学院地理科学与资源研究所 | Sub-pixel drawing method based on vector boundaries |
CN104182953A (en) * | 2013-12-30 | 2014-12-03 | 北京大学 | Normalized multi-end-member decomposition hyperspectral reconstruction method in pixel unmixing inverse process |
CN106157240A (en) * | 2015-04-22 | 2016-11-23 | 南京理工大学 | Remote sensing image super resolution method based on dictionary learning |
CN106548450A (en) * | 2016-09-30 | 2017-03-29 | 西安电子科技大学 | A kind of high spectrum image ultra-resolution method, device and computer-readable medium |
CN106558036A (en) * | 2016-10-27 | 2017-04-05 | 中国科学院光电研究院 | A kind of spaceborne super-resolution imaging design method |
CN107392208A (en) * | 2017-05-23 | 2017-11-24 | 三亚中科遥感研究所 | Object Spectra feature extracting method with purifying is mapped based on spectral space |
CN108921035A (en) * | 2018-06-06 | 2018-11-30 | 中国农业科学院农业资源与农业区划研究所 | Sub-pixed mapping localization method and system based on spatial attraction and pixel concentration class |
CN108985154A (en) * | 2018-06-06 | 2018-12-11 | 中国农业科学院农业资源与农业区划研究所 | Small size atural object sub-pixed mapping localization method and system based on image concentration class |
CN110211042A (en) * | 2019-05-10 | 2019-09-06 | 北京航空航天大学 | The sub-pixed mapping localization method and device of enhanced spectrum image spatial resolution |
CN111899300A (en) * | 2020-07-30 | 2020-11-06 | 北京航空航天大学 | Abundance correction method and device for light field spectral data subpixel positioning |
CN113034637A (en) * | 2021-03-11 | 2021-06-25 | 郑州轻工业大学 | Multi-scale rapid simulated annealing modeling method based on two-dimensional core structure |
CN113409193A (en) * | 2021-06-18 | 2021-09-17 | 北京印刷学院 | Super-resolution reconstruction method and device for hyperspectral image |
CN113624691A (en) * | 2020-05-07 | 2021-11-09 | 南京航空航天大学 | Spectral image super-resolution mapping method based on space-spectrum correlation |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1874499A (en) * | 2006-05-12 | 2006-12-06 | 北京理工大学 | High dynamic equipment for reconstructing image in high resolution |
EP1906357A1 (en) * | 2006-09-28 | 2008-04-02 | Matsushita Electric Industrial Co., Ltd. | Combined optimization of motion estimator and super-resolution |
CN102096920A (en) * | 2010-12-31 | 2011-06-15 | 清华大学 | Target image-based sub-pixel registering method |
CN102136144A (en) * | 2011-04-11 | 2011-07-27 | 北京大学 | Image registration reliability model and reconstruction method of super-resolution image |
-
2012
- 2012-03-01 CN CN201210051365.8A patent/CN102682441B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1874499A (en) * | 2006-05-12 | 2006-12-06 | 北京理工大学 | High dynamic equipment for reconstructing image in high resolution |
EP1906357A1 (en) * | 2006-09-28 | 2008-04-02 | Matsushita Electric Industrial Co., Ltd. | Combined optimization of motion estimator and super-resolution |
CN102096920A (en) * | 2010-12-31 | 2011-06-15 | 清华大学 | Target image-based sub-pixel registering method |
CN102136144A (en) * | 2011-04-11 | 2011-07-27 | 北京大学 | Image registration reliability model and reconstruction method of super-resolution image |
Non-Patent Citations (1)
Title |
---|
ZONGXU PAN: "Super resolution of remote sensing image based on structure similarity in CS frame", 《MULTISPECTRAL IMAGE ACQUISTION PROCESSING AND ANALYSIS 2011》, 8 December 2011 (2011-12-08), pages 1 - 8 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103559506A (en) * | 2013-11-19 | 2014-02-05 | 中国科学院地理科学与资源研究所 | Sub-pixel drawing method based on vector boundaries |
CN103559506B (en) * | 2013-11-19 | 2015-04-15 | 中国科学院地理科学与资源研究所 | Sub-pixel drawing method based on vector boundaries |
CN104182953A (en) * | 2013-12-30 | 2014-12-03 | 北京大学 | Normalized multi-end-member decomposition hyperspectral reconstruction method in pixel unmixing inverse process |
CN104182953B (en) * | 2013-12-30 | 2017-02-01 | 北京大学 | Normalized multi-end-member decomposition hyperspectral reconstruction method in pixel unmixing inverse process |
CN106157240A (en) * | 2015-04-22 | 2016-11-23 | 南京理工大学 | Remote sensing image super resolution method based on dictionary learning |
CN106548450A (en) * | 2016-09-30 | 2017-03-29 | 西安电子科技大学 | A kind of high spectrum image ultra-resolution method, device and computer-readable medium |
CN106558036B (en) * | 2016-10-27 | 2019-08-02 | 中国科学院光电研究院 | A kind of spaceborne super-resolution imaging design method |
CN106558036A (en) * | 2016-10-27 | 2017-04-05 | 中国科学院光电研究院 | A kind of spaceborne super-resolution imaging design method |
CN107392208A (en) * | 2017-05-23 | 2017-11-24 | 三亚中科遥感研究所 | Object Spectra feature extracting method with purifying is mapped based on spectral space |
CN107392208B (en) * | 2017-05-23 | 2020-05-22 | 三亚中科遥感研究所 | Object spectral feature extraction method based on spectral space mapping and purification |
CN108921035B (en) * | 2018-06-06 | 2020-08-11 | 中国农业科学院农业资源与农业区划研究所 | Sub-pixel positioning method and system based on spatial gravitation and pixel aggregation |
CN108985154A (en) * | 2018-06-06 | 2018-12-11 | 中国农业科学院农业资源与农业区划研究所 | Small size atural object sub-pixed mapping localization method and system based on image concentration class |
CN108921035A (en) * | 2018-06-06 | 2018-11-30 | 中国农业科学院农业资源与农业区划研究所 | Sub-pixed mapping localization method and system based on spatial attraction and pixel concentration class |
CN108985154B (en) * | 2018-06-06 | 2020-10-27 | 中国农业科学院农业资源与农业区划研究所 | Small-size ground object sub-pixel positioning method and system based on image concentration |
CN110211042A (en) * | 2019-05-10 | 2019-09-06 | 北京航空航天大学 | The sub-pixed mapping localization method and device of enhanced spectrum image spatial resolution |
CN113624691A (en) * | 2020-05-07 | 2021-11-09 | 南京航空航天大学 | Spectral image super-resolution mapping method based on space-spectrum correlation |
CN113624691B (en) * | 2020-05-07 | 2022-10-04 | 南京航空航天大学 | Spectral image super-resolution mapping method based on space-spectrum correlation |
CN111899300A (en) * | 2020-07-30 | 2020-11-06 | 北京航空航天大学 | Abundance correction method and device for light field spectral data subpixel positioning |
CN111899300B (en) * | 2020-07-30 | 2022-05-31 | 北京航空航天大学 | Abundance correction method and device for light field spectral data subpixel positioning |
CN113034637A (en) * | 2021-03-11 | 2021-06-25 | 郑州轻工业大学 | Multi-scale rapid simulated annealing modeling method based on two-dimensional core structure |
CN113409193A (en) * | 2021-06-18 | 2021-09-17 | 北京印刷学院 | Super-resolution reconstruction method and device for hyperspectral image |
CN113409193B (en) * | 2021-06-18 | 2023-07-04 | 北京印刷学院 | Super-resolution reconstruction method and device for hyperspectral image |
Also Published As
Publication number | Publication date |
---|---|
CN102682441B (en) | 2014-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102682441B (en) | Hyperspectral image super-resolution reconstruction method based on subpixel mapping | |
CN105261000B (en) | A kind of hyperspectral image fusion method mixed based on Endmember extraction and spectrum solution | |
US9317929B2 (en) | Decomposition apparatus and method for refining composition of mixed pixels in remote sensing images | |
WO2024027095A1 (en) | Hyperspectral imaging method and system based on double rgb image fusion, and medium | |
CN103761742B (en) | A kind of high-spectrum remote sensing sparse solution mixing method based on homogeneity index | |
CN105976310A (en) | Partition-based VCA end-member extraction method | |
CN102629374B (en) | Image super resolution (SR) reconstruction method based on subspace projection and neighborhood embedding | |
CN104463223B (en) | Hyperspectral image group sparse unmixing method based on space spectrum information abundance constraint | |
CN111598786B (en) | Hyperspectral image unmixing method based on depth denoising self-coding network | |
CN114694039B (en) | Remote sensing hyperspectral and laser radar image fusion classification method and device | |
CN112712034B (en) | Unmixing method and system for hyperspectral image | |
CN104751181A (en) | High spectral image Deming method based on relative abundance | |
CN102073867B (en) | Sorting method and device for remote sensing images | |
CN102789639A (en) | Method for fusing high-spectrum image and visible light image based on nonnegative matrix decomposition | |
CN106778530B (en) | Hyperspectral image nonlinear unmixing method based on bilinear hybrid model | |
CN103425995A (en) | Hyperspectral image classification method based on area similarity low rank expression dimension reduction | |
CN107590800A (en) | A kind of hyperspectral remote sensing image mixed pixel decomposition method based on object spectrum storehouse | |
CN104200513A (en) | Matrix row-column sampling based multi-light-source rendering method | |
CN103226825A (en) | Low-rank sparse model-based remote sensing image change detection method | |
CN114021603A (en) | Radar signal modulation mode identification method based on model reparameterization | |
CN114612315A (en) | High-resolution image missing region reconstruction method based on multi-task learning | |
CN104318046A (en) | System and method for incrementally converting high dimensional data into low dimensional data | |
CN113421258A (en) | Automatic cloud detection method based on high-resolution remote sensing image | |
CN106971402B (en) | SAR image change detection method based on optical assistance | |
CN102831621A (en) | Video significance processing method based on spectral analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141126 Termination date: 20180301 |
|
CF01 | Termination of patent right due to non-payment of annual fee |