DE10325632A1 - X ray image structure recognition improvement process uses Fourier transform for user selected Gaussian or other weighting filter in frequency domain - Google Patents

X ray image structure recognition improvement process uses Fourier transform for user selected Gaussian or other weighting filter in frequency domain Download PDF

Info

Publication number
DE10325632A1
DE10325632A1 DE2003125632 DE10325632A DE10325632A1 DE 10325632 A1 DE10325632 A1 DE 10325632A1 DE 2003125632 DE2003125632 DE 2003125632 DE 10325632 A DE10325632 A DE 10325632A DE 10325632 A1 DE10325632 A1 DE 10325632A1
Authority
DE
Germany
Prior art keywords
frequency
intensity distribution
structures
signal components
image
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.)
Ceased
Application number
DE2003125632
Other languages
German (de)
Inventor
Michael Dr. Thoms
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.)
Duerr Dental SE
Original Assignee
Duerr Dental SE
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 Duerr Dental SE filed Critical Duerr Dental SE
Priority to DE2003125632 priority Critical patent/DE10325632A1/en
Publication of DE10325632A1 publication Critical patent/DE10325632A1/en
Ceased legal-status Critical Current

Links

Classifications

    • G06T5/94—
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/10—Image enhancement or restoration by non-spatial domain filtering
    • 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/10—Image acquisition modality
    • G06T2207/10116—X-ray image

Abstract

An X ray image structure recognition improvement process (20) stores the image as spatial intensity distribution and uses a Fourier Transform to the frequency domain for filtering with weighting using a control slider (26, 28) selected central value and Gaussian or other profile to enhance average structure dimensions and an inverse transform to return an enhanced image (24). Includes INDEPENDENT CLAIMs for image processing equipment using the procedure.

Description

Die Erfindung betrifft ein Verfahren zur Verbesserung der Erkennbarkeit von unterschiedlichen Strukturen auf Durchstrahlungsbildern sowie eine hierfür geeignete Bildverarbeitungsvorrichtung.The The invention relates to a method for improving the recognizability of different structures on radiographs as well one for this suitable image processing device.

In der medizinischen Diagnostik ist die Auswertung von Durchstrahlungsbildern wie z.B. Röntgenbildern von großer Bedeutung. Knochen, Implantate oder ähnliche Strukturen heben sich im allgemeinen von umgebenden Weichteilgeweben deutlich ab und sind deswegen gut erkennbar. Umgekehrt sind Weichteilgewebestrukturen, etwa Sehnen oder Blutgefäße, in der Regel nur sehr undeutlich auf Durchstrahlungsbildern wiedergegeben. Bei zahlreichen Krankheitsbildern kommt es aber auch oder gerade auf die Erkennbarkeit der Weichteilgewebestrukturen an. Darüber hinaus ist es häufig schwierig, gleichartiges Gewebe voneinander zu unterscheiden. So sind z.B. kleinere Knochen, die auf einem Durchstrahlungsbild über einem größeren Knochen abgebildet sind, oft kaum mit bloßem Auge auszumachen; entsprechendes gilt für Weichteilgewebestrukturen. Ärzte können in derartigen Fällen deswegen allein auf der Grundlage der Durchstrahlungsbilder oft keine oder nur eine recht unzuverlässige Diagnose stellen.In Medical diagnostics is the evaluation of radiographic images such as. X-ray images of great Importance. Bones, implants or similar structures rise up generally clearly and are from surrounding soft tissue therefore easily recognizable. Conversely, soft tissue structures are such as tendons or blood vessels in which Usually only reproduced very clearly on radiographs. With numerous clinical pictures, however, it also or precisely occurs on the recognizability of the soft tissue structures. Furthermore it is common difficult to differentiate between similar tissues. So are e.g. smaller bones on a radiograph over a bigger bones are often depicted with the naked eye; corresponding applies to Soft tissue structures. doctors can in such cases therefore often based solely on the radiographs make no or only a rather unreliable diagnosis.

Eine gewisse Verbesserung hat die Digitalisierung von Durchstrahlungsbildern geschaffen. Durch bekannte Methoden der Bildverarbeitung, etwa der Kontrastverstärkung innerhalb ausgewählter Bildausschnitte, können z.B. Weichteilgewebestrukturen unter Umständen deutlicher hervorgehoben werden. Damit gelingt es aber in der Regel nicht, eine über einem Knochen liegende Sehne erkennbar werden zu lassen. Die kleineren Schwankungen des Signalpegels von Bildsignalanteilen, welche die Sehne repräsentieren, können nämlich vor dem hohen Hintergrund-Signalpegel des Knochens nicht signifikant hervortreten. Ein zur Anzeige verwendeter Bildschirm gibt die geringen Schwankungen des Signalpegels zwar günstigstenfalls noch als Intensitätsschwankungen wieder, doch sind diese meist so gering, daß sie kaum noch für das bloße Auge erkennbar sind.A The digitization of radiographic images has some improvement created. By known methods of image processing, such as the contrast enhancement within selected Image sections, can e.g. Soft tissue structures may be highlighted more clearly become. However, this usually does not succeed, one over one To show the tendon lying on the bone. The smaller fluctuations the signal level of image signal components which represent the chord, can namely before the high background signal level of the bone is not significant emerge. A screen used for the display gives the small ones Fluctuations in the signal level are at best still as intensity fluctuations again, but these are usually so small that they are hardly visible to the naked eye are recognizable.

Aufgabe der Erfindung ist es daher, ein Verfahren und eine Vorrichtung zur Verbesserung der Erkennbarkeit von Strukturen unterschiedlicher Art auf Durchstrahlungsbildern anzugeben.task the invention is therefore a method and an apparatus for Improving the recognizability of different structures Specify type on radiographs.

Gelöst wird diese Aufgabe durch eine Verfahren mit den Merkmalen des Patentanspruchs 1 sowie durch eine Vorrichtung mit den Merkmalen des Patentanspruchs 9.Is solved this task by a method with the features of the claim 1 and by a device with the features of the claim 9th

Der Erfindung liegt die Erkenntnis zugrunde, daß sich in den meisten Fällen die Strukturen, deren Erkennbarkeit verbessert werden soll, hinsichtlich ihrer Größe und Feinheit von den übrigen auf dem Durchstrahlungsbild abgebildeten Strukturen mehr oder weniger deutlich unterscheiden. Da kleinere und feinere Strukturen im Fourierspektrum durch höhere Frequenzen wiedergegeben werden als große grobe Strukturen, kann durch eine Veränderung der Gewichtung zwischen hochfrequenten und niedrigfrequenten Bildsignalanteilen im Fourierspektrum eine Verstärkung des Bildkontrastes entweder für kleine feine oder für große grobe Strukturen erzielt werden. Je nachdem, ob die schwer erkennbaren Strukturen feiner oder gröber sind als die gut erkennbaren Strukturen, wird die Gewichtung der Bildsignalanteile im Frequenzraum entweder zugunsten der hochfrequenten oder der niedrigfrequenten Bildsignalanteile verändert.The Invention is based on the knowledge that in most cases Structures whose recognizability should be improved with regard to their size and delicacy from the rest Structures depicted on the radiograph more or less distinguish clearly. Because smaller and finer structures in the Fourier spectrum by higher Frequencies can be reproduced as large rough structures, by a change the weighting between high-frequency and low-frequency image signal components in the Fourier spectrum an amplification the image contrast either for small fine or for size rough structures can be achieved. Depending on whether the hard to see Structures finer or coarser are as the well recognizable structures, the weighting of the Image signal components in the frequency domain either in favor of the high-frequency or the low-frequency image signal components changed.

