THMMY.gr

εχει ασχοληθει κανεις με τη 2η εργασια?

Λίγο...
Εκεί που πρέπει να βρούμε τη συνάρτηση μεταφοράς θα κόψουμε μάλλον με το χέρι (με το ποντίκι για την ακρίβεια... :P) ένα μέρος από την θορυβοποιημένη εικόνα που μας δίνεται, σωστά;
Αλλιώς είναι περίπλοκο να βρούμε ποιο μέρος είναι κατάλληλο προγραμματιστικά. Ενώ με το μάτι βλέπεις και παίρνεις ότι σου αρέσει πολύ εύκολα...

ναι κατι τετοιο κατεληξα και γω οτι πρεπει να γινει αλλα δεν το χω κανει να δουλεψει ακομα...θα ρθω μαθημα αυριο να ρωτησω...

ειναι κανεις που την ετρεξε και του εμφανισε αντι για 5 εικονες 4  ? χωρις την xinv ..?

Ναι.
Υπάρχει ένα πρόβλημα σε μια εντολή subplot αν θυμάμαι καλά και πέφτει το ένα διάγραμμα πάνω στο άλλο στο grid που δημιουργούν.

Τελικά πως βρίσκουμε την h στην πρώτη συνάρτηση που θα κάνουμε???

μεσω της imcrop θα πας πανω σε ενα αστερι της θολωμενης εικονας .. μικρη διασταση δλδ του στυλ 5χ5 

Η θολωμένη εικόνα είναι εκείνη που μας δίνει στο .tif αρχείο έτσι?

Η επιλογή μέσω της imcrop πως θα γίνει? Δηλαδή πως θα βρω αστέρι?

Θα ψάξω την εικόνα μέσω των τιμών του grayscale?

αν γραψεις imcrop (y) θα σου ανοιξει την εικονα y τη θολωμενη...την y την διαβαζει η συναρτηση σαν ορισμα απο την demo2

δηλαδή εγώ θα την βάλω την εντολή στην demo και κάθε φορά θα επιλέγουμε εμείς την h όταν θα τρέχει η demo

Βασικά οκ για την h.... θα την βρώ την άκρη... ;)

για τον μετασχηματισμό Fourier

θα τον κάνουμε εμείς ή υπάρχει κάποια συνάρτηση?

εγώ βρήκα την fft για fft ( :P )

παίζει να έχει και fft2 για 2Δ σήματα

Α ναι όντως.....

αλλά μετά για τον αντίστροφο μετασχηματισμό βρήκα στο net να το κάνει έτσι --> numeric::invfft(x)


έχετε εσείς κάτι στο μυαλό?

ifft2 νομίζω...

το κάνω...το δέχεται σαν εντολή αλλά εμφανίζει μια άσπρη εικόνα..... τόσο μεγάλο σφάλμα.... :/

Ρε παιδιά....ενώ κάνω το

F=fft2(y);
y1=ifft2(F);
imshow(F)

Μου εμφανίζει μια άσπρη, κάτασπρη εικόνα...και δεν θα έπρεπε να το κάνει
τι παίζει?

Θέλεις να προβάλεις τον ΜΣ Fourier ή την εικόνα y1 και σ'το βγάζει αυτό;


Μια ερώτηση από μένα: το K πώς θα το βρούμε; Λέει με διαδοχικές δοκιμές, ενώ μετά λέει να γράψουμε τη μεθοδολογία υπολογισμού του K. Εννοεί να πούμε ποιες δοκιμές κάναμε (για ποιες τιμές του K);

Βασικά πρέπει να βρεις το Κ που ελαχιστοποιεί το μέσο τετραγωνικό σφάλμα.

Αυτή την απορία την έχω κι εγώ. Καλώ την imcrop και κόβω ένα αστέρι. Από εκεί και πέρα πως παίρνω την h; Απλά με εξομάλυνση της μικρής εικονίτσας που έκοψα;

Ναι, έτσι γίνεται.

Η προθεσμία είναι μέχρι τις 26/05? Γιατί κάτι είχε πει για παράταση...

ναι μέχρι τις 26 είναι

Για την h,  όταν μιλάμε για εξομάλυνση τι ακριβώς εννοεί ο ποιητής ( γιατί σήμερα είμαι και λίγο λιώμα οπότε ούτε το όνομά μου μπορώ να καταλάβω... :P )

