[ΨΕΕ] 2η Εργασία 2013

[Thmmy.gr portal]

  Thmmy.gr portal
   Forum
   Downloads
   Ενεργ. Λογαριασμού
   Επικοινωνία
  
  Χρήσιμα links
   Σελίδα τμήματος
   Βιβλιοθήκη Τμήματος
   Elearning
   Φοιτητικά fora
   Πρόγραμμα Λέσχης
   Πρακτική Άσκηση
   Ηλεκτρονική Εξυπηρέτηση Φοιτητών
   Διανομή Συγγραμμάτων
   Ψηφιακό Καταθετήριο Διπλωματικών
   Πληροφορίες Καθηγητών
   Instagram @thmmy.gr
   mTHMMY
  
  Φοιτητικές Ομάδες
   ACM
   Aristurtle
   ART
   ASAT
   BEAM
   BEST Thessaloniki
   EESTEC LC Thessaloniki
   EΜΒ Auth
   IAESTE Thessaloniki
   IEEE φοιτητικό παράρτημα ΑΠΘ
   SpaceDot
   VROOM
   Panther
  

THMMY.gr Web

Helios_MultiScribbles2SMF Default ThemeSMFone_Blue

[Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads]

Είσαι πρωτοετής;... Καλώς ήρθες! Μπορείς να βρεις πληροφορίες εδώ. Βοήθεια για τους καινούργιους μέσω χάρτη.

Κατεβάστε εδώ το Android Application για εύκολη πρόσβαση στο forum.

Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads με προσοχή στα ονόματα των αρχείων!

[Πληροφορίες Καθηγητών]

Πληροφορίες Καθηγητών

[MARIOS]

1ον) Τρέχω τις εντολές που δίνει για να βάλουμε θόρυβο στην εικόνα και το αποτέλεσμα το βάζω στην imshow και μου δίνει μια πολύ καλή εικόνα.... :/ τι σόι θόρυβος είναι αυτός?

2ον) Πως θα γίνει η εξομάλυνση στην δεύτερη συνάρτηση?

3oν)Το μέτρο της Η    ||Η(ω1,ω2)||^2   δεν είναι η ορίζουσα του πίνακα Η?

[mitsos_dlx]

1) μια πολύ χάλια εικόνα πρέπει να δίνει, σαν αυτή του σχήματος 2 στην εκφώνηση. Την y να προβάλεις, όχι την x.
2) μπορείς να περάσεις μια μάσκα μέσου όρου της γύρω γειτονίας ή να χρησιμοποιήσεις την συνάρτηση medfilt2 που ανέφερε ο ντελόπουλος και κάνει κάτι παρόμοιο. Αυτό για να διώξεις υποτίθεται το θόρυβο. πριν κόψεις το κομμάτι που θέλεις (ή και μετά το κόψιμο αλλά θα παραμορφωθούν λίγο τα άκρα αυτού που έκοψες). Ωστόσο σε εμένα τουλάχιστον τα αποτελέσματα δεν βελτιώνονται καθόλου με αυτό
3) η ορίζουσα από που προκύπτει; νομίζω είναι το abs(H) ή το H.*conj(H) με το πολλαπλασιασμό στοιχείο-στοιχείο.


Τώρα όσο αφορά τα αποτελέσματα είναι απαράδεκτα και προφανώς φταίει η κακή εκτίμηση της h. Δηλαδή από ελάχιστη ως μηδενική βελτίωση σε σχέση με τη y. Επίσης το K σε εμένα ότι και να το βάλω περίπου ίδιο αποτέλεσμα παίρνω ... οπτικά τουλάχιστον, αλλά ίσως έχω κάποιο λάθος.
Και ρωτώ τώρα πως είναι δυνατόν να εκτιμήσουμε καλά την h όταν:
1) πρέπει (ιδανικά) να πάρουμε τη παραμόρφωση που προκαλεί ένα αστέρι του ενός pixel με γύρω-γύρω τελείως μαύρο και ομοιόμορφο φόντο και τέτοιο αστέρι δεν υπάρχει ούτε στη αρχική εικόνα x.
2) Ακόμα και να υπήρχε, στη παραμορφωμένη εικόνα δεν θα δεις τίποτα σε εκείνο το σημείο.
3) Ακόμα και αν μπορούσες να το δεις με κάποιο τρόπο τη παραμορφωμένη έκδοση αυτού του πολύ μικρού αστεριού, για να οφείλεται αυτό μόνο στο συγκεκριμένο άστρο θα έπρεπε η μάσκα παραμόρφωσης να είναι μικρή (και όχι 50x50 για μία εικόνα που κάθε 5-10 pixel έχει και ένα αστεράκι), ώστε η παραμόρφωση να μη επηρεάζεται από γειτονικά αστέρια.

