Hi,
ich sitze hier an einem recht komplexen Thema was meine Physik-Kenntnisseübersteigt. Ich hoffe jemand von euch kann mir dabei helfen.
Es geht dabei darum, die Qualität (Hochkontrast-Auflösung) eines 16bit Graustufenbildes zu ermitteln. Ein Beispiel eines solchen Bildes habe ichangehängt. Der weiße Punkt in dem Bild ist die Aufnahme eines sehr dünnen Drahtes mitder Dichte von 500HU (hounsfield Units) und drum herum ist Wasser mit der Dichte von 0. Der Punkt strahlt jetzt aus und daraus kann man die Qualität des Bildes errechnen.
Das geht, indem man zunächst die Point-Spread-Function ermittelt und daraus dann die Line-Spread-Function. Diese ergibt dann eine Kurve die ungefähr so aussieht wie Bild2 (Nur mal eben von Hand gezeichnet). Im Prinzip also eine ausschwingende Kurve (Überstrahlung des Punktes)
Aus diesen Werten kann man nun den MTF ermitteln, indem man eine FourierTransformation darauf anwendet. Das Ergebnis sollen Werte sein in dem Format lp/cm also LinePairs/cm, was der Auflösung des Bildes und damit der Qualität entspricht.
Da hört es bei mir aber auf. Ich bekomme einfach kein sinnvolles Ergebnis. Das wird daran liegen, dass ich die FFT nie richtig verstanden habe. Ich gehe da mit 2 Werte-Arrays rein. Mit einem Array der reellen Zahlen und einem der imaginären Zahlen und heraus bekomme ich wieder ein Array mit 2*N Werten mit denen ich aber einfach nichts anzufangen weis. Ich weiß auch nichtwas der imaginäre Anteil meiner Pixelwerte eigentlich sein soll. Alles lesen im Internet bringt mich auch nicht wirklich weiter. Ich verstehe zwar prinzipiell was FFT macht, wie man es praktisch anwendet kapier ich aber einfach nicht. Vieelicht stimmt auch schon miene Normierung der LSF Werte nicht und ich gebe der FFT falsche Werte rein?
Hat jemand da eine Idee oder hat sowas eventuell sogar schonmal gemacht oder kann mir einfach mal erklären wie das mit der FFT so richtig gemacht wird?
Danke
Claus
ich sitze hier an einem recht komplexen Thema was meine Physik-Kenntnisseübersteigt. Ich hoffe jemand von euch kann mir dabei helfen.
Es geht dabei darum, die Qualität (Hochkontrast-Auflösung) eines 16bit Graustufenbildes zu ermitteln. Ein Beispiel eines solchen Bildes habe ichangehängt. Der weiße Punkt in dem Bild ist die Aufnahme eines sehr dünnen Drahtes mitder Dichte von 500HU (hounsfield Units) und drum herum ist Wasser mit der Dichte von 0. Der Punkt strahlt jetzt aus und daraus kann man die Qualität des Bildes errechnen.
Das geht, indem man zunächst die Point-Spread-Function ermittelt und daraus dann die Line-Spread-Function. Diese ergibt dann eine Kurve die ungefähr so aussieht wie Bild2 (Nur mal eben von Hand gezeichnet). Im Prinzip also eine ausschwingende Kurve (Überstrahlung des Punktes)
Aus diesen Werten kann man nun den MTF ermitteln, indem man eine FourierTransformation darauf anwendet. Das Ergebnis sollen Werte sein in dem Format lp/cm also LinePairs/cm, was der Auflösung des Bildes und damit der Qualität entspricht.
Da hört es bei mir aber auf. Ich bekomme einfach kein sinnvolles Ergebnis. Das wird daran liegen, dass ich die FFT nie richtig verstanden habe. Ich gehe da mit 2 Werte-Arrays rein. Mit einem Array der reellen Zahlen und einem der imaginären Zahlen und heraus bekomme ich wieder ein Array mit 2*N Werten mit denen ich aber einfach nichts anzufangen weis. Ich weiß auch nichtwas der imaginäre Anteil meiner Pixelwerte eigentlich sein soll. Alles lesen im Internet bringt mich auch nicht wirklich weiter. Ich verstehe zwar prinzipiell was FFT macht, wie man es praktisch anwendet kapier ich aber einfach nicht. Vieelicht stimmt auch schon miene Normierung der LSF Werte nicht und ich gebe der FFT falsche Werte rein?
Hat jemand da eine Idee oder hat sowas eventuell sogar schonmal gemacht oder kann mir einfach mal erklären wie das mit der FFT so richtig gemacht wird?
Danke
Claus
2 Stunden Try & Error erspart 10 Minuten Handbuchlesen.
Pre-Kaffee-Posts sind mit Vorsicht zu geniessen
Pre-Kaffee-Posts sind mit Vorsicht zu geniessen