και για την συνέλιξη για να βρεις το x στην πρώτη συνάρτηση.... έχουμε  y(m,n) και h(k,l) έτσι? Βρήκα στο νετ το εξής
x(m,n)=ΣkΣl[h(k,l)*y(m-k,n-l)]

και στο matlab το έκανα ως εξής:

sum=0;
for (i=1:1:k)
    for (j=1:1:l)
        sum=sum+h(i,j)*y(m-i,n-j);
    end
end

είναι κομπλέ? Πως το βλέπετε?

Βασικά τώρα που το σκέφτομαι είναι μλκ.... έχει κανείς καμιά ιδέα???  :???:

βασικά η y=conv2(A,B) ( τώρα την βρήκα :P ) δουλεύει σωστά???

Γενικά δε θέλει να κάνεις κάτι με συνέλιξη. Τα περνάς στο χώρο συχνοτήτων, τα πολ/ζεις και μετά τα φέρνεις πάλι στον ευθύ χώρο.

Το y = h*x συγκεκριμένα που λες γίνεται με συνέλιξη αλλά σου το δίνει ο ίδιος ο καθηγητής αν δεις πιο κάτω σε έναν κώδικα που έχει γράψει. Και τη συνέλιξη (κυκλική εδώ) την υλοποιεί με τη συνάρτηση imfilter().

Κάτσε περίμενε λίγο....πολλαπλασιάζεις G με Y αλλά δεν ταιριάζουν στις διαστάσεις οι δύο πίνακες....ο ένας έχει τις διαστάσεις της Η και ο άλλος του Υ...

όσο για αυτό που έγραψα κάτω, μπερδεύτικα δεν εννοούσα το h  αλλά το g για να βγάλεις την έξοδο της πρώτης συνάρτησης.

Κάνε zero padding στο G ή στο H. Στο MATLAB μπορείς να κάνεις εύκολα zero padding την στιγμή που κάνεις τον FFT, μπορείς να γράψεις H = fft2(h, m, n) όπου m και n είναι οι διαστάσεις του y (ή και του Y).

Τουλάχιστον εγώ έτσι το έκανα και δουλεύει, αν και δεν βγαίνει και πολύ καλή η εικόνα. Μπορεί να υπάρχει και άλλος τρόπος.

Η ανακατασκευή της εικόνας που παίρνετε είναι καλή? Η δικιά μου όχι, και πιστεύω φταίει το ότι έχει πολλές λεπτομέριες στην αρχική εικόνα, οι οποίες χάνονται.

Επίσης για την xinv και την xinv0 χρησιμοποιούμε το hest?

Βασικά πήγα να δοκιμάσω αυτό που λες, αλλά μου έβγαζε στο τέλος ( αφού το γύριζα από την συχνότητα όλες σχεδόν (αν όχι όλες )) τι τιμές μηδενικές. Όταν πήγα να το κάνω με την conv2() που είπα την έβγαλε κομπλέ ( όπως και την αρχική μάλλον ).

Αυτό το δοκίμασα όχι στην συνάρτηση μέσα αλλά στο command window του MATLAB και χωρίς να έχει θόρυβο η εικόνα... οπότε και δεν ξέρω αν θα λειτουργεί καλά και το φίλτρο....που φαντάζομαι θα λειτουργεί... :/

Ρε παιδια, γιατι λετε οτι δε θελει να κανεισ συνελιξη? απλα αμα παρεισ τη μασκα και το σημα κανεισ ενα conv2 και κανει τη συνεληξη. Σωστα η κανω λαθος? Μια τεραστια απορια που εχω: στην εκτιμηση της παραμορφωσης παιρνω απο ενα μικρο αστερι ,οκ .μετα ομως αμα παρω την εικονα χ και την περασω απο το φιλτρο της παραμορφωσης (με συνεληξη δλδ) μου βγαζει οτι να ναι, καμια σχεση με το y, που υποτιθεται θα πρεπε να μοιαζει!!!

εμένα μου έβγαζε μια άσπρη εικόνα....στον Ντελό που του έστειλα ένα μαιλ, μου είπε ότι μάλλον θα χρειάζεται κανονικοποίηση....δηλαδή διαίρεση του πίνακα που έχεις του τελικού με το max(max(matrix))<-- το μέγιστο στοιχείο του πίνακα