Δεν ξέρω αν μπορεί να γίνει κάτι αλλά θεωρητικά τουλάχιστον η κατάσταση είναι πολύ δύσκολη...

[MARIOS]

1) μια πολύ χάλια εικόνα πρέπει να δίνει, σαν αυτή του σχήματος 2 στην εκφώνηση. Την y να προβάλεις, όχι την x.
2) μπορείς να περάσεις μια μάσκα μέσου όρου της γύρω γειτονίας ή να χρησιμοποιήσεις την συνάρτηση medfilt2 που ανέφερε ο ντελόπουλος και κάνει κάτι παρόμοιο. Αυτό για να διώξεις υποτίθεται το θόρυβο. πριν κόψεις το κομμάτι που θέλεις (ή και μετά το κόψιμο αλλά θα παραμορφωθούν λίγο τα άκρα αυτού που έκοψες). Ωστόσο σε εμένα τουλάχιστον τα αποτελέσματα δεν βελτιώνονται καθόλου με αυτό
3) η ορίζουσα από που προκύπτει; νομίζω είναι το abs(H) ή το H.*conj(H) με το πολλαπλασιασμό στοιχείο-στοιχείο.

1) ok....είχα κάνει ένα λάθος στην imresize όταν την πέρασα στο matlab....
2)θα το κοιτάξω
3)g=h'/(abs(h)^2+1/50)
αυτό όταν πάει να τρέξει μη τετραγωνικό πίνακα μου βγάζει error....

[MARIOS]

thanks και για την απάντηση σου φίλε μου!!!! 

[MARIOS]

1) μια πολύ χάλια εικόνα πρέπει να δίνει, σαν αυτή του σχήματος 2 στην εκφώνηση. Την y να προβάλεις, όχι την x.
2) μπορείς να περάσεις μια μάσκα μέσου όρου της γύρω γειτονίας ή να χρησιμοποιήσεις την συνάρτηση medfilt2 που ανέφερε ο ντελόπουλος και κάνει κάτι παρόμοιο. Αυτό για να διώξεις υποτίθεται το θόρυβο. πριν κόψεις το κομμάτι που θέλεις (ή και μετά το κόψιμο αλλά θα παραμορφωθούν λίγο τα άκρα αυτού που έκοψες). Ωστόσο σε εμένα τουλάχιστον τα αποτελέσματα δεν βελτιώνονται καθόλου με αυτό
3) η ορίζουσα από που προκύπτει; νομίζω είναι το abs(H) ή το H.*conj(H) με το πολλαπλασιασμό στοιχείο-στοιχείο.


Τώρα όσο αφορά τα αποτελέσματα είναι απαράδεκτα και προφανώς φταίει η κακή εκτίμηση της h. Δηλαδή από ελάχιστη ως μηδενική βελτίωση σε σχέση με τη y. Επίσης το K σε εμένα ότι και να το βάλω περίπου ίδιο αποτέλεσμα παίρνω ... οπτικά τουλάχιστον, αλλά ίσως έχω κάποιο λάθος.
Και ρωτώ τώρα πως είναι δυνατόν να εκτιμήσουμε καλά την h όταν:
1) πρέπει (ιδανικά) να πάρουμε τη παραμόρφωση που προκαλεί ένα αστέρι του ενός pixel με γύρω-γύρω τελείως μαύρο και ομοιόμορφο φόντο και τέτοιο αστέρι δεν υπάρχει ούτε στη αρχική εικόνα x.
2) Ακόμα και να υπήρχε, στη παραμορφωμένη εικόνα δεν θα δεις τίποτα σε εκείνο το σημείο.
3) Ακόμα και αν μπορούσες να το δεις με κάποιο τρόπο τη παραμορφωμένη έκδοση αυτού του πολύ μικρού αστεριού, για να οφείλεται αυτό μόνο στο συγκεκριμένο άστρο θα έπρεπε η μάσκα παραμόρφωσης να είναι μικρή (και όχι 50x50 για μία εικόνα που κάθε 5-10 pixel έχει και ένα αστεράκι), ώστε η παραμόρφωση να μη επηρεάζεται από γειτονικά αστέρια.