Besonders deutlich treten die zunächst schwer erkennbaren Strukturen hervor, wenn die zu gewichtenden Bildsignalanteile gemäß Patentanspruch 2 festgelegt werden.Especially they kick clearly at first structures that are difficult to see when the image signal components to be weighted according to claim 2 can be set.

Bei der besonders einfachen Filterung nach Patentanspruch 3 wird die Frequenzraum-Intensitätsverteilung lediglich mit einer Filterfunktion multipliziert.at the particularly simple filtering according to claim 3 is the Frequency-space intensity distribution just multiplied by a filter function.

Durch Verwendung von Frequenz-Zentralwerten und Profilfunktionen zur Festlegung der zu gewichtenden Frequenzbereiche nach den Patentansprüchen 4 und 5 läßt sich mit relativ wenigen Parametern die Filterung gezielt beeinflussen.By Use of frequency central values and profile functions for the definition the frequency ranges to be weighted according to claims 4 and 5 can selectively influence the filtering with relatively few parameters.

Besonders geeignet als Profilfunktion ist gemäß Patentanspruch 6 eine Gaußfunktion, da diese die Eigenschaft hat, auch bei der Fourier-Rücktransformation eine Gaußfunktion zu bleiben. Die Filterung läßt sich dann im Ortsraum als Faltung der Intentsitätsverteilung mit einer Gaußfunktion darstellen. Dadurch wird verhindert, daß es durch die Filterung zu einem Auseinanderfließen von Stellen im Bild kommt, an denen die Intensitätsverteilung sich abrupt ändert und die somit einen besonders hohen Kontrast haben.Especially suitable as a profile function according to claim 6 is a Gaussian function, since this has the property, also with the Fourier inverse transformation a Gaussian function too stay. The filtering can then in space as a convolution of the distribution of intentions with a Gaussian function represent. This prevents it from filtering too a flow apart comes from places in the picture where the intensity distribution changes abruptly and which therefore have a particularly high contrast.

Welche Frequenzen oder Frequenzbereiche in ihrer Gewichtung verändert werden, bestimmt sich nach der mittleren Strukturgröße der Strukturen, deren Erkennbarkeit verbessert werden soll. Die mittlere Strukturgröße oder entsprechende Frequenzbereiche können entweder fest vorgegeben oder gemäß Patentanspruch 7 mit Hilfe von Stellgliedern an der Bildverarbeitungsvorrichtung oder über eine Bedienoberfläche eines übergeordneten Rechners frei wählbar sein. Durch Verändern der maßgeblichen Filterparameter kann ein behandelnder Arzt somit bei unterschiedlichsten Durchleuchtungsbildern gezielt die Erkennbarkeit der ihn interessierenden Strukturen verbessern.Which frequencies or frequency ranges are changed in their weighting is determined by the average structure size of the structures, the recognizability of which is to be improved. The average structure size or corresponding frequency ranges can either be fixed or according to godparents Claim 7 can be freely selected with the help of actuators on the image processing device or via a user interface of a higher-level computer. By changing the relevant filter parameters, a treating doctor can specifically improve the recognizability of the structures that interest him in the most varied of fluoroscopic images.

Darüber hinaus kommt auch eine automatische Bestimmung der Frequenzbereiche im Wege eines adaptiven Verfahrens nach Patentanspruch 8 in Betracht.Furthermore there is also an automatic determination of the frequency ranges in the Way of an adaptive method according to claim 8 into consideration.

Die Auswahl der Strukturen, deren Erkennbarkeit verbessert werden soll, kann dabei z.B. wie in den Patentansprüche 9 oder 10 angeben erfolgen.The Selection of the structures whose recognizability should be improved can e.g. as stated in claims 9 or 10.

Eine zusätzliche Hochfrequenzfilterung nach Patentanspruch 11, z.B. mit einem Gaußfilter nach Patentanspruch 12, führt zu einer Erhöhung des Signal-Rausch-Abstands, da Bildstrukturen wiedergebende Bildsignalanteile gegenüber einem hochfrequenten Untergrundrauschen verstärkt werden. Eine solche Filterung trägt der Tatsache Rechnung, daß bei den in der Praxis häufig darzustellenden Bildern die Fourier-Amplituden mit zunehmender Frequenz f abnehmen.A additional High frequency filtering according to claim 11, e.g. with a Gaussian filter according to claim 12, leads to an increase of the signal-to-noise ratio, since image structures represent image signal components across from high-frequency background noise. Such filtering bears the fact Bill that at often in practice the Fourier amplitudes with increasing frequency f decrease.

Die vorstehend zu dem Verfahren angeführten vorteilhaften Ausgestaltungen und Vorteile gelten sinngemäß auch für die erfindungsgemäße Bildbearbeitungsvorrichtung.The Advantageous embodiments mentioned above for the method and advantages apply mutatis mutandis to the image processing device according to the invention.

Weitere Merkmale und Vorteile der Erfindung ergeben sich aus der nachfolgenden Beschreibung eines Ausführungsbeispiels anhand der Zeichnung. Darin zeigen:Further Features and advantages of the invention result from the following Description of an embodiment based on the drawing. In it show:

1 ein Röntgenbild, auf dem ein Fingerknochen und Weichteilgewebestrukturen erkennbar sind; 1 an x-ray image showing finger bones and soft tissue structures;

2 ein Bockschaltbild einer erfindungsgemäßen Bildverarbeitungsvorrichtung; 2 a block diagram of an image processing device according to the invention;

3 eine eindimensionale periodische Ortraum-Intensitätsverteilung I(x); 3 a one-dimensional periodic spatial intensity distribution I (x);

4 die Frequenzraum-Intensitätsverteilung F(f) der Ortsraum-Intensitätsverteilung I(x) aus 3; 4 the frequency space intensity distribution F (f) of the spatial space intensity distribution I (x) 3 ;

5 eine durch Filterung der in 4 gezeigten Frequenzraum-Intensitätsverteilung F(f) erzeugte gefilterte Frequenzraum-Intensitätsverteilung F' (f) ; 5 one by filtering the in 4 shown frequency space intensity distribution F (f) generated filtered frequency space intensity distribution F '(f);

6 eine durch Fourier-Rücktransformation der in 5 gezeigten gefilterten Frequenzraum-Intensitätsverteilung F'(f) gewonnene Ortsraum-Intensitätsverteilung I'(x); 6 one by Fourier inverse transformation of the in 5 Filtered frequency space intensity distribution F '(f) shown obtained spatial space intensity distribution I'(x);

7 eine Frequenzraum-Intensitätsverteilung F(f) einer eindimensionalen Ortsraum-Intensitätsverteilung mit zwei Profilfunktionen g1(f) und g2(f); 7 a frequency space intensity distribution F (f) of a one-dimensional spatial space intensity distribution with two profile functions g 1 (f) and g 2 (f);

8 eine durch Filterung der in 7 gezeigten Frequenzraum-Intensitätsverteilung F(f) erhaltene gefilterte Frequenzraum-Intensitätsverteilung F'(f). 8th one by filtering the in 7 shown frequency space intensity distribution F (f) obtained filtered frequency space intensity distribution F '(f).