Για να καθαρίσουμε και να δώσουμε σε έξοδο την h θα χρησιμοποιήσουμε την πρώτη συνάρτηση που κάναμε?

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

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

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

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


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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Καλησπέρα

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

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

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

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

εγω εκανα σκετο Y=medfilt2(y), δεν ειναι σωστο? δε βλεπω καμια διαφορα στην εικονα

εννοείς για το κομμάτι που βάζω και το [r p] έτσι το είδα στο documantation στην Matlab.....


αλλά πχ το imcrop το ζητάει και μέσα στην εκφώνηση

παντως το αποτελεσμα με το βινερ μου βγαινει χειροτερο απο το y...

αυτό θα πει βελτίωση ε???
:P

Αν θες βάλε τον κωδικό σου εδώ ή στείλε τον σε πμ μήπως δω τίποτα...

λοιπόν εγώ θέλω να κάνω τετραγωνικό τον h σε περίπτωση που δεν είναι...και χρησιμοποιώ αυτόν τον αλγόριθμο.....αλλά δεν μου βγαίνουν σωστά.... :/

if m~=n
    if m>n
        for i=1:1:n
            z(:,i)=h(:,i);
        end
    end
    if m<n
        for i=1:1:m
            z(i,:)=h(i,:);
        end
    end
    h=z;
end

Ποιον κωδικο εννοεις?

ωπ λάθος..... κώδικα ήθελα να πω.... :P

εσενα σου βγαινει οκ?

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

Να σας πω, έχω κάνει μια συνάρτηση που υπολογίζει το J και το Κ για Jmin
 
function [K min]=minJ(x,y,h)
for K=1:1:10
    J(:,:,K)=(x-wienfilt2(y,h,K))^2;
end
min=MIN(J);
s=0;
i=1;
while s==0
    if min==J(i);
        K=i;
        s=1;
        i=i+1;
    end
end
end

Το ερώτημά μου είναι μάλλον απλό....πως κάνω τώρα πλοτ?

εεε να σας πω....μπορεί κάποιος να μου πει κάποιον τρόπο για να εμφανίσω το J? Ας μην είναι πάνω στην συνάρτηση που έγραψα...αν δεν έχει πρόβλημα ας μου πει κάτι δικό του...απλώς να τελειώνω... please

Εγώ δεν έκανα ξεχωριστή συνάρτηση απλά ένα βρόχο για κάποια k. 
To J είναι διάστασης kx1 και έχει τα σφάλματα για κάθε k. Το διάγραμμα το εμφάνισα απλά με την εντολή:
plot([1:length(J)], J);

Τον κώδικά σου δεν το είδα και πολύ αλλά δε νομίζω να είναι δύσκολο να κάνεις κάτι αντίστοιχο.

Το J=(x-xhat)^2

όπου x, xhat είναι πίνακες....

γιατί είναι kx1?

Βασικά είναι το μέσο τετραγωνικό σφάλμα. Δηλαδή έχεις πχ:
e=(x-xhat).^2
J(k)=mean(e( : ));

Παιδια, στο βινερ, στον παρονομαστη θελει μετρο του πινακα στο τετραγωνο η νορμα? γιατι ειναι τελειως διαφορετικα πραγματα!

Πραγματικά......έχεις έναν καφέ κερασμένο από μένα!

 ^banghead^ ^banghead^ ^banghead^ ^banghead^ ^banghead^

(abs(H)^2   είναι η απόλυτη τιμή των στοιχείων...

πως βρηκατε xinv0 αφού η Η0 έχει αρκετές μηδενικές τιμές άρα δεν μπορώ να πάρω αντιστροφο?

πρόσθεσε eps, ( δηλ γράψε +eps), και έτσι θα προσθέσεις μια πολύ μικρή τιμή οπότε υπάρχει το αντίστροφο.

Θέλει και μπορεί κάποιος να του δώσω την demo2 που έκανα να μου πει αν τρέχει και δεν κρασάρει και σε εκείνον?

Πάντως και εμένα τα αποτελέσματα είναι για τον πούτσο

 το J(k) πως σας ωγαινει, φθινουσα η αυξουσα? εμενα δεν ξερω τι παιζει, αλλαζει...

η νορμα2 ειναι ουσιαστικα το μετρο..