Δεν ξέρω αν μπορεί να γίνει κάτι αλλά θεωρητικά τουλάχιστον η κατάσταση είναι πολύ δύσκολη...

Τώρα για αυτό που λες....εννοείτε ότι έχεις δίκιο και εννοείτε ότι το κάνεις όχι ιδανικά....
εγώ όταν έκανα κάποιες δοκιμές sto command window για να πάρω μια σχετικά μικρή h άρα και πιο κοντά στο ιδανικό, χρησιμοποίησα την imcrop δύο φορές σε σειρά....
h1=imcrop(y);--> διαλέγεις ένα μέρος της εικόνας που θεωρείς ικανοποιητικό
h=imcrop(h1);--> διαλέγεις ένα μέρος ακόμη πιο μικρό που πιο πιθανό να είναι πιο κοντά στο ιδανικό

[MARIOS]

και αυτό που σκέφτηκα τώρα είναι το εξής....αφού επιλέξεις και καθαρίσεις μια εικόνα....που θα έχει θεωρητικά ένα αστέρι....μπορείς να το κάνεις imresize και να το μειώσεις σε έναν μικρότερο πίνακα....και έτσι θα έχεις μεγάλες κεντρικές τιμές και μικρότερες ακριανές

[mitsos_dlx]

Καλή ιδέα... Ισως θα μπορούσαμε να αυξήσουμε την αντίθεση εναλλακτικά. Θα το τσεκάρω το βράδυ.

Για να μη σου βγάζει error μπορείς να κάνεις τις πράξεις με τον H(u1,u2) ανά στοιχείο και όχι ως πίνακες. Νομίζω αυτό είναι σωστό από fourierικής πλευράς... Επίσης θέλει συζυγή στον αριθμητή. Όχι ανάστροφο.

Τελικά το πρόγραμμά μου βγάζει χειρότερα αποτελέσματα... Τσέκαρα το (x-y).^2 και το (x-xhat).^2 και (αν και τα y και xhat φαίνονται τα ίδια σαν εικόνες) το δεύτερο είναι μεγαλύτερο κατά μέσο όρο...

[MARIOS]

Βασικά η h δεν θα πρέπει να είναι τετραγωνική??

Οπότε δεν θα έχω θέμα....

άρα μέσα στην συνάρτηση την δεύτερη που βγάζω την h μετά το crop θα διαγράφω κάποιες στήλες ή γραμμές και θα την κάνω τετραγωνικό πίνακα!!!

[alexx_m87]

εγω σκεφτηκα να κλεψω λιγο και να δω ενα μικρο αστερι στην αρχικη εικονα χ και μετα να κοψω απο την y αυτη την περιοχη, αλλα μαλλον δε θα το κανω.τωρα εχω την εξης απορια: για να φτιαξω το φιλτρο θελω μετ/μο φουριε αλλα στο ματλαμπ κανεις μονο dft, αρα θα προκυψει πινακας και οχι συναρτηση των ω!! οταν κανω ifft2 το προβλημα ειναι οτι ο πινακας που προκυπτει, δηλαδη η μασκα βινερ εχει και αυτη μιγαδικουσ!!!! επισησ αμα δοκιμασω να εφαρμοσω το αντιστροφο φιλτρο απο το h που εκτιμησα, η εικονα ειναι οτι να ναι , καμια σχεση με τη χ

[MARIOS]