1 zeigt ein typisches Röntgenbild 10 eines Fingers 12, auf dem mehrere Fingerknochen 14 sowie diese umgebendes Weichteilgewebe 16 erkennbar sind. Aufgrund ihrer im Verhältnis zu dem Weichteilgewebe 16 hohen Dichte zeichnen sich die Fingerknochen 14 mit hohem Kontrast von diesem ab, während Weichteilgewebestrukturen wie z.B. Sehnen 18 kaum auf dem Röntgenbild 10 zu erkennen sind. Eine Diagnose derartiger Strukturen des Weichteilgewebes 16 anhand des Röntgenbildes 10 ist daher mit größeren Unsicherheiten behaftet. 1 shows a typical x-ray 10 a finger 12 on which several finger bones 14 as well as surrounding soft tissue 16 are recognizable. Because of their relative to the soft tissue 16 The finger bones are characterized by high density 14 with high contrast from this, while soft tissue structures such as tendons 18 hardly on the x-ray 10 are recognizable. A diagnosis of such structures of the soft tissue 16 based on the x-ray image 10 is therefore subject to greater uncertainties.

Selbst wenn das herkömmlich aufgenommene Röntgenbild 10 in einem Scanner digitalisiert und auf einem Bildschirm mit hohem Kontrast dargestellt wird, bleiben die Weichteilgewebestrukturen 18 schwer erkennbar. Der Grund hierfür liegt darin, daß Bildsignalanteile, die die Fingerknochen 14 oder anderes Hartgewebe hoher Dichte wiedergeben, einen sehr hohen Signalpegel haben. Kleinere Schwankungen des Signalpegels, welche die interessierenden Weichteilgewebestrukturen 18 repräsentieren, fallen in Relation zu den hohen Signalpegeln der Fingerknochen 14 kaum ins Gewicht. Ein hochwertiger Bildschirm kann die geringen Schwankungen des Signalpegels zwar unter Umständen noch als Intensitätsschwankungen wiedergeben, doch sind diese dann so gering, daß sie für das bloße Auge kaum erkennbar sind. Entsprechendes gilt im übrigen auch für PSL-Bildplatten (PSL = Photostimulierbare Lumineszenz), die nicht chemisch entwickelt werden, sondern bei denen das latent darin enthaltene Röntgenbild vor dem Betrachten auf einem Bildschirm im Wege eines optomechanischen Scanprozeß gelesen werden muß.Even if the conventionally taken x-ray image 10 digitized in a scanner and displayed on a screen with high contrast, the soft tissue structures remain 18 difficult to see. The reason for this is that the image signal components that the finger bones 14 or other high density hard tissue, have a very high signal level. Smaller fluctuations in the signal level, which the soft tissue structures of interest 18 represent, fall in relation to the high signal levels of the finger bones 14 hardly any weight. A high-quality screen may reflect the slight fluctuations in the signal level as fluctuations in intensity, but these are then so small that they are barely visible to the naked eye. The same also applies to PSL image plates (PSL = photo-stimulable luminescence), which are not developed chemically, but in which the latent X-ray image contained in the way before viewing on a screen an optomechanical scanning process must be read.

Um die Erkennbarkeit der Weichteilgewebestrukturen 18 zu verbessern, wird das digitalisierte Röntgenbild 10 deswegen in einer Bildverarbeitungsvorrichtung 20 aufbereitet, deren Aufbau in 2 gezeigt ist. Die Bildverarbeitungsvorrichtung 20 umfaßt einen Speicher MEM, in welchem die in einem Scanner SCAN erzeugten digitalen Bilddaten gespeichert werden können. Die digitalen Bilddaten geben für jeden Bildpunkt P = (x, y) Intensitäten I(x, y) wieder, wobei z.B. jeder Bildpunkt P mit 16 Bit codiert sein kann, so daß mehr als 65.000 Helligkeitswerte unterschieden werden können.To make the soft tissue structures recognizable 18 will improve the digitized x-ray image 10 therefore in an image processing device 20 processed, their structure in 2 is shown. The image processing device 20 comprises a memory MEM in which the digital image data generated in a scanner SCAN can be stored. The digital image data represent intensities I (x, y) for each pixel P = (x, y), it being possible, for example, for each pixel P to be coded with 16 bits, so that more than 65,000 brightness values can be distinguished.

Der Speicher MEM ist mit einer Fourier-Transformationseinheit FT verbunden, mit der sich aus dem Speicher MEM ausgelesene digitale Bilddaten einer Fourier-Transformation unterziehen lassen. Die von der Fourier-Transformationseinheit FT erzeugte Frequenzraum-Intensitätsverteilung F(fx, fy) ist eine komplexe Funktion über dem durch die Koordinaten fx und fy aufgespannten Frequenzraum und hat anschaulich die Bedeutung eines Amplitudendichtespektrums.The memory MEM is connected to a Fourier transformation unit FT, with which digital image data read out from the memory MEM can be subjected to a Fourier transformation. The frequency space intensity distribution F (f x , f y ) generated by the Fourier transformation unit FT is a complex function over the frequency space spanned by the coordinates f x and f y and clearly has the meaning of an amplitude density spectrum.

Die Bildverarbeitungsvorrichtung 20 umfaßt außerdem eine Filtereinheit FIL, mit der die Frequenzraum-Intensitätsverteilung F(fx, fy) so gefiltert werden kann, daß die Gewichtung unterschiedlicher Frequenzbereiche verändert wird. Dies wird weiter unten mit Bezug auf die 3 bis 6 näher erläutert.The image processing device 20 also includes a filter unit FIL, with which the frequency space intensity distribution F (f x , f y ) can be filtered so that the weighting of different frequency ranges is changed. This is discussed below with reference to the 3 to 6 explained in more detail.

Schließlich weist die Bildverarbeitungsvorrichtung 20 eine Fourier-Rück-Transformationseinheit FT–1 auf, welche die durch die Filtereinheit FIL gefilterte Frequenzraum-Intensitätsverteilung F'(fx, fy) zurück in den Ortsraum transformiert, wodurch eine modifizierte Ortsraum-Intensitätsverteilung I'(x, y) entsteht. An einen Ausgang 22 der Bildverarbeitungsvorrichtung 20 kann ein Bildschirm 24 angeschlossen werden, auf dem die modifizierte Ortsraum-Intensitätsverteilung I'(x, y) zur Anzeige gebracht werden kann.Finally, the image processing device 20 a Fourier-back transformation unit FT -1 , which transforms the frequency space intensity distribution F '(f x , f y ) filtered by the filter unit FIL back into the spatial space, whereby a modified spatial space intensity distribution I' (x, y) arises , To an exit 22 the image processing device 20 can be a screen 24 be connected, on which the modified spatial intensity distribution I '(x, y) can be displayed.

Die Filterung des Frequenzraum-Intensitätsverteilung F(fx, fy) in der Filtereinheit wird im folgenden anhand der 2 bis 6 näher erläutert.The filtering of the frequency space intensity distribution F (f x , f y ) in the filter unit is described below using the 2 to 6 explained in more detail.

3 zeigt für eine Bildkoordinate x eine Intensitätsverteilung I(x) im Ortsraum, wobei der Einfachheit halber eine periodische Verteilung angenommen ist. Die Intensitätsverteilung I(x) stellt eine Überlagerung einer großräumigen kosinusförmigen Intensitätsverteilung mit der Periode P1 mit einer kleinräumigen kosinusförmigen Intensitätsverteilung mit der Periode P2 dar. Die großräumige Intensitätsverteilung soll in diesem vereinfachten Beispiel die – hier kosinusförmig angenommene – Form von Knochen wiedergeben, während die kleinräumige Intensitätsverteilung die – ebenfalls kosinusförmig angenommene – Form von sich in Durchleuchtungsrichtung darüber angeordneten kosinusförmigen Weichteilgewebestrukturen repräsentiert, deren charakteristische Abmessungen deutlich kleiner als die der Knochen sind. Eine halbe Periodenlänge, d.h. ein Wellenberg der Kosinusfunktion, entspricht dabei jeweils der charakteristischen Strukturgröße. 3 shows an intensity distribution I (x) in the spatial domain for an image coordinate x, a periodic distribution being assumed for the sake of simplicity. The intensity distribution I (x) represents a superimposition of a large-scale cosine-shaped intensity distribution with the period P 1 with a small-scale cosine-shaped intensity distribution with the period P 2 the small-scale intensity distribution represents the - also assumed cosine-shaped - form of cosine-shaped soft tissue structures arranged above it in the direction of fluoroscopy, the characteristic dimensions of which are significantly smaller than those of the bones. Half a period length, ie a wave crest of the cosine function, corresponds to the characteristic structure size.

Wie in 3 erkennbar ist, tritt vor dem relativ hohen Signalpegel der dem Knochen zugeordneten großräumigen Intensitätsverteilung die den Weichteilgewebestrukturen zugeordnete kleinräumige Intensitätsverteilung so in den Hintergrund, daß die relativ kleinen Schwankungen der kleinräumigen Intensitätsverteilung kaum erkennbar werden, wenn die in 3 gezeigte Gesamtintensitätsverteilung I(x, y) auf einem Bildschirm angezeigt würde.As in 3 is recognizable, in front of the relatively high signal level of the large-area intensity distribution assigned to the bone, the small-area intensity distribution assigned to the soft tissue structures takes a back seat so that the relatively small fluctuations in the small-area intensity distribution are hardly recognizable when the in 3 total intensity distribution I (x, y) shown would be displayed on a screen.

4 zeigt die Frequenzraum-Intensitätsverteilung F(f) für die in 3 gezeigte Ortsraum-Intensitätsverteilung I(x). Die Frequenzraum-Intensitätsverteilung F(f) weist neben einem einen Gleichanteil der Intensitätsverteilung I(x) wiedergebenden Beitrag bei der Frequenz f = 0 einen Beitrag mit dem Betrag f1 und einen Beitrag mit dem Betrag f2 auf, wobei die kleinere Frequenz f1 den großräumigen Anteil der Ortsraum-Intensitätsverteilung und die größere Frequenz f2 den kleinräumigen Anteil der Ortsraum-Intensitätsverteilung repräsentiert. 4 shows the frequency space intensity distribution F (f) for the in 3 shown spatial area intensity distribution I (x). The frequency space intensity distribution F (f) has, in addition to a contribution reflecting a constant component of the intensity distribution I (x) at the frequency f = 0, a contribution with the amount f 1 and a contribution with the amount f 2 , the smaller frequency f 1 represents the large-scale portion of the spatial intensity distribution and the larger frequency f 2 represents the small-scale portion of the spatial intensity distribution.

Die Filterung der Frequenzraum-Intensitätsverteilung F(x) wird nun so durchgeführt, daß die Amplituden der Beitrage mit dem Frequenzbetrag f1 verringert und die Amplituden der Beiträge mit dem Frequenzbetrag f2 vergrößert werden. Dies kann z.B. durch die folgenden Operationen erzielt werden:

Figure 00100001
wobei r1 und r2 Verstärkungsfaktoren mit r1 > 1 und r2 < 1 sind. Die durch die Filterung erhaltene gefilterte Frequenzraum-Intensitätsverteilung F'(f) ist in 5 gezeigt. Die Gewichtung der Bildsignalanteile mit der Frequenz f1 und der Bildsignalanteile mit der Frequenz f2 ist dabei, wie ein Vergleich mit 4 ergibt, zugunsten des hochfrequenteren Bildsignalanteils der Frequenz f2 verändert worden.The filtering of the frequency space intensity distribution F (x) is now carried out in such a way that the amplitudes of the contributions are reduced with the frequency amount f 1 and the amplitudes of the contributions are increased with the frequency amount f 2 . This can be achieved, for example, by the following operations:
Figure 00100001
where r 1 and r 2 are amplification factors with r 1 > 1 and r 2 <1. The filtered frequency space intensity distribution F '(f) obtained by the filtering is in 5 shown. The weighting of the image signal components with the Fre quenz f 1 and the image signal components with the frequency f 2 is like a comparison with 4 results, have been changed in favor of the higher-frequency image signal component of the frequency f 2 .

6 zeigt die modifizierte Intensitätsverteilung I'(x), die durch Fourier-Rücktransformation aus der gefilterten Frequenzraum-Intensitätsverteilung F'(f) gewonnen ist. In der Darstellung ist gut zu erkennen, daß die kleinräumigen Schwankungen der Intensität, die die Weichteilgewebestrukturen wiedergeben sollen, jetzt eine größere Amplitude als vor der Filterung haben und sich deswegen erheblich besser von den großräumigen Schwankungen abheben, die die Knochen darstellen. 6 shows the modified intensity distribution I '(x), which is obtained by Fourier inverse transformation from the filtered frequency space intensity distribution F' (f). It can be clearly seen in the illustration that the small-scale fluctuations in the intensity which the soft tissue structures are intended to reproduce are now larger in amplitude than before the filtering and are therefore much better contrasted with the large-scale fluctuations which represent the bones.

Das vorstehend geschilderte Beispiel stellt durch seine Beschränkung auf kosinusförmige Strukturen in nur einer Dimension eine sehr grobe Vereinfachung dar, aus der das Wesen der Filterung aber besonders anschaulich hervorgeht. In realen Bildern haben die abgebildeten Strukturen hingegen einen in weiten Grenzen beliebigen Verlauf, weswegen die durch Fourier-Transformation erhaltene Frequenzraum-Intensitätsverteilung eine in der Frequenz kontinuierliche Funktion darstellt. Würden dann lediglich die Amplituden einzelner Frequenzen angehoben oder abgesenkt, wie dies in dem vorstehend geschilderten Beispiel der Fall ist, so wirkte sich dies auf das resultierende gefilterte Bild nur unmerklich aus.The The example described above is based on its limitation cosine Structures in only one dimension are a very rough simplification from which the essence of filtering is particularly vivid evident. In real pictures, the structures shown have on the other hand an arbitrary course within wide limits, which is why the Frequency space intensity distribution obtained by Fourier transform represents a continuous function in frequency. Would then only the amplitudes of individual frequencies are raised or lowered, as is the case in the example described above, so this had an imperceptible effect on the resulting filtered image out.

Aus diesem Grunde erfolgt die Gewichtung der Bildsignalanteile vorzugsweise nicht nur für einzelne diskrete Frequen zen, sondern für Frequenzbänder. Jedes Frequenzband, das in seiner Gewichtung verändert werden soll, wird mit Hilfe einer Profilfunktion festgelegt. Besonders geeignet als Profilfunktion ist eine Gaußfunktion, da diese die Eigenschaft hat, auch nach der Fourier-Rücktransformation die Gestalt einer Gaußfunktion beizubehalten. Eine Gewichtung der Bildsignalanteile durch Multiplikation der Frequenzraum-Intensitätsverteilung mit einer Gaußfunktion entspricht somit im Ortsraum einer Faltung der Intensitätsverteilung I(x, y) mit einer Gaußfunktion. Dies hat wiederum zur Folge, daß Stellen, an denen die Intensitätsverteilung sich abrupt ändert und die somit einen besonders hohen Kontrast haben, nach der Filterung nicht räumlich auseinanderzufließen scheinen.Out for this reason, the weighting of the image signal components is preferably carried out not only for individual discrete frequencies, but for frequency bands. Any frequency band that its weighting can be changed is determined using a profile function. Especially A Gaussian function is suitable as a profile function, since this is the property has, even after the Fourier inverse transformation the shape of a Gaussian function maintain. A weighting of the image signal components by multiplication the frequency space intensity distribution with a Gaussian function corresponds to a convolution of the intensity distribution in the local space I (x, y) with a Gaussian function. This in turn means that positions where the intensity distribution changes abruptly and thus have a particularly high contrast after filtering not spatially auseinanderzufließen seem to be.

7 zeigt die Frequenzraum-Intensitätsverteilung F(f) für eine beliebige, d.h. nicht aus kosinusförmigen Verteilungen zusammengesetzte, eindimensionale Ortsraum-Intensitätsverteilung. Gestrichelt eingezeichnet sind eine erste und eine zweite Profilfunktion g1(f) bzw. g2(f), bei denen es sich in beiden Fällen um durch die Gleichung gj(f) = exp(–(fZj – f)2/2wj 2) (2)gegebene Gaußfunktionen mit dem Zentralwert fZ1 bzw. fZ1 und dem Abstand w1 bzw. w2 zwischen dem Zentralwert und dem Wendepunkt handelt. Die Abstände w1 und w2 sind ein Maß für die Breite der glockenförmigen Profilfunktionen g1(f) und g2(f). Das Filter wirkt in diesem Beispiel derart, daß die Fourier-Amplituden von Frequenzen, die innerhalb der um den Zentralwert fZ1 liegenden Profilkurve g1(f) liegen, abgesenkt werden. Fourier-Amplituden von Frequenzen, die innerhalb der um den Zentralwert fZ2 liegenden Profilkurve g2(f) liegen, werden hingegen angehoben. 7 shows the frequency space intensity distribution F (f) for an arbitrary one-dimensional spatial space intensity distribution, that is to say that it is not composed of cosine-shaped distributions. A first and a second profile function g 1 (f) and g 2 (f), which are in both cases by the equation, are shown in dashed lines G j (f) = exp (- (f zj - f) 2 / 2w j 2 ) (2) given Gaussian functions with the central value f Z1 or f Z1 and the distance w 1 or w 2 between the central value and the turning point. The distances w 1 and w 2 are a measure of the width of the bell-shaped profile functions g 1 (f) and g 2 (f). In this example, the filter acts in such a way that the Fourier amplitudes of frequencies which lie within the profile curve g 1 (f) lying around the central value f Z1 are reduced. On the other hand, Fourier amplitudes of frequencies which lie within the profile curve g 2 (f) around the central value f Z2 are increased.

Im Einzelnen erfolgt die Filterung der Frequenzraum-Intensitätsverteilung F(f) dabei gemäß der Gleichung F'(f) = F(f)·TF(f), (3)wobei TF(f) eine Filterfunktion ist, die gegeben ist durch TF(f) = (1 + r1·g1(f))(1 + r2·g2(f)). (4) In particular, the frequency space intensity distribution F (f) is filtered according to the equation F '(f) = F (f) * T F (f), (3) where T F (f) is a filter function given by T F (f) = (1 + r 1 ·G 1 (f)) (1 + r 2 ·G 2 (F)). (4)

Die Verstärkungskoeffizienten r1 und r2 geben dabei an, wie stark die Fourier-Amplituden innerhalb der durch die Profilfunktionen vorgegebenen Frequenzbereiche verändert werden sollen. Im dargestellten Beispiel ist r1 > 0, da die Fourier-Amplituden um die kleinere Frequenz fZ1 herum angehoben werden sollen. Für den Verstärkungskoeffizienten r2 gilt hingegeben r2 < 0, was zu einer Verringerung der Fourier-Amplituden führt.The gain coefficients r 1 and r 2 indicate how strongly the Fourier amplitudes are to be changed within the frequency ranges specified by the profile functions. In the example shown, r 1 > 0, since the Fourier amplitudes are to be increased around the smaller frequency f Z1 . For the gain coefficient r 2 there applies r 2 <0, which leads to a reduction in the Fourier amplitudes.

Die Auswirkungen der in 7 dargestellten Filterung auf die Frequenzraum-Intensitätsverteilung F(f) sind in 8 gezeigt. Fourier-Amplituden für Frequenzen, die innerhalb der Filterprofile g1(f) und g2(f) liegen, sind durch die Filterung mit der Filterfunktion TF(f) abgesenkt bzw. ange hoben. Auf diese Weise treten großräumige Strukturen der Fingerknochen 14 auf dem im Bildschirm 24 gezeigten Bild gegenüber kleinräumigen Strukturen wie den Sehnen 18 zurück, so daß letztere besser für einen Arzt erkennbar werden.The impact of in 7 Filtering shown on the frequency space intensity distribution F (f) are in 8th shown. Fourier amplitudes for frequencies that lie within the filter profiles g 1 (f) and g 2 (f) are reduced or raised by the filtering with the filter function T F (f). In this way, large-scale structures of the finger bones occur 14 on the in the screen 24 shown image compared to small-scale structures such as the tendons 18 back, so that the latter can be better recognized by a doctor.

Aus den Gleichungen (2) bis (4) geht hervor, daß die Filterung der Frequenzraum-Intensitätsverteilungn F(f) mit Hilfe der Profilfunktionen in dem dargestellten Beispiel durch die Wertepaare (fZj, wj) bestimmt wird. Die Zentralwerte fZj sind dabei vorzugsweise wie auch in dem oben anhand der 3 bis 6 erläuterten Beispiel so zu wählen, daß die diesen Frequenzen fZj entsprechenden halben Periodenlängen etwa in der Größenordnung der Abmessungen derjenigen Strukturen liegen, die durch die Filterung bei der Bilddarstellung hervorgehoben oder abgeschwächt werden sollen. Sind diese typischen Strukturgrößen für alle denkbaren Anwendungen gleich, so können die Zentralwerte fZj und auch die Profilbreiten wj fest in der Bildverarbeitungsvorrichtung 10 vorgegeben sein. Vorzugsweise sind diese Parameter jedoch von dem Arzt mit Hilfe von an der Bildverarbeitungsvorrichtung 10 vorgesehenen Bedienelementen 26, 28 frei wählbar, um die Erkennbarkeit der Weichteilgewebestrukturen zu verbessern. Je nach Vergrößerung, Art und Anordnung der Knochen einerseits und der Art der zu untersuchenden Weichteilgewebestrukturen andererseits können die jeweils charakteristischen Abmessungen recht unterschiedlich sein, so daß eine Justierung während der Betrachtung des Bildschirms 24 in der Regel zu den besten Ergebnissen führen wird.From equations (2) to (4) it can be seen that the filtering of the frequency space intensity distribution F (f) is determined with the help of the profile functions in the example shown by the value pairs (f Zj , w j ). The central values f Zj are preferably, as in that above, based on the 3 to 6 to select the example explained so that the half-period lengths corresponding to these frequencies f Zj are approximately in the order of magnitude of the dimensions of those structures which are to be emphasized or weakened by the filtering in the image display. If these typical structure sizes are the same for all conceivable applications, the central values f Zj and also the profile widths w j can be fixed in the image processing device 10 be given. However, these parameters are preferably provided by the doctor with the aid of on the image processing device 10 provided controls 26 . 28 freely selectable to improve the recognizability of the soft tissue structures. Depending on the magnification, type and arrangement of the bones on the one hand and the type of soft tissue structures to be examined on the other hand, the respective characteristic dimensions can be quite different, so that an adjustment while viewing the screen 24 will usually lead to the best results.

