CN114511473A - Hyperspectral remote sensing image denoising method based on unsupervised adaptive learning - Google Patents

Hyperspectral remote sensing image denoising method based on unsupervised adaptive learning Download PDF

Info

Publication number
CN114511473A
CN114511473A CN202210409199.8A CN202210409199A CN114511473A CN 114511473 A CN114511473 A CN 114511473A CN 202210409199 A CN202210409199 A CN 202210409199A CN 114511473 A CN114511473 A CN 114511473A
Authority
CN
China
Prior art keywords
denoising
image
network
discriminator
loss
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
Application number
CN202210409199.8A
Other languages
Chinese (zh)
Other versions
CN114511473B (en
Inventor
王心宇
罗朝之
钟燕飞
张良培
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202210409199.8A priority Critical patent/CN114511473B/en
Publication of CN114511473A publication Critical patent/CN114511473A/en
Application granted granted Critical
Publication of CN114511473B publication Critical patent/CN114511473B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A40/00Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
    • Y02A40/10Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Software Systems (AREA)
  • Molecular Biology (AREA)
  • Computational Linguistics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Mathematical Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

The invention relates to a hyperspectral remote sensing image denoising method based on unsupervised adaptive learning. Aiming at the problem that the generalization of the model is reduced by the degradation difference between the simulated image and the real image, the unsupervised self-adaptive learning strategy is provided, the pre-training is carried out on the high-quality ground image, the discriminator is designed to carry out modeling on noise, the discriminator is used for carrying out fine adjustment on the denoising parameter when the real image is processed, and the generalization of the model on the real image is improved. According to the invention, a depth denoising network based on space-spectrum residual errors and a discriminator based on global information modeling are designed in a model so as to fully excavate hyperspectral depth prior. The method can solve the problem of simulation training data in the hyperspectral remote sensing image deep learning denoising, reduces the dependence of a deep learning model on the simulation training data, and effectively improves the applicability and precision of hyperspectral denoising.

Description

Hyperspectral remote sensing image denoising method based on unsupervised adaptive learning
Technical Field
The invention relates to the technical field of remote sensing image processing, in particular to a hyperspectral remote sensing image denoising method based on unsupervised adaptive learning.
Background
With the increase of the spectral resolution of the sensor, the hyperspectral remote sensing is rapidly developed. The hyperspectral satellite can quickly acquire images of hundreds of wave bands in the same area, and has the advantage of integrating maps. However, in the process of acquiring data by the hyperspectral remote sensing satellite, due to the bad atmospheric conditions, the imperfect correction process and the defects of the sensor, the hyperspectral remote sensing image is inevitably polluted by various noises, so that the degradation problems of stripe noise, random noise, image contrast reduction and the like are generated, and the usability of the image data is severely limited. In the image data containing mixed noise, the image information consists of original radiation information reflecting ground real characteristics and noise information, and the noise seriously reduces the image quality and precision. Therefore, it is vital to remove mixed noise in the hyperspectral remote sensing satellite image and truly restore image radiation information, and the removal of the mixed noise can greatly improve the usability of the hyperspectral remote sensing satellite image data and is beneficial to the subsequent processing of the image, such as resource exploration, change detection, land utilization type investigation and the like.
At present, the existing hyperspectral mixed noise removal method has great limitation in practical application. The filtering-based method is limited by a fixed transform domain, the geometric characteristics of hyperspectral data are difficult to fully mine, and a proper cut-off frequency is difficult to determine in practical application; the regularization model-based method is based on artificial assumptions on noise characteristics and potential noise-free image characteristics, accurate modeling of complex noise in massive hyperspectral data is difficult, and when more constraint terms are added, the solving process is very complex, the operation efficiency is low, and the practical application requirements are difficult to meet; the deep learning-based method has high-efficiency operation efficiency, but the training process of the deep learning-based method depends on a large number of noise-free image pairs, because airborne and spaceborne hyperspectral remote sensing platforms can not obtain real image pairs generally, deep learning models are mostly trained on simulated image pairs, and the generalization of the models on real data is reduced. Therefore, the deep learning model capable of improving the generalization of the model and efficiently completing the task of removing the mixed noise of the hyperspectral remote sensing image in a robust mode is researched, the quality of the existing satellite data is improved, and beneficial guidance is provided for efficient preprocessing of the hyperspectral remote sensing satellite.
Disclosure of Invention
The invention aims to provide a hyperspectral remote sensing image denoising method based on unsupervised adaptive learning.
The hyperspectral remote sensing image denoising method based on unsupervised adaptive learning provided by the invention provides an unsupervised adaptive learning strategy aiming at the problem that the model generalization is reduced by the degradation difference between the simulation data and the real data, pre-trains on a high-quality ground image, designs a discriminator to model noise, and finely adjusts denoising parameters by the discriminator when the real data is processed, thereby improving the model generalization. According to the invention, a depth denoising network based on space-spectrum residual errors and a discriminator based on global information modeling are designed in a model so as to fully excavate hyperspectral depth prior.
The hyperspectral remote sensing image denoising method based on unsupervised adaptive learning provided by the invention has the following three remarkable characteristics. Firstly, a training strategy based on unsupervised adaptive learning enables a model to learn a common image prior of a simulated image and can finely adjust denoising parameters according to a real image, so that the generalization of the model is improved; secondly, a discriminator is innovatively introduced to learn the noise mode, unsupervised learning on a real image is realized, and the global space-spectrum modeling capability of the model is effectively improved through global pooling; and thirdly, aiming at the characteristics of the hyperspectral denoising task, a space-spectrum residual volume block is designed to fully excavate the bottom layer characteristics in the model, so that the spatial detail information in the denoising result is more accurately recovered.
The invention provides a hyperspectral remote sensing image denoising method based on unsupervised adaptive learning, which comprises the following implementation steps of:
step 1, inputting a ground hyperspectral image
Figure 975323DEST_PATH_IMAGE001
Normalizing the image data to make the pixel values distributed in the range of 0-1, and then locating the pixel valuesForming a block of image of a certain size, adding analog noise to the block to generate a degraded image
Figure 325752DEST_PATH_IMAGE002
Wherein, in the step (A),HWBthe number of image rows, columns and wave bands respectively;
step 2, the initial parameters of the de-noising network are adjusted
Figure 308751DEST_PATH_IMAGE003
And discriminator initial parameter
Figure 235119DEST_PATH_IMAGE004
Carrying out initialization;
step 3, training a denoising network and a discriminator based on an algorithm of adaptive moment estimation, and setting the current iteration times ask
Step 3.1, simulating the degraded image
Figure 190437DEST_PATH_IMAGE005
Inputting the data into a denoising network, and outputting a denoising result by the denoising network through calculation
Figure 216162DEST_PATH_IMAGE006
Wherein, in the step (A),
Figure 748774DEST_PATH_IMAGE007
is as followskThe mapping function of the de-noising network at the time of the sub-iteration,
Figure 852734DEST_PATH_IMAGE008
is as followskDenoising parameters of the network during the secondary iteration;
step 3.2, respectively using the noiseless images
Figure 52771DEST_PATH_IMAGE009
And denoising results
Figure 249397DEST_PATH_IMAGE010
Inputting the image into a discriminator, and calculating the probability that the current input is a noise-free image
Figure 410251DEST_PATH_IMAGE009
Corresponding probability value of
Figure 412843DEST_PATH_IMAGE011
De-noising result
Figure 342752DEST_PATH_IMAGE010
Corresponding probability value of
Figure 477324DEST_PATH_IMAGE012
Wherein, in the step (A),
Figure 984529DEST_PATH_IMAGE013
is as followskThe mapping function of the discriminator at the time of the sub-iteration,
Figure 462914DEST_PATH_IMAGE014
is as followskParameters of the discriminator during the secondary iteration;
step 3.3, calculating the loss of the denoising network, including MSE loss
Figure 247331DEST_PATH_IMAGE015
And fight against loss
Figure 113656DEST_PATH_IMAGE016
Step 3.4, calculating the discrimination loss of the discriminator
Figure 45840DEST_PATH_IMAGE017
Figure 29714DEST_PATH_IMAGE018
(3)
WhereinnFor the set batch size,
Figure 793271DEST_PATH_IMAGE019
represents the L2 norm;
step 3.5, loss according to MSE
Figure 768180DEST_PATH_IMAGE015
And fight against loss
Figure 125343DEST_PATH_IMAGE016
Calculating gradient, updating denoising network parameters by back propagation, and outputtingk + Parameters of denoised network at 1 iteration
Figure 273427DEST_PATH_IMAGE020
According to the judgment of the loss
Figure 32436DEST_PATH_IMAGE021
Calculating gradient, back-propagating and updating discriminator parameters, and outputtingk + Parameters of the discriminator at 1 iteration
Figure 673852DEST_PATH_IMAGE022
Step 3.6, judging whether the current iteration times exceed a certain number, if not, performing step 3.1, and if so, performing step 4;
step 4, inputting a real hyperspectral image
Figure 377366DEST_PATH_IMAGE023
Normalizing the image data to make the pixel value distribution in the range of 0 to 1,
Figure 204508DEST_PATH_IMAGE024
Figure 880340DEST_PATH_IMAGE025
Figure 524948DEST_PATH_IMAGE026
the number of real image lines, columns and wave bands respectively;
step 5, carrying out algorithm based on adaptive moment estimation on the trained denoising network
Figure 591124DEST_PATH_IMAGE027
Fine-tuning is performed, wherein,
Figure 517229DEST_PATH_IMAGE028
setting the current iteration times as
Figure 844305DEST_PATH_IMAGE029
Step 6, the real degraded image is processed
Figure 331918DEST_PATH_IMAGE030
Input to post-trim denoising network
Figure 150970DEST_PATH_IMAGE031
In the method, the denoising network outputs a final denoising result through calculation
Figure 444548DEST_PATH_IMAGE032
Wherein, in the step (A),
Figure 32655DEST_PATH_IMAGE033
and 5, outputting the denoised network parameters output in the step 5.
Further, the analog noise in step 1 includes gaussian noise
Figure 927055DEST_PATH_IMAGE034
A strip noise and an impulse noise,
Figure 92457DEST_PATH_IMAGE035
is the standard deviation of gaussian noise.
Further, in the step 2, a He initialization-based mode is adopted to carry out initial parameters on the denoising network
Figure 861830DEST_PATH_IMAGE036
And discriminator initial parameter
Figure 304444DEST_PATH_IMAGE037
Initialization is performed.
Further, in step 2, the layer 1 of the denoising network is composed of a convolution layer Conv and a Parametric ReLU activation function, the layers 2-13 are space-spectrum residual convolution blocks D-block, the layer 14 is composed of a convolution layer Conv and a batch normalization layer BN, a feature map of a space-spectrum residual convolution module is output, the layers 15-16 are composed of a convolution layer Conv, a batch normalization layer BN and a Parametric ReLU activation function, finally, a denoising result is output through convolution, and the space-spectrum residual convolution block and the input and the output of the denoising network are connected by using a skip layer; the 1 st layer of the discriminator consists of convolution layer Conv and Leaky ReLU activation functions, the 2 nd to 8 th layers are LeakyConv and consist of convolution layer Conv, batch normalization layer BN and Leaky ReLU activation functions, the step length of the LeakyConv of the even layer is set to be 2, and the 9 th layer outputs a probability matrix through the convolution layer Conv and Sigmoid activation functions;
the space-spectrum residual error convolution block is formed by two layers of multi-channel two-dimensional convolution and is followed by a batch normalization layer BN, the first layer BN is followed by a Parametric ReLU activation function, and meanwhile, a feature graph output by the D-block and an input are added by using layer skipping connection.
Further, loss of MSE
Figure 868281DEST_PATH_IMAGE038
And fight against loss
Figure 786558DEST_PATH_IMAGE039
The specific calculation formula of (A) is as follows;
Figure 795840DEST_PATH_IMAGE040
Figure 217594DEST_PATH_IMAGE041
wherein the content of the first and second substances,nis the set batch size.
Further, the specific implementation manner of step 5 is as follows;
step 5.1, real degraded image
Figure 952332DEST_PATH_IMAGE042
Inputting the data into a denoising network, and outputting a denoising result by the denoising network through calculation
Figure 764430DEST_PATH_IMAGE043
Wherein, in the process,
Figure 78868DEST_PATH_IMAGE044
is as follows
Figure 620708DEST_PATH_IMAGE045
Denoising parameters of the network during the secondary iteration;
step 5.2, denoising the noise-free image
Figure 310969DEST_PATH_IMAGE046
Inputting the noise-free image into a discriminator, and calculating the probability that the current denoising result is a noiseless image by the discriminator
Figure 548046DEST_PATH_IMAGE047
Wherein, in the step (A),
Figure 462913DEST_PATH_IMAGE048
the discriminator parameter outputted in the step 3;
step 5.3, calculating the loss of the de-noising network, including the consistency constraint loss
Figure 593680DEST_PATH_IMAGE049
And unsupervised adaptive loss
Figure 670220DEST_PATH_IMAGE050
Step 5.4, loss constraint according to consistency
Figure 158708DEST_PATH_IMAGE049
And unsupervised adaptive loss
Figure 611686DEST_PATH_IMAGE051
Calculating gradient, back-propagating to update de-noising network parameters, and outputting
Figure 862539DEST_PATH_IMAGE052
Parameters of denoised network at sub-iteration
Figure 578822DEST_PATH_IMAGE053
And 5.5, judging whether the current iteration times exceed a certain number, if not, performing the step 5.1, and if so, performing the step 6.
Further, consistency constraint loss
Figure 180705DEST_PATH_IMAGE054
And unsupervised adaptive loss
Figure 437374DEST_PATH_IMAGE051
The specific calculation formula of (A) is as follows;
Figure 919564DEST_PATH_IMAGE055
Figure 931383DEST_PATH_IMAGE056
the method of the invention has the following remarkable effects: (1) after the analog image general image prior is learned, the real image prior can be further learned, and the model generalization is improved; (2) a discriminator is designed to learn the noise mode in the image, unsupervised learning on the real degraded image is realized, and the global information modeling capability of the model is further improved by using global pooling; (3) and designing a space-spectrum residual convolution block to ensure that the network utilizes the bottom layer characteristics and accurately recovers the spatial detail information.
Drawings
FIG. 1 is an overall flow diagram of an embodiment of the present invention.
Fig. 2 is a structural diagram of a denoising network, a discriminator and a spatio-spectral residual convolution block in step 2 according to the embodiment of the present invention.
Fig. 3 is a final hyperspectral image denoising result output in step 6 in the embodiment of the present invention.
Detailed Description
The technical scheme of the invention is further specifically described by the following embodiments and the accompanying drawings.
As shown in FIG. 1, the hyperspectral remote sensing image denoising method based on unsupervised adaptive learning provided by the invention comprises the following steps:
step 1, inputting 100 ground hyperspectral images from an ICVL data set
Figure 427086DEST_PATH_IMAGE057
Normalizing the image data to make the pixel value distribution in the range of 0-1, processing the pixel value distribution into 40 x 10 image blocks, and adding analog noise to the image blocks to generate a degraded image
Figure 425129DEST_PATH_IMAGE058
The analog noise includes Gaussian noise
Figure 650574DEST_PATH_IMAGE059
A strip noise and an impulse noise, wherein,HWBthe number of image rows, columns and bands,
Figure 505397DEST_PATH_IMAGE060
is the standard deviation of gaussian noise; considering the difference of the noise intensity of each wave band of the real hyperspectral image, randomly selecting the simulated noise intensity of each wave band, wherein the standard deviation of Gaussian noise
Figure 190194DEST_PATH_IMAGE061
Standard deviation of band noise
Figure 116562DEST_PATH_IMAGE062
Percentage of impulse noise
Figure 868617DEST_PATH_IMAGE063
Step 2, carrying out initial parameter pair on the denoising network based on He initialization mode
Figure 97605DEST_PATH_IMAGE064
And discriminator initial parameter
Figure 895796DEST_PATH_IMAGE065
Initializing, wherein a denoising network, a discriminator and a space-spectrum residual volume block (D-block) in the denoising network are shown in FIG. 2; the empty-spectrum residual volume block is composed of two layers of multichannel two-dimensional convolution (Conv), a Batch Normalization Layer (BN) is arranged behind the empty-spectrum residual volume block, a Parametric ReLU (PReLU) activation function is arranged behind the first Layer BN, and meanwhile, the feature graph of the D-block output and the input are added by using Layer skip connection; the layer 1 of the denoising network is composed of Conv and a PReLU activation function, the layers 2-13 are D-blocks, the layer 14 is composed of Conv and BN, a characteristic diagram of an empty-spectrum residual convolution module is output, the layers 15-16 are composed of Conv, BN and a PReLU activation function, finally, a denoising result is output through convolution, and the empty-spectrum residual convolution block and the input and the output of the denoising network are connected by using skip layers; the 1 st layer of the discriminator consists of Conv and Leaky ReLU activation functions, the 2 nd to 8 th layers are LeakyConv which consists of Conv, BN and Leaky ReLU activation functions, the step size of the LeakyConv at the even layer is set to be 2, and the 9 th layer outputs a probability matrix through the Conv and Sigmoid activation functions.
Step 3, training a denoising network and a discriminator based on an algorithm of adaptive moment estimation, setting the batch size to be 128, and setting the current iteration times to bek
Step 3.1, simulating the degraded image
Figure 501221DEST_PATH_IMAGE005
Inputting the data into a denoising network, and outputting a denoising result by the denoising network through calculation
Figure 603388DEST_PATH_IMAGE006
Wherein, in the step (A),
Figure 127911DEST_PATH_IMAGE007
is as followskThe mapping function of the de-noising network at the time of the sub-iteration,
Figure 85502DEST_PATH_IMAGE008
is as followskDenoising parameters of the network during the secondary iteration;
step 3.2, respectively using the noiseless images
Figure 963460DEST_PATH_IMAGE009
And denoising results
Figure 690107DEST_PATH_IMAGE010
Inputting the image into a discriminator, and calculating the probability that the current input is a noise-free image
Figure 651110DEST_PATH_IMAGE009
Corresponding probability value of
Figure 266637DEST_PATH_IMAGE011
De-noising result
Figure 338498DEST_PATH_IMAGE010
Corresponding probability value of
Figure 185231DEST_PATH_IMAGE012
Wherein, in the step (A),
Figure 926922DEST_PATH_IMAGE013
is as followskThe mapping function of the discriminator at the time of the sub-iteration,
Figure 921423DEST_PATH_IMAGE014
is as followskParameters of the discriminator during the secondary iteration;
step 3.3, calculating the loss of the denoising network, including MSE loss
Figure 203500DEST_PATH_IMAGE015
And fight against loss
Figure 343888DEST_PATH_IMAGE016
Wherein, in the step (A),nis the set batch size;
Figure 53218DEST_PATH_IMAGE040
Figure 535015DEST_PATH_IMAGE041
step 3.4, calculating the discrimination loss of the discriminator
Figure 620782DEST_PATH_IMAGE066
Figure 379791DEST_PATH_IMAGE018
(3)
Step 3.5, loss according to MSE
Figure 587918DEST_PATH_IMAGE015
And fight against loss
Figure 727650DEST_PATH_IMAGE016
Calculating gradient, back-propagating to update de-noising network parameters, and outputtingk + Parameters of denoised network at 1 iteration
Figure 554792DEST_PATH_IMAGE020
According to the judgment of the loss
Figure 292941DEST_PATH_IMAGE021
Calculating gradient, back-propagating and updating discriminator parameters, and outputtingk + Parameters of the arbiter at 1 iteration
Figure 344074DEST_PATH_IMAGE022
And 3.6, judging whether the current iteration times exceed 500 times, if not, performing the step 3.1, and if so, performing the step 4.
Step 4, inputting the hyperspectral image of the WHU-Hi-Baoxie unmanned aerial vehicle
Figure 675829DEST_PATH_IMAGE067
Normalizing the image data to make the pixel value distribution in the range of 0 to 1, wherein,
Figure 165716DEST_PATH_IMAGE024
Figure 191660DEST_PATH_IMAGE025
Figure 882536DEST_PATH_IMAGE026
the number of real image rows, columns and bands.
Step 5, carrying out algorithm based on adaptive moment estimation on the trained denoising network
Figure 560642DEST_PATH_IMAGE027
Fine-tuning is performed, wherein,
Figure 526324DEST_PATH_IMAGE028
setting the current iteration times as
Figure 114431DEST_PATH_IMAGE029
Step 5.1, real degraded image
Figure 835262DEST_PATH_IMAGE042
Inputting the data into a denoising network, and outputting a denoising result by the denoising network through calculation
Figure 436883DEST_PATH_IMAGE043
Wherein, in the step (A),
Figure 143939DEST_PATH_IMAGE068
is as follows
Figure 711186DEST_PATH_IMAGE045
Denoising parameters of the network during the secondary iteration;
step 5.2, denoising the noise-free image
Figure 540602DEST_PATH_IMAGE046
Inputting the noise-free image into a discriminator, and calculating the probability that the current denoising result is a noiseless image by the discriminator
Figure 334246DEST_PATH_IMAGE047
Wherein, in the step (A),
Figure 704047DEST_PATH_IMAGE048
the discriminator parameter outputted in the step 3;
step 5.3, calculating the loss of the denoising network, including consistency constraint loss
Figure 564949DEST_PATH_IMAGE049
And unsupervised adaptive loss
Figure 502949DEST_PATH_IMAGE050
Figure 642943DEST_PATH_IMAGE055
Figure 19698DEST_PATH_IMAGE056
Step 5.4, loss constraint according to consistency
Figure 171325DEST_PATH_IMAGE054
And unsupervised adaptive loss
Figure 670439DEST_PATH_IMAGE051
Calculating gradient, updating denoising network parameters by back propagation, and outputting
Figure 468369DEST_PATH_IMAGE052
Parameters of de-noising network at sub-iteration
Figure 320918DEST_PATH_IMAGE053
And 5.5, judging whether the current iteration times exceed 10 times, if not, performing the step 5.1, and if so, performing the step 6.
Step 6, the real degraded image is processed
Figure 982844DEST_PATH_IMAGE069
Input to post-trim denoising network
Figure 262647DEST_PATH_IMAGE031
In the method, the denoising network outputs a final denoising result of the hyperspectral image of the WHU-Hi-Baoxie unmanned aerial vehicle through calculation
Figure 377233DEST_PATH_IMAGE032
As shown in fig. 3, in which,
Figure 830211DEST_PATH_IMAGE033
and 5, outputting the denoised network parameters output in the step 5.
The specific embodiments described herein are merely illustrative of the spirit of the invention. Various modifications or additions may be made to the described embodiments or alternatives may be employed by those skilled in the art without departing from the spirit or ambit of the invention as defined in the appended claims.

Claims (7)

1. A hyperspectral remote sensing image denoising method based on unsupervised adaptive learning is characterized by comprising the following steps:
step 1, inputting a ground hyperspectral image
Figure 506677DEST_PATH_IMAGE001
Normalizing the image data to make the pixel values distributed in the range of 0-1, processing the pixel values into image blocks with a certain size, and adding analog noise to the image blocks to generate a degraded image
Figure 82015DEST_PATH_IMAGE002
Wherein, in the step (A),HWBthe number of image rows, columns and wave bands respectively;
step 2, carrying out initial parameters on the denoising network
Figure 559264DEST_PATH_IMAGE003
And discriminator initial parameter
Figure 143829DEST_PATH_IMAGE004
Carrying out initialization;
step 3, training a denoising network and a discriminator based on an algorithm of adaptive moment estimation, and setting the current iteration times ask
Step 3.1, simulating the degraded image
Figure 186871DEST_PATH_IMAGE005
Inputting the data into a denoising network, and outputting a denoising result by the denoising network through calculation
Figure 74056DEST_PATH_IMAGE006
Wherein, in the process,
Figure 632076DEST_PATH_IMAGE007
is as followskThe mapping function of the de-noising network at the time of the sub-iteration,
Figure 394234DEST_PATH_IMAGE008
is a firstkDenoising parameters of the network during the secondary iteration;
step 3.2, respectively using the noiseless images
Figure 619679DEST_PATH_IMAGE009
And denoising results
Figure 474502DEST_PATH_IMAGE010
Inputting the image into a discriminator, and calculating the probability that the current input is a noise-free image
Figure 660764DEST_PATH_IMAGE009
Corresponding probability value of
Figure 587132DEST_PATH_IMAGE011
De-noising result
Figure 542449DEST_PATH_IMAGE010
Corresponding probability value of
Figure 335218DEST_PATH_IMAGE012
Wherein, in the step (A),
Figure 867831DEST_PATH_IMAGE013
is as followskThe mapping function of the discriminator at the time of the sub-iteration,
Figure 207676DEST_PATH_IMAGE014
is as followskParameters of the discriminator during the secondary iteration;
step 3.3, calculating the loss of the denoising network, including MSE loss
Figure 407714DEST_PATH_IMAGE015
And to combat the loss
Figure 73181DEST_PATH_IMAGE016
Step 3.4, calculating the discrimination loss of the discriminator
Figure 765194DEST_PATH_IMAGE017
Figure 298943DEST_PATH_IMAGE018
(3)
WhereinnFor the set batch size to be set,
Figure 461809DEST_PATH_IMAGE019
represents the L2 norm;
step 3.5, loss according to MSE
Figure 829337DEST_PATH_IMAGE015
And fight against loss
Figure 336541DEST_PATH_IMAGE016
Calculating gradient, back-propagating to update de-noising network parameters, and outputtingk + De-noising at 1 iterationParameters of a network
Figure 80506DEST_PATH_IMAGE020
According to the judgment of the loss
Figure 864923DEST_PATH_IMAGE021
Calculating gradient, back-propagating and updating discriminator parameters, and outputtingk + Parameters of the arbiter at 1 iteration
Figure 731248DEST_PATH_IMAGE022
Step 3.6, judging whether the current iteration times exceed a certain number, if not, performing step 3.1, and if so, performing step 4;
step 4, inputting a real hyperspectral image
Figure 651316DEST_PATH_IMAGE023
Normalizing the image data to make the pixel value distribution in the range of 0 to 1, wherein,
Figure 402234DEST_PATH_IMAGE024
Figure 431370DEST_PATH_IMAGE025
Figure 140700DEST_PATH_IMAGE026
the number of real image lines, columns and wave bands respectively;
step 5, carrying out algorithm based on adaptive moment estimation on the trained denoising network
Figure 497863DEST_PATH_IMAGE027
Fine-tuning is performed, wherein,
Figure 911527DEST_PATH_IMAGE028
setting the current iteration times as
Figure 965809DEST_PATH_IMAGE029
Step 6, the real degraded image is processed
Figure 783723DEST_PATH_IMAGE030
Input to post-trim denoising network
Figure 752816DEST_PATH_IMAGE031
In the method, the denoising network outputs a final denoising result through calculation
Figure 642275DEST_PATH_IMAGE032
Wherein, in the step (A),
Figure 990210DEST_PATH_IMAGE033
and 5, outputting the denoised network parameters output in the step 5.
2. The hyperspectral remote sensing image denoising method based on unsupervised adaptive learning of claim 1, wherein: the analog noise in step 1 includes gaussian noise
Figure 634818DEST_PATH_IMAGE034
A strip noise and an impulse noise,
Figure 264776DEST_PATH_IMAGE035
is the standard deviation of gaussian noise.
3. The hyperspectral remote sensing image denoising method based on unsupervised adaptive learning of claim 1, wherein: step 2, adopting a He-based initialization mode to carry out initial parameters on the denoising network
Figure 895609DEST_PATH_IMAGE036
And discriminator initial parameter
Figure 488264DEST_PATH_IMAGE037
Initialization is performed.
4. The hyperspectral remote sensing image denoising method based on unsupervised adaptive learning of claim 1, wherein: in the step 2, the layer 1 of the denoising network consists of a convolution layer Conv and a Parametric ReLU activation function, the layers 2-13 are space-spectrum residual convolution blocks D-block, the layer 14 consists of the convolution layer Conv and a batch normalization layer BN, a characteristic diagram of a space-spectrum residual convolution module is output, the layers 15-16 consist of the convolution layer Conv, the batch normalization layer BN and the Parametric ReLU activation function, finally, a denoising result is output through convolution, and the input and the output of the space-spectrum residual convolution blocks and the denoising network are connected by using skip layers; the 1 st layer of the discriminator consists of convolution layer Conv and Leaky ReLU activation functions, the 2 nd to 8 th layers are LeakyConv and consist of convolution layer Conv, batch normalization layer BN and Leaky ReLU activation functions, the step length of the LeakyConv of the even layer is set to be 2, and the 9 th layer outputs a probability matrix through the convolution layer Conv and Sigmoid activation functions;
the space-spectrum residual error convolution block is formed by two layers of multi-channel two-dimensional convolution and is followed by a batch normalization layer BN, the first layer BN is followed by a Parametric ReLU activation function, and meanwhile, a feature graph output by the D-block and an input are added by using layer skipping connection.
5. The hyperspectral remote sensing image denoising method based on unsupervised adaptive learning of claim 1, wherein: loss of MSE
Figure 179140DEST_PATH_IMAGE038
And fight against loss
Figure 857246DEST_PATH_IMAGE039
The specific calculation formula of (A) is as follows;
Figure 291769DEST_PATH_IMAGE040
Figure 175149DEST_PATH_IMAGE041
wherein the content of the first and second substances,nis the set batch size.
6. The hyperspectral remote sensing image denoising method based on unsupervised adaptive learning of claim 1, wherein: the specific implementation manner of the step 5 is as follows;
step 5.1, real degraded image
Figure 895981DEST_PATH_IMAGE042
Inputting the data into a denoising network, and outputting a denoising result by the denoising network through calculation
Figure 999066DEST_PATH_IMAGE043
Wherein, in the step (A),
Figure 971701DEST_PATH_IMAGE044
is as follows
Figure 273369DEST_PATH_IMAGE045
Denoising parameters of the network during the secondary iteration;
step 5.2, denoising the noise-free image
Figure 306048DEST_PATH_IMAGE046
Inputting the noise-free image into a discriminator, and calculating the probability that the current denoising result is a noiseless image by the discriminator
Figure 958746DEST_PATH_IMAGE047
Wherein, in the step (A),
Figure 761836DEST_PATH_IMAGE048
the discriminator parameter outputted in the step 3;
step 5.3, calculating the loss of the de-noising network, including the consistency constraint loss
Figure 324535DEST_PATH_IMAGE049
And unsupervised adaptive loss
Figure 652748DEST_PATH_IMAGE050
Step 5.4, loss constraint according to consistency
Figure 668109DEST_PATH_IMAGE049
And unsupervised adaptive loss
Figure 779284DEST_PATH_IMAGE051
Calculating gradient, back-propagating to update de-noising network parameters, and outputting
Figure 321124DEST_PATH_IMAGE052
Parameters of de-noising network at sub-iteration
Figure 194140DEST_PATH_IMAGE053
And 5.5, judging whether the current iteration times exceed a certain number, if not, performing the step 5.1, and if so, performing the step 6.
7. The hyperspectral remote sensing image denoising method based on unsupervised adaptive learning of claim 6, wherein: loss of consistency constraint
Figure 821431DEST_PATH_IMAGE054
And unsupervised adaptive loss
Figure 673980DEST_PATH_IMAGE051
The specific calculation formula of (A) is as follows;
Figure 335905DEST_PATH_IMAGE055
Figure 678025DEST_PATH_IMAGE056
CN202210409199.8A 2022-04-19 2022-04-19 Hyperspectral remote sensing image denoising method based on unsupervised adaptive learning Active CN114511473B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210409199.8A CN114511473B (en) 2022-04-19 2022-04-19 Hyperspectral remote sensing image denoising method based on unsupervised adaptive learning

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210409199.8A CN114511473B (en) 2022-04-19 2022-04-19 Hyperspectral remote sensing image denoising method based on unsupervised adaptive learning

Publications (2)

Publication Number Publication Date
CN114511473A true CN114511473A (en) 2022-05-17
CN114511473B CN114511473B (en) 2022-07-05

Family

ID=81555535

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210409199.8A Active CN114511473B (en) 2022-04-19 2022-04-19 Hyperspectral remote sensing image denoising method based on unsupervised adaptive learning

Country Status (1)

Country Link
CN (1) CN114511473B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115546653A (en) * 2022-12-02 2022-12-30 耕宇牧星(北京)空间科技有限公司 Remote sensing image denoising method based on depth enhancement network
CN116071265A (en) * 2023-02-21 2023-05-05 哈尔滨工业大学 Image denoising method based on space self-adaptive self-supervision learning
CN116861167A (en) * 2023-06-12 2023-10-10 河北工程大学 FBG spectrum cyclic denoising method based on deep learning
CN117876692A (en) * 2024-03-11 2024-04-12 中国石油大学(华东) Feature weighted connection guided single-image remote sensing image denoising method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110320162A (en) * 2019-05-20 2019-10-11 广东省智能制造研究所 A kind of semi-supervised high-spectral data quantitative analysis method based on generation confrontation network
CN110991636A (en) * 2019-11-14 2020-04-10 东软医疗系统股份有限公司 Training method and device of generative confrontation network, image enhancement method and equipment
US20210304364A1 (en) * 2020-03-26 2021-09-30 Bloomberg Finance L.P. Method and system for removing noise in documents for image processing
CN113808042A (en) * 2021-09-16 2021-12-17 沈阳工业大学 SAR image denoising method based on wavelet transformation and generation countermeasure network
WO2022041521A1 (en) * 2020-08-31 2022-03-03 浙江大学 Low-dose sinogram denoising and pet image reconstruction method based on teacher-student generators

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110320162A (en) * 2019-05-20 2019-10-11 广东省智能制造研究所 A kind of semi-supervised high-spectral data quantitative analysis method based on generation confrontation network
CN110991636A (en) * 2019-11-14 2020-04-10 东软医疗系统股份有限公司 Training method and device of generative confrontation network, image enhancement method and equipment
US20210304364A1 (en) * 2020-03-26 2021-09-30 Bloomberg Finance L.P. Method and system for removing noise in documents for image processing
WO2022041521A1 (en) * 2020-08-31 2022-03-03 浙江大学 Low-dose sinogram denoising and pet image reconstruction method based on teacher-student generators
CN113808042A (en) * 2021-09-16 2021-12-17 沈阳工业大学 SAR image denoising method based on wavelet transformation and generation countermeasure network

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
X. WANG,ET AL.: "A Self-Supervised Denoising Network for Satellite-Airborne-Ground Hyperspectral Imagery", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
冯旭斌: "基于深度学习的光学遥感图像去噪与超分辨率重建算法研究", 《CNKI博士论文》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115546653A (en) * 2022-12-02 2022-12-30 耕宇牧星(北京)空间科技有限公司 Remote sensing image denoising method based on depth enhancement network
CN116071265A (en) * 2023-02-21 2023-05-05 哈尔滨工业大学 Image denoising method based on space self-adaptive self-supervision learning
CN116071265B (en) * 2023-02-21 2023-08-08 哈尔滨工业大学 Image denoising method based on space self-adaptive self-supervision learning
CN116861167A (en) * 2023-06-12 2023-10-10 河北工程大学 FBG spectrum cyclic denoising method based on deep learning
CN116861167B (en) * 2023-06-12 2023-12-15 河北工程大学 FBG spectrum cyclic denoising method based on deep learning
CN117876692A (en) * 2024-03-11 2024-04-12 中国石油大学(华东) Feature weighted connection guided single-image remote sensing image denoising method
CN117876692B (en) * 2024-03-11 2024-05-17 中国石油大学(华东) Feature weighted connection guided single-image remote sensing image denoising method

Also Published As

Publication number Publication date
CN114511473B (en) 2022-07-05

Similar Documents

Publication Publication Date Title
CN114511473B (en) Hyperspectral remote sensing image denoising method based on unsupervised adaptive learning
CN109035142B (en) Satellite image super-resolution method combining countermeasure network with aerial image prior
Vasamsetti et al. Wavelet based perspective on variational enhancement technique for underwater imagery
CN111028177A (en) Edge-based deep learning image motion blur removing method
Sonogashira et al. High-resolution bathymetry by deep-learning-based image superresolution
CN111008936B (en) Multispectral image panchromatic sharpening method
CN110648292A (en) High-noise image denoising method based on deep convolutional network
CN111179196B (en) Multi-resolution depth network image highlight removing method based on divide-and-conquer
CN112767267B (en) Image defogging method based on simulation polarization fog-carrying scene data set
CN111738954B (en) Single-frame turbulence degradation image distortion removal method based on double-layer cavity U-Net model
CN112330569A (en) Model training method, text denoising method, device, equipment and storage medium
CN110807428B (en) Coal sample identification method, device, server and storage medium
Feng et al. Turbugan: An adversarial learning approach to spatially-varying multiframe blind deconvolution with applications to imaging through turbulence
Hai et al. Advanced retinexnet: a fully convolutional network for low-light image enhancement
US20220067882A1 (en) Image processing device, computer readable recording medium, and method of processing image
Li et al. Efficient burst raw denoising with variance stabilization and multi-frequency denoising network
CN115760641B (en) Remote sensing image cloud and fog removing method and equipment based on multiscale characteristic attention network
Huang et al. Attention-based for multiscale fusion underwater image enhancement
CN115829870A (en) Image denoising method based on variable scale filtering
CN114758030A (en) Underwater polarization imaging method integrating physical model and deep learning
Su et al. Deconvolution of defocused image with multivariate local polynomial regression and iterative wiener filtering in DWT domain
CN113450267A (en) Transfer learning method capable of rapidly acquiring multiple natural degradation image restoration models
He et al. An Unsupervised Dehazing Network with Hybrid Prior Constraints for Hyperspectral Image
Sangeetha et al. Performance analysis of exemplar based image inpainting algorithms for natural scene image completion
CN114359103B (en) Hyperspectral image defogging method, hyperspectral image defogging device, computer product and storage medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant