THMMY.gr

Στο demo3a εκει που λέει Ir1=deconv2naive(Io1,Id1,Id1); νομιζω οτι πρέπει το τελευταίο ορισμα να ειναι Id2. Aντιστοιχα και στα άλλα scripts

επειδη την δουλεψα την εργασια μου φαινονται ολα σωστα,
επισης Ιr1 προκυπτει απο Ιd1 ... το χει σαν test κατα πασα πιθανοτητα

και μιας και ανοιξε το θεμα :
Παραμετρική εκτίμηση του μηχανισμού παραμόρφωσης στο πεδίο του χώρου
τι διαδικασια ακουλουθησατε για να βρειτε το sigma??

γιατι αυτο μου φαινεται λιγο κουλο:

while(norm(y-yhat)>e)
    sigma=sigma+1;
    h=fspecial('gaussian',k,sigma);
    yhat=imfilter(x,h,'full','conv');
end

Για να βρεις το sigma, τρέξε την gaussian filter της προηγούμενης εργασίας για κάποιο εύρος s, κάνε συνέλιξη με την εικόνα εισόδου και μετά αφαίρεσε αυτό που σου βγάζει με την εικόνα που προκύπτει απ το fspecial. Ο πινακας διαφοράς που έχει το μικρότερο άθροισμα στοιχειων ειναι ο σωστός(δεν ξέρω αμα κανει το ιδιο το scriptaki που πόσταρες  :P )

ναι ρε αυτο κανω ,λιγο κουλο μου φαινεται ομως
γιατι δινει και 3 μοναδες! και ειναι απλο

Τι να σε πω, εμενα δεν μ ερχεται κατι αλλο και βγάζει και σωστό αποτελεσμα. Ισως τν επιασαν καλοσυνες :P

Είναι απλό γιατί είδαμε ότι σ=5  και κάθε φορά ελέγχουμε για ένα σύνολο ακεραίω σ.. αν ήταν 0.256 θα εψαχνες διαφορετικα ! Ο μεταπτυχιακος πάντως στο email που εστειλα ειπε οτι ειναι οκ ο τροπος που λετε και εκανα κι εγω ! Στο demo3c, στην δευτερη εικονα , σας βγαζει κατι ασπρομαυρο απροσδιοριστο?

ναι σαν χιονια,ετσι πρεπει να βγει γιατι ενισχυει τον προσθετικο θορυβο το αντιστροφο φιλτρο

demo3e σας βγαζει best K=20?

και επιπλεον χρειαζεται να αφαιρεσουμε μεσους ορους των εικονων για ναι ειναι
unbiased οι εκτιμησεις?

Εμένα το καλύτερο αποτέλεσμα μου το βγάζει για Κ=3.
Δεν είμαι σίγουρος τι εννοείς με το 'μεσους ορους εικονων'.

Για την μείωση του λευκού θορύβου θα χρησιμοποιήσουμε low-pass filter;

τι error βγαζεις? σου βγαζει της ταξεως 2.4551?
βασικα ακυρο εβαζα λαθος διαφορα Ir2-Id2

Εμενα μου το βγαζει για Κ=4. Την H την υπολογισατε με την estFrDom στην συχνότητα κατευθειαν ή με την estSpDom στον χώρο και μετα στην συχνότητα;

Και εγώ Κ=4 με σφάλμα 0.0025 περίπου, χρησιμοποιώντας την estFrDom.... Στο δεύτερο ερώτημα αυτά που λέει στο τέλος ότι οι εικόνες πρέπει να είναι ευθυγραμμισμένες τί εννοεί? είναι κάτι που πρέπει να το ελέγξουμε?

Στην ανακατασκευή της εικόνας χωρίς θόρυβο, βγάζετε ακριβως την αρχική εικόνα?

Που αναφέρεσαι; Υπάρχουν αρκετές περιπτώσεις που δεν έχουμε θόρυβο

Στο "Αντίστροφο φίλτρο" (deconv2naive), που εκτελεί το demo3a(), αναφέρεται.

Στο ερώτημα "Αντιστροφή φίλτρου όταν είναι γνωστή η αρχική και η παραμορφωμένη εικόνα"(όπου δεν υπάρχει θόρυβος),  η ανακατασκευασμένη εικόνα που βρισκετε, ποσο απέχει απο την αρχικη?

Ιδιες μου βγαινουν ,τι εκανες ν σου πω τι μπορει να κανεις λαθος η εναλλακτικα
pm και στειλε κωδικα

Ερωτηση και γω, με το story βγαλατε ακρη ?

Οταν χρησιμοποιώ το estFrDom κανω τον DFT των Ιο1 και Ιd1, μετα zeropad την Ιο1 και τις διαιρώ για να βρω το H. Μετα διαιρώ το Id2 με το Η για να βρω την εικόνα. Εδω σκέφτομαι μήπως πρέπει να κάνω καποια ευθυγράμιση, οπως περιγράφεται στις διαφανειες του κεφαλιου 7. Αλλα οταν κάνω το πράγμα χειροτερεύει  :P Οταν χρησιμοποιώ το estSpDom ενώ βρίσκω σωστά το gaussian φιλτρο, στη συνέχεια παιρνω το DFT του, αποκοπτω καποιες πολύ χαμηλές τιμες και το αντιστρέφω, αρα έχω το 1/Η. Στη συνέχεια, κόβω απο εκει καποιες ψηλες τιμες και το ξαναφέρνω στο χρόνο οπου το συνελίσσω με το Id2(η 1η μέθοδος που υπάρχει εδώ συγκεκριμενα http://www.owlnet.rice.edu/~elec539/Projects99/BACH/proj2/inverse.html ). Παλι δεν έχω πολύ καλα αποτελεσματα. Γραφω και τα 2 γιατι δεν έχω αποφασισει ποια απ τις 2 συναρτήσεις θα χρησιμοποιήσω. Αυτα :P

δοκιμασε αυτο
 H=estFrDom(x,y)
X=fft2(x,m,n);
Y=fft2(y);
Η=Υ./Χ;
m,n μεγεθος του y

επισης αυτο που αποκοπτεις χαμηλες και υψηλες δεν χρειαζεται
γιατι η αρχικη εικονα ειναι noiseless.Αυτο που λεει στο tutorial το κανει για
μην ενισχυσει τον θορυβο της αρχικης.Αλλα για να εισαι απολυτα σωστος τσεκαρεις για Η αν εχει μηδενικα για να μην βγαζει 1./Η inf, ιδιο και για Χ

Και μετά απλά φέρνω και την Id2 στη συχνότητα και διαιρώ Id2./H?

ναι

Οντως δουλεύει, να σαι καλα, με διευκόλυνες παρα πολυ  :) :) :)

+1 καποιος κανα hint?

Θεωρητικα είναι ιδια περιπτωση με το προηγουμενο απ αυτήν ερώτημα(gaussian φιλτρο και λευκός θόρυβος). Αρα αμα χρησιμοποιησουμε την deconv2noise θα βγει. Θεωρητικα μονο βέβαια γιατί εγώ δεν το καταφερα η αληθεια ειναι. Ισως χρειαζονται καποια απο τα φιλτρα της εργασιας 2.

Καλημέρα!
Στο  2ο ερώτημα που ζητάει οι εικόνες  x και y να είναι μεταξύ τους ευθυγραμμισμένες, τι εννοεί ?

Έχει καταφέρει κανείς να κάνει την ιστορία; Γιατί ό,τι και αν έχω προσπαθήσει βελτιώνει κάπως το αποτέλεσμα, αλλά δεν καταλήγω σε εικόνα που να διαβάζεται.

+1 αν το έχει καταλάβει κανείς ας βοηθήσει

Δοκιμασε να δουλέψεις χωρις να το λάβεις υποψιν σου και αμα δουλεύει αυτο που κάνεις,φανταζομαι ειναι οκ :P Εγω τουλαχιστον οταν χρησιμοποιησα την estSpDom για το deconv2naive2, εβγαλα καλο αποτελεσμα χωρίς να κάνω κάτι για τον ευθυγραμισμο των εικόνων

Το έχω βγάλει για την απλή περίπτωση και όντως δουλεύει. Αλλά μάλλον ζητάει και κάτι παραπάνω :P

'Εβγαλε τελικά κανείς καμιά άκρη με αυτό?

Όχι και γενικά με το imfilter('full') μας έχει πηδήξει λίγο γιατί το padding είναι μονό. Έχει καταλάβει κανείς ποιό είναι το σωστό. Ειδικά στο demo3e που ψάχνεις το MSR χρειάζεται να τα χεις σωστά.

Στο demo3d με το Wiener, παίρνετε καλά αποτελέσματα? Γιατί εμένα η ανακατασκευασμένη δεν είναι πολύ καλύτερη από τη distorted, και δεν ξέρω αν αυτό είναι το αναμενόμενο :P

Και σε μένα έτσι είναι, νομίζω είναι αναμενόμενο από τη στιγμή που δε γνωρίζεις τα χαρακτηριστικά του θορύβου και παίρνεις τον προσεγγιστικό τύπο με το Κ..

+1

Στην deconv2noise τι προεπεξεργασία κάνετε για την Id1? δοκίμασα διάφορα, αλλά χωρίς αποτέλεσμα :(

Εγω εκανα τα εξης:πήρα την h απο την EstSpDom, εκανα μια συνάρτηση παρόμοια με την bestK για να βρω το βέλτιστο K για to φιλτρο wiener και τελος απλα εφαρμόζω αυτό το φίλτρο στην Id2. Δεν βγαζει τελειο αποτελεσμα, ειναι βελτιωμενο παντως

Thnx, θα το δοκιμάσω. Δεν κάνεις κάτι στην Id1 δλδ?

Μπα οχι, μονο αυτα π έγραψα παραπανω

Ε ναι όμως δεν λέει να κάνουμε pre-processing την Ιd1?

Ναι preprocess στην Id1 κάνεις περνώντας την από το h που βρίσκεις από την estSpDom. Το φίλτρο που βρίσκεις τελικά και μετά το wiener το εφαρμόζεις στην Id2 ...

Τελικά με το story έβγαλε κανείς άκρη;