Die in ihrer Gewichtung zu verändernden Bildsignalanteile können alternativ hierzu auch von der Bildverarbeitungsvorrichtung 20 selbständig im Wege eines adaptiven Verfahrens festgelegt werden. Hierzu ist durch einen zu behandelnden Arzt festzulegen, welche Weichteilgewebestrukturen besser erkennbar dargestellt werden sollen. Hierzu wählt der Arzt vorzugsweise einen Randbereich der betreffenden Weichteilgewebestruktur aus und markiert diesen. In 1 ist beispielhaft eine mit 30 bezeichnete Markierung gezeigt, die der Arzt z.B. mit einem Cursor auf einem zur Anzeige des Röntgenbildes 10 verwendeten Bildschirm erzeugen kann und in dem dargestellten Beispiel zwei Punkte und eine diese verbindende Linie 32 umfaßt.As an alternative to this, the weighting of the image signal components to be changed can also be obtained from the image processing device 20 be determined independently by means of an adaptive procedure. For this purpose, a doctor to be treated must determine which soft tissue structures should be shown more clearly. For this purpose, the doctor preferably selects and marks an edge area of the soft tissue structure in question. In 1 is an example with 30 designated marking shown that the doctor, for example, with a cursor on a to display the X-ray image 10 can use the screen used and in the example shown two points and a line connecting them 32 includes.

Die Bildverarbeitungseinrichtung 20 führt nun für eine Vielzahl von Frequenzbereichen die oben erläuterte Filterung aus und überprüft jeweils, inwieweit hierdurch der Kontrast entlang der Linie 32 zwischen den Punkten der Markierung 30 verbessert wird. Zur Anzeige gebracht wird dann die aus derjenigen Filterung gewonnene modifizierte Intensitätsverteilung, bei der der höchste Kontrast erzielt wurde.The image processing device 20 now carries out the filtering described above for a large number of frequency ranges and checks to what extent this results in the contrast along the line 32 between the points of the marking 30 is improved. The modified intensity distribution obtained from the filtering in which the highest contrast was achieved is then displayed.

Bei den bislang beschriebenen Beispielen wurde davon ausgegangen, daß durch die Filterung in lediglich einem Frequenzbereich die Fourier-Amplituden vergrößert und in lediglich einem Frequenzbereich die Fourier-Amplituden verkleinert werden. Zur Verbesserung der Erkennbarkeit kommt es jedoch nur auf das Amplitudenverhältnis an, so daß prin zipiell auch eine der genannten Maßnahmen zur Verbesserung der Erkennbarkeit genügt. Andererseits kann es auch zweckmäßig sein, in mehr als zwei Frequenzbereichen Fourier-Amplituden zu verändern, um die gewünschte Verbesserung der Darstellung zu erzielen. Für den Zähler j in den Gleichungen (2), (3) und (4) bedeutet dies, daß er nicht nur die Werte 1 und 2, sondern auch größere Werte annehmen kann.at the examples described so far have been assumed that by the filtering in only one frequency range the Fourier amplitudes enlarged and the Fourier amplitudes are reduced in only one frequency range become. However, it only improves the recognizability the amplitude ratio so that in principle also one of the measures mentioned suffice to improve recognizability. On the other hand, it can be appropriate change Fourier amplitudes in more than two frequency ranges to achieve the desired improvement to achieve the representation. For the counter j in equations (2), (3) and (4) this means that it is not can only assume the values 1 and 2, but also larger values.

Das vorstehend im Zusammenhang mit den 7 und 8 erläuterte Beispiel ist zur Vereinfachung der Darstellung auf eine eindimensionale Intensitätsverteilung I(x) beschränkt. Zur Bildverarbeitung müssen die oben angegebenen Gleichungen auf zwei Dimensionen erweitert werden. Die Gleichung (3) läßt sich dann unter Berücksichtung der Gleichungen (2) und (4) schreiben als

Figure 00160001
The above in connection with the 7 and 8th The example explained is restricted to a one-dimensional intensity distribution I (x) in order to simplify the illustration. For image processing, the equations given above must be expanded to two dimensions. Equation (3) can then be written taking into account equations (2) and (4) as
Figure 00160001

Die Filterfunktion Tf(fx, fy) kann außerdem noch mit einer weiteren Profilfunktion multipliziert werden, um den Signal-Rausch-Abstand zu verbessern. Handelt es sich bei dieser Profilfunktion z.B. um eine Gaußfunktion mit der Zentralfrequenz fZ = 0 und einer Breite w, die der Frequenz entspricht, bei der der hochfrequente Bildsignalanteil im Rauschen untergeht, so werden die Bildsignalanteile gegenüber einem Untergrundrauschen verstärkt. Diese Wahl der Profilfunktion trägt der Tatsache Rechnung, daß bei den in der Praxis häufig darzustellenden Bildern die Fourier- Amplituden mit zunehmender Frequenz f abnehmen, so daß bei hohen Frequenzen meist ein stets vorhandenes Rauschsignal überwiegt.The filter function T f (f x , f y ) can also be multiplied by another profile function in order to improve the signal-to-noise ratio. If this profile function is, for example, a Gaussian function with the central frequency f Z = 0 and a width w which corresponds to the frequency at which the high-frequency image signal component is lost in noise, the image signal components are amplified against background noise. This choice of the profile function takes into account the fact that the Fourier amplitudes decrease with increasing frequency f in the images that are frequently to be displayed in practice, so that at high frequencies an always present noise signal predominates.

Es versteht sich, daß die vorstehenden Ausführungen und Erläuterungen eines Ausführungsbeispiels nur beispielhaft und insbesondere nicht auf die Verbesserung der Erkennbarkeit von Weichteilgewebestrukturen beschränkt sind. Wie eingangs erläutert, können in der vorstehend beschriebenen Weise nicht nur Weichteilgewebestrukturen von Hartgewebestrukturen wie etwa Knochen, sondern auch gleichartige Gewebestrukturen besser voneinander unterschieden werden, sofern sie sich in ihrer Größe voneinander unterscheiden.It it is understood that the above statements and explanations one embodiment only exemplary and in particular not on the improvement of recognizability are limited by soft tissue structures. As explained at the beginning, in in the manner described above, not just soft tissue structures of hard tissue structures such as bones, but also similar ones Tissue structures can be better distinguished from one another, provided they differ in size differ.

Claims (14)