εγω σκεφτηκα να κλεψω λιγο και να δω ενα μικρο αστερι στην αρχικη εικονα χ και μετα να κοψω απο την y αυτη την περιοχη, αλλα μαλλον δε θα το κανω.τωρα εχω την εξης απορια: για να φτιαξω το φιλτρο θελω μετ/μο φουριε αλλα στο ματλαμπ κανεις μονο dft, αρα θα προκυψει πινακας και οχι συναρτηση των ω!! οταν κανω ifft2 το προβλημα ειναι οτι ο πινακας που προκυπτει, δηλαδη η μασκα βινερ εχει και αυτη μιγαδικουσ!!!! επισησ αμα δοκιμασω να εφαρμοσω το αντιστροφο φιλτρο απο το h που εκτιμησα, η εικονα ειναι οτι να ναι , καμια σχεση με τη χ

βασικά γιατί έχει μιγαδικούς?

 G=....
g=ifft2(G);
και δεν έχει μηδενικά....έκανα μερικά παραδείγματα και δεν το είδα αυτό....

[alexx_m87]

δεν ξερω ρε γμτ...ισως εχει προβλημα το ματλαμπ.γιατι μια μου το βγαζει ετσι, μια κανονικα, ενω ειμαι σιγοθροσ οτι σε ifft2 δεν εχει μιγαδικουσ1!!! δοκιμασε κανεισ το αντιστροφο φιλτρο? θα το κανουμε ετσι???????:
H=fft2(h)
Hinv=inv(H)
hinv=ifft2(Hinv)

εεεεεε???????

[MARIOS]

εγώ να ρωτήσω για την 2η συνάρτηση....

για την εξομάλυνση....αφού κάνεις το crop και σου αρέσει κτλ κτλ...

απλώς κάνεις
h=medfilt2(h,[m n]);
και τελειώνεις την συνάρτηση?

αλέξη....λογικά έτσι είναι.....

μετά από αυτό την περνάς πάλι την hinv στην wienfilt έτσι?
για την xinv0 τι ακριβώς ζητάει να κάνουμε γιατί έχω μπερδευτεί εκεί....

[alexx_m87]

ο ντελο στο μθμ ειπε να κανουμε εξομαλυνση την εικονα  y πριν κανουμε κροπ για τη μασκα καλυτερα.οποτε κανεισ Y=medfilt2(y) και κοβεις απο την Υ. ειπε να το δοθμε και χωρις την εξομαλυνση να δουμε πως θα ναι, αλλα καμια διαφορα δε βλεπω εγω. ρε παιδια, δεν ετυχε σε κανεναν να κανει ιφφτ2 και να του βγει πινακας με μιγαδικους? τι παιζει???επισησ το φιλτρο βινερ μπορω να πω χειροτερο και απο το y το κανει το σημα.

[Groucho]

δεν ξερω ρε γμτ...ισως εχει προβλημα το ματλαμπ.γιατι μια μου το βγαζει ετσι, μια κανονικα, ενω ειμαι σιγοθροσ οτι σε ifft2 δεν εχει μιγαδικουσ1!!! δοκιμασε κανεισ το αντιστροφο φιλτρο? θα το κανουμε ετσι???????:
H=fft2(h)
Hinv=inv(H)
hinv=ifft2(Hinv)

εεεεεε???????

ifft2(H,'symmetric') για να μην εχει μιγαδικους

[MARIOS]

Καλησπέρα

Να ρωτήσω κάτι? Στην filtest αυτό που ουσιαστικά βάζουμε είναι αυτές οι δύο εντολές έτσι?

[r p]=size(y);
y=medfilt2(y,[r p]);--> εξομάλυνση
h=imcrop(y);--> επιλογή

και από κει ότι καταλαβαίνει ο καθένας για την συνέχεια του....πχ εγώ έκανα και ένα Imresize(H,1/8) για να πάρω μια πιο "ιδανική" h.
και στην συνέχεια εάν δεν ήταν τετραγωνικός ο πίνακας h το κάνω τετραγωνικό αφαιρώντας κάποιες γραμμές ή στήλες.

Αλλά στην ουσία είναι οι τρεις πιο πάνω εντολές....