THMMY.gr

στην 3η εργασια ποιο φιλτρο θα παρουμε για h

To φίλτρο Haar.

Όχι το haar ... Αφού δίνει αρχείο με πίνακα. Αυτό πρέπει να χρησιμοποιήσουμε.
Τώρα αν αυτό που δίνει είναι το h0 η το g0 δεν είναι και πολύ σαφές, αλλά έχω μια υποψία ότι το ίδιο μας κάνει αν υπολογίσεις τα υπόλοιπα με τις σχέσεις που δίνει.

Το λέει ξεκάθαρα ότι είναι το h0.
Εκεί που κολλάω εγώ είναι το γεγονός ότι στο τέλος η εικόνα θα έχει και αρνητικές τιμές. Τι τις κάνουμε αυτές?

Το L που το χρησιμοποιούμε στο τέλος στην εκτύπωση των αποτελεσμάτων λέει ότι είναι η μέγιστη τιμή φωτεινότητας.
Δηλαδή είναι το μέγιστο του x ή παίρνουμε το μέγιστο κάθε φορά από το y00, y01, y10, y11;

Επίσης η εικόνα στο τέλος (το xhat) μου βγαίνει πολύ σκοτεινή. Πρέπει να την κανονικοποιήσω για να γίνει κανονικά φωτεινή ή κάπου έχω κάνει λάθος;

Η εξίσωση 14 δεν ισχύει αν ξεκινήσουμε με τους συντελεστές που δίνονται στο πίνακα db10. Για αυτό η εικόνα βγαίνει σκοτεινή.
Αν τους προσαρμόσετε ώστε να ισχύει η εξίσωση, η εικόνα θα βγαίνει με ίδια φωτεινότητα. Δείτε πόσο λάθος βγαίνει η 14 και θα καταλάβετε τι πρέπει να γίνει... ;)

 :???:
Κανονικοποίηση θα κάνουμε; Δεν μπορώ να καταλάβω τι ακριβώς πρέπει να αλλάξω.

Ρε παιδιά. Κατάλαβε κανείς πως βγαίνουν τα h0, h1, g0, g1 από τους τύπους? Έκανα την μετατροπή στον χρόνο αλλά οι δ[] που προέκυψαν είναι πολύ μακρυά από το 1~20 (από τις διαστάσεις του db10). Και τα μηδενικά που θα μπουν?

Thanks.

for daubechies filter

http://www.mathworks.com/help/wavelet/ref/dbaux.html (http://www.mathworks.com/help/wavelet/ref/dbaux.html)

Εγώ για τον υπολογισμό χρησιμοποίησα από το βιβλίο του Gonzalez τους τύπους 7.1-23 σελ 358. Δεν προκύπτουν δ[].

Το πρόβλημα για μένα είναι αυτό που λέει ο mitsos_dlx πιο πάνω. Τι αλλαγές πρέπει να κάνω στο db10 πριν το χρησιμοποιήσω;

Να σας πω.... πρώτον πως ανοίγουμε το αρχείο mat μέσα στο πρόγραμμα.... δηλαδή δεν μπορούμε να κάνουμε χειροκίνητα open....
δεύτερον.... οκ το ανοίγουμε και αποθηκεύουμε τους συντελεστές σε έναν πίνακα....από κει και έπειτα πως προχωράμε?

δηλαδή
open-->db10
y=db10;

μετά τι?

Αφού τους ανοίξεις μετά κάνεις
[h0,h1,g0,g1] = orthonormalFitlerBanks(y);
και αφού έχεις τα φίλτρα κάνεις ανάλυση και σύνθεση...

Όσο για το αρχείο db10 εμένα μια χαρά μου βγήκε. Στο τέλος πρέπει να κάνετε:   imshow(xhat/max(xhat( : )) ) (το οποίο είναι λογική κίνηση)
Αν δεν διορθώνεται με αυτό, τότε έχετε λάθος στον κώδικα...

ναι όμως για να ανοίξουν δεν πρέπει να έχεις κάποια εντολή, fopen(db10) πχ

load(db10). Ελεος πια

παιδια, η σχεση 15 στο χρονο πως ειναι??

άμα μου βρεις την 16....

εγώ θα σου πω τις άλλες δύο...;)

παλικάρι......calm down

\m/

και να σε καλά......εγώ θεωρούσα ότι κάτι σε fopen ή open θα είναι.....οπότε πραγματικά χαλάρωσε.....

ελα ρε, απο τις ιδιοτητεσ του μετ/μου Ζ, βγαινει οτι h(n)=g(-n)

μα εγώ δεν έχω διαβάσει μετασχηματισμό z.... :P

g1(n) = ((-l)^(n))*g0(2K - 1 - n)
hi(n) = gi(2K - 1 - n), i = {0,1}

Ψήνετε να ζητήσουμε καμιά παράταση;

Ναι γιατί όχι....


Λοιπόν από τις ιδιότητες από ότι είδα είναι
x(-n)=X(z^(-1))
και επειδή είναι αιτιατό το σήμα είναι
x(-n)=x(n)?
δηλαδή στην προκειμένη θα έχω
h0(n)=g0(n)?


και ισχύει: h0(n)=db10;

Το αιτιατό σήμα πού χρησιμεύει;
Η ιδιότητα είναι όντως x(-n) <-> X(z^-1)  (όχι =)
Άρα αν έχεις ένα πίνακα [1 2 3 4] για g(n), η h(n) που θα προκύψει από την Η(z)=G(z^-1) θα είναι η h(n)=g(-n).
Δηλαδή αν το πρώτο στοιχείο της h(n) ήταν στη θέση n=0 τότε (θεωρητικά) θα πρέπει να κάνεις ένα πίνακα με θέσεις -3 έως 0 και να βάλεις [4 3 2 1].
Βέβαια επειδή δεν μπορείς να πεις στο matlab h(-3)=g(3), ... , h(0)=g(0) γιατί θα σου πει ότι οι πίνακες του ξεκινάνε από 1 και αυξάνουν θα πρέπει να κάνεις κάτι άλλο. Κάπου εδώ κολλάνε τα μηδενικά αλλά και να μη τα βάλεις πρακτικά δεν θα υπάρξει πρόβλημα. Θεωρητικά μπορεί να είναι λάθος να πεις h=[4 3 2 1] αλλά θα δώσει ίδιο αποτέλεσμα σε συνελίξεις με το αν της έβαζες μηδενικά γύρω-γύρω (όχι ίδιο βασικά αλλά για την δουλειά που κάνουμε θα είναι ok).

ναι μεχρι την κυριακη αν ειναι !!!

Του είπα σήμερα να ανοίξει τις ομάδες και έβαλε μόνος του ημερομηνία λήξης 17/6  :) άρα ως τη Δευτέρα

Μπράβο, Ντελόπουλε!

Να σας πω, πως θα δείξουμε ότι ικανοποιείτε η σχέση 14?

θα πει καποιος ρε παιδια ποσο βγαινει η σχεση 15 στο χρονο?

g1(n) = ((-l)^(n))*g0(2K - 1 - n)

Από τις εξισώσεις 7.23 που έχει στο βιβλίο των Gonzalez-Wood

Μετατρέποντας τη στον χρόνο. Θα βρεις με τις ιδιότητες σε τι αντιστοιχούν στο χρόνο τα G(z^-1), G(-z) κλπ. Όπως βρήκες και τα αλλά φίλτρα δηλαδή από τις υπόλοιπες εξισώσεις. Μετά οι πολλαπλασιασμοί γίνονται συνέλιξη, μια πρόσθεση και έτοιμο. Το 2 στο δεξί μέλος αντιστοιχεί σε 2δ(n)

Μην το ψάχνετε τόσο απλά είναι τα πράγματα. Εγώ δεν είδα καν τις σχέσεις για να βγάλω τα φίλτρα. Δέστε εδ (http://www.mathworks.com/help/wavelet/ref/orthfilt.html) εκεί που έχει τις εικόνες των φίλτρων μπορείτε εύκολα να καταλάβετε πως βγαίνει το καθένα. Επίσης αν θέλετε να είστε σίγουροι ότι δουλεύει ο κώδικάς σας τρέξτε τον με h = [1/sqrt(2) 1/sqrt(2)] πρέπει να βγάλει το ίδιο αποτέλεσμα και εκεί.

Δηλαδή μπορούμε στην ψύχρα να χρησιμοποιήσουμε την συνάρτηση αυτή...???

[Lo_D,Hi_D,Lo_R,Hi_R] = orthfilt(w);

 :P

΄Πάντως επειδή έχω μπερδευτεί ( Βλέπε τα χτεσίνα ξύδια.... ) στην 2Δ analysis2d(ho,h1,x)

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

Λοιπόν και η απορία μου είναι η εξής....
Έχω τον αρχικό MxN πίνακα και την [y0 y1]=analysis(h0,h1,x)

Στην θέση του x θα ξεχωρίζω κάθε φορά μια στήλη.... τα h0,h1 κομπλέ...θα το γράψω έτσι

for i=1:1:dis(2)  //dis(2)-->αριθμός στηλών
     [y0(i) y1(i)]=analysis(h0,h1,x(:,i));
end

και μετά?

Πως ξανασχηματίζω τους πίνακες???

*το [y0(i) y1(i)] φαντάζομαι είναι σωστό σαν έκφραση....
* μήπως κάτι τέτοιο θα έλυνε το πρόβλημά μου?

f0=x;
f1=x;
for i=1:1:dis(2)  //dis(2)-->αριθμός στηλών
     [f0(:,i) f1(:,i)]=analysis(h0,h1,x(:,i));
end

Real, i need some help here.......... please....