DE102005032287A1 - Two-dimensional initial image processing method for diagnostic application, involves setting specific factor on value above minimal valve in case of affirmation and on minimal value in case of negation - Google Patents
Two-dimensional initial image processing method for diagnostic application, involves setting specific factor on value above minimal valve in case of affirmation and on minimal value in case of negation Download PDFInfo
- Publication number
- DE102005032287A1 DE102005032287A1 DE200510032287 DE102005032287A DE102005032287A1 DE 102005032287 A1 DE102005032287 A1 DE 102005032287A1 DE 200510032287 DE200510032287 DE 200510032287 DE 102005032287 A DE102005032287 A DE 102005032287A DE 102005032287 A1 DE102005032287 A1 DE 102005032287A1
- Authority
- DE
- Germany
- Prior art keywords
- computer
- image
- processing method
- pixel
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration by the use of local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
Description
Die vorliegende Erfindung betrifft ein Bearbeitungsverfahren für ein zweidimensionales Anfangsbild. Die vorliegende Erfindung betrifft weiterhin einen Datenträger mit einem auf dem Datenträger gespeicherten Computerprogramm für einen Rechner zur Durchführung eines derartigen Bearbeitungsverfahrens. Schließlich betrifft die vorliegende Erfindung noch einen Rechner mit einem solchen Datenträger.The The present invention relates to a processing method for a two-dimensional one Initial image. The present invention further relates to a disk with a stored on the disk Computer program for a computer to carry out such a processing method. Finally, the present concerns Invention still a computer with such a disk.
Anfangsbilder im Sinne der vorliegenden Erfindung können z.B. zweidimensionale Bildsignale sein, die von Röntgen-Durchstrahlungsanlagen geliefert werden. Es handelt sich hierbei in der Regel um diskrete Bildsignale. Die Bildsignale können Rohsignale oder aber bereits teilweise aufbereitete Signale sein. Sie enthalten oftmals medizinische Bedeutungsinhalte und werden insbesondere zur Diagnose und teilweise auch im Rahmen der Therapie bzw. Intervention verwendet. Sie weisen in der Regel eine hohe Dynamik auf. Weiterhin beinhalten die Bilder in vielen Fällen feine Details, die relevante Bildinformation enthalten, z.B. kleine Gefäße, die mit Kontrastmittel gefüllt sind, oder Führungsdrähte. In vielen Fällen kann der behandelnde Arzt jedoch nicht oder nur schlecht Abbildungsparameter und Bildaufbereitungsparameter derart einstellen, dass das Bild gut aufbereitet wird. Dennoch soll vor al-lem in der szenischen Bildgebung, also wenn eine Sequenz von Anfangsbildern vorliegt, das dargestellte Bild auch auf Anhieb alle relevanten Details zeigen.initial images For the purposes of the present invention, e.g. two-dimensional Image signals are from X-ray transmission systems to be delivered. These are usually discrete Image signals. The image signals can Be raw signals or already partially processed signals. They often contain medical meaning content and become in particular for diagnosis and partly also in the context of therapy or intervention used. They usually have high dynamics on. Furthermore, in many cases the images contain fine details, the relevant ones Image information, e.g. small vessels containing contrast media filled are, or guidewires. In many cases However, the attending physician may not or only poorly mapping parameters and adjust rendering parameters such that the image is treated well. Nevertheless, especially in scenic imaging, ie if there is a sequence of initial images, the one shown Picture also immediately show all relevant details.
Um eine allzu große Dynamik in den dargestellten Bildern zu vermeiden, werden im Stand der Technik oftmals Harmonisierungen mit einem vorgegebenen festen Filterkern vorgenommen, der im Wesentlichen einem Hochpassfilter entspricht. Ein wesentlicher Nachteil von derartigen Hochpassfiltern ist jedoch die Ausbildung von Überschwingern an Kanten, die den Arzt im Ex tremfall sogar zu Fehldiagnosen veranlassen können, in der Regel aber zumindest einen ungewohnten Bildeindruck bewirken.Around too big To avoid dynamics in the images shown, are in the state The technique often harmonizes with a given fixed Filter core made essentially a high-pass filter equivalent. A major disadvantage of such high-pass filters is however the education of overshooters On edges that cause the doctor in the extreme case even to misdiagnosis can, but usually cause at least an unfamiliar image impression.
Weiterhin ist der Rauscheindruck des Bildes für den Betrachter von Bedeutung. Dies gilt vor allem in der röntgenbasierten Bildgebung. Denn wenn das Nutzsignal und das Rauschen nicht voneinander getrennt werden können, ist die Wahl für den richtigen Kontrast sehr schwierig. Ist der Kontrast zu groß, wird das Auge des Betrachters zu stark durch das Rauschen irritiert. Ist der Kontrast zu gering, ist das Signal kaum noch oder sogar gar nicht mehr zu erkennen.Farther the noise impression of the picture is of importance to the viewer. This is especially true in the X-ray based Imaging. Because if the useful signal and the noise are not from each other can be separated is the choice for the right contrast very difficult. If the contrast is too big, it will Eye of the viewer too much irritated by the noise. is the contrast is too low, the signal is barely or even even no longer recognizable.
Zur
Bildaufbereitung ist – beispielsweise
aus der
- – Der Rechner zerlegt das Anfangsbild in Teilbilder und ein Restbild. Jedes Teilbild enthält dabei diejenigen Anteile des zweidimensionalen Anfangsbildes, die mit Frequenzen örtlich variieren, die in einem vorbestimmten, für das jeweilige Teilbild charakteristischen, die Frequenz Null nicht enthaltenden Teilbildfrequenzbereich liegen. Das Restbild enthält zumindest einen Gleichanteil des zweidimensionalen Anfangsbildes, der örtlich invariabel ist.
- – Für jedes Pixel jedes Teilbildes ermittelt der Rechner einen für das jeweilige Pixel des jeweiligen Teilbildes spezifischen Faktor, mit dem er dieses Pixel multipliziert. Auf diese Weise ermittelt der Rechner bearbeitete Teilbilder.
- – Weiterhin ermittelt der Rechner anhand des Restbildes ein bearbeitetes Restbild.
- – Die bearbeiteten Teilbilder und das bearbeitete Restbild kombiniert der Rechner zu einem Endbild.
- – Wenn für das Pixel eine vorbestimmte Bedingung erfüllt ist, setzt der Rechner den Faktor auf einen Wert oberhalb eines Minimalwertes. Anderenfalls setzt er den Faktor auf den Minimalwert. Der Minimalwert ist dabei für alle Pixel des jeweiligen Teilbildes der gleiche.
- - The computer decomposes the initial image into partial images and a residual image. In this case, each partial image contains those portions of the two-dimensional initial image that vary locally with frequencies which are in a predetermined partial image frequency range that is characteristic of the respective partial image and does not include the frequency zero. The residual image contains at least a DC component of the two-dimensional initial image, which is locally invariable.
- For each pixel of each partial image, the computer determines a specific factor for the respective pixel of the respective partial image with which it multiplies this pixel. In this way, the computer determines processed partial images.
- - Furthermore, the computer determines a processed residual image based on the residual image.
- - The edited partial images and the processed residual image are combined by the calculator into a final image.
- If the pixel meets a predetermined condition, the calculator sets the factor to a value above a minimum value. Otherwise it sets the factor to the minimum value. The minimum value is the same for all pixels of the respective sub-image.
Bei
der
Die
Vorgehensweise der
Vom Erfinder der vorliegenden Erfindung wurde erkannt, dass die Veränderung des Rauscheindrucks bei dem Verfahren des Standes der Technik ihre Ursache darin hat, dass bei der Ermittlung der bearbeiteten Teilbilder jeweils versucht wird, das Nutzsignal unbeeinflusst zu lassen, das Rauschen dagegen zu dämpfen.from Inventors of the present invention recognized that the change the noise impression in the process of the prior art their The reason for this is that when determining the edited drawing files each attempt is made to leave the useful signal unaffected, the To dampen noise, on the other hand.
Die Aufgabe der vorliegenden Erfindung besteht darin, ein Bearbeitungsverfahren der eingangs genannten Art und die entsprechenden Gegenstände derart auszugestalten, dass die Verbesserung der Erkennbarkeit der relevanten Bildinformation erhalten bleibt, der Rauscheindruck aber als subjektiv erheblich angenehmer empfunden wird.The object of the present invention is to provide a machining method of the initially mentioned type and the corresponding objects in such a way that the improvement of the recognizability of the relevant image information is retained, the noise impression but perceived as subjectively much more pleasant.
Die Aufgabe wird dadurch gelöst, dass der Rechner zum Ermitteln des für das jeweilige Pixel des jeweiligen Teilbildes spezifischen Faktors ermittelt, ob er in dem jeweiligen Teilbild innerhalb eines zweidimensionalen Filterkerns um das jeweilige Pixel herum eine Struktur erkennt, und den Faktor bejahendenfalls auf einen Wert oberhalb des Minimalwertes und verneinendenfalls auf den Minimalwert setzt.The Task is solved by that the calculator for determining the for each pixel of the respective Partial image specific factor determines whether he is in the respective Partial image within a two-dimensional filter core around the respective one Pixels around recognizes a structure, and the factor if so to a value above the minimum value and in the negative case set to the minimum value.
Hiermit korrespondierend wird die Aufgabe durch einen Datenträger gelöst, auf dem ein Computerprogramm für einen Rechner zur Durchführung eines derartigen Bearbeitungsverfahrens gespeichert ist. Auch wird die Aufgabe durch einen Rechner mit einem solchen Datenträger gelöst.Herewith Correspondingly, the task is solved by a disk on a computer program for a computer to carry out of such a processing method is stored. Also will solved the problem by a computer with such a disk.
Zum Erkennen einer Struktur und dem damit korrespondierenden Ermitteln des Faktors ist es insbesondere möglich, dass der Rechner
- – innerhalb des Filterkerns um das jeweilige Pixel des jeweiligen Teilbildes herum die statistischen Werte von mehreren Pixelgruppen bildet, wobei die Pixel jeder Pixelgruppe zumindest im Wesentlichen auf einer Geraden liegen und als Zentralpixel das jeweilige Pixel enthalten, und
- – anhand eines für die Schwankung der statistischen Werte charakteristischen Wertes den Faktor für das jeweilige Pixel des jeweiligen Teilbildes ermittelt.
- The statistical values of a plurality of pixel groups are formed within the filter kernel around the respective pixel of the respective partial image, the pixels of each pixel group lying at least essentially on a straight line and containing the respective pixel as the central pixel, and
- - Determines the factor for the respective pixel of the respective field based on a value characteristic of the fluctuation of the statistical values.
Das Ermitteln des für die Schwankung der statistischen Werte charakteristischen Wertes kann der Rechner beispielsweise anhand mindestens zweier der Werte Maximum, Minimum und Durchschnittswert der statistischen Werte vornehmen.The Determine the for the fluctuation of statistical values of characteristic value For example, the calculator can use at least two of the values Make the maximum, minimum and average of the statistical values.
Als geeignete statistische Werte der Pixelgruppen kommen insbesondere deren Mittelwerte und deren Varianzen in Betracht.When suitable statistical values of the pixel groups come in particular their mean values and their variances are considered.
Vorzugsweise setzt der Rechner den Faktor dann und nur dann auf den Minimalwert, wenn der für die Schwankung der statistischen Werte charakteristische Wert unterhalb eines Schwellwerts liegt. Weiterhin hebt der Rechner den Faktor vorzugsweise linear mit dem für die Schwankung der statistischen Werte charakteristischen Wert an, wenn der für die Schwankung der statistischen Werte charakteristische Wert oberhalb des Schwellwerts liegt.Preferably then the calculator sets the factor to the minimum value if that for the Fluctuation of statistical values characteristic value below a threshold. Furthermore, the calculator raises the factor preferably linear with the for the fluctuation of statistical values characteristic value, if the for the fluctuation of statistical values characteristic value above of the threshold.
Der Schwellwert ist vorzugsweise für alle Teilbilder der gleiche. Er kann vom Rechner anhand des Anfangsbildes und/oder anhand einer Anwendereingabe ermittelt werden. Bei spielsweise kann ein Prozentsatz des durchschnittlichen Helligkeitswertes des Anfangsbildes zur Berechnung des Schwellwertes dienen, wobei der Prozentsatz vom Anwender vorgegeben werden kann.Of the Threshold is preferably for all partial images the same. He can from the computer based on the initial image and / or determined by a user input. For example can be a percentage of the average brightness value of the Initial image used to calculate the threshold, the Percentage can be specified by the user.
Analog zum Schwellwert ist vorzugsweise auch der Minimalwert für alle Teilbilder der gleiche. Alternativ ist aber auch möglich, dass der Minimalwert individuell für das jeweilige Teilbild ist. Auch der Minimalwert bzw. die Minimalwerte kann (können) vorzugsweise dem Rechner von einem Anwender vorgegeben werden.Analogous the threshold is preferably also the minimum value for all fields the same. Alternatively, it is also possible that the minimum value individually for the respective partial image is. Also the minimum value or the minimum values can (preferably) The computer can be specified by a user.
Es ist möglich, dass das Restbild ausschließlich den Gleichanteil enthält. Es ist aber auch möglich, dass das Restbild zusätzlich zum Gleichanteil diejenigen Anteile des zweidimensionalen Anfangsbildes enthält, die mit Frequenzen örtlich variieren, die kleiner als die Frequenzen des niedrigsten Teilbildfrequenzbereichs sind.It is possible, that the residual image is exclusive contains the DC component. But it is also possible that the residual image in addition to the DC component those parts of the two-dimensional initial image contains those with frequencies locally which are smaller than the frequencies of the lowest field frequency range are.
Der Rechner ermittelt das bearbeitete Restbild vorzugsweise dadurch, dass er das Restbild wie ein Teilbild behandelt. Alternativ ist es aber auch möglich, dass der Rechner zum Ermitteln des bearbeiteten Restbildes die Pixel des Restbildes mit einem ortsunabhängigen Restbildfaktor multipliziert.Of the Computer determines the processed residual image preferably by that he treats the rest of the picture like a partial picture. Alternatively it is but it is also possible that the computer for determining the edited residual image, the pixels of the residual image multiplied by a location-independent residual image factor.
Der Filterkern weist eine Filterkerngröße und mindestens einen weiteren Parameter auf. Zumindest der weitere Parameter sollte vorzugsweise für alle Teilbilder der gleiche sein.Of the Filter core has a filter core size and at least one other Parameters. At least the other parameter should preferably for all Partial images will be the same.
Weitere Vorteile und Einzelheiten ergeben sich aus der nachfolgenden Beschreibung eines Ausführungsbeispiels in Verbindung mit den Zeichnungen. Dabei zeigen in Prinzipdarstellung:Further Advantages and details will become apparent from the following description an embodiment in conjunction with the drawings. Here are a schematic diagram:
Gemäß
Das
Computerprogramm
Bei
Aufruf des Computerprogramms
Gemäß
Die
Schritte S1 bis S3 sind – einzeln
oder zusammen – nur
optional. Sie sind also nicht zwingend erforderlich. Aus diesem
Grund sind sie in
In
einem Schritt S4 nimmt der Rechner über die Schnittstelle
In einem Schritt S5 ermittelt der Rechner sodann anhand des Grundfaktors F1, eines Zusatzfaktors F2 und des Anfangsbildes B einen Schwellwert t. Er ermittelt nämlich den durchschnittlichen Helligkeitswert aller Pixel des Anfangsbildes B und multipliziert diesen Durchschnittswert mit dem Grundfaktor F1 und dem Zusatzfaktor F2.In a step S5, the calculator then determines based on the basic factor F1, an additional factor F2 and the initial image B a threshold value t. He ascertains the average brightness value of all pixels of the initial image B and multiply this average value by the basic factor F1 and the additional factor F2.
In
einem Schritt S6 zerlegt der Rechner das Anfangsbild B in Teilbilder
Ti (i = 0 ... n) und ein Restbild R. Die
Zerlegung kann dabei auf beliebige Weise erfolgen. Beispielsweise
kann eine Zerlegung mittels einer Gauß-Pyramide erfolgen. Gemäß Ausführungsbeispiel – siehe
Der
Aufbau und die Wirkungsweise der Laplace-Pyramide ist Fachleuten
allgemein bekannt: Die Reduktionsblöcke
Im Ergebnis enthält somit das Teilbild T0 diejenigen Anteile des Anfangsbildes B, die mit Frequenzen örtlich variieren, die in einem für das Teilbild T0 charakteristischen höchsten Teilbildfrequenzbereich liegen. Das Teilbild T1 enthält diejenigen Anteile des Anfangsbildes B, die mit Frequenzen örtlich variieren, die in einem korrespondierenden zweithöchsten Teilbildfrequenzbereich liegen. Analoge Ausführungen gelten für die anderen Teilbilder Ti. Für die korrespondierenden Teilbildfrequenzbereiche gilt, dass die jeweiligen Frequenzen immer kleiner werden.As a result, the partial image T 0 thus contains those portions of the initial image B which vary locally with frequencies which lie in a highest partial image frequency range that is characteristic of the partial image T 0 . The partial image T 1 contains those portions of the initial image B that vary locally with frequencies that lie in a corresponding second highest field frequency range. Analogous explanations apply to the other partial images T i . For the corresponding field frequency ranges, the respective frequencies are becoming smaller and smaller.
Das Restbild R enthält zumindest einen Gleichanteil des zweidimensionalen Anfangsbildes B. Es enthält also den Anteil des zweidimensionalen Anfangsbildes B, der örtlich nicht variiert. Wenn die Laplace-Pyramide hoch genug ist, enthält das Restbild R sogar ausschließlich den Gleichanteil. Wenn die Laplace-Pyramide hingegen nicht hoch genug ist, enthält das Restbild R zusätzlich zum Gleichanteil diejenigen Anteile des zweidimensionalen Anfangsbildes B, die mit Frequenzen örtlich variieren, die kleiner als die Frequenzen des niedrigsten Teilbildfrequenzbereichs sind.The Residual image R contains at least a DC component of the two-dimensional initial image B. It contains So the proportion of the two-dimensional initial image B, which is not locally varied. If the Laplace pyramid is high enough, contains the residual image R even exclusively the DC component. If the Laplace pyramid, however, not high is enough contains the residual image R in addition to the DC component those parts of the two-dimensional initial image B, which with frequencies locally which are smaller than the frequencies of the lowest field frequency range are.
In einem Schritt S7 prüft der Rechner als nächstes, ob er das Restbild R wie ein Teilbild behandeln soll. Wenn dies der Fall ist, stuft der Rechner in einem Schritt S8 das Restbild R als (zusätzliches) Teilbild ein. Anderenfalls multipliziert der Rechner in einem Schritt S9 die Pixel des Restbildes R mit einem ortsunabhängigen Restbildfaktor F3 und ermittelt so ein bearbeitetes Restbild R'. Der Restbildfaktor F3 wird vom Rechner dabei derart bestimmt, dass die Pixel des bearbeiteten Restbildes R' den Helligkeitswert H aufweisen.In a step S7, the computer next checks whether it should treat the residual image R as a partial image. If this is the case, in a step S8 the computer classifies the residual image R as an (additional) partial image. Otherwise, the calculator multiplies in a step S9, the pixels of the residual image R with a location-independent residual image factor F3 and thus determines a processed residual image R '. The residual image factor F3 is determined by the computer in such a way that the pixels of the processed residual image R 'have the brightness value H.
In einem Schritt S10 selektiert der Rechner eines der Teilbilder Ti und in einem Schritt S11 ein Pixel des selektierten Teilbildes Ti.In a step S10, the computer selects one of the partial images T i and in a step S11 a pixel of the selected partial image T i .
In
einem Schritt S12 prüft
der Rechner sodann, ob er innerhalb eines zweidimensionalen Filterkerns
um das im Schritt S11 selektierte Pixel herum eine Struktur erkennen
kann. In Abhängigkeit
davon berechnet der Rechner dann einen Faktor g für dieses
Pixel dieses Teilbildes Ti und speichert
diesen Faktor g. Der Schritt S12 wird später in Verbindung mit
In
einem Schritt S13 prüft
der Rechner, ob er den Schritt S12 bereits für alle Pixel des selektierten Teilbildes
Ti ausgeführt hat. Wenn dies nicht der
Fall ist, geht der Rechner zu einem Schritt S14 über, in dem er ein anderes,
bisher noch nicht behandeltes Pixel des selektierten Teilbildes
Ti selektiert. Sodann geht er zum Schritt
Wenn der Rechner hingegen bereits für alle Pixel des selektierten Teilbildes Ti den korrespondierenden pixelspezifischen Faktor g ermittelt hat, geht der Rechner zu einem Schritt S15 über. Im Schritt S15 multipliziert der Rechner jedes Pixel des selektierten Teilbildes Ti mit dem korrespondierenden pixelspezifischen Faktor g und ermittelt so ein bearbeitetes Teilbild Ti'.By contrast, if the computer has already determined the corresponding pixel-specific factor g for all pixels of the selected partial image T i , the computer transfers to a step S 15. In step S15, the computer multiplies each pixel of the selected partial image T i by the corresponding pixel-specific factor g and thus determines a processed partial image T i '.
Gemäß
In
einem Schritt S16 prüft
der Rechner, ob er bereits alle Teilbilder Ti (ggf.
einschließlich
des Restbildes R, sofern dieses im Schritt S8 als Teilbild eingestuft
wurde) bearbeitet hat. Wenn dies nicht der Fall ist, selektiert
der Rechner in einem Schritt S17 ein anderes, bisher noch nicht
bearbeitetes Teilbild Ti (ggf. auch das
Restbild R) und geht zum Schritt S11 zurück. Anderenfalls kombiniert
der Rechner in einem Schritt S18 die bearbeiteten Teilbilder Ti' und das
bearbeitete Restbild R' zu
einem Endbild B'.
Das Kombinieren erfolgt – siehe
Wenn
das Restbild R im Schritt S8 als Teilbild eingestuft wurde, wird
weiterhin ein Schritt S19 ausgeführt.
Im Schritt S19 ermittelt der Rechner den durchschnittlichen Helligkeitswert
der Pixel des Endbildes B'.
Sodann ermittelt er einen ortsunabhängigen Endbildfaktor F4 und
multipliziert das Endbild B' mit
diesem Endbildfaktor F4. Der Endbildfaktor F4 ist dabei derart bestimmt,
dass der durchschnittliche Helligkeitswert der Pixel des Endbildes
B' auf den Helligkeitswert
H skaliert wird. Dieser Schritt S19 wird gemäß
In
einem Schritt S20 wird das Endbild B' vom Rechner schließlich über eine geeignete Ausgabeeinrichtung
Als letztes prüft der Rechner in einem Schritt S21, ob ein weiteres Anfangsbild B bearbeitet werden soll. Wenn dies der Fall ist, wird zum Schritt S4 zurückgegangen. Anderenfalls ist die Ausführung des Bearbeitungsverfahrens beendet.When last checks the computer in a step S21, whether another initial image B to be edited. If this is the case, it becomes the step S4 declined. Otherwise, the execution of the Processing procedure ended.
Zum
Ermitteln, ob innerhalb des zweidimensionalen Filterkerns um ein
Pixel herum in einem der Teilbilder Ti eine
Struktur erkannt wird, also zur Implementierung des Schrittes 512,
wird gemäß
Zunächst
bildet der Rechner innerhalb des Filterkerns um das jeweilige Pixel
des jeweiligen Teilbildes Ti herum Mittelwerte μj von mehreren
Pixelgruppen. Die Pixel jeder Pixelgruppe liegen dabei zumindest im
Wesentlichen auf einer Geraden und enthalten das jeweilige selektierte
Pixel als Zentralpixel. Wenn der zweidimensionale Filterkern beispielsweise – siehe
First, the computer forms average values μj of several pixel groups within the filter kernel around the respective pixel of the respective partial image T i . The pixels of each pixel group lie at least substantially on a straight line and contain the respective selected pixel as a central pixel. For example, if the two-dimensional filter kernel - see
In
einem Schritt S32 ermittelt der Rechner sodann einen Durchschnittswert μ der Mittelwerte μj und den
Minimalwert μ' der Mittelwerte μj. Weiterhin ermittelt
er im Schritt S32 die Differenz
In
einem Schritt S33 prüft
der Rechner, ob die Differenz
Der Rechner setzt also den Faktor dann und nur dann auf den Minimalwert M bzw. Mi, wenn der für die Schwankung der Mittelwerte μj charakteristische Wert δ unterhalb des Schwellwerts t liegt. Liegt der charakteristische Wert δ hingegen oberhalb des Schwellwerts t, hebt der Rechner den Faktor linear mit diesem Wert δ an. In jedem Fall aber wird der Faktor für das jeweilige Pixel des jeweiligen Teilbildes Ti anhand des wertes δ ermittelt.The calculator then sets the factor to the minimum value M or M i only when the value δ characteristic of the fluctuation of the average values μj is below the threshold value t. On the other hand, if the characteristic value δ is above the threshold value t, the calculator increases the factor linearly with this value δ. In any case, however, the factor for the respective pixel of the respective partial image T i is determined on the basis of the value δ.
Die
Mittelwerte μj
sind statistische Werte der Pixelgruppen. Es können aber auch andere statistische
Werte herangezogen werden, insbesondere die Varianzen σ der Pixelgruppen.
In diesem Fall wird – alternativ
zur Vorgehensweise von
Zunächst
ermittelt der Rechner in einem Schritt S41 innerhalb des Filterkerns
um das selektierte Pixel des jeweiligen Teilbildes Ti herum
Varianzen σj
von mehreren Pixelgruppen. Die Pixelgruppen sind dabei dieselben
Pixelgruppen, die bereits in Verbindung mit den
First, the computer determines in a step S41 within the filter kernel around the selected pixel of the respective subpicture T i variances σj of several pixel groups. The pixel groups are the same pixel groups already in connection with the
In einem Schritt S42 ermittelt der Rechner sodann das Maximum σ und das Minimum σ' der Varianzen σj sowie die Differenz δ dieser beiden Werte σ, σ'. Diese Differenz δ ist wieder charakteristisch für die Schwankung der Varianzen σj. Sie zeigt insbesondere an, ob das betrachtete Pixel Bestandteil einer Bildkante ist.In In a step S42, the computer then determines the maximum σ and the Minimum σ 'of the variances σj and the Difference δ this both values σ, σ '. This difference δ is again characteristic of the Variation of variances σj. In particular, it indicates whether the considered pixel is part of it an image edge is.
Die
Schritte S43 bis S45 entsprechen den Schritten S33 bis S35 von
Aus
den obigen Ausführungen
zu den
Wie bereits erwähnt, kann der Minimalwert M bzw. Mi alternativ für alle Teilbilder Ti der gleiche sein oder aber individuell für das jeweilige Teilbild Ti. Wenn der Minimalwert Mi individuell für das jeweilige Teilbild Ti ist, kann er insbesondere dem Gewichtungsfaktor der kurz vor der vorliegenden Erfindung eingereichten Patentanmeldung mit dem Titel „Bearbeitungsverfahren für ein zweidimensionales Anfangsbild und hiermit korrespondierende Gegenstände", Erfinder Philipp Bernhardt und Markus Lendl, Anmeldetag beim Deutschen Patentamt 22.06.2005, amtliches Aktenzeichen 10 2005 028 892.8, entsprechen. Unabhängig davon, welche Vorgehensweise ergriffen wird, ist der Minimalwert M bzw. Mi aber für alle Pixel des jeweiligen Teilbildes Ti der gleiche.As already mentioned, the minimum value M or M i can alternatively be the same for all partial images T i or else individually for the respective partial image T i . If the minimum value M i is individual to the respective subpicture T i , it may in particular be the weighting factor of the patent application filed shortly before the present invention titled "Processing method for a two-dimensional initial image and objects corresponding thereto", inventors Philipp Bernhardt and Markus Lendl, filing date the German Patent Office 22.06.2005, official file number 10 2005 028 892.8. Regardless of which procedure is taken, the minimum value M or M i is the same for all pixels of the respective subpicture T i .
Mittels
der erfindungsgemäßen Vorgehensweise
ist es insbesondere möglich,
den natürlichen Rauscheindruck
zu erhalten und dennoch ein gutes Rauschreduktionsverfahren anzuwenden.
Letztendlich wird dadurch die intuitive Erkennbarkeit von Nutzsignalen
durch den Anwender
Claims (20)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE200510032287 DE102005032287B4 (en) | 2005-07-11 | 2005-07-11 | Processing method for a two-dimensional initial image and objects corresponding thereto |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE200510032287 DE102005032287B4 (en) | 2005-07-11 | 2005-07-11 | Processing method for a two-dimensional initial image and objects corresponding thereto |
Publications (2)
Publication Number | Publication Date |
---|---|
DE102005032287A1 true DE102005032287A1 (en) | 2007-01-25 |
DE102005032287B4 DE102005032287B4 (en) | 2010-11-11 |
Family
ID=37575386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE200510032287 Expired - Fee Related DE102005032287B4 (en) | 2005-07-11 | 2005-07-11 | Processing method for a two-dimensional initial image and objects corresponding thereto |
Country Status (1)
Country | Link |
---|---|
DE (1) | DE102005032287B4 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102009032441A1 (en) | 2009-07-09 | 2011-01-27 | Siemens Aktiengesellschaft | Two-dimensional X-ray image processing method for X-ray system, involves realizing detector-independent visual impressions by detector-specific balancing factors, and adding balanced partial images adjacent to final image |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5526446A (en) * | 1991-09-24 | 1996-06-11 | Massachusetts Institute Of Technology | Noise reduction system |
US5561724A (en) * | 1993-06-11 | 1996-10-01 | Hitachi Medical Corporation | Method conducting smoothing processing in a linear direction having a minimized magnitude of density change |
US5602934A (en) * | 1993-09-08 | 1997-02-11 | The Regents Of The University Of California | Adaptive digital image signal filtering |
US6252931B1 (en) * | 1998-10-24 | 2001-06-26 | U.S. Philips Corporation | Processing method for an original image |
DE10146582A1 (en) * | 2001-09-21 | 2003-04-24 | Micronas Munich Gmbh | Device and method for the subband decomposition of image signals |
-
2005
- 2005-07-11 DE DE200510032287 patent/DE102005032287B4/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5526446A (en) * | 1991-09-24 | 1996-06-11 | Massachusetts Institute Of Technology | Noise reduction system |
US5561724A (en) * | 1993-06-11 | 1996-10-01 | Hitachi Medical Corporation | Method conducting smoothing processing in a linear direction having a minimized magnitude of density change |
US5602934A (en) * | 1993-09-08 | 1997-02-11 | The Regents Of The University Of California | Adaptive digital image signal filtering |
US6252931B1 (en) * | 1998-10-24 | 2001-06-26 | U.S. Philips Corporation | Processing method for an original image |
DE10146582A1 (en) * | 2001-09-21 | 2003-04-24 | Micronas Munich Gmbh | Device and method for the subband decomposition of image signals |
Non-Patent Citations (2)
Title |
---|
GONZALEZ, R.C. ; WOODS, R.E.: Digital Image Proces-sing - second edition, Addison-Wesley, 2002, S. 239-241 * |
GONZALEZ, R.C.; WOODS, R.E.: Digital Image Proces- sing - second edition, Addison-Wesley, 2002, S. 239-241 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102009032441A1 (en) | 2009-07-09 | 2011-01-27 | Siemens Aktiengesellschaft | Two-dimensional X-ray image processing method for X-ray system, involves realizing detector-independent visual impressions by detector-specific balancing factors, and adding balanced partial images adjacent to final image |
DE102009032441B4 (en) * | 2009-07-09 | 2016-06-09 | Siemens Healthcare Gmbh | Method for automatically processing an X-ray image, X-ray system and computer program product |
Also Published As
Publication number | Publication date |
---|---|
DE102005032287B4 (en) | 2010-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE102013021991B4 (en) | Adaptive filter mechanism for the removal of coding image errors in video data | |
DE19916354B4 (en) | Method and device for analyzing image structures | |
DE102007003714A1 (en) | Method and system for generating object image layers | |
DE102006005803A1 (en) | Method for noise reduction in imaging methods | |
DE102007013570A1 (en) | Method for noise reduction in digital images with locally different and directional noise | |
DE102007028226A1 (en) | Evaluation method for a temporal sequence of X-ray images and objects corresponding thereto | |
WO2017153355A1 (en) | Method and device for carrying out eye gaze mapping | |
DE102006020430A1 (en) | Method and apparatus for noise reduction in a sequence of fluoroscopic images | |
EP1197067A1 (en) | Method and circuit system for improving image contrast | |
DE102005032287B4 (en) | Processing method for a two-dimensional initial image and objects corresponding thereto | |
WO2004110276A1 (en) | Method and unit for refining x-ray images | |
DE102005028892A1 (en) | Processing a two-dimensional initial image by dismantling initial image into partial images and residual image, assigning one of the partial images as pilot image and assigning pilot frequency, and determining and summing weighting factors | |
DE10324908A1 (en) | Self learning medical X ray image processing procedure modifies image data using parameter set updated from standard based on copies previous parameters | |
EP1489842B1 (en) | Motion vector based method and device for the interpolation of picture elements | |
DE102006044189A1 (en) | A method for providing a variety of processed image data, methods for processing a variety of image data and an X-ray image system | |
EP3637362A1 (en) | Method for adjusting an image impression | |
DE19604050B4 (en) | Edge pixel smoothing of block-coded images - has edge pixels processed by edge filter having low=pass stages appropriate to motion vector ranges | |
DE69823998T2 (en) | SYSTEM AND METHOD FOR LOCAL IMAGE PROCESSING | |
DE102005046203B3 (en) | Diagnostic picture data selectively representing method, involves determining three-dimensional picture information, and subdividing determined information into representative and non-representative unit data values | |
DE102006013472A1 (en) | Method and device for processing image data for automatic detection and / or segmentation of anatomical features from computed tomography images | |
EP3304169B1 (en) | Method for operating an augmented reality glasses device and augmented reality glasses device | |
DE10301341A1 (en) | Real time processing of video data, e.g. using a restoration algorithm, whereby in a first playback step only a section of the image is played back and processing parameters then selected for processing the whole image sequence | |
DE102019202878A1 (en) | Method for noise reduction in computer tomography data, X-ray equipment, computer programs and electronically readable data carriers | |
EP3664036B1 (en) | Method and system for calculating an output of a scrollable stack of images provided by a tomographer | |
EP3689005B1 (en) | Method and apparatus for computer-assisted processing of audiosignals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
OP8 | Request for examination as to paragraph 44 patent law | ||
8120 | Willingness to grant licences paragraph 23 | ||
8364 | No opposition during term of opposition | ||
R020 | Patent grant now final |
Effective date: 20110211 |
|
R081 | Change of applicant/patentee |
Owner name: SIEMENS HEALTHCARE GMBH, DE Free format text: FORMER OWNER: SIEMENS AKTIENGESELLSCHAFT, 80333 MUENCHEN, DE |
|
R119 | Application deemed withdrawn, or ip right lapsed, due to non-payment of renewal fee |