Verfahren zur Verbesserung der Erkennbarkeit von unterschiedlichen Strukturen (18) auf Durchstrahlungsbildern (10) unter Verwendung einer Bildverarbeitungsvorrichtung (20), umfassend die folgenden Schritte: a) Abspeichern eines in elektronischer Form vorliegenden Durchstrahlungsbildes (10) als Ortsraum-Intensitätsverteilung; b) Durchführen einer Fourier-Transformation zur Ermittlung einer Frequenzraum-Intensitätsverteilung; c) Filtern der Frequenzraum-Intensitätsverteilung durch Verändern der Gewichtung zwischen hochfrequenten und niedrigfrequenten Bildsignalanteilen, wobei die Festlegung der stärker zu gewichtenden Bildsignalanteile unter Berücksichtigung einer mittleren Strukturgröße der Strukturen erfolgt, deren Erkennbarkeit verbessert werden soll; d) Durchführen einer Rück-Fourier-Transformation der gefilterten Frequenzraum-Intensitätsverteilung, wodurch eine modifizierte Ortsraum-Intensitätsverteilung erhalten wird, in der diese Strukturen besser erkennbar sind.Methods to improve the recognizability of different structures ( 18 ) on radiographs ( 10 ) using an image processing device ( 20 ), comprising the following steps: a) storing a radiographic image available in electronic form ( 10 ) as spatial intensity distribution; b) performing a Fourier transformation to determine a frequency space intensity distribution; c) filtering the frequency space intensity distribution by changing the weighting between high-frequency and low-frequency image signal components, the definition of the image signal components to be weighted taking into account an average structure size of the structures, the recognizability of which is to be improved; d) performing a back Fourier transformation of the filtered frequency-space intensity distribution, as a result of which a modified spatial-space intensity distribution is obtained in which these structures can be better recognized. Verfahren nach Anspruch 1, bei dem die Festlegung der stärker zu gewichtenden Bildsignalanteile derart erfolgt, daß die diesen Frequenzen entsprechenden Periodenlängen etwa annähernd doppelt so groß sind wie die mittlere Strukturgröße der Strukturen erfolgt, deren Erkennbarkeit verbessert werden soll.The method of claim 1, wherein the determining the stronger to be weighted image signal components in such a way that these Periods corresponding to frequencies approximately twice as long are so big like the average structure size of the structures the recognizability of which is to be improved. Verfahren nach Anspruch 1 oder 2, bei dem die gefilterte Frequenzraum-Intensitätsverteilung F'(fx, fy) gegeben ist durch F'(fx, fy) = TF(fx, fy)·F(fx, fy),wobei F(fx, fy) die Frequenzraum-Intensitätsverteilung der Ortsraum-Intensitätsverteilung, fx, fy Frequenzen im zweidimensionalen Frequenzraum und TF(fx, fy) eine Filterfunktion zur Gewichtung von Bildsignalanteilen sind.Method according to Claim 1 or 2, in which the filtered frequency space intensity distribution F '(f x , f y ) is given by F '(f x , f y ) = T F (f x , f y ) · F (f x , f y ) where F (f x , fy) the frequency space intensity distribution of the spatial spatial intensity distribution, f x , f y frequencies in the two-dimensional frequency space and T F (f x , f y ) are a filter function for weighting image signal components. Verfahren nach einem der vorhergehenden Ansprüche, bei dem die Festlegung der Filterfunktion durch wenigstens einen Frequenz-Zentralwert und wenigstens eine Profilfunktion erfolgt, die die Gewichtung der Bildsignalanteile in Abhängigkeit von dem Abstand zu dem Frequenz-Zentralwert verändert.Method according to one of the preceding claims, which the definition of the filter function by at least one frequency central value and at least one profile function takes place, which weights the image signal components dependent on from the distance to the frequency center value changed. Verfahren nach den Ansprüchen 3 und 4, bei dem die Filterfunktion TF(fx, fy) gegeben ist durch
Figure 00200001
wobei g(fx, fy, fZj, wj) die Profilfunktion, rj ein Gewichtungsfaktor für den Frequenz-Zentralwert fZj und wj ein Breitenparameter ist, der ein Maß für die Breite der Profilfunktion ist.Method according to claims 3 and 4, in which the filter function T F (f x , f y ) is given by
Figure 00200001
where g (f x , f y , f Zj , w j ) is the profile function, r j is a weighting factor for the frequency central value f Zj and w j is a width parameter which is a measure of the width of the profile function. Verfahren nach Anspruch 5, bei dem die Profilfunktion g(fx, fy, fZj, wj) eine durch
Figure 00200002
gegebene Gaußfunktion ist.A method according to claim 5, wherein the profile function g (f x , f y , f Zj , w j ) one by
Figure 00200002
given Gaussian function. Verfahren nach einem der Ansprüche 4 bis 6, bei dem die Festlegung des wenigstens einen Frequenz-Zentralwerts und der wenigstens einen Profilfunktion mit Hilfe von an der Bildverarbeitungseinrichtung (20) vorgesehenen Stellgliedern (26, 28) veränderbar ist.Method according to one of claims 4 to 6, in which the determination of the at least one frequency central value and the at least one profile function with the aid of on the image processing device ( 20 ) provided actuators ( 26 . 28 ) is changeable. Verfahren nach einem der Ansprüche 4 bis 6, bei dem nach einer von einer Bedienperson durchgeführten Auswahl einer Struktur, deren Erkennbarkeit verbessert werden soll, auf dem Durchleuchtungsbild der wenigstens eine Frequenz-Zentralwert und die wenigstens eine Profilfunktion so adaptiv festgelegt wird, daß der Kontrast dieser Struktur erhöht wird.Method according to one of claims 4 to 6, in which according to performed by an operator Choosing a structure whose visibility should be improved the at least one central frequency value on the fluoroscopic image and the at least one profile function is so adaptively determined that the contrast this structure increases becomes. Verfahren nach Anspruch 8, bei dem die Auswahl der Struktur durch Vorgabe eines Punkts auf der Begrenzung der Struktur und einer Richtung erfolgt, entlang derer der Kontrast erhöht werden soll.The method of claim 8, wherein the selection of the Structure by specifying a point on the boundary of the structure and in a direction along which the contrast is increased should. Verfahren nach Anspruch 8, bei dem die Auswahl der Struktur durch Vorgabe zweier Punkte innerhalb der Struktur erfolgt, zwischen denen der Kontrast erhöht werden soll.The method of claim 8, wherein the selection of the Structure is made by specifying two points within the structure, between which the contrast increases shall be. Verfahren nach einem der vorhergehenden Ansprüche, bei dem die Frequenzraum-Intentsitätsverteilung F(fx, fy) zusätzlich einer Hochfrequenzfilterung unterzogen wird.Method according to one of the preceding claims, in which the frequency space intensity distribution F (f x , f y ) is additionally subjected to high-frequency filtering. Verfahren nach Anspruch 11, bei dem das Hochfrequenzfilter durch ein Gaußfilter mit dem Frequenz-Zentralwert 0 gegeben ist.The method of claim 11, wherein the high frequency filter through a Gaussian filter with the frequency central value 0 is given. Bildverarbeitungsvorrichtung zur Verbesserung der Erkennbarkeit von unterschiedlichen Strukturen (18) auf Durchstrahlungsbildern (10), mit: a) einem Speicher (MEM) zum Abspeichern eines in elektronischer Form vorliegenden Durchstrahlungsbildes (10) als Ortsraum-Intensitätsverteilung; b) einer Fourier-Transformationseinheit (FT) zum Durchführen einer Fourier-Transformation zur Ermittlung einer Frequenzraum-Intensitätsverteilung; c) einem Filter (FIL) zum Filtern der Frequenzraum-Intensitätsverteilung durch Verändern der Gewichtung zwischen hochfrequenten und niedrigfrequenten Bildsignalanteilen, wobei die Festlegung der stärker zu gewichtenden Bildsignalanteile unter Berücksichtigung einer mittleren Strukturgröße der Strukturen erfolgt, deren Erkennbarkeit verbessert werden soll; d) einer Fourier-Rück-Transformationseinheit (FT–1) zum Durchführen einer Rück-Fourier-Transformation der gefilterten Frequenzraum-Intensitätsverteilung, wodurch eine modifizierte Ortsraum-Intensitätsverteilung erhalten wird, in der diese Strukturen besser erkennbar sind.Image processing device to improve the recognizability of different structures ( 18 ) on radiographs ( 10 ), with: a) a memory (MEM) for storing an electronically transmitted image ( 10 ) as spatial intensity distribution; b) a Fourier transformation unit (FT) for performing a Fourier transformation for determining a frequency space intensity distribution; c) a filter (FIL) for filtering the frequency space intensity distribution by changing the weighting between high-frequency and low-frequency image signal components, the definition of the image signal components to be weighted taking into account an average structure size of the structures whose detectability is to be improved; d) a Fourier-back transformation unit (FT- 1 ) for performing a back-Fourier transformation of the filtered frequency-space intensity distribution, whereby a modified spatial-space intensity distribution is obtained in which these structures can be better recognized. Bildverarbeitungsvorrichtung nach Anspruch 13, bei der die Festlegung der stärker zu gewichtenden Bildsignalanteile derart erfolgt, daß die diesen Frequenzen entsprechenden Periodenlängen zumindest annähernd doppelt so groß sind wie die mittlere Strukturgröße der Strukturen, deren Erkennbarkeit verbessert werden soll.An image processing apparatus according to claim 13, in who's laying down the stronger to be weighted image signal components in such a way that these Frequencies corresponding to period lengths at least approximately twice are so big like the average structure size of the structures, whose recognizability should be improved.
DE2003125632 2003-06-06 2003-06-06 X ray image structure recognition improvement process uses Fourier transform for user selected Gaussian or other weighting filter in frequency domain Ceased DE10325632A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
DE2003125632 DE10325632A1 (en) 2003-06-06 2003-06-06 X ray image structure recognition improvement process uses Fourier transform for user selected Gaussian or other weighting filter in frequency domain

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE2003125632 DE10325632A1 (en) 2003-06-06 2003-06-06 X ray image structure recognition improvement process uses Fourier transform for user selected Gaussian or other weighting filter in frequency domain

