1 / 32

Ψηφιακή Επεξεργασία και Ανάλυση Εικόνας Ακαδημαϊκό Έτος 2007-8

Ψηφιακή Επεξεργασία και Ανάλυση Εικόνας Ακαδημαϊκό Έτος 2007-8. Παρουσίαση Νο. 3 Δισδιάστατα σήματα και συστήματα #2. Projection-Slice Theorem (1/ 4 ). Χρησιμοποιείται στην τομογραφία Απεικόνιση του εσωτερικού στερεών σωμάτων (της τομής f(x,y) ) με μη επεμβατικό τρόπο

huong
Download Presentation

Ψηφιακή Επεξεργασία και Ανάλυση Εικόνας Ακαδημαϊκό Έτος 2007-8

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Ψηφιακή Επεξεργασία και Ανάλυση ΕικόναςΑκαδημαϊκό Έτος 2007-8 Παρουσίαση Νο. 3 Δισδιάστατα σήματα και συστήματα #2

  2. Projection-Slice Theorem (1/4) • Χρησιμοποιείται στην τομογραφία • Απεικόνιση του εσωτερικού στερεών σωμάτων (της τομής f(x,y) )με μη επεμβατικό τρόπο • Θυμίζουμε ότι για την συνεχή φυσική σκηνή f(x,y) το ζεύγος Fourier είναι ΤΜΗΥΠ / ΕΕΣΤ

  3. Projection-Slice Theorem (2/4) • Κατασκευή της 1-D συνάρτησης pθ(t) • Κάθε τιμή της pθ(t)δίνεται από το άθροισμα των τιμών της τομής κατά μήκος των γραμμών που φαίνονται στο σχήμα ΤΜΗΥΠ / ΕΕΣΤ

  4. Projection-Slice Theorem (3/4) • Μετασχηματισμός Radon (προβολή της f(x,y) υπό γωνία θ ) Αλλαγή συντεταγμένων από το σύστημα t,u στο x,y • Ησυνάρτηση f(x,y) εξαρτάται από το υλικό της τομήςαφού κάθε υλικό παρουσιάζει διαφορετικές ιδίοτητες (π.χ., απορρόφηση σε κάποια ακτινοβολία) ΤΜΗΥΠ / ΕΕΣΤ

  5. Projection-Slice Theorem (4/4) • Θεώρημα Projection-Slice ΤΜΗΥΠ / ΕΕΣΤ

  6. Ανακατασκευή από προβολές 1/2 • Για την πλήρη ανακατασκευή του στερεού σώματος • Εκτιμάται η F(Ω1,Ω2) και κατ’ επέκταση η f(x,y) από τον υπολογισμό της pθ(t) για διάφορες γωνίες στο διάστημα [0 π] • Εφαρμόζεται η παραπάνω διαδικασία για διάφορα επίπεδα (τομές) και στη συνέχεια γίνεται η κατάλληλη αντιστοίχησή τους ΤΜΗΥΠ / ΕΕΣΤ

  7. Ανακατασκευή από προβολές 2/2 • Ns : αριθμός εκπεμπομένων φωτονίων ανά χρονική μονάδα • Nr : αριθμός λαμβανομένων φωτονίων ανά χρονική μονάδα • Σύμφωνα με το νόμο Lambert-Beer : • Όπου η α(d)είναι ησταθερά εκθετικής απόσβεσης κατά μήκος της διαδρομής d και δίνεται ως Αν πάρουμε φυσικό λογάριθμο και στις δύο πλευρές του Νόμου Lambert-Beer βρίσκουμε την απορροφητικότητα δια μέσου της συγκεκριμένης διαδρομής: ΤΜΗΥΠ / ΕΕΣΤ

  8. Σχέση Χώρου – Συχνότητας • F(Ω1,Ω2)=FT[f(t1,t2)] f(t1,t2): απόλυτα ολοκληρώσιμη • Εάν f(t1,t2):ζωνοπεριορισμένη  f(n1,n2) • F(ω1,ω2)=“DS”FT[f(n1,n2)] (discrete space) • Εάν f(n2,n2) περιορισμένης χωρικής επέκτασης  F(ω1,ω2) (δειγματοληψία του DSFT) • F(ω1,ω2) = DFT[f(n1,n2)] ΤΜΗΥΠ / ΕΕΣΤ

  9. Περιγραφή συχνοτικού περιεχομένου (1/5) • To ζεύγος DFT για 2D σήματα είναι (με συμβολισμό πιο συνήθη σε ΨΕΑΕ): • Εμφάνιση συχνοτήτων σε μονοδιάστατα σήματα ΤΜΗΥΠ / ΕΕΣΤ

  10. Περιγραφή συχνοτικού περιεχομένου (2/5) • Όμοια για τα δισδιάστατα σήματα • Η ολίσθηση είναι της μορφής F(u,v)F(u-N/2,v-N/2) • Χαμηλές συχνότητες • Ομαλές περιοχές στην εικόνα • Υψηλές συχνότητες • Περιγράμματα και άκρα εικόνας λόγω παραθύρωσης ΤΜΗΥΠ / ΕΕΣΤ

  11. Περιγραφή συχνοτικού περιεχομένου (3/5) v u • Μετασχηματισμός Fourier (στο σύστημα αξόνων u-v), p0(t) (κατά μήκος του άξονα v), pπ/2(t) (κατά μήκος του άξονα u) • Η υψηλή ενέργεια που εμφανίζεται στις υψηλές συχνότητες κατά μήκος των αξόνων δικαιολογείται από τις έντονες διακυμάνσεις και τυχόν ασυνέχειες κατά μήκος του άξονα t των καμπυλών pθ(t) ΤΜΗΥΠ / ΕΕΣΤ

  12. Περιγραφή συχνοτικού περιεχομένου (4/5) • Μετασχηματισμός Fourier, p3π/4(t) (κατά μήκος της v=-u), pπ/4(t) (κατά μήκος της v=u) • Ενώ στο άξονα v=-u υπάρχει ενέργεια στις υψηλές συχνότητες (η αντίστοιχη καμπύλη μεταβάλλεται γρήγορα) δεν συμβαίνει το ίδιο με τον άξονα v=u (η αντίστοιχη καμπύλη είναι σαφώς ομαλότερη) ΤΜΗΥΠ / ΕΕΣΤ

  13. Περιγραφή συχνοτικού περιεχομένου (5/5) • Άλλα παραδείγματα: ΤΜΗΥΠ / ΕΕΣΤ

  14. Πληροφορία πλάτους-φάσης (1/3) • Ο μετασχηματισμός Fourier διακριτού χρόνου είναι μιγαδική ποσότητα και έτσι μπορεί να γραφεί σε πολική μορφή • Η περισσότερη πληροφορία υπάρχει στην φάση Φx() ΤΜΗΥΠ / ΕΕΣΤ

  15. Πληροφορία πλάτους-φάσης (2/3) • Έστω ηεικόνα x με F{x}=MxejΦxκαι δημιουργώ τις εικόνες xp, xmώστε xp=F-1{cejΦx}, xm=F-1{Μxejf}, όπου c, f σταθερές • Εικόνα x, εικόνα xpμε c=1, εικόνα xmμε f=0,4π ΤΜΗΥΠ / ΕΕΣΤ

  16. Πληροφορία πλάτους-φάσης (3/3) • Έστω οι εικόνες x,y με F{x}=MxejΦx, F{y}=MyejΦyκαι δημιουργώ τις εικόνες z, w ώστε z=F-1{MyejΦx}, w=F-1{ΜxejΦy} Εικόνα x Εικόνα yΕικόνα zΕικόνα w ΤΜΗΥΠ / ΕΕΣΤ

  17. Το πρόβλημα της ανάκτησης πλάτους (1/4) • Σύμφωνα με αυτό γίνεται προσπάθεια ανακατασκευής ενός σήματος x(n1,n2)μόνο από την φάση Φxτου μετασχηματισμού Fourier • Το σήμα x'(n1,n2)που προκύπτει ισούται με αx(n1,n2) • Περιορισμοί • Το x(n1,n2) πραγματικό και πεπερασμένο • Ο μετασχηματισμός Fourier δεν παραγοντοποιείται σε γινόμενο πολυωνύμων των ejω1και ejω2 • Στις εικόνες ισχύουν οι παραπάνω περιορισμοί ΤΜΗΥΠ / ΕΕΣΤ

  18. Το πρόβλημα της ανάκτησης πλάτους (2/4) • Δυο λύσεις για το πρόβλημα • Μέσω γραμμικού συστήματος • Μέσω επαναληπτικής μεθόδου • Γραμμικό σύστημα • Αν Ν2 άγνωστοι τότε δειγματοληπτούμε σεΝ2 σημεία (ω1,ω2) • Το σύστημα που προκύπτει είναι της μορφής Ax=0 οπότε η λύση είναι της μορφής αx(n1,n2) ΤΜΗΥΠ / ΕΕΣΤ

  19. Το πρόβλημα της ανάκτησης πλάτους (3/4) • Η μέθοδος γραμμικού συστήματος είναι κατάλληλη για μικρές εικόνες αφού για παράδειγμα αν έχουμε εικόνα 512x512 τότε απαιτούνται 262144 εξισώσεις!!! • Επαναληπτική μέθοδος • Θέτει περιορισμούς • Στον χωρικό πεδίο: πραγματικές τιμές και τα εκτός εικόνας εικονοστοιχεία μηδέν • Στο συχνοτικό πεδίο: ο μετασχηματισμός Fourier δεν πρέπει να παραγοντοποιείται • Κριτήριο σταματήματος: Η φάση του αποτελέσματος πρέπει • να είναι ίση (ή ~ = ) με αυτή της αρχικής εικόνας ΤΜΗΥΠ / ΕΕΣΤ

  20. Το πρόβλημα της ανάκτησης πλάτους (4/4) • Επαναληπτική διαδικασία ΤΜΗΥΠ / ΕΕΣΤ

  21. Δειγματοληψία (1/5) • Η εικόνα ψηφιοποιείται για να γίνει η επεξεργασία της από υπολογιστή • Η δειγματοληψία αναφέρεται στην διακριτοποίηση των αξόνων x,y της εικόνας (δηλ. του πεδίου ορισμού) • Από το σήμα xα(t1,t2) λαμβάνεται το σήμα x(n1,n2) • Ομοιόμορφη δειγματοληψία (εύκολη υλοποίηση) • Ανομοιόμορφη δειγματοληψία (π.χ., περισσότερα δείγματα στα περιγράμματα) • Δυο ενδιαφέροντα προβλήματα • Σχέση φασμάτων συνεχούς και διακριτού σήματος • Ανακατασκευή συνεχούς από διακριτό σήμα ΤΜΗΥΠ / ΕΕΣΤ

  22. Δειγματοληψία (2/5) • Πλέγμα δειγματοληψίας • t=Vn, t=[t1,t2]T, n=[n1, n2]T, V=[vij]2x2 • x(n)=xα(Vn) • Ορθογώνιο πλέγμα • Εξαγωνικό πλέγμα (βέλτιστο για κυκλικά φάσματα) ΤΜΗΥΠ / ΕΕΣΤ

  23. Δειγματοληψία (3/5) • Σχέση φασμάτων (FT – DSFT) • x(n)=xα(Vn) • Σχέση φασμάτωνμε χρήση ορθογώνιου πλέγματος • x(n1,n2)=xα(n1T1,n2T2) , T=T1xT2 (χωρική περίοδος) ΤΜΗΥΠ / ΕΕΣΤ

  24. Δειγματοληψία (4/5) • Παράδειγμα συνεχούς και διακριτού φάσματος ΤΜΗΥΠ / ΕΕΣΤ

  25. Δειγματοληψία (5/5) • Ανακατασκευή αναλογικού σήματος • Εφαρμογή ορίου Nyquist σε κάθε συνιστώσα κατά την δειγματοληψία για να αποφευχθεί η αναδίπλωση συχνότητας ΤΜΗΥΠ / ΕΕΣΤ

  26. Κβαντισμός (1/3) • Ο κβαντισμός αναφέρεται στη διακριτοποίηση της τιμής f του κάθε εικονοστοιχείου (δηλαδή του πεδίου τιμών) • Διαδικασία κβαντισμού • To πεδίο τιμών της fχωρίζεται σε L υποδιαστήματα • Οι ακραίες τιμές κάθε υποδιαστήματος είναι τα όρια απόφασης di-1, di και μια τιμή μεταξύ αυτών ονομάζεται επίπεδο κβάντισης ri • Αν di-1< f≤ di τότε το fq=Q(f)=ri ΤΜΗΥΠ / ΕΕΣΤ

  27. Κβαντισμός (2/3) • Θόρυβος κβαντισμού: f=fq+eq • To eqεξαρτάται από το σήμακαι το L • Ομοιόμορφος κβαντισμός • di-di-1=Δ,ri=(di+di-1)/2, 1≤i≤L, Δ=(dL-d0)/2(# ψηφίων) ΤΜΗΥΠ / ΕΕΣΤ

  28. Παράδειγμα κβαντισμού (3/3) • Επανακβαντισμός εύρους τιμών [0:255] στα πέντε επίπεδα [0:51:255] και αντιστοίχιση στην κεντρική τιμή • Παρατηρήστε ότι στην νέα εικόνα εμφανίζονται μόνο πέντε επίπεδα του γκρι ΤΜΗΥΠ / ΕΕΣΤ

  29. Διανυσματικός κβαντισμός (1/4) • Σε αυτή την περίπτωση οι ποσότητες που κβαντίζονται είναι Ν-διάστατα διανύσματα που παράγονται είτε από μια διανυσματική πηγή είτε μετά από ομαδοποίηση βαθμωτών ποσοτήτων • Δημιουργείται ένας Ν-διάστατος χώρος • Διαδικασία κβαντισμού (δοθέντος του διαν. κβαντιστή) • Ο χώρος είναι διαχωρισμένος σε L υποπεριοχές Ri • Τα όρια αυτών των υποπεριοχών είναι τα όρια απόφασηςκαι επιλέγεται κάποιο διάνυσμα εντός των ορίων που ονομάζεται επίπεδο κβάντισης ri • Αν το διάνυσμα g ανήκει στην περιοχή Ri τότε το gq==Q(g)=riκαι g=gq+eq ΤΜΗΥΠ / ΕΕΣΤ

  30. Παράδειγμα διανυσματικού κβαντισμού (2/4) • Έστω g=[g1 g2]T, 0≤g1,g2≤1 και L=4 • Οι περιοχές μπορούν να έχουν οποιοδήποτε σχήμα και τα διανύσματα ri μπορεί να μην είναι στο κέντρο κάθε περιοχής ΤΜΗΥΠ / ΕΕΣΤ

  31. Διανυσματικός κβαντισμός (3/4) • Επιλογή Ri, ri • Βέλτιστη επιλογή με βάση κάποιο κριτήριο (π.χ. Ευκλείδιο απόσταση) • D=E[d(g,gq)], d(g,gq)=(gq-g)T(gq-g) • Για τον βέλτιστο διανυσματικό κβαντιστή πρέπει να ισχύουν οι εξής ιδιότητες • (1)gq=Q(g)=ri, αν και μόνο αν d(g,ri)≤d(g,rj), j≠i, 1≤j≤L • (2) To riπρέπει να είναι κεντροειδές, δηλαδή ΤΜΗΥΠ / ΕΕΣΤ

  32. Διανυσματικός κβαντισμός (4/4) • Ο καθορισμός των ri μπορεί να γίνει με τον αλγόριθμο Κ-means • (1) Επιλέγονται τυχαία M (Μ>>L) διανύσματα gi (διανύσματα εκπαίδευσης) • (2) Από αυτά επιλέγονται L και θεωρούνται τα κεντροειδή διανύσματα. • (3) Χρησιμοποιώντας το MSE D (MSE Distance) κατηγοριοποιούνται τα διανύσματα στις L υποπεριοχές(ιδιότητα 1 του βέλτιστου κβαντιστή) • (4) Για κάθε υποπεριοχή υπολογίζονται τα νέα κεντροειδή • Όταν χρησιμοποιείται το MSE D τότε τα κεντροειδή υπολογίζονται από τον μέσο όρο των διανυσμάτων κάθε περιοχής • (5) Επιστροφή στο βήμα 3 χρησιμοποιώντας να νέα κεντροειδή – αυτή η διαδικασία συνεχίζεται έως ότου τα κεντροειδή παραμείνουν τα ίδια ΤΜΗΥΠ / ΕΕΣΤ

More Related