DE102021103293A1 - Ein Bayes'sches Bildentrauschungsverfahren, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert - Google Patents

Ein Bayes'sches Bildentrauschungsverfahren, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert Download PDF

Info

Publication number
DE102021103293A1
DE102021103293A1 DE102021103293.8A DE102021103293A DE102021103293A1 DE 102021103293 A1 DE102021103293 A1 DE 102021103293A1 DE 102021103293 A DE102021103293 A DE 102021103293A DE 102021103293 A1 DE102021103293 A1 DE 102021103293A1
Authority
DE
Germany
Prior art keywords
distribution
noise
function
distance
noisy
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.)
Pending
Application number
DE102021103293.8A
Other languages
English (en)
Inventor
Yuxiang Xing
Li Zhang
Hewei Gao
Zhi Deng
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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Publication of DE102021103293A1 publication Critical patent/DE102021103293A1/de
Pending legal-status Critical Current

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
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/22Matching criteria, e.g. proximity measures
    • 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
    • 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/047Probabilistic or stochastic networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/60Image enhancement or restoration using machine learning, e.g. neural networks
    • 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/048Activation functions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20076Probabilistic image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Computing Systems (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Health & Medical Sciences (AREA)
  • Computational Linguistics (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Operations Research (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)

Abstract

Die vorliegende Erfindung bezieht sich auf Gebiet für digitale Bildverarbeitung und Computer Vision, insbesondere auf ein Bayes'sches Bildentrauschungsverfahren, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert. In der Bayes'schen posterioren Wahrscheinlichkeitstheorie hängt die Schätzung der rauschfreien Bilder aus den rauschhaltigen Bilder von der Modellierung der vorherigen Verteilung der rauschfreien Bilder ab. Die vorliegende Erfindung schlägt zunächst ein Verfahren zum Erhalten der Verteilung von rauschfreien Bildern durch Lernen aus rauschhaltigen Bildproben unter Kenntnis von Verteilungsmodell der additiven Räusche vor. Das vorliegende Verfahren transformiert die Beschränkung der vorherigen Verteilung der rauschfreien Bilder in die Beschränkung der vorherigen Verteilung der rauschhaltigen Bilder zum Bayes'schen Entrauschen. Die vorliegende Erfindung schlägt ferner ein Bayes'sches Entrauschungsrealisierungsverfahren zum überwachungsfreien Trainieren eines bildentrauschenden neuronalen Netzwerks vor. Durch das vorliegende Verfahren wird ermöglicht, rauschhaltige Bildproben vollständig auszunutzen, um die impliziten Verteilungseigenschaften von rauschfreien Bildern genau zu lernen, wodurch ein effizientes Entrauschen von Bildern realisiert wird.

Description

  • Technisches Gebiet
  • Die vorliegende Erfindung bezieht sich auf Gebiet für digitale Bildverarbeitung und Computer Vision, insbesondere auf ein Bayes'sches Bildentrauschungsverfahren, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert.
  • Stand der Technik
  • Bei der Bildaufnahme von Fotografie, Überwachung der öffentlichen Sicherheit, medizinischer Bildgebung, mikroskopischer Bildgebung usw. enthalten die aufgenommenen Bilder aufgrund von Faktoren wie zufälligen Eigenschaften physikalischer Signale, Umgebungsstörungen und Zustandsänderungen des Systems selbst usw. unweigerliche Räusche, die erhebliche Auswirkungen auf Bildqualität haben. Das Bildentrauschungsverfahren zielt darauf ab, die zufälligen Rauschkomponenten im Bild zu reduzieren und die effektive Information der Bilder wiederherzustellen. Es ist ein wichtiges Verfahren und eine wichtige Technologie bei der Bildverarbeitung.
  • Klassische Bildentrauschungsverfahren umfassen viele Filtermethoden wie Mittelwertfilterung, Medianfilterung, Gaußsche Tiefpassfilterung usw.. Diese Filter begrenzen die rauschbezogenen Komponenten in Bilder, um den Effekt der Geräuschentfernung zu erzielen. Das Filterverfahren kann sowohl in der Bilddomäne als auch in der Transformationsdomäne nach einer spezifischen Transformation der Bilder, wie beispielsweise einer Frequenzdomäne, einer Waveletdomäne usw., implementiert werden. Filterverfahren haben jedoch inhärente Nachteile bei der Aufrechterhaltung von Bilddetails. Zusätzlich zu den analytischen Filtermethoden ist die iterative Bildoptimierungsmethode, die auf der Bayes'schen posterioren Wahrscheinlichkeitsverteilung basiert, eine andere Art der Bildentrauschungsmethode. Bei dieser Art von Methode sind zusätzliche Strafbegriffe wie die Minimierung der Gesamtvariation und die Minimierung der potentiellen Energiefunktion des Markov-Zufallsfeldes eingeführt, und somit erreicht den Zweck der Rauschentfernung. Durch das Bayes'sche posteriore Schätzverfahren kann eine genauere rauschfreie Bildschätzung erzeugt werden, aber sein Ergebnis hängt wesentlich von der auferlegten Strafe ab.
  • In den letzten Jahren hat die Anwendung von Deep-Learning-Methoden zur Bildentrauschung mit der Entwicklung künstlicher neuronaler Netze und der Deep-Learning-Technologie den Effekt der Bildentrauschung erheblich verbessert. Unter dem Entfernen von Räusche können Bilddetails genauer wiederhergestellt werden als mit klassischen Methoden. Obwohl verschiedene neuronale Netze in vielen Anwendungen leistungsstarke Funktionen zum Entrauschen und Wiederherstellen von Bildern aufweisen, basieren sie alle auf einer großen Anzahl von Trainingsdatensätzen. Überwachtes Lernen erfordert gepaarte rauschhaltige Bilder und rauschfreie Bilder als Grundwahrheit, um das Netzwerk gemeinsam zu trainieren. Beim generativen kontroversen Lernen können die Verteilungsregeln und Verteilungseigenschaften von rauschfreien Bildern aus ungepaarten rauschfreien Bilddatensätzen extrahiert und somit kann eine Bildentrauschung durchgeführt werden, wodurch die Notwendigkeit einer Datenpaarung entfällt. In den meisten Anwendungen ist es jedoch oft schwierig oder teuer, echte rauschfreie Daten zu erhalten, was die Anwendung von Deep Learning beim Entrauschen von Bildern erschwert.
  • Darstellung der Erfindung
  • Der Zweck der vorliegenden Erfindung besteht darin, ein Bayes'sches Bildentrauschungsverfahren vorzuschlagen, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert, das den praktisch leicht verfügbaren rauschhaltigen Datensatz und das Rauschmodell verwendet, um auf die Verteilungseigenschaften der rauschfreien Bilder zu schließen, und dann die Verteilung der rauschfreien Bilder verwendet, um die Schätzung der rauschfreien Bilder und das Entrauschen der rauschhaltigen Bilder zu erreichen.
  • Das durch die vorliegende Erfindung vorgeschlagenene Bayes'sche Bildentrauschungsverfahren, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert, umfasst folgende Schritte:
    • (1) Berechnen einer Entfernungskostenfunktion einer Verteilung von rauschfreien Bilder gemäß einer Verteilung von rauschhaltigen Bilder;
    • (2) Durchführen einer Bayes'schen Entrauschung an den rauschhaltigen Bilder gemäß der Entfernungskostenfunktion der Verteilung der rauschfreien Bilder in Schritt (1).
  • Im obigen Verfahrensschritt (1) wird die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder gemäß der Verteilung der rauschhaltigen Bilder berechnet. Für Bilder, deren Verteilung genau oder annähernd analytisch ausgedrückt werden kann, wird ein auf einer charakteristischen Funktion basierende Verfahren verwendet, um die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder zu schätzen, einschließlich der folgenden Schritte:
    • (1) Aufzeichnen eines vom Bildgebungssystem erzeugten rauschhaltigen Bild-Zufallsvektors als Z, Schätzen einer Wahrscheinlichkeitsdichtefunktion fz des rauschhaltigen Bild-Zufallsvektors Z in analytischer Form, und Aufzeichnen eines additiven Rausch-Zufallsvektors, der in dem rauschhaltigen Bild-Zufallsvektor Z enthalten ist, als N , Schätzen einer Wahrscheinlichkeitsdichtefunktion fn des additiven Rausch-Zufallsvektors N in analytischer Form;
    • (2) Durchführen jeweils einer Fourier-Transformation der Wahrscheinlichkeitsdichtefunktionen fz und fN in Schritt (1), um eine charakteristische Funktion Fz des rauschhaltigen Bild-Zufallsvektors und eine charakteristische Funktion FN des additiven Rausch-Zufallsvektors zu erhalten;
    • (3) Aufzeichnen eines vom Bildgebungssystem erzeugten rauschfreien Bild-Zufallsvektors als s, und Dividieren der charakteristischen Funktion Fz des rauschhaltigen Bild-Zufallsvektors durch die charakteristische Funktion FN des additiven Rausch-Zufallsvektors in Schritt (2), um eine charakteristische Funktion Fs des rauschfreien Bild-Zufallsvektors S zu erhalten: F S = F Z / F N ,
      Figure DE102021103293A1_0001
      oder: Durch ein Zieloptimierungsverfahren wird ein geschätzter Wert Fs* der charakteristischen Funktion Fs des Zufallsvektors der rauschfreien Bilder gemäß einer Konsistenz von Fs mit Fz/FN erhalten: F S * = argmin F S Φ ( F S , F Z / F N ) ,
      Figure DE102021103293A1_0002
      wobei Φ(g,g) ein Abstandsmaß zwischen den charakteristischen Funktionen zweier Zufallsvektoren ist, und argmin einen Parameter darstellt, bei dem eine Zielfunktion einen Minimalwert annimmt;
    • (4) Durchführen einer inversen Fourier-Transformation für die charakteristische Funktion Fs oder den geschätzten Wert Fs* der charakteristischen Funktion Fs des Zufallsvektors der rauschfreien Bilder in Schritt (3), um eine Wahrscheinlichkeitsdichtefunktion fs des Zufallsvektors der rauschfreien Bilder zu erhalten, d.h. eine Verteilung des Zufallsvektors S der rauschfreien Bilder zu erhalten;
    • (5) Durchführen einer analytischen Filterentrauschung oder einer iterativen Bayes'schen Entrauschung an einer Gruppe von rauschhaltigen Bildproben z1, z2, z3,...zn, um einen Anfangswert von geschätzten rauschfreien Bildproben Ŝ1,2,3,...n zu erhalten, und Aufzeichnen eines Satzes von ŝ1,ŝ2,3,...,ŝn als ŝ;
    • (6) Berechnen eines Ausdrucks einer Verteilungsentfernung D(ŝ,ŝ) zwischen dem Satz ŝ der geschätzten rauschfreien Bildproben Ŝ1,2,3,...n in Schritt (5) und dem in Schritt (4) erhaltenen Zufallsvektor s der rauschfreien Bilder gemäß einer Verteilungsentfernungsfunktion D, d.h. Entfernungskostenfunktion der Verteilung der rauschfreien Bilder, wobei die Methode zur Berechnung der Verteilungsentfernung D(ŝ,s) ein Verfahren zum Identifizieren der Funktion der Verteilungsentfernung des neuronalen Netzes ist, einschließlich der folgenden Schritte:
      • (6-1) Einstellen der Verteilungsentfernungsfunktion D so, dass sie ein diskriminantes neuronales Netzwerk ΩD als Verteilungsentfernungsoperator enthält, und eine Parametergruppe des diskriminanten neuronalen Netzwerks ΩD ist θD;
      • (6-2) Eingeben des in Schritt (5) von Anspruch 2 beschriebenen Satzes ŝ der geschätzten rauschfreien Bildproben ŝ1,2,3,...ŝn und des in Schritt (4) von Anspruch 2 erhaltenen Zufallsvektors s der rauschfreien Bilder in die Verteilungsentfernungsfunktion D, und Verwenden des diskriminanten neuronalen Netzwerks ΩD in Schritt (1) zur Berechnung eines Kreuzentropie-Abstandes zwischen dem Satz ŝ und dem Zufallsvektor s der rauschfreien Bilder: D ( S ^ , S ) = i = l : n log ( Ω D ( s i ; θ D ) ) + log ( 1 Ω D ( s ^ i ; θ D ) )
        Figure DE102021103293A1_0003
        oder eines Wasserstein-Abstandes: D ( S ^ , S ) = i = l : n | Ω D ( s i ; θ D ) Ω D ( s ^ i ; θ D ) | ;
        Figure DE102021103293A1_0004
    • Die obigen S1,S2,S3,....,Sn sind n Abtastproben des Zufallsvektors S der rauschfreien Bilder; oder Verwenden des diskriminanten neuronalen Netzwerks ΩD in Schritt (1), um eine Diskriminanznetzwerk-Verteilungsentfernungskostenfunktion von generativen gegnerischen Netzwerken Methoden im Fachbereich zu berechnen. Die genannten generativen gegnerischen Netzwerken Methoden umfassen eines von SGAN, RSGAN, RaSGAN, LSGAN, RaLSGAN, HingeGAN, WGAN-GP, RSGAN-GP oder RaSGAN-GP
  • Im obigen Verfahrensschritt (1) wird die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder gemäß der Verteilung der rauschhaltigen Bilder berechnet. Für den Fall, dass die Verteilung schwer analytisch auszudrücken ist, oder für Bilder, deren Verteilung bei Berechnung zu kompliziert ist, wird ein indirektes Verfahren zur Beschränkung einer Verteilungsentfernung verwendet, um die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder zu schätzen, einschließlich der folgenden Schritte:
    • (1) Durchführen einer analytischen Filterentrauschung oder einer iterativen Bayes'schen Entrauschung an einer Gruppe von rauschhaltigen Bildproben z1,z2,z3,... , zn, um einen Anfangswert von geschätzten rauschfreien Bildproben Ŝ1,2,3,...n zu erhalten;
    • (2) Einstellen eines Zufallsvektors von in den rauschhaltigen Bildproben Z1,Z2,Z3,... , Zn in Schritt (1) enthaltenen additiven Räusche als N , Schätzen der Wahrscheinlichkeitsdichtefunktion ƒN des Zufallsvektors N für additive Räusche in analytischer Form, und Aufnahme von Rauschproben n 1 * , n 2 * , n 3 * , , n n *
      Figure DE102021103293A1_0005
      mit einer gleichen Anzahl von der rauschhaltigen Bildproben aus der Verteilung des Zufallsvektors N für additive Räusche; oder: Aufnahme von Rauschproben des Bildgebungssystems, wenn kein Signal zu erhalten ist, um n 1 * , n 2 * , n 3 * , , n n *
      Figure DE102021103293A1_0006
      zu erhalten;
    • (3) Addieren der Rauschproben n 1 * , n 2 * , n 3 * , , n n *
      Figure DE102021103293A1_0007
      in Schritt (2) mit der geschätzten rauschfreien Bildproben Ŝ1,2,3,...n in Schritt (1), um geschätzte rauschhaltige Bildproben ẑ1,z2,z3,...ẑn zu erhalten, wobei z ^ i = s ^ i + n i * ,
      Figure DE102021103293A1_0008
      i=1:n;
    • (4) Aufzeichnen eines Satzes der rauschhaltigen Bildproben z1,z2,z3,..., zn in Schritt (1) als Z und eines Satzes der geschätzten rauschhaltigen Bildproben ẑ1,ẑ2,3,...,ẑn in Schritt (3) als Ẑ, Berechnen einer Verteilungsentfernung D(Ẑ,Z) zwischen dem Satz Ẑ der diskreten geschätzten rauschhaltigen Bildproben und dem Satz Z der wahren diskreten rauschhaltigen Bildproben gemäß einer Verteilungsentfernungsfunktion D, wobei ein Ausdruck der Verteilungsentfernung D(Ẑ,Z) eine Entfernungskostenfunktion der Verteilung der rauschfreien Bilder ist, wobei die Methode zur Berechnung der Verteilungsentfernung D(Ẑ,Z) ein Verfahren zum Identifizieren der Funktion der Verteilungsentfernung des neuronalen Netzes ist, einschließlich der folgenden Schritte:
      • (4-1) Einstellen der Verteilungsentfernungsfunktion D so, dass sie ein diskriminantes neuronales Netzwerk ΩD als Verteilungsentfernungsoperator enthält, und eine Parametergruppe des diskriminanten neuronalen Netzwerks ΩD ist θD;
      • (4-2) Eingeben des Satzes Ẑ der diskreten geschätzten rauschhaltigen Bildproben und dem Satz Z der wahren diskreten rauschhaltigen Bildproben in die Verteilungsentfernungsfunktion D, und Verwenden des diskriminanten neuronalen Netzwerks ΩD in Schritt (1) zur Berechnung eines Kreuzentropie-Abstandes zwischen dem Satz Ẑ der diskreten geschätzten rauschhaltigen Bildproben und dem Satz Z der wahren diskreten rauschhaltigen Bildproben: D ( Z ^ , Z ) = i = l : n log ( Ω D ( z i ; θ D ) ) + log ( 1 Ω D ( z ^ i ; θ D )
        Figure DE102021103293A1_0009
        oder eines Wasserstein-Abstandes: D ( Z ^ , Z ) = i = l : n | Ω D ( z i ; θ D ) Ω D ( z ^ i ; θ D ) | .
        Figure DE102021103293A1_0010
        Oder wird ein diskriminantes neuronales Netzwerk ΩD in Schritt (1) verwendet, um eine Diskriminanznetzwerk-Verteilungsentfernungskostenfunktion von generativen gegnerischen Netzwerken Methoden zu berechnen. Die Diskriminanznetzwerk-Verteilungsentfernungskostenfunktion von der genannten generativen gegnerischen Netzwerken Methoden ist eines von SGAN, RSGAN, RaSGAN, LSGAN, RaLSGAN, HingeGAN, WGAN-GP, RSGAN-GP oder RaSGAN-GP
  • Die Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N in analytischer Form ist dabei derart geschätzt: Gemäß den physikalischen Gesetzen des Bildgebungssystems wird eine Verteilungsform des additiven Rausch-Zufallsvektors N erhalten, und dann werden Verteilungsparameter additiven Rausch-Zufallsvektors N durch Simulation oder tatsächliche Experimente gemessen, so dass die Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N erhalten wird. Oder wird eine gemischte Gaußsche Verteilung zum Anpassen an eine unbekannte Verteilung verwendet. Nehmen zunächst an, dass der Zufallsvektor N für additive Räusche der gemischten Gaußschen Verteilung entspricht. Und dann werden die gemischten Gaußschen Verteilungsparameter durch Simulation oder tatsächliche Experimente gemessen, so dass die Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N erhalten wird.
  • Das Verfahren zum Schätzen der Wahrscheinlichkeitsdichtefunktion fN des additiven Rausch-Zufallsvektors N in analytischer Form lautet: Für nichtadditive Rauschverteilung kann sie durch eine spezifische mathematische Transformation in ungefähre additive Räusche umgewandelt werden, und dann wird die Wahrscheinlichkeitsdichtefunktion fN des additiven Rausch-Zufallsvektors N erhalten.
  • Bei dem obigen Verfahren lautet das Verfahren zum Schätzen der Wahrscheinlichkeitsdichtefunktion fz des rauschhaltigen Bild-Zufallsvektors z in analytischer Form: Die Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors z wird durch Zählen von Frequenzen einer großen Anzahl von rauschhaltigen Bildproben erhalten. Oder lautet das Verfahren zum Schätzen der Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors z in analytischer Form: Verteilungsmodellannahmen im Fachbereich werden verwendet, um eine Verteilungsform des rauschhaltigen Bild-Zufallsvektors z zu erhalten, wobei Verteilungsparameter des rauschhaltigen Bild-Zufallsvektors z durch Simulation oder tatsächliche Experimente gemessen werden, so dass die Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors z geschätzt wird.
  • „Durchführen einer Bayes'schen Entrauschung an den rauschhaltigen Bilder gemäß der Entfernungskostenfunktion der Verteilung der rauschfreien Bilder“ im obigen Verfahrensschritt (2) umfasst folgende Schritte:
    • (1) Für rauschhaltige Bildproben z1,z2,z3,..., Zn und geschätzte rauschfreie Bildproben ŝ1,ŝ2,ŝ3,...,ŝn wird eine Bayes'sche posteriore Verteilungskostenfunktion Ψ konstruiert: ψ ( s ^ 1 , s ^ 2 , s ^ 3 , , s ^ n ) = i = l : n φ f ( z i s ^ i ) + λ D ,
      Figure DE102021103293A1_0011
      wobei Ψ die Bayes'sche posteriore Verteilungskostenfunktion ist, wobei bei der Bayes'schen posterioren Verteilungskostenfunktion der erste Term φf ein Datentreue-Term zwischen z1,z2,z3,..., zn und ŝ1,2,3,...ŝn ist, und der zweite Term D ein Verteilungsabstand D(ŝ,s) oder ein Verteilungsabstand D(Ẑ,Z) ist, der λ̂ ein Hyperparameter ist, durch den eine Intensität des Verteilungsabstands D(Ŝ,S) eingestellt wird;
    • (2) Die Bayes'sche posteriore Verteilungskostenfunktion Ψ von Schritt (1) wird minimiert und optimiert, und die endgültigen geschätzten rauschfreien Bilder ŝ1,2,3,...ŝn werden erhalten, um eine Bildentrauschung zu realisieren. Das Verfahren zur Minimierung und Optimierung der Bayes'schen posterioren Verteilungskostenfunktion Ψ ist ein übliches klassisches mathematisches Optimierungsverfahren, wie ein Lagrange-Multiplikatorverfahren, ein alternierendes Iterationsverfahren, eine alternative konjugierte Methode der kleinsten Fehler-Quadrate oder ein Gradientenoptimierungsalgorithmus.
  • Das Verfahren zur Minimierung und Optimierung der Bayes'schen posterioren Verteilungskostenfunktion Ψ kann auch eine Optimierungsmethode für neuronale Netze sein, die folgende Schritte umfasst:
    • (1) Ein bildentrauschendes neuronales Netzwerk ΩF wird als Bildentrauschungsoperator ŝiF(ziE), i=1:n konstruiert, wobei θE eine Parametergruppe des bildentrauschenden neuronalen Netzwerks ist.
    • (2) Die Bayes'sche posteriore Verteilungskostenfunktion Ψ wird als Trainingskostenfunktion des bildentrauschenden neuronalen Netzwerks verwendet, und eine große Anzahl von rauschhaltigen Bildprobendatensätzen wird verwendet, um das bildentrauschende neuronale Netzwerk zu trainieren.
    • (3) Die rauschhaltigen Bildproben werden als Eingabe des bildrauschenden neuronalen Netzwerks verwendet, und die geschätzten rauschfreien Bilder werden ausgegeben und erhalten. Diese rauschfreien Bilder sind also die Bilder, die nach der Optimierung und der Minimierung des neuronalen Netzwerks erhalten werden.
  • Die Verteilungsentfernungsfunktion D in dem obigen Verfahren ist eine analytische dominante Verteilungsentfernungsfunktion, wobei die analytische dominante Verteilungsentfernungsfunktion eine K-L-Divergenzdominanzverteilungsentfernungsfunktion oder eine J-S-Divergenzdominanzverteilungsentfernungsfunktion umfasst.
  • Das in der vorliegenden Erfindung vorgeschlagene Bayes'sche Bildentrauschungsverfahren, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert, hat die folgenden Vorteile:
    1. 1. Bei dem erfindungsgemäßen Verfahren zum Lernen der Verteilung von rauschfreien Bildern basierend auf der Verteilung von rauschhaltigen Bildern kann die Verteilung der rauschfreien Bilder indirekt durch Schätzung der Verteilung der rauschhaltigen Bilder geschätzt werden, wenn die vom Signal unabhängige additive Rauschverteilung bekannt ist. Diese Transformation der Verteilungsbeschränkungen vermeidet die Notwendigkeit rauschfreier Bilddatensätze für die traditionelle rauschfreie Bildverteilungsmodellierung, wie z. B. Wörterbuchlernmethoden, und erleichtert die Implementierung der rauschfreien Bildverteilungsmodellierung.
    2. 2. Die a priori Bildbeschränkungen, die in dem erfindungsgemäßen Verfahren betroffen sind, werden aus den Datensätze abgeleitet und nicht künstlich festgelegt, so dass es Flexibilität und Datenspezifität gibt. Verglichen mit den künstlich auferlegten Einschränkungen bei dem traditionellen Bayes'schen Schätzverfahren sind die a priori Bildbeschränkungen im erfindungsgemäßen Verfahren genauer und die rauschfreien Bilder können genauer geschätzt werden.
    3. 3. Das Verfahren in der vorliegenden Erfindung vernetzt die Verteilungsentfernungsbeschränkung und den Bayes'schen Schätzprozess und nutzt die leistungsfähigen Bildverarbeitungs- und Charakterisierungsfähigkeiten des neuronalen Netzwerks vollständig aus, so dass die Beschränkung der Bildverteilung leicht zu implementieren ist. Gleichzeitig kann das Netzwerk das Bild nach Abschluss des Trainings von Ende zu Ende auf analytische Weise entrauschen, wodurch der Rechenaufwand und die Berechnungszeit im iterativen Prozess der Bayes'schen Schätzmethode erheblich reduziert werden.
  • Figurenliste
    • 1 ist ein schematisches Diagramm des bildentrauschenden neuronalen Netzwerks und seiner Kostenfunktion in der zweiten Ausführungsform der vorliegenden Erfindung.
    • 2 ist ein schematisches Diagramm der Implementierung spezifischer neuronaler Netzwerkstrukturen, die in der zweiten und dritten Ausführungsform der vorliegenden Erfindung betroffen sind, wobei (a) ein bildentrauschendes neuronales Netzwerk ist und (b) ein diskriminierendes neuronales Netzwerk ist.
    • 3 ist ein schematisches Diagramm des bildentrauschenden neuronalen Netzwerks, des diskriminierenden neuronalen Netzwerks und der Kostenfunktionen der beiden in der dritten Ausführungsform der vorliegenden Erfindung.
  • Detaillierte Beschreibung
  • Bei dem durch die vorliegende Erfindung vorgeschlagenenen Bayes'schen Bildentrauschungsverfahren, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert:
    • (1) Berechnen einer Entfernungskostenfunktion einer Verteilung von rauschfreien Bilder gemäß einer Verteilung von rauschhaltigen Bilder;
    • (2) Durchführen einer Bayes'schen Entrauschung an den rauschhaltigen Bilder gemäß der Entfernungskostenfunktion der Verteilung der rauschfreien Bilder in Schritt (1).
  • Im obigen Verfahrensschritt (1) wird die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder gemäß der Verteilung der rauschhaltigen Bilder berechnet. Für Bilder, deren Verteilung genau oder annähernd analytisch ausgedrückt werden kann, wird ein auf einer charakteristischen Funktion basierende Verfahren verwendet, um die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder zu schätzen, einschließlich der folgenden Schritte:
    • (1) Aufzeichnen eines vom bei diesem Verfahren verwendeten Bildgebungssystem erzeugten rauschhaltigen Bild-Zufallsvektors als Z , Schätzen einer Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors Z in analytischer Form, und Aufzeichnen eines additiven Rausch-Zufallsvektors, der in dem rauschhaltigen Bild-Zufallsvektor Z enthalten ist, als N , Schätzen einer Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N in analytischer Form;
    • (2) Durchführen jeweils einer Fourier-Transformation der Wahrscheinlichkeitsdichtefunktionen ƒz und ƒN in Schritt (1), um eine charakteristische Funktion Fz des rauschhaltigen Bild-Zufallsvektors und eine charakteristische Funktion FN des additiven Rausch-Zufallsvektors zu erhalten;
    • (3) Aufzeichnen eines vom bei diesem Verfahren verwendeten Bildgebungssystem erzeugten rauschfreien Bild-Zufallsvektors als s , und Dividieren der charakteristischen Funktion Fz des rauschhaltigen Bild-Zufallsvektors durch die charakteristische Funktion FN des additiven Rausch-Zufallsvektors in Schritt (2), um eine charakteristische Funktion Fs des rauschfreien Bild-Zufallsvektors s zu erhalten: F S = F Z / F N ,
      Figure DE102021103293A1_0012
      oder: Durch ein Zieloptimierungsverfahren wird ein geschätzter Wert Fs* der charakteristischen Funktion Fs des Zufallsvektors der rauschfreien Bilder gemäß einer Konsistenz von Fs mit Fz/FN erhalten: F S * = argmin F S Φ ( F S , F Z / F N ) ,
      Figure DE102021103293A1_0013
      wobei Φ(g,g) ein Abstandsmaß zwischen den charakteristischen Funktionen zweier Zufallsvektoren ist, und argmin einen Parameter darstellt, bei dem eine Zielfunktion einen Minimalwert annimmt;
    • (4) Durchführen einer inversen Fourier-Transformation für die charakteristische Funktion Fs oder den geschätzten Wert Fs* der charakteristischen Funktion Fs des Zufallsvektors der rauschfreien Bilder in Schritt (3), um eine Wahrscheinlichkeitsdichtefunktion ƒs des Zufallsvektors der rauschfreien Bilder zu erhalten, d.h. eine Verteilung des Zufallsvektors S der rauschfreien Bilder zu erhalten;
    • (5) Durchführen einer analytischen Filterentrauschung oder einer iterativen Bayes'schen Entrauschung an einer Gruppe von rauschhaltigen Bildproben z1,z2,z3,...,zn, die mit diesem Verfahren entrauscht werden sollen, um einen Anfangswert von geschätzten rauschfreien Bildproben Ŝ1,2,3,...n zu erhalten, und Aufzeichnen eines Satzes von ŝ1,23,...,ŝn als ŝ;
    • (6) Berechnen eines Ausdrucks einer Verteilungsentfernung D(ŝ,ŝ) zwischen dem Satz ŝ der geschätzten rauschfreien Bildproben Ŝ1,2,3,...n in Schritt (5) und dem in Schritt (4) erhaltenen Zufallsvektor \s der rauschfreien Bilder gemäß einer Verteilungsentfernungsfunktion D, d.h. Entfernungskostenfunktion der Verteilung der rauschfreien Bilder, wobei die Methode zur Berechnung der Verteilungsentfernung D(Ŝ,S) ein Verfahren zum Identifizieren der Funktion der Verteilungsentfernung des neuronalen Netzes ist, einschließlich der folgenden Schritte:
    • (6-1) Einstellen der Verteilungsentfernungsfunktion D so, dass sie ein diskriminantes neuronales Netzwerk ΩD als Verteilungsentfernungsoperator enthält, und eine Parametergruppe des diskriminanten neuronalen Netzwerks ΩD ist θD;
    • (6-2) Eingeben des in Schritt (5) von Anspruch 2 beschriebenen Satzes ŝ der geschätzten rauschfreien Bildproben ŝ1,2,3,...ŝn und des in Schritt (4) von Anspruch 2 erhaltenen Zufallsvektors s der rauschfreien Bilder in die Verteilungsentfernungsfunktion D, und Verwenden des diskriminanten neuronalen Netzwerks ΩD in Schritt (1) zur Berechnung eines Kreuzentropie-Abstandes zwischen dem Satz ŝ und dem Zufallsvektor s der rauschfreien Bilder: D ( S ^ , S ) = i = l : n log ( Ω D ( s i ; θ D ) ) + log ( 1 Ω D ( s ^ i ; θ D )
      Figure DE102021103293A1_0014
      oder eines Wasserstein-Abstandes: D ( S ^ , S ) = i = l : n | Ω D ( s i ; θ D ) Ω D ( s ^ i ; θ D ) | ;
      Figure DE102021103293A1_0015
      Die obigen s1,s2,s3,...,sn sind n Abtastproben des Zufallsvektors s der rauschfreien Bilder; oder Verwenden des diskriminanten neuronalen Netzwerks ΩD in Schritt (1), um eine Diskriminanznetzwerk-Verteilungsentfernungskostenfunktion von anderen üblichen generativen gegnerischen Netzwerken Methoden im Fachbereich zu berechnen. Die Diskriminanznetzwerk-Verteilungsentfernungskostenfunktion von der genannten generativen gegnerischen Netzwerken Methoden ist eines von SGAN, RSGAN, RaSGAN, LSGAN, RaLSGAN, HingeGAN, WGAN-GP, RSGAN-GP oder RaSGAN-GP
  • Im obigen Verfahrensschritt (1) wird die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder gemäß der Verteilung der rauschhaltigen Bilder berechnet. Für den Fall, dass die Verteilung schwer analytisch auszudrücken ist, oder für Bilder, deren Verteilung bei Berechnung zu kompliziert ist, wird ein indirektes Verfahren zur Beschränkung einer Verteilungsentfernung verwendet, um die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder zu schätzen, einschließlich der folgenden Schritte:
    • (1) Durchführen einer analytischen Filterentrauschung oder einer iterativen Bayes'schen Entrauschung an einer Gruppe von rauschhaltigen Bildproben z1,z2,z3,...,zn, die mit diesem Verfahren entrauscht werden sollen, um einen Anfangswert von geschätzten rauschfreien Bildproben Ŝ1,2,3,...n zu erhalten;
    • (2) Einstellen eines Zufallsvektors von in den rauschhaltigen Bildproben z1,z2,z3,..., zn in Schritt (1) enthaltenen additiven Räusche als N , Schätzen der Wahrscheinlichkeitsdichtefunktion ƒN des Zufallsvektors N für additive Räusche in analytischer Form, und Aufnahme von Simulationsrauschproben n 1 * , n 2 * , n 3 * , , n n *
      Figure DE102021103293A1_0016
      mit einer gleichen Anzahl von der rauschhaltigen Bildproben aus der Verteilung des Zufallsvektors N für additive Räusche; oder: Aufnahme von Rauschproben des Bildgebungssystems, wenn kein Signal zu erhalten ist, um n 1 * , n 2 * , n 3 * , , n n *
      Figure DE102021103293A1_0017
      zu erhalten;
    • (3) Addieren der Rauschproben n 1 * , n 2 * , n 3 * , , n n *
      Figure DE102021103293A1_0018
      in Schritt (2) mit der geschätzten rauschfreien Bildproben Ŝ1,2,3,...n in Schritt (1), um geschätzte rauschhaltige Bildproben ẑ1,z2,z3,...ẑn zu erhalten, wobei z ^ i = s ^ i + n i * ,
      Figure DE102021103293A1_0019
      i=1:n;
    • (4) Aufzeichnen eines Satzes der rauschhaltigen Bildproben z1,z2,z3,..., zn in Schritt (1) als Z und eines Satzes der geschätzten rauschhaltigen Bildproben ẑ1,z2,z3,...ẑn in Schritt (3) als Ẑ, Berechnen einer Verteilungsentfernung D(Ẑ,Z) zwischen dem Satz Ẑ der diskreten geschätzten rauschhaltigen Bildproben und dem Satz Z der wahren diskreten rauschhaltigen Bildproben gemäß einer Verteilungsentfernungsfunktion D, wobei ein Ausdruck der Verteilungsentfernung D(Ẑ,Z) eine Entfernungskostenfunktion der Verteilung der rauschfreien Bilder ist, wobei die Methode zur Berechnung der Verteilungsentfernung D(Ẑ,Z) ein Verfahren zum Identifizieren der Funktion der Verteilungsentfernung des neuronalen Netzes ist, einschließlich der folgenden Schritte:
      • (4-1) Einstellen der Verteilungsentfernungsfunktion D so, dass sie ein diskriminantes neuronales Netzwerk ΩD als Verteilungsentfernungsoperator enthält, und eine Parametergruppe des diskriminanten neuronalen Netzwerks ΩD ist θD;
      • (4-2) Eingeben des Satzes Ẑ der diskreten geschätzten rauschhaltigen Bildproben und dem Satz Z der wahren diskreten rauschhaltigen Bildproben in die Verteilungsentfernungsfunktion D, und Verwenden des diskriminanten neuronalen Netzwerks ΩD in Schritt (1) zur Berechnung eines Kreuzentropie-Abstandes zwischen dem Satz Ẑ der diskreten geschätzten rauschhaltigen Bildproben und dem Satz Z der wahren diskreten rauschhaltigen Bildproben: D ( Z ^ , Z ) = i = l : n log ( Ω D ( z i ; θ D ) ) + log ( 1 Ω D ( z ^ i ; θ D )
        Figure DE102021103293A1_0020
        oder eines Wasserstein-Abstandes: D ( Z ^ , Z ) = i = l : n | Ω D ( z i ; θ D ) Ω D ( z ^ i ; θ D ) | .
        Figure DE102021103293A1_0021
        Oder wird ein diskriminantes neuronales Netzwerk ΩD in Schritt (1) verwendet, um eine Diskriminanznetzwerk-Verteilungsentfernungskostenfunktion von generativen gegnerischen Netzwerken Methoden im Fachbereich zu berechnen. Die genannten generativen gegnerischen Netzwerken Methoden umfassen eines von SGAN, RSGAN, RaSGAN, LSGAN, RaLSGAN, HingeGAN, WGAN-GP, RSGAN-GP oder RaSGAN-GP
  • Die Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N in analytischer Form ist dabei derart geschätzt: Gemäß den physikalischen Gesetzen des bei diesem Verfahren verwendeten Bildgebungssystems wird eine Verteilungsform des additiven Rausch-Zufallsvektors N erhalten, und dann werden Verteilungsparameter additiven Rausch-Zufallsvektors N durch Simulation oder tatsächliche Experimente gemessen, so dass die Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N erhalten wird. Oder wird eine gemischte Gaußsche Verteilung zum Anpassen an eine unbekannte Verteilung verwendet. Nehmen zunächst an, dass der Zufallsvektor N für additive Räusche der gemischten Gaußschen Verteilung entspricht. Und dann werden die gemischten Gaußschen Verteilungsparameter durch Simulation oder tatsächliche Experimente gemessen, so dass die Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N erhalten wird.
  • Das Verfahren zum Schätzen der Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N in analytischer Form lautet: Für nichtadditive Rauschverteilung kann sie durch eine spezifische mathematische Transformation in ungefähre additive Räusche umgewandelt werden, und dann wird die Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N erhalten.
  • Bei dem obigen Verfahren lautet das Verfahren zum Schätzen der Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors z in analytischer Form: Die Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors z wird durch Zählen von Frequenzen einer großen Anzahl von rauschhaltigen Bildproben erhalten. Oder lautet das Verfahren zum Schätzen der Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors z in analytischer Form: Verteilungsmodellannahmen im Fachbereich werden verwendet, um die Verteilungsform des rauschhaltigen Bild-Zufallsvektors z zu erhalten, wobei die Verteilungsparameter des rauschhaltigen Bild-Zufallsvektors z durch Simulation oder tatsächliche Experimente gemessen werden, so dass die Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors z geschätzt wird.
  • „Durchführen einer Bayes'schen Entrauschung an den rauschhaltigen Bilder gemäß der Entfernungskostenfunktion der Verteilung der rauschfreien Bilder“ im obigen Verfahrensschritt (2) umfasst folgende Schritte:
    • (1) Für eine Gruppe von rauschhaltigen Bildproben z1,z2,z3,..., zn, die mit diesem Verfahren entrauscht werden sollen, und geschätzte rauschfreie Bildproben ŝ1,2,3,...,ŝn wird eine Bayes'sche posteriore Verteilungskostenfunktion Ψ konstruiert: ψ ( s ^ 1 , s ^ 2 , s ^ 3 , , s ^ n ) = i = l : n φ f ( z i s ^ i ) + λ D ,
      Figure DE102021103293A1_0022
      wobei Ψ die Bayes'sche posteriore Verteilungskostenfunktion ist, wobei bei der Bayes'schen posterioren Verteilungskostenfunktion der erste Term φf ein Datentreue-Term zwischen z1,z2,z3,.., zn und ŝ1,ŝ2,ŝ3,...ŝn ist, und der zweite Term D ein Verteilungsabstand D(ŝ,s) oder ein Verteilungsabstand D(Ẑ,Z) ist, der λ̂ ein Hyperparameter ist, durch den eine Intensität des Verteilungsabstands D(Ŝ,S) eingestellt wird;
    • (2) Die Bayes'sche posteriore Verteilungskostenfunktion Ψ von Schritt (1) wird minimiert und optimiert, und die endgültigen geschätzten rauschfreien Bilder ŝ1,ŝ2,ŝ3,...,ŝn werden erhalten, um eine Bildentrauschung zu realisieren. Das Verfahren zur Minimierung und Optimierung der Bayes'schen posterioren Verteilungskostenfunktion Ψ ist ein übliches klassisches mathematisches Optimierungsverfahren, wie ein Lagrange-Multiplikatorverfahren, ein alternierendes Iterationsverfahren, eine alternative konjugierte Methode der kleinsten Fehler-Quadrate oder ein Gradientenoptimierungsalgorithmus.
  • Das Verfahren zur Minimierung und Optimierung der Bayes'schen posterioren Verteilungskostenfunktion Ψ kann auch eine Optimierungsmethode für neuronale Netze sein, die folgende Schritte umfasst:
    • (1) Ein bildentrauschendes neuronales Netzwerk ΩF wird als Bildentrauschungsoperator ŝiF(ziE), i=1:n konstruiert, wobei θE eine Parametergruppe des bildentrauschenden neuronalen Netzwerks ist.
    • (2) Die Bayes'sche posteriore Verteilungskostenfunktion Ψ wird als Trainingskostenfunktion des bildentrauschenden neuronalen Netzwerks verwendet, und eine große Anzahl von rauschhaltigen Bildprobendatensätzen wird verwendet, um das bildentrauschende neuronale Netzwerk zu trainieren.
    • (3) Die rauschhaltigen Bildproben werden als Eingabe des bildrauschenden neuronalen Netzwerks verwendet, und die geschätzten rauschfreien Bilder werden ausgegeben und erhalten. Diese rauschfreien Bilder sind also die Bilder, die nach der Optimierung und der Minimierung des neuronalen Netzwerks erhalten werden.
  • Die Verteilungsentfernungsfunktion D in dem obigen Verfahren ist eine analytische dominante Verteilungsentfernungsfunktion, wobei die analytische dominante Verteilungsentfernungsfunktion eine K-L-Divergenzdominanzverteilungsentfernungsfunktion oder eine J-S-Divergenzdominanzverteilungsentfernungsfunktion umfasst.
  • Deep Learning hat in Bildentrauschungsanwendungen eine gute Leistung gezeigt. Im Vergleich zu klassischen Filter- und Entrauschungsmethoden hat Deep Learning eine signifikante Verbesserung der Wirkung erzielt, Deep Learning-Methoden stützen sich jedoch stark auf rauschfreie Bilddatensätze. Die strengen Anforderungen an rauschfreie Datensätze begrenzen den Implementierungseffekt und den Anwendungsbereich von Deep Learning bei Problemen mit der Entrauschung realer Bilder. In der Bayes'schen posterioren Wahrscheinlichkeitstheorie hängt die Schätzung der rauschfreien Bilder aus den rauschhaltigen Bilder von der Modellierung der vorherigen Verteilung der rauschfreien Bilder ab. Die vorliegende Erfindung schlägt zunächst ein Verfahren zum Erhalten der Verteilung von rauschfreien Bildern durch Lernen aus rauschhaltigen Bildproben unter Kenntnis von Verteilungsmodell der additiven Räusche vor. Die vorliegende Erfindung transformiert ferner die Beschränkung der vorherigen Verteilung der rauschfreien Bilder in die Beschränkung der vorherigen Verteilung der rauschhaltigen Bilder zum Bayes'schen Entrauschen. Die vorliegende Erfindung schlägt ferner ein Bayes'sches Entrauschungsrealisierungsverfahren zum überwachungsfreien Trainieren eines bildentrauschenden neuronalen Netzwerks vor. Die vorliegende Erfindungsidee besteht darin, rauschhaltige Bildproben vollständig auszunutzen, um die impliziten Verteilungseigenschaften von rauschfreien Bildern genau zu lernen, wodurch ein effizientes Entrauschen von Bildern realisiert wird.
  • Das Bild in dem Verfahren der vorliegenden Erfindung bezieht sich auf ein allgemeines räumliches zweidimensionales oder dreidimensionales Signal, einschließlich allgemeiner Bilder mit sichtbarem Licht, klinisch-medizinischen Bildern, mikroskopischen Bildern usw.. Dieses Verfahren kann auf viele Bereiche wie Fotografie, Bildgebung der öffentlichen Sicherheit, medizinische Bildgebung und mikroskopische Bildgebung usw. angewendet werden.
  • Das Folgende beschreibt Ausführungsformen des Verfahrens der vorliegenden Erfindung:
    • Beispiel 1: Entrauschen von rauschhaltigen handgeschriebenen digitalen Bildern
  • Das bei diesem Verfahren verwendete Bildgebungssystem ist: Handgeschriebenes digitales Erfassungssystem.
  • Die rauschhaltigen Bilder sind: rauschhaltige handgeschriebene digitale Bilder, die Bildauflösung beträgt 28*28.
  • Die rauschfreien Bilder sind: rauschfreie handgeschriebene digitale Bilder, die Bildauflösung beträgt 28*28.
    • (1) Berechnen der Entfernungskostenfunktion der Verteilung der rauschfreien handgeschriebenen digitalen Bilder gemäß der Verteilung der rauschhaltigen handgeschriebenen digitalen Bilder;
    • (1-1) Aufzeichnen des vom handgeschriebenen digitalen Erfassungssystem erzeugten Zufallsvektors der rauschhaltigen handgeschriebenen digitalen Bilder als z, Schätzen der Wahrscheinlichkeitsdichtefunktion ƒz des Zufallsvektors Z der rauschhaltigen handgeschriebenen digitalen Bilder in analytischer Form,
    • (1-2) Vom handgeschriebenen digitalen Erfassungssystem sind 70,000 rauschhaltige handgeschriebene digitale Bilder gesammelt. Die Wahrscheinlichkeitsdichtefunktion ƒz des Zufallsvektors Z der rauschhaltigen handgeschriebenen digitalen Bilder wird durch Zählen von Frequenzen der 70,000 rauschhaltigen Bildproben erhalten.
    • (1-3) Aufzeichnen des additiven Rausch-Zufallsvektors, der in dem Zufallsvektor Z der rauschhaltigen handgeschriebenen digitalen Bilder enthalten ist, als N, Schätzen der Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N in analytischer Form;
    • (1-4) Eine gemischte Gaußsche Verteilung wird zum Anpassen an eine additive Rauschverteilung verwendet. Nehmen zunächst an, dass jedes Pixel des Zufallsvektors N für additive Räusche der fünfdimensionalen gemischten Gauß-Verteilung entspricht, wobei Rauschproben aus dem signalfreien Teil der rauschhaltigen handgeschriebenen digitalen Bilder erhalten werden und die Parameter der fünfdimensionalen Gaußschen Mischungsverteilung mittels dieser Rauschproben durch eine Maximum-Likelihood-Methode geschätzt werden, so dass die Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N erhalten wird.
    • (1-5) Durchführen jeweils einer Fourier-Transformation der Wahrscheinlichkeitsdichtefunktionen ƒz und ƒN, um die charakteristische Funktion Fz des Zufallsvektors der rauschhaltigen handgeschriebenen digitalen Bilder und die charakteristische Funktion FN des additiven Rausch-Zufallsvektors zu erhalten;
    • (1-6) Aufzeichnen des vom handgeschriebenen digitalen Erfassungssystem erzeugten Zufallsvektors der rauschfreien handgeschriebenen digitalen Bilder als S , und Dividieren der charakteristischen Funktion Fz des Zufallsvektors der rauschhaltigen handgeschriebenen digitalen Bilder durch die charakteristische Funktion FN des additiven Rausch-Zufallsvektors, um die charakteristische Funktion Fs des Zufallsvektors S der rauschfreien handgeschriebenen digitalen Bilder zu erhalten: F S = F Z / F N ;
      Figure DE102021103293A1_0023
    • (2) Durchführen einer inversen Fourier-Transformation für die charakteristische Funktion Fs des Zufallsvektors der rauschfreien handgeschriebenen digitalen Bilder, um die Wahrscheinlichkeitsdichtefunktion ƒs des Zufallsvektors der rauschfreien handgeschriebenen digitalen Bilder zu erhalten, d.h. die Verteilung des Zufallsvektors s der rauschfreien handgeschriebenen digitalen Bilder zu erhalten, einschließlich der folgenden Schritte:
      • (2-1) Durchführen einer analytischen Filterentrauschung nun an einer Gruppe von rauschhaltigen handgeschriebenen digitalen Bildproben z1,z2,z3,..., zn, die mit diesem Verfahren entrauscht werden sollen, um den Anfangswert der geschätzten rauschfreien handgeschriebenen digitalen Bildproben Ŝ1,2,3,...n zu erhalten, und Aufzeichnen eines Satzes von ŝ1,ŝ2,ŝ3,...,nals ŝ;
      • (2-2) Berechnen eines Ausdrucks der Verteilungsentfernung D(ŝ,ŝ) zwischen dem Satz ŝ von geschätzten rauschfreien handgeschriebenen digitalen Bildproben ŝ1,2,3,...,ŝn und dem erhaltenen Zufallsvektor s der rauschfreien handgeschriebenen digitalen Bilder gemäß einer K-L-Verteilungsentfernungsfunktion D : D ( S ^ , S ) = i = 1 n 1 n ( ln 1 n ƒ s ( s ^ i ) )
        Figure DE102021103293A1_0024
    • (3) Ein auf der charakteristischen Funktion basierendes Verfahren wird verwendet, um die Entfernungskostenfunktion der Verteilung der rauschfreien handgeschriebenen digitalen Bilder zu schätzen, wobei eine Bayes'sche Entrauschung an den rauschhaltigen handgeschriebenen digitalen Bilder durchgeführt wird, einschließlich der folgenden Schritte:
      • (3-1) Für die Gruppe von rauschhaltigen handgeschriebenen digitalen Bildproben z1,z2,z3,.., zn, die mit diesem Verfahren entrauscht werden sollen, und entsprechende geschätzte rauschfreie handgeschriebenen digitalen Bildproben ŝ1,2,3,...ŝn wird eine Bayes'sche posteriore Verteilungskostenfunktion Ψ konstruiert: ψ ( s ^ 1 , s ^ 2 , s ^ 3 , , s ^ n ) = i = l : n φ f ( z i s ^ i ) + λ D ,
        Figure DE102021103293A1_0025
        wobei Ψ die Bayes'sche posteriore Verteilungskostenfunktion ist, wobei bei der Bayes'schen posterioren Verteilungskostenfunktion der erste Term φr ein Datentreue-Term zwischen z1,z2,z3,..., zn und ŝ1,2,3,...ŝn ist, und der zweite Term D also D(ŝ,s) ist, der λ̂ ein Hyperparameter ist, durch den eine Intensität des Verteilungsabstands D(Ŝ,S) eingestellt wird;
      • (3-2) Die Bayes'sche posteriore Verteilungskostenfunktion Ψ wird durch ein Gradientenoptimierungsalgorithmus minimiert und optimiert, und die endgültigen geschätzten rauschfreien handgeschriebenen digitalen Bilder ŝ1,2,3,...,ŝn werden erhalten, um eine Bildentrauschung zu realisieren.
  • Beispiel 2: Entrauschen von nachtfotografischen Bildern einer Kamera
  • Das bei diesem Verfahren verwendete Bildgebungssystem ist: Kamera.
  • Die rauschhaltigen Bilder sind: rauschhaltige fotografische Bildblöcke, die Bildauflösung beträgt 512*512.
  • Die rauschfreien Bilder sind: rauschfreie fotografische Bildblöcke, die Bildauflösung beträgt 512*512.
  • Die Auflösung der ursprünglichen rauschhaltigen fotografischen Bilder ist normalerweise größer als 512*512. Aus ursprünglichen rauschhaltigen fotografischen Bilder können rauschhaltige fotografische Bildblöcke von 512*512 allmählich abgefangen werden und wird jeder rauschhaltiger fotografischer Bildblock entrauscht, um geschätzte rauschfreie fotografische Bildblöcke zu erhalten. Die rauschfreien fotografischen Bildblöcke werden zusammengefügt, um geschätzte rauschfreie fotografische Bilder zu erhalten.
    • (1) Berechnen der Entfernungskostenfunktion der Verteilung der rauschfreien fotografischen Bildblöcke gemäß der Verteilung der rauschhaltigen fotografischen Bildblöcke; Dabei wird das indirekte Verteilungsentfernungsbeschränkungsverfahren verwendet, um die Entfernungskostenfunktion der Verteilung der rauschfreien fotografischen Bildblöcke zu schätzen, einschließlich der folgenden Schritte:
      • (1-1) Durchführen einer analytischen Filterentrauschung an einer Gruppe von rauschhaltigen fotografischen Bildblockproben z1,z2,z3,..., zn, die mit diesem Verfahren entrauscht werden sollen, um den Anfangswert der geschätzten rauschfreien Bildblockproben Ŝ1,2,3,...n zu erhalten;
      • (1-2) Einstellen eines Zufallsvektors von in den rauschhaltigen fotografischen Bildblockproben z1,z2,z3,..., zn enthaltenen additiven Räusche als N , Schätzen der Wahrscheinlichkeitsdichtefunktion ƒN des Zufallsvektors N für additive Räusche in analytischer Form;
      • (1-3) Eine gemischte Gaußsche Verteilung wird zum Anpassen an eine additive Rauschverteilung verwendet. Nehmen zunächst an, dass jedes Pixel des Zufallsvektors N für additive Räusche der fünfdimensionalen gemischten Gauß-Verteilung entspricht, wobei Rauschproben gemessen werden, wenn das Kameraobjektiv vollständig blockiert ist, und die Parameter der fünfdimensionalen Gaußschen Mischungsverteilung mittels dieser Rauschproben durch eine Maximum-Likelihood-Methode geschätzt werden, so dass die Wahrscheinlichkeitsdichtefunktion des additiven Rausch-Zufallsvektors ƒN erhalten wird.
      • (1-4) Und Aufnahme von Simulationsrauschproben n 1 * , n 2 * , n 3 * , , n n *
        Figure DE102021103293A1_0026
        mit einer gleichen Anzahl von der rauschhaltigen Bildproben aus der Verteilung des Zufallsvektors N für additive Räusche;
  • Addieren der Rauschproben n 1 * , n 2 * , n 3 * , , n n *
    Figure DE102021103293A1_0027
    mit der geschätzten rauschfreien fotografischen Bildblockproben ŝ1,ŝ2,ŝ3,...n, um geschätzte rauschhaltige fotografischen Bildblockproben ẑ1,z2,z3,...ẑn zu erhalten,
    wobei z ^ i = s ^ i + n i * ,
    Figure DE102021103293A1_0028
    i=1:n;
  • Aufzeichnen eines Satzes von rauschhaltigen fotografischen Bildblockproben z1,z2,z3,.., zn als Z und eines Satzes von geschätzten rauschhaltigen Bildproben ẑ1,ẑ2,ẑ3,...,ẑn in Schritt (3) als Ẑ . Berechnen eines Ausdrucks einer Verteilungsentfernung D(Ẑ,Z) zwischen dem Satz Ẑ von diskreten geschätzten rauschhaltigen fotografischen Bildblockproben und dem Satz Z der wahren diskreten rauschhaltigen fotografischen Bildblockproben mittels einer J-S-Verteilungsentfernungsfunktion D:
    • Schätzen der Wahrscheinlichkeitsdichtefunktion ƒz des Zufallsvektors Z der rauschhaltigen fotografischen Bildblöcke zunächst in analytischer Form: Verwenden einer Modellannahme der gemeinsamen Gaußschen Verteilung im Fachbereich, um die Verteilungsform des Zufallsvektors Z der rauschhaltigen fotografischen Bildblöcke zu erhalten. Der Zufallsvektor Z der rauschhaltigen fotografischen Bildblöcke ist von fünfhundertdimensionaler unabhängiger Standard-Gauß-Verteilung durch lineare Transformation erzeugt. Durch eine Maximum-Likelihood-Methode werden Verteilungsparameter des Zufallsvektors Z der rauschhaltigen fotografischen Bildblöcke, d.h. ein linearer Transformationskoeffizient, ermittelt, um die Wahrscheinlichkeitsdichtefunktion ƒz des Zufallsvektors Z der rauschhaltigen Bilder zu erhalten.
  • Schätzen der Wahrscheinlichkeitsdichtefunktion ƒ des Zufallsvektors ẑ der rauschhaltigen fotografischen Bildblöcke auf die gleiche Weise in analytischer Form: D ( Z ^ , Z ) = [ i = 1 n 1 n ( ln 1 n ƒ z ( z ^ i ) ) + i = 1 n 1 n ( ln 1 n ƒ z ^ ( z i ) ) ] / 2,
    Figure DE102021103293A1_0029
    d.h. die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder.
    • (2) Durchführen einer Bayes'schen Entrauschung an den rauschhaltigen fotografischen Bildblöcke gemäß der Entfernungskostenfunktion der Verteilung der rauschfreien fotografischen Bildblöcke:
      • (2-1) Für die Gruppe von rauschhaltigen fotografischen Bildblockproben z1,z2,z3,..., zn, die mit diesem Verfahren entrauscht werden sollen, und entsprechende geschätzte rauschfreie fotografischen Bildblockproben ŝ1,2,3,...ŝn wird eine Bayes'sche posteriore Verteilungskostenfunktion Ψ konstruiert: ψ ( s ^ 1 , s ^ 2 , s ^ 3 , , s ^ n ) = i = l : n φ f ( z i s ^ i ) + λ D ,
        Figure DE102021103293A1_0030
        wobei Ψ die Bayes'sche posteriore Verteilungskostenfunktion ist, wobei bei der Bayes'schen posterioren Verteilungskostenfunktion der erste Term φr ein Datentreue-Term zwischen z1,z2,z3,..., zn und ŝ1,2,3,...ŝn ist, und der zweite Term D also D(Ẑ,Z) ist, der λ̂ ein Hyperparameter ist, durch den eine Intensität des Verteilungsabstands D(Ẑ,Z) eingestellt wird;
      • (2-2) Bei der Bayes'schen posterioren Verteilungskostenfunktion Ψ wird eine Optimierungsmethode für neuronale Netze verwendet, und die endgültigen geschätzten rauschfreien fotografischen Bildblöcke ŝ1,ŝ2,ŝ3,...ŝn werden erhalten, um eine Bildentrauschung zu realisieren.
  • Ein bildentrauschendes neuronales Netzwerk ΩF wird als Bildentrauschungsoperator ŝiF(ziE), i=1:n konstruiert, wobei θE eine Parametergruppe des bildentrauschenden neuronalen Netzwerks ist; Die Bayes'sche posteriore Verteilungskostenfunktion Ψ wird als Trainingskostenfunktion des bildentrauschenden neuronalen Netzwerks verwendet, und eine große Anzahl von rauschhaltigen Bildprobendatensätzen wird verwendet, um das bildentrauschende neuronale Netzwerk zu trainieren: min θ E i = l : n φ f ( z i s ^ i ) + λ D .
    Figure DE102021103293A1_0031
    Die rauschhaltigen fotografischen Bildblöcke werden als Eingabe des bildrauschenden neuronalen Netzwerks verwendet, und die geschätzten rauschfreien fotografischen Bildblöcke werden ausgegeben und erhalten. Diese rauschfreien fotografischen Bildblöcke sind also die Bilder, die nach der Optimierung und der Minimierung des neuronalen Netzwerks erhalten werden.
  • Das bildentrauschende neuronale Netzwerk und seine Kostenfunktion sind wie in 1 gezeigt, wobei eine spezifische Struktur des bildentrauschenden neuronalen Netzwerks wie in 2 (a) gezeigt ist.
  • Beispiel 3: Entrauschen von Röntgen-CT-Projektionsbildern
  • Das bei diesem Verfahren verwendete Bildgebungssystem ist: Röntgen-CT.
  • Die rauschhaltigen Bilder sind: rauschhaltige Projektionsbilder, die Bildauflösung beträgt 1024*720.
  • Die rauschfreien Bilder sind: rauschfreie Projektionsbilder, die Bildauflösung beträgt 1024*720.
    • (1) Berechnen der Entfernungskostenfunktion der Verteilung der rauschfreien Projektionsbilder gemäß der Verteilung der rauschhaltigen Projektionsbilder; Dabei wird das indirekte Verteilungsentfernungsbeschränkungsverfahren verwendet, um die Entfernungskostenfunktion der Verteilung der rauschfreien Projektionsbilder zu schätzen, wobei konkrete Schritte wie folgend sind:
      • (1-1) Durchführen einer iterativen Bayes'schen Entrauschung an einer Gruppe von rauschhaltigen Projektionsbildproben z1,z2,z3,..., zn, die mit diesem Verfahren entrauscht werden sollen, um den Anfangswert der geschätzten rauschfreien Projektionsbildproben Ŝ1,2,3,...n zu erhalten;
      • (1-2) Die in rauschhaltigen Projektionsbilder enthaltenen Räusche sind nicht additiv. Eine exponentielle Transformation und eine verallgemeinerte Anscombe-Transformation für rauschhaltige Projektionsbilder z1,z2,z3,..., zn sind durchgeführt, um rauschhaltige Anscombe-Projektionsbildproben zA1,zA2,zA3,..., zAn zu erhalten: z A i = 2 γ × γ I 0 exp ( z i ) + 0.375 γ 2 + σ e 2 , i = 1 : n ,
        Figure DE102021103293A1_0032
        wobei I0 ein Röntgen-CT-Signalwert ist, wenn kein Signal vorhanden ist, wobei γ ein Beitragswert jedes vom Röntgen-CT-Detektor empfangenen Photons zum elektrischen Signal ist, und σ e 2
        Figure DE102021103293A1_0033
        eine elektronische Rauschvarianz ist. I0, γ und σ e 2
        Figure DE102021103293A1_0034
        sind jeweils durch tatsächliche Experimente gemessen.
      • (1-3) Einstellen des Zufallsvektors der in den rauschhaltigen Anscombe-Projektionsbildproben zA1,zA2,zA3,..., zAn enthaltenen additiven Räusche als N . Jedes Pixel von N entspricht der Standard-Gauß-Verteilung gemäß den Gesetzen der Physik, so dass die Wahrscheinlichkeitsdichtefunktion ƒN eine Wahrscheinlichkeitsdichtefunktion von 1024*720-dimensionaler unabhängiger Standard-Gauß-Verteilungs ist. Und Aufnahme von Simulationsrauschproben n 1 * , n 2 * , n 3 * , , n n *
        Figure DE102021103293A1_0035
        mit der gleichen Anzahl von der rauschhaltigen Bildproben aus der Verteilung des Zufallsvektors N für additive Räusche;
      • (1-4) Eine exponentielle Transformation und eine verallgemeinerte Anscombe-Transformation für geschätzten rauschfreien Projektionsbildblockproben ŝ1,2,3,...ŝn sind durchgeführt, um rauschfreie Anscombe-Projektionsbildproben ŝA1,ŝA2,ŝA3,...,ŝ̂An zu erhalten. Addieren der Rauschproben n 1 * , n 2 * , n 3 * , , n n *
        Figure DE102021103293A1_0036
        mit der geschätzten rauschfreien Anscombe-Projektionsbildproben ŝA1,ŝA2,ŝA3,..., ŝAn , um geschätzte rauschhaltige Anscombe-Projektionsbildproben ẑA1,A2,A3,...Anzu erhalten, wobei ẑAi =ŝAi +ni ,i=1:n;
      • (1-5) Aufzeichnen eines Satzes von rauschhaltigen Anscombe-Projektionsbildproben zA1' zA2, zA3,..., zAn als zA und eines Satzes von geschätzten rauschhaltigen Anscombe-Projektionsbildproben ẑA1,A2,ẑA3,...,ẑAn als ẑA . Berechnen eines Ausdrucks einer Verteilungsentfernung D(ẐA,ẐA) zwischen dem Satz ẑA von diskreten geschätzten rauschhaltigen Anscombe-Projektionsbildproben und dem Satz ẑA der wahren diskreten rauschhaltigen Anscombe-Projektionsbildproben mittels einer Entfernungsfunktion D für ein diskriminierendes neuronales Netzwerk:
        • Einstellen der Verteilungsentfernungsfunktion D so, dass sie ein diskriminantes neuronales Netzwerk ΩD als Verteilungsentfernungsoperator enthält, und eine Parametergruppe des diskriminanten neuronalen Netzwerks ΩD ist θD;
        • Eingeben des Satzes ẑA von rauschhaltigen Anscombe-Projektionsbildproben zA1,zA2,zA3,..., zAn und des Satzes ẑA von geschätzten rauschhaltigen Anscombe-Projektionsbildproben ẑA1,ẑA2,ẑA3,...An in die Verteilungsentfernungsfunktion D, und Verwenden des diskriminanten neuronalen Netzwerks ΩD in Schritt (1) zur Berechnung eines Kreuzentropie-Abstandes zwischen dem Satz ZA und dem Satz ẐA : D ( Z ^ A , Z A ) = i = 1 n [ log ( Ω D ( z Ai ; θ D ) ) + log ( 1 Ω D ( z ^ Ai ; θ D ) ] .
          Figure DE102021103293A1_0037
    • (2) Durchführen einer Bayes'schen Entrauschung an den rauschhaltigen Projektionsbilder gemäß der Entfernungskostenfunktion der Verteilung der rauschfreien Projektionsbilder:
      • (2-1) Für die Gruppe von rauschhaltigen Projektionsbildproben z1-z2,z3,..., zn, die mit diesem Verfahren entrauscht werden sollen, und entsprechende geschätzte rauschfreie Projektionsbildproben ŝ1,ŝ2,ŝ3...,ŝn. wird eine Bayes'sche posteriore Verteilungskostenfunktion Ψ konstruiert: ψ ( s ^ 1 , s ^ 2 , s ^ 3 , , s ^ n ) = i = l : n φ f ( z i s ^ i ) + λ D ,
        Figure DE102021103293A1_0038
        wobei Ψ die Bayes'sche posteriore Verteilungskostenfunktion ist, wobei bei der Bayes'schen posterioren Verteilungskostenfunktion der erste Term φr ein Datentreue-Term zwischen z1,z2,z3,..., zn und ŝ1,2,3,...ŝn ist, und der zweite Term D also D(ẐA,ZA) ist, der λ̂ ein Hyperparameter ist, durch den eine Intensität des Verteilungsabstands D(ẐA,ZA) eingestellt wird;
      • (2-2) Bei der Bayes'schen posterioren Verteilungskostenfunktion Ψ wird eine Optimierungsmethode für neuronale Netze verwendet, und die endgültigen geschätzten rauschfreien Projektionsbildproben ŝ1,2,3,...ŝn werden erhalten, um eine Bildentrauschung zu realisieren.
      • (2-3) Ein bildentrauschendes neuronales Netzwerk ΩF wird als Bildentrauschungsoperator ŝiF(ziE), i = 1: n konstruiert, wobei θE eine Parametergruppe des bildentrauschenden neuronalen Netzwerks ist; Die Bayes'sche posteriore Verteilungskostenfunktion Ψ wird als Trainingskostenfunktion des bildentrauschenden neuronalen Netzwerks verwendet, und eine große Anzahl von rauschhaltigen Projektionsbildprobendatensätzen wird verwendet, um alternierend das bildentrauschende neuronale Netzwerk ΩF und das diskriminante neuronale Netzwerk ΩD in der Entfernungskostenfunktion der Verteilung der rauschfreien Projektionsbilder zu trainieren: min θ E max θ D i = l : n φ f ( z i s ^ i ) + λ D .
        Figure DE102021103293A1_0039
        Die rauschhaltigen Projektionsbildproben werden als Eingabe des bildrauschenden neuronalen Netzwerks verwendet, und die geschätzten rauschfreien Projektionsbilder werden ausgegeben und erhalten. Diese rauschfreien Projektionsbilder sind also die Bilder, die nach der Optimierung und der Minimierung des neuronalen Netzwerks erhalten werden. Das bildentrauschende neuronale Netzwerk, das diskriminante neuronale Netzwerk und Kostenfunktionen der beiden sind wie in 3 gezeigt, wobei eine spezifische Struktur des bildentrauschenden neuronalen Netzwerks wie in 2 (a) und eine spezifische Struktur des diskriminanten neuronalen Netzwerks wie in 2 (b) gezeigt ist.

Claims (19)

  1. Ein Bayes'sches Bildentrauschungsverfahren, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert, dadurch gekennzeichnet, daß das Verfahren folgende Schritte umfasst: (1) Berechnen einer Entfernungskostenfunktion einer Verteilung von rauschfreien Bilder gemäß einer Verteilung von rauschhaltigen Bilder; (2) Durchführen einer Bayes'schen Entrauschung an den rauschhaltigen Bilder gemäß der Entfernungskostenfunktion der Verteilung der rauschfreien Bilder in Schritt (1).
  2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daß in Schritt (1) die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder gemäß der Verteilung der rauschhaltigen Bilder berechnet wird, umfassend folgende Schritte: Für Bilder, deren Verteilung genau oder annähernd analytisch ausgedrückt werden kann, wird ein auf einer charakteristischen Funktion basierende Verfahren verwendet, um die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder zu schätzen, einschließlich der folgenden Schritte: (1) Aufzeichnen eines vom Bildgebungssystem erzeugten rauschhaltigen Bild-Zufallsvektors als Z, Schätzen einer Wahrscheinlichkeitsdichtefunktion fz des rauschhaltigen Bild-Zufallsvektors Z in analytischer Form, und Aufzeichnen eines additiven Rausch-Zufallsvektors, der in dem rauschhaltigen Bild-Zufallsvektor Z enthalten ist, als N , Schätzen einer Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N in analytischer Form; (2) Durchführen jeweils einer Fourier-Transformation der Wahrscheinlichkeitsdichtefunktionen ƒz und ƒN in Schritt (1), um eine charakteristische Funktion Fz des rauschhaltigen Bild-Zufallsvektors und eine charakteristische Funktion FN des additiven Rausch-Zufallsvektors zu erhalten; (3) Aufzeichnen eines vom Bildgebungssystem erzeugten rauschfreien Bild-Zufallsvektors als s, und Dividieren der charakteristischen Funktion Fz des rauschhaltigen Bild-Zufallsvektors durch die charakteristische Funktion FN des additiven Rausch-Zufallsvektors in Schritt (2), um eine charakteristische Funktion Fs des rauschfreien Bild-Zufallsvektors s zu erhalten: F S = F Z / F N ,
    Figure DE102021103293A1_0040
    oder: Durch ein Zieloptimierungsverfahren wird ein geschätzter Wert Fs* der charakteristischen Funktion Fs des Zufallsvektors der rauschfreien Bilder gemäß einer Konsistenz von Fs mit FZ/FN erhalten: F S * = argmin F S Φ ( F S , F Z / F N ) ,
    Figure DE102021103293A1_0041
    wobei Φ(g,g) ein Abstandsmaß zwischen den charakteristischen Funktionen zweier Zufallsvektoren ist, und argmin einen Parameter darstellt, bei dem eine Zielfunktion einen Minimalwert annimmt; (4) Durchführen einer inversen Fourier-Transformation für die charakteristische Funktion Fs oder den geschätzten Wert Fs* der charakteristischen Funktion Fs des Zufallsvektors der rauschfreien Bilder in Schritt (3), um eine Wahrscheinlichkeitsdichtefunktion ƒs des Zufallsvektors der rauschfreien Bilder zu erhalten, d.h. eine Verteilung des Zufallsvektors S der rauschfreien Bilder zu erhalten; (5) Durchführen einer analytischen Filterentrauschung oder einer iterativen Bayes'schen Entrauschung an einer Gruppe von rauschhaltigen Bildproben z1,z2,z3,..., zn, um einen Anfangswert von geschätzten rauschfreien Bildproben Ŝ1,2,3,...n zu erhalten, und Aufzeichnen eines Satzes von ŝ1,2,3,...n als ŝ; (6) Berechnen einer Verteilungsentfernung D(Ŝ,S) zwischen dem Satz ŝ der geschätzten rauschfreien Bildproben Ŝ1,2,3,...n in Schritt (5) und dem in Schritt (4) erhaltenen Zufallsvektor s der rauschfreien Bilder gemäß einer Verteilungsentfernungsfunktion D, wobei ein Ausdruck der Verteilungsentfernung D(Ŝ,S) die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder ist.
  3. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daß in Schritt (1) die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder gemäß der Verteilung der rauschhaltigen Bilder berechnet wird, umfassend folgende Schritte: Für den Fall, dass die Verteilung schwer analytisch auszudrücken ist, oder für Bilder, deren Verteilung bei Berechnung zu kompliziert ist, wird ein indirektes Verfahren zur Beschränkung einer Verteilungsentfernung verwendet, um die Entfernungskostenfunktion der Verteilung der rauschfreien Bilder zu schätzen, einschließlich der folgenden Schritte: (1) Durchführen einer analytischen Filterentrauschung oder einer iterativen Bayes'schen Entrauschung an einer Gruppe von rauschhaltigen Bildproben Z1,Z2,Z3,..., Zn, um einen Anfangswert von geschätzten rauschfreien Bildproben Ŝ1,2,3,...n zu erhalten; (2) Einstellen eines Zufallsvektors von in den rauschhaltigen Bildproben z1,z2,z3,... zn in Schritt (1) enthaltenen additiven Räusche als N , Schätzen der Wahrscheinlichkeitsdichtefunktion ƒN des Zufallsvektors N für additive Räusche in analytischer Form, und Aufnahme von Rauschproben n 1 * , n 2 * , n 3 * , , n n *
    Figure DE102021103293A1_0042
    mit einer gleichen Anzahl von der rauschhaltigen Bildproben aus der Verteilung des Zufallsvektors N für additive Räusche; oder: Aufnahme von Rauschproben des Bildgebungssystems, wenn kein Signal zu erhalten ist, um n 1 * , n 2 * , n 3 * , , n n *
    Figure DE102021103293A1_0043
    zu erhalten; (3) Addieren der Rauschproben n 1 * , n 2 * , n 3 * , , n n *
    Figure DE102021103293A1_0044
    in Schritt (2) mit der geschätzten rauschfreien Bildproben Ŝ1,2,3,...n in Schritt (1), um geschätzte rauschhaltige Bildproben ẑ1,2,ẑ3,...ẑn zu erhalten, wobei zi=si+ni, i=1:n; (4) Aufzeichnen eines Satzes der rauschhaltigen Bildproben z1,z2,z3,..., Zn in Schritt (1) als Z und eines Satzes der geschätzten rauschhaltigen Bildproben ẑ1,z2,z3,...ẑn in Schritt (3) als Ẑ, Berechnen einer Verteilungsentfernung D(Ẑ,Z) zwischen dem Satz Ẑ von diskreten geschätzten rauschhaltigen Bildproben und dem Satz Z der wahren diskreten rauschhaltigen Bildproben gemäß einer Verteilungsentfernungsfunktion D, wobei ein Ausdruck der Verteilungsentfernung D(Ẑ,Z) eine Entfernungskostenfunktion der Verteilung der rauschfreien Bilder ist.
  4. Verfahren nach Anspruch 2 oder 3, dadurch gekennzeichnet, daß die Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N in analytischer Form derart geschätzt ist: Gemäß den physikalischen Gesetzen des Bildgebungssystems wird eine Verteilungsform des additiven Rausch-Zufallsvektors N erhalten, und dann werden Verteilungsparameter additiven Rausch-Zufallsvektors N durch Simulation oder tatsächliche Experimente gemessen, so dass die Wahrscheinlichkeitsdichtefunktion fN des additiven Rausch-Zufallsvektors N erhalten wird.
  5. Verfahren nach Anspruch 2 oder 3, dadurch gekennzeichnet, daß die Wahrscheinlichkeitsdichtefunktion fN des additiven Rausch-Zufallsvektors N in analytischer Form derart geschätzt ist: Eine gemischte Gaußsche Verteilung zum Anpassen an eine unbekannte Verteilung wird verwendet, Nehmen zunächst an, dass der Zufallsvektor N für additive Räusche der gemischten Gaußschen Verteilung entspricht, und dann die gemischten Gaußschen Verteilungsparameter durch Simulation oder tatsächliche Experimente gemessen werden, so dass die Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N erhalten wird.
  6. Verfahren nach Anspruch 2 oder 3, dadurch gekennzeichnet, daß das Verfahren zum Schätzen der Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N in analytischer Form lautet: Für nichtadditive Rauschverteilung kann sie durch eine spezifische mathematische Transformation in ungefähre additive Räusche umgewandelt werden, und dann die Wahrscheinlichkeitsdichtefunktion ƒN des additiven Rausch-Zufallsvektors N gemäß dem Verfahren nach Anspruch 4 oder 5 erhalten wird.
  7. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daß das Verfahren zum Schätzen der Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors z in analytischer Form lautet: die Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors z durch Zählen von Frequenzen einer großen Anzahl von rauschhaltigen Bildproben erhalten wird.
  8. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daß das Verfahren zum Schätzen der Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors z in analytischer Form lautet: Verteilungsmodellannahmen im Fachbereich verwendet werden, um eine Verteilungsform des rauschhaltigen Bild-Zufallsvektors z zu erhalten, wobei Verteilungsparameter des rauschhaltigen Bild-Zufallsvektors z durch Simulation oder tatsächliche Experimente gemessen werden, so dass die Wahrscheinlichkeitsdichtefunktion ƒz des rauschhaltigen Bild-Zufallsvektors z geschätzt wird.
  9. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daß „Durchführen einer Bayes'schen Entrauschung an den rauschhaltigen Bilder gemäß der Entfernungskostenfunktion der Verteilung der rauschfreien Bilder“ folgende Schritte umfasst: (1) Für rauschhaltige Bildproben z1,z2,z3,..., Zn und geschätzte rauschfreie Bildproben ŝ1,ŝ2,ŝ3,...,ŝn wird eine Bayes'sche posteriore Verteilungskostenfunktion Ψkonstruiert: ψ ( s ^ 1 , s ^ 2 , s ^ 3 , , s ^ n ) = i = l : n φ f ( z i s ^ i ) + λ D ,
    Figure DE102021103293A1_0045
    wobei Ψ die Bayes'sche posteriore Verteilungskostenfunktion ist, wobei bei der Bayes'schen posterioren Verteilungskostenfunktion der erste Term φr ein Datentreue-Term zwischen z1 ,z2,z3,..., zn und ŝ1,2,3,...ŝn ist, und der zweite Term D ein Verteilungsabstand D(ŝ,s) oder ein Verteilungsabstand D(Ẑ,Z) ist, der λ̂ ein Hyperparameter ist, durch den eine Intensität des Verteilungsabstands D(Ŝ,S) eingestellt wird; (2) Die Bayes'sche posteriore Verteilungskostenfunktion Ψ von Schritt (1) wird minimiert und optimiert, und die endgültigen geschätzten rauschfreien Bilder ŝ1,ŝ2,ŝ3,...,ŝn werden erhalten, um eine Bildentrauschung zu realisieren.
  10. Verfahren nach Anspruch 9, dadurch gekennzeichnet, daß das Verfahren zur Minimierung und Optimierung der Bayes'schen posterioren Verteilungskostenfunktion Ψ in Schritt (2) ein übliches klassisches mathematisches Optimierungsverfahren ist.
  11. Verfahren nach Anspruch 10, dadurch gekennzeichnet, daß das Verfahren zur Minimierung und Optimierung der Bayes'schen posterioren Verteilungskostenfunktion Ψ ein übliches klassisches mathematisches Optimierungsverfahren ist, und ein von Lagrange-Multiplikatorverfahren, alternierendem Iterationsverfahren, alternativer konjugierter Methode der kleinsten Fehler-Quadrate oder Gradientenoptimierungsalgorithmus umfasst.
  12. Verfahren nach Anspruch 9, dadurch gekennzeichnet, daß das Verfahren zur Minimierung und Optimierung der Bayes'schen posterioren Verteilungskostenfunktion Ψ eine Optimierungsmethode für neuronale Netze ist, die folgende Schritte umfasst: (1) Ein bildentrauschendes neuronales Netzwerk ΩF, wird als Bildentrauschungsoperator ŝiF(ziE), i=1:n konstruiert, wobei θE eine Parametergruppe des bildentrauschenden neuronalen Netzwerks ist; (2) Die Bayes'sche posteriore Verteilungskostenfunktion Ψ wird als Trainingskostenfunktion des bildentrauschenden neuronalen Netzwerks verwendet, und eine große Anzahl von rauschhaltigen Bildprobendatensätzen wird verwendet, um das bildentrauschende neuronale Netzwerk zu trainieren; (3) Die rauschhaltigen Bildproben werden als Eingabe des bildrauschenden neuronalen Netzwerks verwendet, und die geschätzten rauschfreien Bilder werden ausgegeben und erhalten, wobei diese rauschfreien Bilder also die Bilder sind, die nach der Optimierung und der Minimierung des neuronalen Netzwerks erhalten werden.
  13. Verfahren nach Anspruch 2 oder 3, dadurch gekennzeichnet, daß die Verteilungsentfernungsfunktion D in Schritt (6) eine analytische dominante Verteilungsentfernungsfunktion ist.
  14. Verfahren nach Anspruch 13, dadurch gekennzeichnet, daß die analytische dominante Verteilungsentfernungsfunktion eine K-L-Divergenzdominanzverteilungsentfernungsfunktion oder eine J-S-Divergenzdominanzverteilungsentfernungsfunktion umfasst.
  15. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daß das Verfahren zum Berechnen und Erhalten einer Funktion einer Verteilungsentfernung D(ŝ,ŝ) in Schritt (5) gemäß einer Verteilungsentfernungsfunktion D ein Verfahren zum Identifizieren der Funktion der Verteilungsentfernung des neuronalen Netzes ist, einschließlich der folgenden Schritte: (1) Einstellen der Verteilungsentfernungsfunktion D so, dass sie ein diskriminantes neuronales Netzwerk ΩD als Verteilungsentfernungsoperator enthält, und eine Parametergruppe des diskriminanten neuronalen Netzwerks ΩD ist θD; (2) Eingeben des in Schritt (5) von Anspruch 2 beschriebenen Satzes ŝ der geschätzten rauschfreien Bildproben ŝ1,ŝ2,ŝ3,...ŝn und des in Schritt (4) von Anspruch 2 erhaltenen Zufallsvektors s der rauschfreien Bilder in die Verteilungsentfernungsfunktion D, und Verwenden des diskriminanten neuronalen Netzwerks ΩD in Schritt (1) zur Berechnung einer Funktion eines Kreuzentropie-Abstandes zwischen dem Satz ŝ und dem Zufallsvektor s der rauschfreien Bilder: D ( S ^ , S ) = i = l : n log ( Ω D ( s i ; θ D ) ) + log ( 1 Ω D ( s i ; θ D )
    Figure DE102021103293A1_0046
    oder einer Funktion eines Wasserstein-Abstandes: D ( S ^ , S ) = i = l : n | Ω D ( s i ; θ D ) Ω D ( s i ; θ D ) | ;
    Figure DE102021103293A1_0047
    wobei die obigen s1,s2,s3,....,sn n Abtastproben des Zufallsvektors s der rauschfreien Bilder sind.
  16. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daß das Verfahren zum Berechnen und Erhalten einer Funktion einer Verteilungsentfernung D(Ŝ,S) in Schritt (5) gemäß einer Verteilungsentfernungsfunktion D ein Verfahren zum Identifizieren der Funktion der Verteilungsentfernung des neuronalen Netzes ist, einschließlich der folgenden Schritte: (1) Einstellen der Verteilungsentfernungsfunktion D so, dass sie ein diskriminantes neuronales Netzwerk ΩD als Verteilungsentfernungsoperator enthält, und eine Parametergruppe des diskriminanten neuronalen Netzwerks ΩD ist θD; (2) Eingeben des in Schritt (5) von Anspruch 2 beschriebenen Satzes ŝ der geschätzten rauschfreien Bildproben ŝ1,2,3,...ŝn und des in Schritt (4) von Anspruch 2 erhaltenen Zufallsvektors s der rauschfreien Bilder in die Verteilungsentfernungsfunktion D, und Verwenden des diskriminanten neuronalen Netzwerks ΩD in Schritt (1), um eine Diskriminanznetzwerk-Verteilungsentfernungskostenfunktion von generativen gegnerischen Netzwerken Methoden im Fachbereich zu berechnen.
  17. Verfahren nach Anspruch 3, dadurch gekennzeichnet, daß das Verfahren zum Berechnen einer Funktion einer Verteilungsentfernung D(Ẑ,Z) in Schritt (6) gemäß einer Verteilungsentfernungsfunktion D ein Verfahren zum Identifizieren der Funktion der Verteilungsentfernung des neuronalen Netzes ist, einschließlich der folgenden Schritte: (1) Einstellen der Verteilungsentfernungsfunktion D so, dass sie ein diskriminantes neuronales Netzwerk \ΩD als Verteilungsentfernungsoperator enthält, und eine Parametergruppe des diskriminanten neuronalen Netzwerks ΩD ist θD; (2) Eingeben des Satzes Ẑ der diskreten geschätzten rauschhaltigen Bildproben und dem Satz Z der wahren diskreten rauschhaltigen Bildproben in die Verteilungsentfernungsfunktion D, und Verwenden des diskriminanten neuronalen Netzwerks ΩD in Schritt (1) zur Berechnung einer Funktion eines Kreuzentropie-Abstandes zwischen dem Satz Ẑ der diskreten geschätzten rauschhaltigen Bildproben und dem Satz Z der wahren diskreten rauschhaltigen Bildproben: D ( Z ^ , Z ) = i = 1 : n log ( Ω D ( z i ; θ D ) ) + log ( 1 Ω D ( z ^ i ; θ D )
    Figure DE102021103293A1_0048
    oder einer Funktion eines Wasserstein-Abstandes: D ( Z ^ , Z ) = i = l : n | Ω D ( z i ; θ D ) Ω D ( z ^ i ; θ D ) | .
    Figure DE102021103293A1_0049
  18. Verfahren nach Anspruch 3, dadurch gekennzeichnet, daß das Verfahren zum Berechnen einer Verteilungsentfernung D(Ẑ,Z) in Schritt (6) gemäß einer Verteilungsentfernungsfunktion D ein Verfahren zum Identifizieren der Funktion der Verteilungsentfernung des neuronalen Netzes ist, einschließlich der folgenden Schritte: (1) Einstellen der Verteilungsentfernungsfunktion D so, dass sie ein diskriminantes neuronales Netzwerk ΩD als Verteilungsentfernungsoperator enthält, und eine Parametergruppe des diskriminanten neuronalen Netzwerks ΩD ist θD; (2) Eingeben des Satzes Ẑ der diskreten geschätzten rauschhaltigen Bildproben und dem Satz Z der wahren diskreten rauschhaltigen Bildproben in die Verteilungsentfernungsfunktion D, und Verwenden des diskriminanten neuronalen Netzwerks ΩD in Schritt (1), um eine Diskriminanznetzwerk-Verteilungsentfernungskostenfunktion von generativen gegnerischen Netzwerken Methoden im Fachbereich zu berechnen.
  19. Verfahren nach Anspruch 16 oder 18, dadurch gekennzeichnet, daß die genannten generativen gegnerischen Netzwerken Methoden im Fachbereich in Schritt (2) eines von SGAN, RSGAN, RaSGAN, LSGAN, RaLSGAN, HingeGAN, WGAN-GP, RSGAN-GP oder RaSGAN-GP sind.
DE102021103293.8A 2020-07-20 2021-02-11 Ein Bayes'sches Bildentrauschungsverfahren, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert Pending DE102021103293A1 (de)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202010696129.6A CN111932468A (zh) 2020-07-20 2020-07-20 一种基于含噪声图像分布约束的贝叶斯图像去噪方法
CN202010696129.6 2020-07-20

Publications (1)

Publication Number Publication Date
DE102021103293A1 true DE102021103293A1 (de) 2022-01-20

Family

ID=73312864

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102021103293.8A Pending DE102021103293A1 (de) 2020-07-20 2021-02-11 Ein Bayes'sches Bildentrauschungsverfahren, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert

Country Status (3)

Country Link
US (1) US11880903B2 (de)
CN (1) CN111932468A (de)
DE (1) DE102021103293A1 (de)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114529480A (zh) * 2022-03-02 2022-05-24 北京大学 一种基于随机重配对的无监督一锅式多帧图像去噪方法
CN115308215A (zh) * 2022-10-11 2022-11-08 南通市怡天时纺织有限公司 基于激光束的织物织造疵点检测方法
CN117669335A (zh) * 2024-02-01 2024-03-08 北京国星创图科技有限公司 一种基于太空环境数据挖掘方法及系统

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112767271B (zh) * 2021-01-19 2023-09-22 西安交通大学 一种基于三维变分网络的高光谱图像深度降噪的方法
CN115482162B (zh) * 2022-09-02 2023-10-24 华南理工大学 一种基于随机重排和无标签模型的隐式图像盲去噪方法
CN115859007B (zh) * 2023-02-17 2023-05-12 广东石油化工学院 一种石化仪表采样数据滑窗约束容错滤波降噪方法和装置
CN116339128B (zh) * 2023-05-30 2023-07-28 北京国电富通科技发展有限责任公司 一种输渣系统落渣雨量在线监控方法及系统
CN116451598B (zh) * 2023-06-20 2023-09-05 杭州经纬信息技术股份有限公司 基于去噪扩散概率模型的太阳辐照度预测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10726525B2 (en) * 2017-09-26 2020-07-28 Samsung Electronics Co., Ltd. Image denoising neural network architecture and method of training the same

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114529480A (zh) * 2022-03-02 2022-05-24 北京大学 一种基于随机重配对的无监督一锅式多帧图像去噪方法
CN115308215A (zh) * 2022-10-11 2022-11-08 南通市怡天时纺织有限公司 基于激光束的织物织造疵点检测方法
CN117669335A (zh) * 2024-02-01 2024-03-08 北京国星创图科技有限公司 一种基于太空环境数据挖掘方法及系统
CN117669335B (zh) * 2024-02-01 2024-05-07 北京国星创图科技有限公司 一种太空环境数据挖掘方法及系统

Also Published As

Publication number Publication date
US11880903B2 (en) 2024-01-23
CN111932468A (zh) 2020-11-13
US20220020123A1 (en) 2022-01-20

Similar Documents

Publication Publication Date Title
DE102021103293A1 (de) Ein Bayes'sches Bildentrauschungsverfahren, das auf Verteilungsbeschränkungen von rauschhaltigen Bildern basiert
DE112017002799B4 (de) Verfahren und system zum generieren multimodaler digitaler bilder
EP3336800B1 (de) Bestimmen einer trainingsfunktion zum generieren von annotierten trainingsbildern
DE112009004059B4 (de) Verfahren zum Entfernen von Unschärfe von einem Bild und Aufzeichnungsmedium, auf dem das Verfahren aufgezeichnet ist
DE60020795T2 (de) Echtzeitmodellierung des Verhaltens menschlicher Sehvermögenssysteme
DE102020125207A1 (de) Unbeaufsichtigtes training von neuronalem netz zur bildkompression mit hohem dynamikumfang
DE102019207238A1 (de) Bereitstellen eines Ergebnisbilddatensatzes
DE112018002822T5 (de) Klassifizieren neuronaler netze
DE19827034A1 (de) Iteratives Filtersystem für medizinische Bilder
DE112016005905T5 (de) Verfahren und System zum Verschmelzen erfasster Messungen
DE602004001366T2 (de) Gewichtete, überbestimmte Rauschunterdrückung
DE112019007393T5 (de) Verfahren und System zum Trainieren eines Modells zur Bilderzeugung
DE102020134755A1 (de) Verarbeiten von videoframes über neuronales faltungsnetzwerk durch verwenden von statistiken vorhergehender frames
DE60115056T2 (de) Bildverbesserung
DE102022113243A1 (de) Score-basierte generative Modellierung im latenten Raum
Raj et al. Denoising of medical images using total variational method
DE102005045179A1 (de) Bildverbindung auf der Grundlage unabhängiger Rauschbedingungen
DE602005005348T2 (de) Konforme segmentierung von organen auf medizinischen bildern
Wang et al. A modified higher-order singular value decomposition framework with adaptive multilinear tensor rank approximation for three-dimensional magnetic resonance Rician noise removal
Chen et al. Enhancement and denoising method for low-quality MRI, CT images via the sequence decomposition Retinex model, and haze removal algorithm
Xu et al. Multiplicative decomposition based image contrast enhancement method using PCNN factoring model
Radhika et al. An adaptive optimum weighted mean filter and bilateral filter for noise removal in cardiac MRI images
DE102021124428A1 (de) Trainieren energiebasierter variational autoencoder
DE102021114064A1 (de) Rauschadaptive Nicht-Blind-Bildschärfung
DE102020210775A1 (de) Magnetresonanztomographie-Rekonstruktion durch Verwenden von maschinellem Lernen und Bewegungskompensation

Legal Events

Date Code Title Description
R012 Request for examination validly filed