Publications (1)

Publication Number Publication Date
DE10325632A1 true DE10325632A1 (en) 2004-12-23

Family

ID=33482610

Family Applications (1)

Application Number Title Priority Date Filing Date
DE2003125632 Ceased DE10325632A1 (en) 2003-06-06 2003-06-06 X ray image structure recognition improvement process uses Fourier transform for user selected Gaussian or other weighting filter in frequency domain

Country Status (1)

Country Link
DE (1) DE10325632A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006010372A1 (en) * 2004-07-27 2006-02-02 Dürr Dental GmbH & Co. KG Method and device for improving perceptibility different structures on radiographs
US20100040301A1 (en) * 2006-07-20 2010-02-18 Bam Bundesanstalt Fuer Material Forschung Und-Pruefung Image processing method without parameters
CN101430787B (en) * 2008-11-14 2012-07-04 西安交通大学 Digital sternum heart calcification image intensification method
WO2017089215A1 (en) * 2015-11-26 2017-06-01 Koninklijke Philips N.V. Apparatus having a user interface for enhancing medical images

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006010372A1 (en) * 2004-07-27 2006-02-02 Dürr Dental GmbH & Co. KG Method and device for improving perceptibility different structures on radiographs
EA011562B1 (en) * 2004-07-27 2009-04-28 Дюрр Денталь Гмбх & Ко.Кг Method and advice for improving perceptibility different structures on radiographs
US8244019B2 (en) 2004-07-27 2012-08-14 Duerr Dental Gmbh & Co. Kg Method and device for improving perceptibility different structures on radiographs
US20100040301A1 (en) * 2006-07-20 2010-02-18 Bam Bundesanstalt Fuer Material Forschung Und-Pruefung Image processing method without parameters
CN101430787B (en) * 2008-11-14 2012-07-04 西安交通大学 Digital sternum heart calcification image intensification method
WO2017089215A1 (en) * 2015-11-26 2017-06-01 Koninklijke Philips N.V. Apparatus having a user interface for enhancing medical images
US10657632B2 (en) 2015-11-26 2020-05-19 Koninklijke Philips N.V. Apparatus having a user interface for enhancing medical images

Similar Documents

Publication Publication Date Title
DE60300462T2 (en) Method for sharpening a digital image with signal noise evaluation
DE69936036T2 (en) Method and apparatus for automatically compensating for time gain and / or page intensity gain in B-mode ultrasound imaging
DE60038382T2 (en) INTRAVASCULAR ULTRASONIC IMAGE ANALYSIS USING AN ACTIVE CONTOUR METHOD
DE69832357T2 (en) NOISE REDUCTION IN ONE IMAGE
DE69832412T2 (en) METHOD AND APPARATUS FOR DETECTING BLUTSPECKLE IN AN INTRAVASCULAR ULTRASOUND PICTURE SYSTEM
DE10312018B4 (en) System and method for ultrasound image processing
DE69937284T2 (en) Automatic regulation of speed and pulse repetition frequency for Doppler ultrasound spectrograms
DE2952422C3 (en) Method and device for processing an X-ray image in an X-ray image copying system
DE60300097T2 (en) Method for sharpening a digital image without amplification noise
DE19857526A1 (en) Removing stripe artifacts in medical imaging system
DE102006005803A1 (en) Method for noise reduction in imaging methods
EP0938063B1 (en) Method of two-dimensional imaging of structures for medical diagnosis
DE102007057013A1 (en) Digital mammogram e.g. raw mammogram, processing device for breast, involves optimizing contrast of pixel value of brightness optimized mammogram and/or sharpness of structure contained in displayed object
DE102010043975B4 (en) Procedure for reducing the radiation dose used as part of an imaging X-ray examination and computer system
DE102004026996A1 (en) Method for generating magnetic resonance recordings and control device for a magnetic resonance tomograph
DE102007013570A1 (en) Method for noise reduction in digital images with locally different and directional noise
DE102008023915A1 (en) Filtering control variable setting method for noise reduction in medical image in X-ray diagnostic facility, involves comparing standard noise deviation and standard structure deviation and setting control variable as function of comparison
DE10163215B4 (en) System and method with automatically optimized image generation
DE60202588T2 (en) Method for noise reduction
DE102007058498A1 (en) Method and apparatus for noise reduction in medical images
DE102004056589A1 (en) Method and device for performing segmentation-based image operations
WO2006010372A1 (en) Method and device for improving perceptibility different structures on radiographs
DE10325632A1 (en) X ray image structure recognition improvement process uses Fourier transform for user selected Gaussian or other weighting filter in frequency domain
DE60214967T2 (en) Method for improving the contrast of an image
DE10214114A1 (en) CT or X-ray medical image filtering method in which image integral SNRs are calculated and compared with calculated limit SNR values while carrying out iterative filtering and, if necessary, changing the filter dose

Legal Events

Date Code Title Description
8127 New person/name/address of the applicant

Owner name: DUERR DENTAL AG, 74321 BIETIGHEIM-BISSINGEN, DE

8110 Request for examination paragraph 44
R016 Response to examination communication
R002 Refusal decision in examination/registration proceedings
R003 Refusal decision now final