THMMY.gr

Μαθήματα Κύκλου Ηλεκτρονικής & Υπολογιστών => Ψηφιακή Επεξεργασία Εικόνας => Topic started by: ΚΗΜΜΥ on March 26, 2012, 13:27:52 pm



Title: [ΨΕΕ] Εργασιες 2012
Post by: ΚΗΜΜΥ on March 26, 2012, 13:27:52 pm
Ανεβηκε η πρωτη εργασια στο ετημμυ.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: il capitano on March 26, 2012, 17:58:26 pm
Το πορτοκάλι μ' άρεσε πολύ! :P


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: ΚΗΜΜΥ on March 26, 2012, 22:35:03 pm
Ρε καμμενε...


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on April 04, 2012, 16:24:59 pm
κανει κανεις την εργασια με την interp2? Γιατι δεν πολυβολευει αυτη η συναρτηση... :-\


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: BabylonX on April 05, 2012, 00:18:30 am
Εγώ την έκανα χωρίς την interp2.

Βασικά ας δώσει κάποιος κανένα hint για τη demo... Έχω κάνει τη συνάρτηση bayer2rgb και δεν έχω ιδέα πως να κάνω τη demo. Δηλαδή τι πρέπει να κάνει αυτή η demo; Να δημιουργεί όλα τα δεδομένα και μετά να τρέχει τη συνάρτηση; Γίνεται αυτό; Και πως θα ξεχωρίζουν οι χρήσεις για type=1 και type=2;

Ο Ντελόπουλος είπε "Θέλω να κάνω διπλό κλικ το demo και να τρέχει"....

Dafuq?  :o


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Salvation on April 05, 2012, 00:23:00 am
Απλά θα κάνεις ένα script που θα το ονομασεις demo1, μέσα στο οποίο θα διαβαζεις την εικονα με το πορτοκαλι, imread, και μετα θα καλεις την συναρτηση bayer2rgb 2 φορές, μία με όρισμα με type=1 και μια με type 2 και θα σου δειχνει 2 figures με τα αποτελεμσατα για nearest neighbor και bilinear. (Αυτα με την προυποθεση ότι μέσα στη bayer2rgb έχει βαλει imshow).



Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Vicariously,I on April 05, 2012, 01:10:59 am
πες'τε κάτι για τη μεθοδο της παρεμβολης.. ειπε κατι στο μαθημα γιατι δεν εχω ιδεα τι να κανω.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on April 05, 2012, 01:46:29 am
πες'τε κάτι για τη μεθοδο της παρεμβολης.. ειπε κατι στο μαθημα γιατι δεν εχω ιδεα τι να κανω.

Αν googlάρεις bayer bilinear interpolation θα βρεις αρκετά σχετικά με το θέμα. Στο μάθημα δεν είπε κάτι το ιδιαίτερο πέραν της εκφώνησης.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: sok on April 05, 2012, 13:58:07 pm
Από που διαβάζεται για την εργασία ρε παιδιά;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Vicariously,I on April 06, 2012, 01:28:59 am
εχει κανεις πληροφοριες για το πως γινεται η nearest neighbour interpolation γιατι ετσι οπως το κανω μου βγαινει μπλε! το πορτοκαλι :P


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on April 06, 2012, 01:43:50 am
εχει κανεις πληροφοριες για το πως γινεται η nearest neighbour interpolation γιατι ετσι οπως το κανω μου βγαινει μπλε! το πορτοκαλι :P

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


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Vicariously,I on April 06, 2012, 02:19:40 am
ισχυει απλα στο bilinear τα παρινω κανονικα και βγαινει....ο,τι ναναι


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on April 06, 2012, 10:02:42 am
Οχι δεν ισχυει, αναποδα πρεπει να τα παρεις και στις δυο περιπτωσεις. Τα μπλε ειναι τα κοκκινα και τα κοκκινα ειναι τα μπλε γενικα.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Tius on April 07, 2012, 04:05:55 am
Diladi pos einai to sxima?   


G R G R
B G B G
G R G R

? kai pos to katalabainei kaneis afto?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Neal on April 07, 2012, 12:45:56 pm
Ναι, έτσι είναι το σχήμα. Δε νομίζω να υπάρχει τρόπος να καταλάβει κανείς το bayer pattern απ'τον RAW πίνακα.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on April 07, 2012, 13:37:38 pm
Αυτο το καταλαβαινεις αν τυπωσεις το πορτοκαλι οπως ειναι το πλεγμα στο σχημα τοτε βγαινει μπλε αντι για πορτοκαλι. Αρα το μπλε ειναι στη θεση του κοκκινου. Δεν ειναι τιποτα, απλα οπου κοκκινο βαλε στις πραξεις το μπλε και το αντιστροφο και καθαρισες! ;)


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: atalanta on April 09, 2012, 12:32:53 pm
γιατί στην bilinear interpolation μου βγαίνει έτσι...που κάνω λάθος γμτ?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on April 09, 2012, 12:59:17 pm
γιατί στην bilinear interpolation μου βγαίνει έτσι...που κάνω λάθος γμτ?

ُΕπειδή κάθε pixel έχει μόνο ένα χρώμα. Πρέπει για κάθε pixel να βρεις και τις 3 συνιστώσες rgb

Υ.Γ Βασικά η εικονα φαίνεται από κάτω, μήπως τυπώνεις κάτι επιπλέον από πάνω που δε χρειάζεται?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: atalanta on April 09, 2012, 13:31:55 pm
γιατί στην bilinear interpolation μου βγαίνει έτσι...που κάνω λάθος γμτ?

ُΕπειδή κάθε pixel έχει μόνο ένα χρώμα. Πρέπει για κάθε pixel να βρεις και τις 3 συνιστώσες rgb

Υ.Γ Βασικά η εικονα φαίνεται από κάτω, μήπως τυπώνεις κάτι επιπλέον από πάνω που δε χρειάζεται?

τις εχω βρει και τις 3 συνιστωσες rgb νομιζω..στην μεθοδο του κοντινοτερου γειτονα μου βγαινει κανονικη η εικονα..


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: BabylonX on April 09, 2012, 13:33:09 pm
Την είχα πατήσει και εγώ όπως η atalanta. Δε πρέπει να στείλεις την πληροφορία του χρώματος μόνο στα γύρω γύρω pixel. Πρέπει να πάει και στο pixel που είσαι. Μάλλον δηλαδή ξεχνάς τις μισές αντιγραφές.Πχ το κόκκινο από το [2,1] πρέπει να πάει στα [1,1], [1,2], [2,2] αλλά και στο [2,1]


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: ΚΗΜΜΥ on April 09, 2012, 20:20:31 pm
Δευτερη εργασια στο ετημμυ. Με γεια μας


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on April 11, 2012, 19:56:40 pm
Τι μαλακια ειναι αυτο για εργασια... >:(


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: ΚΗΜΜΥ on April 11, 2012, 23:10:02 pm
Τι μαλακια ειναι αυτο για εργασια... >:(

++


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: SolidSNK on April 12, 2012, 06:42:40 am
Εμένα με χαλάει που ο Ντελό δε μας έδωσε κάτι προσθετικό  :-\


Τον αγαπάμε ακόμα όμως  :P


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: averelgr on April 16, 2012, 22:22:44 pm
Παιδιά έχετε καταλάβει τι παίζει με την πολυπλοκότητα? Κάποια βοήθεια κανείς?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Faros on April 17, 2012, 12:20:43 pm
Μπορει κανεις να ανεβάσει την εκφώνηση της 2ης εργασίας γιατι δε μπαίνει στο e-thmmy?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: averelgr on April 17, 2012, 12:26:47 pm
Ορίστε!


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: phen3x on April 21, 2012, 16:08:14 pm
Την πολυπλοκότητα θα την βγάλεις ανάλογα με τον κώδικα που έχεις γράψει ρε συ.. Λογικά θα είναι κοντά στο Ν^2logN. Με την απόδειξη βγάλατε καμία άκρη; Εγώ σκέφτομαι μόνο το προφανές αλλά δεν γράφω κάτι αυστηρά μαθηματικό και πιστεύω δεν είναι αρκετό...


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: gerasimos on April 22, 2012, 16:13:34 pm
Κανείς καποια βοήθεια για την απόδειξη??


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Εσύ on April 22, 2012, 21:14:09 pm
utr.
Επίσης,
αν πολαπλασιάσουμε σε μια γενική περίπτωση τους πίνακες   Χ και Υ στο πεδίο της διακριτής συχνότητας κατά τον τρόπο που προτείνεται στην εκφώνηση (με N1+M1-1) έχουμε ασυμφωνία διαστάσεων. Εχει κάνει κανείς κάτι πάνω σ αυτό;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Salvation on April 22, 2012, 21:48:14 pm
Δεν ειναι ο κλασσικος πολλαπλασιασμός πινάκων. Είναι πολλαπλασιασμός στοιχείο προς στοιχείο.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Panos_Kanel on April 23, 2012, 10:38:36 am
utr.
Επίσης,
αν πολαπλασιάσουμε σε μια γενική περίπτωση τους πίνακες   Χ και Υ στο πεδίο της διακριτής συχνότητας κατά τον τρόπο που προτείνεται στην εκφώνηση (με N1+M1-1) έχουμε ασυμφωνία διαστάσεων. Εχει κάνει κανείς κάτι πάνω σ αυτό;

Δοκίμασε να κάνεις debug και πήγαινε πάνω στις μεταβλητές να δεις τι διαστάσεις πολλαπλασιάζεις.Εγώ έτσι το σωσα


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Εσύ on April 23, 2012, 11:50:43 am
Δεν ειναι ο κλασσικος πολλαπλασιασμός πινάκων. Είναι πολλαπλασιασμός στοιχείο προς στοιχείο.
Να σαι καλά!

Από που το βρήκες;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: tomshare3 on April 23, 2012, 20:45:02 pm
στην απόδειξη, για το z(n1,n2) θα ισχύει η εξίσωση (1) ή η (5)?? έχει προχωρήσει κανείς??


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: ΚΗΜΜΥ on April 24, 2012, 12:51:00 pm
Εργασία #2: παράταση προθεσμίας παράδοσης

24 Απρ 2012 11:48 πμ
Ντελόπουλος

Θα έχετε τη δυνατότητα να υποβάλετε την δεύτερη εργασία μέχρι το βράδυ της Πέμπτης 26.4.2012



Δεχομαστε συγχαρητηρια. Επισης καλη επιτυχια για την προοδο της γραφικης :D


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: SolidSNK on April 24, 2012, 14:22:47 pm
Πότε είναι η πρόοδος μάθαμε?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: ΚΗΜΜΥ on April 24, 2012, 14:25:11 pm
Θα πουμε σε λιγο στο μαθημα ;)


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: johnnysp on April 24, 2012, 19:23:30 pm
Στο κομμάτι της επαλήθευσης των τυχαίων εικόνων, στην ουσία δημιουργούμε 2 τυχαίες εικόνες, τις δίνουμε ως όρισμα στη συνάρτηση που φτιάξαμε παίρνοντας το αποτέλεσμα και στην συνέχεια υπολογίζουμε τον χρόνο αυτού του αποτελέσματος;;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on April 24, 2012, 19:25:38 pm
Στο κομμάτι της επαλήθευσης των τυχαίων εικόνων, στην ουσία δημιουργούμε 2 τυχαίες εικόνες, τις δίνουμε ως όρισμα στη συνάρτηση που φτιάξαμε παίρνοντας το αποτέλεσμα και στην συνέχεια υπολογίζουμε τον χρόνο αυτού του αποτελέσματος;;

Ναι

Πότε είναι η πρόοδος μάθαμε?

Τελικα δεν θα γινει.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: johnnysp on April 25, 2012, 19:21:11 pm
Σχετικά με τον χρόνο εκτέλεσης...υπολογίζουμε την πολυπλοκότητα της συνάρτησης μας και στην συνέχεια για να βρούμε τον χρόνο πρέπει θεωρητικά να διαιρέσουμε με την συχνότητα του ρολογιού. Έτσι δεν είναι;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on April 25, 2012, 20:32:01 pm
Σχετικά με τον χρόνο εκτέλεσης...υπολογίζουμε την πολυπλοκότητα της συνάρτησης μας και στην συνέχεια για να βρούμε τον χρόνο πρέπει θεωρητικά να διαιρέσουμε με την συχνότητα του ρολογιού. Έτσι δεν είναι;

Σου λέει να χρησιμοποιήσεις έτοιμη συνάρτηση (την etime)


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: johnnysp on April 25, 2012, 20:37:29 pm
Ναι οκ! Απλά μετά σου λέει να επαληθεύσεις με την etime την πολυπλοκότητα που υπολογίστηκε παραπάνω. Αυτή η πολυπλοκότητα δεν πρέπει κάπως να μετατραπεί σε χρόνο;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on April 25, 2012, 20:45:17 pm
Ναι οκ! Απλά μετά σου λέει να επαληθεύσεις με την etime την πολυπλοκότητα που υπολογίστηκε παραπάνω. Αυτή η πολυπλοκότητα δεν πρέπει κάπως να μετατραπεί σε χρόνο;

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

εδιτ : Τα παραπάνω με μια επιφύλαξη πάντα.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: miaou on May 19, 2012, 17:37:54 pm


3η Εργασία is up!
"Πιθανότατα στο διάστημα αυτό θα ανατεθεί και η τέταρτη εργασία με προθεσμία λίγο μετά την παραπάνω ημερομηνία."


hell yeah.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: SolidSNK on May 28, 2012, 17:13:49 pm
Για την 3η εργασία, ο αστερίσκος δίπλα από το H αναπαριστά το συζυγές της ή κάτι άλλο που έχασα?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: tomshare3 on May 28, 2012, 17:26:48 pm
ναι ο συζυγης ειναι!

εβγαλε κανεις αποτελέσματα?? ακομα και για βελτιστο Κ, η xhat μου βγαίνει χειροτερη απο την y..  :-\
Με fft2 και ifft2 παμε στις συχνοτητες και πισω? ή θελει κατι αλλο??


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: SolidSNK on May 28, 2012, 17:31:41 pm
Ευχαριστώ :)


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on May 28, 2012, 19:08:35 pm
ναι ο συζυγης ειναι!

εβγαλε κανεις αποτελέσματα?? ακομα και για βελτιστο Κ, η xhat μου βγαίνει χειροτερη απο την y..  :-\
Με fft2 και ifft2 παμε στις συχνοτητες και πισω? ή θελει κατι αλλο??

Ναι, αλλά θέλει να κάνεις zero padding στο h για να έχει ίδιο μέγεθος με το y.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: neon on May 29, 2012, 14:15:35 pm
1) "xinv ως το αποτέλεσμα του φιλτραρίσματος της y με το αντίστροφο φίλτρο H(ω1,ω2)^−1 "

πώς υπολογίζεται το Η^-1 ?

2) "xinv0 ως το αποτέλεσμα του φιλτραρίσματος της καθαρής από θόρυβο y0"

Σίγουρα ζητάει το y0 εδώ; Μήπως θέλει το θόρυβο v αντί του y0 ?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on May 29, 2012, 14:33:11 pm
1) "xinv ως το αποτέλεσμα του φιλτραρίσματος της y με το αντίστροφο φίλτρο H(ω1,ω2)^−1 "

πώς υπολογίζεται το Η^-1 ?

2) "xinv0 ως το αποτέλεσμα του φιλτραρίσματος της καθαρής από θόρυβο y0"

Σίγουρα ζητάει το y0 εδώ; Μήπως θέλει το θόρυβο v αντί του y0 ?

1)  Θα βρεις τα H ,Y με μετασχηματισμό fourier, μετά λες XINV = Y./H  και μετά θα βρεις το xinv με αντίστροφο fourier

2) To y0 θέλει. Με y που ζητάει παραπάνω θα βγάλει χάλια αποτελέσματα. Με y0 θα βγάλει καλά.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: neon on May 29, 2012, 14:44:00 pm

1)  Θα βρεις τα H ,Y με μετασχηματισμό fourier, μετά λες XINV = Y./H  και μετά θα βρεις το xinv με αντίστροφο fourier


what about imfilter(y,x/y,..) ?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on May 29, 2012, 14:47:46 pm

1)  Θα βρεις τα H ,Y με μετασχηματισμό fourier, μετά λες XINV = Y./H  και μετά θα βρεις το xinv με αντίστροφο fourier


what about imfilter(y,x/y,..) ?

Τι είναι το x/y?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: neon on May 29, 2012, 14:51:58 pm

1)  Θα βρεις τα H ,Y με μετασχηματισμό fourier, μετά λες XINV = Y./H  και μετά θα βρεις το xinv με αντίστροφο fourier


what about imfilter(y,x/y,..) ?

Τι είναι το x/y?

μλκία.. το σωστό που εννοούσα ήταν ifft2(fft2(x)/fft2(y))


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: neon on May 30, 2012, 11:37:17 am
ναι ο συζυγης ειναι!

εβγαλε κανεις αποτελέσματα?? ακομα και για βελτιστο Κ, η xhat μου βγαίνει χειροτερη απο την y..  :-\
Με fft2 και ifft2 παμε στις συχνοτητες και πισω? ή θελει κατι αλλο??

Ναι, αλλά θέλει να κάνεις zero padding στο h για να έχει ίδιο μέγεθος με το y.

Έχει διαφορά αν κάνεις zero padding την h ή την H ; Και αντίστοιχα στην h ή στην G ;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: *** on May 30, 2012, 13:43:13 pm
ναι έχει. αν κάνεις zero padding σε μια ακολουθία, δε παντάρεται αντίστοιχα ο DFT της. μπορεί να μην έχει και κανένα μηδενικό στοιχείο.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: neon on May 30, 2012, 14:25:07 pm
Τι τιμές Κ έχετε γενικώς; Βασικά με τη τιμή ξεκινάτε τις δοκιμές;

Εμένα στον υπολογισμό της G μου βγάζει "Warning: Matrix is close to singular or badly scaled.". Είναι λόγω λάθος Κ;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Selportion on May 30, 2012, 18:13:31 pm
Αντί να κάνεις zero padding στο fft2(h), και μετά συνέλιξη, δεν μπορείς να μην το κάνεις αλλά να εφαρμόσεις το φίλτρο h με την εντολή imfilter(y,h,'circular') όπως το κάνει για να εφαρμόσει τις παραμορφώσεις; Ξέρει κανείς τι διαφορές έχει; Κυκλική συνέλιξη κάνει;

Και δεν διευκρινίζεται καλά: θέλει xinv= φιλτράρισμα της y με το H-1,
αλλά xinv0= φιλτράρισμα της y0 με το H-1 ή με το wiener φίλτρο; Από το όνομα "inv" υποθέτω το 1ο...


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on May 31, 2012, 10:45:24 am
Εγω να ρωτησω, αφου η συνελιξη θα δινει διαστασεις (Μ+Κ-1)x(N+L-1) πως ειναι δυνατον η xhat να εχει διαστασεις MxN?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on May 31, 2012, 10:51:46 am
Εγω να ρωτησω, αφου η συνελιξη θα δινει διαστασεις (Μ+Κ-1)x(N+L-1) πως ειναι δυνατον η xhat να εχει διαστασεις MxN?


Όντως, σύμφωνα με την προηγούμενη εργασία ήταν (Μ+Κ-1)x(N+L-1), αλλά τώρα έκανα μόνο padding στο h ώστε να έχει ίδιο μέγεθος με το y και το xhat.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on May 31, 2012, 11:08:39 am
Εγω να ρωτησω, αφου η συνελιξη θα δινει διαστασεις (Μ+Κ-1)x(N+L-1) πως ειναι δυνατον η xhat να εχει διαστασεις MxN?


Όντως, σύμφωνα με την προηγούμενη εργασία ήταν (Μ+Κ-1)x(N+L-1), αλλά τώρα έκανα μόνο padding στο h ώστε να έχει ίδιο μέγεθος με το y και το xhat.

Αυτο ομως σημαινει οτι κανεις κυκλικη συνελιξη, ομως στην εκφωνηση δινει xhat=g*y που ειναι γραμμικη συνελιξη. Οποτε για να κανεις τον αντιστοιχο πολλαπλασιαμο στη συχνοτητα θελει το κολπακι της δευτερης εργασιας. Εγω παλι δεν κανω πολ/σμο αλλα συνελιξη στο χωρο με conv2 για να μην μπλεκω με αυτα ( :P ) αλλα μεχρι στιγμης παιρνω λαθος αποτελεσματα...   :(


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: johnnysp on May 31, 2012, 11:29:19 am
Η xhat θα πρέπει να είναι ακριβώς ίδια με την x; Ή απλά να την προσεγγίζει;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on May 31, 2012, 12:10:05 pm
Η xhat θα πρέπει να είναι ακριβώς ίδια με την x; Ή απλά να την προσεγγίζει;

Να προσεγγίζει

Εγω να ρωτησω, αφου η συνελιξη θα δινει διαστασεις (Μ+Κ-1)x(N+L-1) πως ειναι δυνατον η xhat να εχει διαστασεις MxN?


Όντως, σύμφωνα με την προηγούμενη εργασία ήταν (Μ+Κ-1)x(N+L-1), αλλά τώρα έκανα μόνο padding στο h ώστε να έχει ίδιο μέγεθος με το y και το xhat.

Αυτο ομως σημαινει οτι κανεις κυκλικη συνελιξη, ομως στην εκφωνηση δινει xhat=g*y που ειναι γραμμικη συνελιξη. Οποτε για να κανεις τον αντιστοιχο πολλαπλασιαμο στη συχνοτητα θελει το κολπακι της δευτερης εργασιας. Εγω παλι δεν κανω πολ/σμο αλλα συνελιξη στο χωρο με conv2 για να μην μπλεκω με αυτα ( :P ) αλλα μεχρι στιγμης παιρνω λαθος αποτελεσματα...   :(

και για το y λέει y = h*x + v , αλλά το h το εφαρμόζει με imfilter με παράμετρο circular , που λογικά δεν είναι κυκλική συνέλιξη?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on May 31, 2012, 12:43:12 pm
και για το y λέει y = h*x + v , αλλά το h το εφαρμόζει με imfilter με παράμετρο circular , που λογικά δεν είναι κυκλική συνέλιξη?

Χμ ενα δικιο το εχεις. Το zero padding το θελει σιγουρα, απλα το μεγεθος ειναι το θεμα, και η πλακα ειναι οτι μου βγαζει χειροτερα αποτελεσματα οταν κανω γραμμικη απ οτι κυκλικη... :-\


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: neon on May 31, 2012, 19:35:37 pm
y: Grayscale εικόνα διάστασης M × N, στο πεδίο του χώρου,
xhat είναι Grayscale εικόνα διάστασης M × N, στο πεδίο του χώρου

γιατί το τονίζει; πρέπει να κάνουμε κάποια ιδιαίτερη μεταχείριση που μου διαφεύγει;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on May 31, 2012, 19:39:33 pm
οχι απλα μαλλον τονιζει οτι πρεπει να παρεις το φουριε τους για να βγει το g και μετα η συνελιξη/πολσμος. Και τελικα να επιστρεφει το σημα στο χωρο και οχι στη συχνοτητα (δηλαδη να γινει ifft)


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: tomshare3 on June 01, 2012, 16:10:57 pm
Τι τιμές Κ έχετε γενικώς; Βασικά με τη τιμή ξεκινάτε τις δοκιμές;


κοντα στο 1250 μου βγαινει το βελτιστο.. αλλα απο 1000 μεχρι 1500 περιπου δεν βλεπω οπτικη διαφορα.. δεν ξερω αν ειναι σωστο αλλα δεν το ψαχνω αλλο!


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Salvation on June 01, 2012, 22:45:35 pm
Στο αποτέλεσμα της Wiener τη σκιακέρια την βλέπετε;

Εγώ παιρνω αυτο:

http://s11.postimage.org/te2ptsh5v/untitled.jpg


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: miaou on June 02, 2012, 12:14:37 pm
το γεγονός ότι τα γκρι κουτάκια είναι μετατοπισμένα έχει σχέση με το πως κάνουμε zero padding την h?

γτ και σε μένα έτσι βγαίνει


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: johnnysp on June 02, 2012, 16:31:57 pm
Σχετικά με την γραφική παράσταση της J ως προς Κ..Σας βγαίνει αύξουσα;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Leo_Paranoid on June 02, 2012, 16:32:59 pm
Έχει καταλάβει κανείς τι θέλει για το xinv0????


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: miaou on June 02, 2012, 16:36:07 pm
Έχει καταλάβει κανείς τι θέλει για το xinv0????

να βάλεις στο φίλτρο σου, όχι την Υ, αλλά την Υ0, τα υπόλοιπα ορίσματα ίδια.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on June 02, 2012, 16:46:42 pm
Σχετικά με την γραφική παράσταση της J ως προς Κ..Σας βγαίνει αύξουσα;

έχει κάποιο ελάχιστο κάπου και μετά ανεβαίνει πάλι


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Leo_Paranoid on June 02, 2012, 16:53:17 pm
Thx man


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: miaou on June 02, 2012, 16:59:04 pm
Thx man
lady plz   :P


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Leo_Paranoid on June 02, 2012, 17:31:55 pm
Σχετικά με την γραφική παράσταση της J ως προς Κ..Σας βγαίνει αύξουσα;

Ναι και εμένα αύξουσα βγαίνει αλλά δεν νομίζω ότι θα έπρεπε.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: johnnysp on June 02, 2012, 18:15:07 pm
Κι εγώ πιστεύω θα έπρεπε να βγει φθίνουσα...


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Leo_Paranoid on June 02, 2012, 19:00:17 pm
Πάρτε το μέτρο του (xhat-xinv0).^2 και θα έχετε φθίνουσα.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: *** on June 02, 2012, 19:05:27 pm
δεν πρέπει να βγαίνει γνησίως μονότονη. Αλλάζει μονοτονία, έχει ελάχιστο.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: johnnysp on June 02, 2012, 19:50:49 pm
Ακόμα και το μέτρο να πάρω βγαίνει αύξουσα..


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on June 02, 2012, 20:16:20 pm
το μεσο τετραγωνικο σφαλμα πως το υπολογιζετε? Με μεσο ορο και προς τις δυο κατευθυνσεις?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: johnnysp on June 02, 2012, 20:37:23 pm
Πήρα τον μέσο όρο όλων..


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: tomshare3 on June 13, 2012, 16:54:23 pm
finally.... ανεβηκε η 4η εργασια!!  προθεσμία 10/7


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: miaou on June 13, 2012, 19:06:23 pm
μόνο σε μένα ανεβαίνει το αίμα στο κεφάλι όταν βλέπω ημερομηνία παραδόσης και αριθμός σελίδων εργασίας?!?!?!?!?!?!?!?!?!?!?!?!?!??!



Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: AgentCain on June 13, 2012, 19:10:22 pm
Εχεις 1 βδομαδα μετα την εξεταστικη  :P


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on June 13, 2012, 23:55:02 pm
Μια απορία που μου ήρθε μετά από γρήγορη ματιά ,

Η συνάρτηση orthonormalFilterBanks παίρνει ως είσοδο h και θέλει να υπολογιστούν τα h0, h1, g0, g1 (στο πεδίο του χρόνου) σύμφωνα με τις εξισώσεις 14-17. Οι εξισώσεις όμως αυτές είναι στο πεδίο του Z. Θέλει δηλαδή να πάρουμε τον αντίστροφο μετασχηματισμό Z των εξισώσεων για να τα υπολογίσουμε?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: averelgr on June 16, 2012, 11:04:09 am
Υπάρχει κανείς που να έχει πρόβλημα και να μην μπορεί να υλοποιήση την εργασία? Δηλαδή να μην μπορεί να ασχοληθεί μετά την εξεταστική και να υποχρεούται έτσι να ασχοληθεί μέσα στην εξεταστική?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: neon on June 16, 2012, 11:57:40 am
Υπάρχει κανείς που να έχει πρόβλημα και να μην μπορεί να υλοποιήση την εργασία? Δηλαδή να μην μπορεί να ασχοληθεί μετά την εξεταστική και να υποχρεούται έτσι να ασχοληθεί μέσα στην εξεταστική?

not really


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: BabylonX on June 17, 2012, 23:45:15 pm
Υπάρχει κανείς που να έχει πρόβλημα και να μην μπορεί να υλοποιήση την εργασία? Δηλαδή να μην μπορεί να ασχοληθεί μετά την εξεταστική και να υποχρεούται έτσι να ασχοληθεί μέσα στην εξεταστική?

Στην τελική πιάσε τον και ζήτα του να την παραδώσεις πιο μετά και να σου περάσει το βαθμό το Σεπτέμβρη. Εγώ αυτό θα έκανα στη θέση σου. Είναι πολύ συνεννοήσιμος άνθρωπος. Εάν μπορεί και δεν έχει κάποιο κώλυμα ο ίδιος για κάποιο λόγο πιστεύω θα σε βοηθήσει με ότι τρόπο μπορεί. Για το Ντελόπουλο μιλάμε.

Ξέρει κανείς πως δίνουμε; Ανοικτά/Κλειστά βιβλία;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on June 17, 2012, 23:49:54 pm
Ξέρει κανείς πως δίνουμε; Ανοικτά/Κλειστά βιβλία;
Μέχρι τώρα όλα τα μαθήματα του Ντελόπουλου ήταν με ανοιχτά. Υποθέτω θα ισχύει και σε αυτό το μάθημα.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: neon on June 18, 2012, 11:57:06 am
Ξέρει κανείς πως δίνουμε; Ανοικτά/Κλειστά βιβλία;
Μέχρι τώρα όλα τα μαθήματα του Ντελόπουλου ήταν με ανοιχτά. Υποθέτω θα ισχύει και σε αυτό το μάθημα.

Ναι ανοιχτά, το είπε. Λάθος thread που τα λέμε αυτά βέβαια..


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Tempus Fugit on June 20, 2012, 01:01:51 am
Μια απορία που μου ήρθε μετά από γρήγορη ματιά ,

Η συνάρτηση orthonormalFilterBanks παίρνει ως είσοδο h και θέλει να υπολογιστούν τα h0, h1, g0, g1 (στο πεδίο του χρόνου) σύμφωνα με τις εξισώσεις 14-17. Οι εξισώσεις όμως αυτές είναι στο πεδίο του Z. Θέλει δηλαδή να πάρουμε τον αντίστροφο μετασχηματισμό Z των εξισώσεων για να τα υπολογίσουμε?

γενικά ανέφερε στην τάξη, οτι αν H(-z) τοτε στο χρονο ειναι (-1)^n * h(n) και αντιστοιχα -H(-z) ειναι (-1)^(n+1) * h(n). Αυτο ομως αναφερονταν στις σχεσεις των QMF. Αν δεις και wikipedia (http://en.wikipedia.org/wiki/Z-transform#Properties) εχει αντιστοιχες ιδιοτητες. Το θεμα μου ειναι πως υλοποιεις το H(z^-1) οταν εχεις πεπερασμενο φιλτρο, δεδομενου οτι το h(-n) δεν εχει νοημα. Αν ασχοληθει κανεις και το βρει ας πει :Ρ


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on June 20, 2012, 11:56:37 am
Μια απορία που μου ήρθε μετά από γρήγορη ματιά ,

Η συνάρτηση orthonormalFilterBanks παίρνει ως είσοδο h και θέλει να υπολογιστούν τα h0, h1, g0, g1 (στο πεδίο του χρόνου) σύμφωνα με τις εξισώσεις 14-17. Οι εξισώσεις όμως αυτές είναι στο πεδίο του Z. Θέλει δηλαδή να πάρουμε τον αντίστροφο μετασχηματισμό Z των εξισώσεων για να τα υπολογίσουμε?

γενικά ανέφερε στην τάξη, οτι αν H(-z^1) τοτε στο χρονο ειναι (-1)^n * h(n) και αντιστοιχα -H(-z) ειναι (-1)^(n+1) * h(n). Αυτο ομως αναφερονταν στις σχεσεις των QMF. Αν δεις και wikipedia (http://en.wikipedia.org/wiki/Z-transform#Properties) εχει αντιστοιχες ιδιοτητες. Το θεμα μου ειναι πως υλοποιεις το H(z^-1) οταν εχεις πεπερασμενο φιλτρο, δεδομενου οτι το h(-n) δεν εχει νοημα. Αν ασχοληθει κανεις και το βρει ας πει :Ρ

Ναι το βρήκα τελικά

Για τις περιπτώσεις που λέει Η(z^-1) πήρα στο χρόνο το φλιπαρισμένο h. Για το g1 πήρα τον τύπο 5 σε αυτό το site http://cnx.org/content/m11064/latest/

Αν δεις και αυτό το pdf που είναι με βάση το βιβλίο του Gonzalez σου λέει και ακριβώς τις τιμές των h0,h1,g0,g1

http://www.cis.rit.edu/class/simg782/lectures/lecture_17/lec782_04_16.pdf

Επίσης το matlab σου βγάζει τους συντελεστές με αυτή τη συνάρτηση
[h0 h1 g0 g1]=wfilters('haar')

αλλά μάλλον δε θέλει αυτό ο Ντελόπουλος  :P


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Tempus Fugit on June 20, 2012, 12:10:36 pm
ωραιος. Κι εγω φλιπαρισμενο σκεφτομουν αλλα δεν μπορουσα να το αποδειξω καπως :Ρ

Θα δουμε, αν δουλεψει...


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: thessalos on June 21, 2012, 23:55:30 pm
Οι εργασιες μετρανε για σεπτεμβρη? Μπορω να κανω καποια που δεν ειχα κανει τωρα?
Ξερει κανεις?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: neon on June 22, 2012, 13:05:40 pm
Στην Εργασία όταν λέει χρησιμοποιήστε την οικογένεια φίλτρων που παράγονται από το φίλτρο Haar, εννοεί ότι στις σχέσεις 14-17 το Ho(z) είναι το (http://latex.codecogs.com/gif.latex?\frac{1}{\sqrt{2}}+\frac{1}{\sqrt{2}}&space;z^-1) ?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Neal on June 22, 2012, 14:07:54 pm
Οι εργασιες μετρανε για σεπτεμβρη? Μπορω να κανω καποια που δεν ειχα κανει τωρα?
Ξερει κανεις?

Στείλε mail στον Ντελό να ρωτήσεις και ενημέρωσε και εμάς.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on June 22, 2012, 14:08:50 pm
Στην Εργασία όταν λέει χρησιμοποιήστε την οικογένεια φίλτρων που παράγονται από το φίλτρο Haar, εννοεί ότι στις σχέσεις 14-17 το Ho(z) είναι το (http://latex.codecogs.com/gif.latex?\frac{1}{\sqrt{2}}+\frac{1}{\sqrt{2}}&space;z^-1) ?

Ναι


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: bard on June 26, 2012, 19:38:03 pm
Ερωτήσεις και (ίσως)Απαντήσεις....
Στις σχέσεις 14-17 λέει ότι Κ είναι το μήκος της κρουστικής απόκρισης h0. Με βάση τον αντίστροφο μετασχηματισμό z, ισχύει:
               g0(n)=h0(-n)   (1)
               g1(n)=h1(-n)   (2)
 και ότι   g1(n)=(-1)2k-n * g0(2k-1-n)=(-1)-n*g0(2k-1-n)  (3)

Λέει δλδ ότι g1(2k-1)=g0(0),g1(2k-2)=-g0(1) κλπ άρα είναι σαν έχουμε φλιπάρει και μετακινήσει την g0, απλά πρέπει να προσέχουμε πότε είναι -1 και πότε +1. Πέρνοντας τις (1) και (2) μπορούμε να βρούμε τις υπόλοιπες....

Το ερώτημα τώρα... είναι σίγουρα Κ το μήκος; μήπως είναι 2Κ; οπότε είναι πάντα σαν να φλιπάρουμε την h0 με χρήση της wrev() της matlab.

Επίσης, πρέπει να κάνουμε κάποια κανονικοποιήση πριν την εφαρμογή;

Εγώ δοκίμασα την παραπάνω ιδέα στο φίλτρο Haar και οι κρουστικές συμπίπτουν με το αποτέλεσμα της wfilters('haar'). Απλά δν ξέρω πως να το δικαιολογίσω.

Τέλος, υπάρχει και η orthfilt() η οποία μάλλον κάνει αυτό που μας ζητείται, έτσι για έλεγχο.... 


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Salvation on June 26, 2012, 19:48:25 pm
Στο βιβλίο του gonzalez λεει οτι το 2Κ είναι το μήκος. Έγώ έτσι θα το πάρω. Δεν εμπιστεύομαι τον Ντελο :D


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: bard on June 29, 2012, 16:33:33 pm
Ωραία, γιατί έτσι έχει περισσότερο νόημα....
Επίσης, χρειάζεται κάποια κανονικοποίηση; Δηλαδή πρέπει οι όροι της κρουστικής να έχουν κάποιο μέτρο ή άθροισμα;....


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Neal on July 02, 2012, 19:45:05 pm
Τα X(z^(1/2)) και X(z^(-1/2)) πως τα υπολογίζουμε στην εξ. (4);

...μπας και πάμε επιτέλους διακοπές ήσυχοι


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: lostakias on July 02, 2012, 20:38:02 pm
Εγώ έχω την εξής απορία: Πώς βρίσκουμε την τιμή του L; Εγώ για να το βρω έκανα δοκιμές με διάφορες τιμές και κατέληξα στην τιμή που μου έδωσε το καλύτερο αποτέλεσμα.Είναι σωστό αυτό ή θέλει κάτι συγκεκριμένο να βάλουμε;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Neal on July 03, 2012, 08:11:58 am
Ας πει ο καθένας την απορία του τότε χωρίς να απαντάει στις προηγούμενες που μπορεί. Θα έχει πολύ νόημα το topic.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on July 03, 2012, 11:29:09 am
Το Κ είναι μάλλον το μήκος του h
Το 2Κ λογικά προκύπτει αν προσθέσεις τα μηδενικά που λέει ο Ντελόπουλος ότι χρειάζονται.

Αλλά επειδή λειτουργεί και χωρίς να βάλεις τα μηδενικά, εγώ έγραψα έναν πιο μπακάλικο, απλοποιημένο τύπο(χωρίς Κ ή L) ώστε να βγαίνουν οι συντελεστές όπως θέλει.

Και με το μάτι να το δεις, από εδώ πχ http://www.cis.rit.edu/class/simg782/lectures/lecture_17/lec782_04_16.pdf  σελ 8
καταλαβαίνεις πως βγαίνει το g1 από το h0



Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: johnnysp on July 04, 2012, 16:40:02 pm
Έβγαλε κανείς άκρη; Το αποτέλεσμά μου είναι 4 μαύρα κουτάκια...  :(


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Salvation on July 04, 2012, 17:01:43 pm
Πριν δώσεις την εικόνα σου ως όρισμα στη συνάρτηση analysis2d γράψε:

farm=double(farm);
farm=farm/255;
imshow(farm);

Το imshow είναι απλά για να επιβεβαιώσεις ότι φαίνεται η εικόνα.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: ktsourap on July 04, 2012, 17:32:22 pm
Μήπως το -1 στον τύπο 15 είναι λάθος ρε παιδιά?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on July 04, 2012, 18:05:36 pm
Πριν δώσεις την εικόνα σου ως όρισμα στη συνάρτηση analysis2d γράψε:

farm=double(farm);
farm=farm/255;

imshow(farm);

Το imshow είναι απλά για να επιβεβαιώσεις ότι φαίνεται η εικόνα.

υπάρχει και η εντολή im2double που κάνει το ίδιο


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: lostakias on July 04, 2012, 20:19:51 pm
Μήπως το -1 στον τύπο 15 είναι λάθος ρε παιδιά?

Είναι δεν είναι λάθος,δεν μας αφορά και τόσο, αφού στο βιβλίο του Gonzalez μας δίνονται έτοιμες οι σχέσεις αυτές στο πεδίο του χρόνου(βλέπε σελ.358, τύπος 7.1-23).Υλοποιείς τον κώδικά σου μ'αυτές τις σχέσεις κι είσαι εντάξει


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: averelgr on July 04, 2012, 20:28:30 pm
Μάλλον το πρόβλημα σου είναι στην imshow() τσέκαρε ότι η είσοδος που παίρνει είναι είτε uint8 δηλαδή από 0-255 είτε double από 0 εως 1..


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: bard on July 05, 2012, 15:33:09 pm
Μήπως το -1 στον τύπο 15 είναι λάθος ρε παιδιά?

Είναι δεν είναι λάθος,δεν μας αφορά και τόσο, αφού στο βιβλίο του Gonzalez μας δίνονται έτοιμες οι σχέσεις αυτές στο πεδίο του χρόνου(βλέπε σελ.358, τύπος 7.1-23).Υλοποιείς τον κώδικά σου μ'αυτές τις σχέσεις κι είσαι εντάξει

Επίσης στους τύπου 7.1-10 και 7.1-11 φαίνεται ότι ούτε το ένα είναι λάθος ούτε το άλλο, αρκεί να ισχύει μία από τις δύο σχέσεις. Με το σκεπτικό αυτό δεν πειράζει να έχεις διαφορετικά πρόσημα σε σχέση με τους συντελεστές που βγαίνουν από το Matlab για το φίλτρο Haar.
Δλδ
h0=[1 1]*0,707         g0=[1 1]*0,707
h1=[-1 1]*0,707        g1=[1 -1]*0,707


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Tempus Fugit on July 06, 2012, 01:45:13 am
Ωραία, γιατί έτσι έχει περισσότερο νόημα....
Επίσης, χρειάζεται κάποια κανονικοποίηση; Δηλαδή πρέπει οι όροι της κρουστικής να έχουν κάποιο μέτρο ή άθροισμα;....
Δε θέλει, γι' αυτό και είναι ορθοκανονικά φίτλρα. Αν το πάρεις σαν διάνυσμα 1x2 το h0=[1,1] τότε ρίζα 2 είναι το μέτρο
του ουσιαστικά, και έχει ήδη διαιρεθεί μ' αυτό.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: phen3x on July 06, 2012, 21:53:04 pm
Το z πώς το ορίζετε ρε σεις; Με syms δεν γίνεται γιατί το βλέπει σαν σύμβολο μετά και τον αντίστροφο μετασχηματισμό. Το όρισμα που δίνουμε στη συνάρτηση πώς χρησιμοποιείται;  :-\ :-\


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: ilithios on July 08, 2012, 12:49:42 pm
Δεν το ορίζουμε... Δες πιο πάνω τα replies, δουλεύουμε στο πεδίο του χώρου


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Salvation on July 08, 2012, 16:39:10 pm
Η εικόνα που θα πάρουμε θα είναι ακριβώς η ίδια με την αρχικη; Δλδ αν παρετε τη διαφορα της farm με τη xhat θα είναι όλο μηδενικα;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Selportion on July 08, 2012, 20:06:44 pm
Υποθέτω πως θα έπρεπε, δηλαδή αφού μιλάμε για πλήρη ανακατασκευή.. Αλλά εφόσον μπορεί να στρογγυλοποιούνται αριθμοί ή να πετιούνται φανταστικά μέρη στους fft/ifft κτλ, μπορεί να έχει απόκλιση.

Εγώ αντιμετωπίζω πρόβλημα στο ότι ενώ η ανάλυση γίνεται κανονικά, στη σύνθεση λόγω της υπερδειγματοληψίας μία στις 2 στήλες και μία στις δύο γραμμές είναι μαύρη ακόμη και μετά το φιλτράρισμα με τα g0,1. Πώς εφαρμόσατε τα φίλτρα; Με filter, με conv ή με zeropadding και πολλαπλασιασμό στοιχείο επί στοιχείο στο πεδίο συχνότητας;

Και μια άλλη σκέψη, όλα αυτά τα βρίσκουμε θεωρώντας δεδομένο πως h0=[1 1] και κανονικοποιημένο. Αν πούμε πώς h0=[1 1 1 1 1]; Ή με ακόμη περισσότερους όρους; Το φίλτρο haar είναι εξορισμού μήκους 2;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on July 08, 2012, 20:41:15 pm
Υποθέτω πως θα έπρεπε, δηλαδή αφού μιλάμε για πλήρη ανακατασκευή.. Αλλά εφόσον μπορεί να στρογγυλοποιούνται αριθμοί ή να πετιούνται φανταστικά μέρη στους fft/ifft κτλ, μπορεί να έχει απόκλιση.

Εγώ αντιμετωπίζω πρόβλημα στο ότι ενώ η ανάλυση γίνεται κανονικά, στη σύνθεση λόγω της υπερδειγματοληψίας μία στις 2 στήλες και μία στις δύο γραμμές είναι μαύρη ακόμη και μετά το φιλτράρισμα με τα g0,1. Πώς εφαρμόσατε τα φίλτρα; Με filter, με conv ή με zeropadding και πολλαπλασιασμό στοιχείο επί στοιχείο στο πεδίο συχνότητας;

Και μια άλλη σκέψη, όλα αυτά τα βρίσκουμε θεωρώντας δεδομένο πως h0=[1 1] και κανονικοποιημένο. Αν πούμε πώς h0=[1 1 1 1 1]; Ή με ακόμη περισσότερους όρους; Το φίλτρο haar είναι εξορισμού μήκους 2;

Το h0 είναι [1/riza2 1/riza2] και είναι η κρουστική απόκριση του φίλτρου Haar H0. Κανονικά έχει άπειρο μήκος
[1/riza2 1/riza2 0 0 0 0 0 0 0 ...], αλλά μόνοι οι δύο πρώτοι όροι είναι διάφοροι του μηδενός (λογικό, γιατί αλλιώς θα ήταν ασταθές φίλτρο ). Γι'αυτό είναι μήκους 2.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Merlin on July 08, 2012, 21:04:45 pm
Το zero padding γινεται στην αρχη,στο τελος η ειναι αδιαφορο?Επισης για τα 1Δ φιλτρα h0,h1,g0,g1 χρησιμοποιουμε fft και ifft ενω για τις x,y0,y1 και y00 ,y01, y10, y11 fft2 και ifft2 σωστα?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Selportion on July 08, 2012, 21:50:11 pm
Ναι σωστά, αλλά αν υλοποιήσεις τις 2d συναρτήσεις καλώντας πολλές φορές τις 1d, νομίζω σε όλη την εργασία δεν θα χρειαστεί fft2 ifft2. To zeropadding εξαρτάται σε πιο σημείο και γιατί το κάνεις. Μπορεί να γίνει αν δεν κάνω λάθος απλά και με όρισμα fft(pinakas, diastasi1, diastasi2), ώστε προσθέτει μηδενικά μέχρι να έχει το ζητούμενο μέγεθος.
Γιαυτό επανέρχομαι αν κάποιος γνωρίζει:
Quote
αντιμετωπίζω πρόβλημα στο ότι ενώ η ανάλυση γίνεται κανονικά, στη σύνθεση λόγω της υπερδειγματοληψίας μία στις 2 στήλες και μία στις δύο γραμμές είναι μαύρη ακόμη και μετά το φιλτράρισμα με τα g0,1. Πώς εφαρμόσατε τα φίλτρα; Με filter, με conv ή με zeropadding και πολλαπλασιασμό στοιχείο επί στοιχείο στο πεδίο συχνότητας;


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: tati on July 08, 2012, 22:49:28 pm
εγω το έκανα με filter και δεν ειχα τέτοιο πρόβλημα ..


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Tempus Fugit on July 09, 2012, 15:09:35 pm
Ναι σωστά, αλλά αν υλοποιήσεις τις 2d συναρτήσεις καλώντας πολλές φορές τις 1d, νομίζω σε όλη την εργασία δεν θα χρειαστεί fft2 ifft2. To zeropadding εξαρτάται σε πιο σημείο και γιατί το κάνεις. Μπορεί να γίνει αν δεν κάνω λάθος απλά και με όρισμα fft(pinakas, diastasi1, diastasi2), ώστε προσθέτει μηδενικά μέχρι να έχει το ζητούμενο μέγεθος.
Γιαυτό επανέρχομαι αν κάποιος γνωρίζει:
Quote
αντιμετωπίζω πρόβλημα στο ότι ενώ η ανάλυση γίνεται κανονικά, στη σύνθεση λόγω της υπερδειγματοληψίας μία στις 2 στήλες και μία στις δύο γραμμές είναι μαύρη ακόμη και μετά το φιλτράρισμα με τα g0,1. Πώς εφαρμόσατε τα φίλτρα; Με filter, με conv ή με zeropadding και πολλαπλασιασμό στοιχείο επί στοιχείο στο πεδίο συχνότητας;

πρώτα προσθέτεις μηδενικά ανά μια θέση, μετά φιλτράρεις και τα δύο κομμάτια,
και μετά απλή πρόσθεση. Την πρόσθεση κάνω με filter(g,1,w)
Προσπάθησε να ελέγξεις αν η τετράδα φίλτρων που υπολογίζεις είναι η σωστή, εκτυπώνοντας τις τιμές.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: lynx on July 09, 2012, 17:34:01 pm
τελικά γιατί δουλεύετε στο πεδίο του χώρου και όχι του z;

τα ορίσματα των συναρτήσεων είναι στο πεδίο του χώρου;

αν ναι, στην orthonormalFilterBanks(h) τι δίνω ως είσοδο; (για το φίλτρο Haar)


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: ktsourap on July 10, 2012, 01:07:35 am
τελικά πόσα μηδενικά θέλει να προσθέσουμε στην orthonormalFilterBanks(h) και που??
εμένα χωρις μηδενικά δουλεύει κομπλέ!! αν βάλω όπου και να τα βάλω μου δημιουργεί αντίστοιχες πρώτες μηδενικές στήλες και γραμμές στην τελική ανακατασκευασμένη εικόνα!!


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: ktsourap on July 10, 2012, 01:08:46 am
τελικά γιατί δουλεύετε στο πεδίο του χώρου και όχι του z;

τα ορίσματα των συναρτήσεων είναι στο πεδίο του χώρου;

αν ναι, στην orthonormalFilterBanks(h) τι δίνω ως είσοδο; (για το φίλτρο Haar)

h=[0.7071    0.7071] απλά και όλα στο πεδίο του χρόνου


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: miaou on July 10, 2012, 10:34:45 am
Λευτεριά σε όσους γύρισαν από Ποσείδι για την εργασία..




Άλλη μια παράτασ12.7.2012


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on July 10, 2012, 15:05:43 pm
παιδια να ρωτησω στην orthonormalFilterBanks το g1 ειναι αρνητικο στα μονα και θετικο στα ζυγα λογω του (-1)-n ετσι δν ειναι? Οταν ομως τρεχω τη συναρτηση του matlab για επαληθευση μου βγαζει οτι στο g1 τα ζυγα ειναι τα αρνητικα και τα μονα τα θετικα :-\


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Tempus Fugit on July 10, 2012, 15:29:39 pm
παιδια να ρωτησω στην orthonormalFilterBanks το g1 ειναι αρνητικο στα μονα και θετικο στα ζυγα λογω του (-1)-n ετσι δν ειναι? Οταν ομως τρεχω τη συναρτηση του matlab για επαληθευση μου βγαζει οτι στο g1 τα ζυγα ειναι τα αρνητικα και τα μονα τα θετικα :-\

γιατι εις την -n και οχι στην n;
Οι ασσοι που χρησιμοποιησα εγω ηταν [1 -1 1 -1 ...]

1. Κανεις invert να βρεις το G0(z-1)
2. Πολλαπλασιαζεις με τους ασσους να βρεις το G0(-z-1)
3. Αλλάζεις πρόσημο για να βρεις το -G0(-z-1)

Στο χρονο ολα αυτα δλδ.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on July 10, 2012, 15:35:40 pm
το -n βγαινει απο το οτι ειναι αρχικα 2k-1-n. Ομως υπαρχει και το ολικο - (αυτο που λες στο βημα 3) αρα γινεται 2k-1-n+1=2k-n. Το 2k βγαινει παντα ζυγος αρα δν μας ενδιαφερει ετσι μενει το -n το οποιο βεβαια ειναι το ιδιο με το n σε αυτη την περιπτωση. Με αυτο που λες ομως και συ βγαζεις αρνητικα τα μονα και θετικα τα ζυγα (αν αντιστρεψεις τον πινακα των ασσων) ομως με την wfilters βγαινει το αντιθετο.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on July 10, 2012, 16:47:03 pm
βγαινει και σε κανεναν αλλον το y00 να εχει την περισσοτερη πληροφορια και τα υπολοιπα τρια να ειναι σχεδον μαυρα με μερικο ασπρο ιχνος?


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Salvation on July 10, 2012, 16:53:40 pm
Ναι το y00 έχει την περισσότερη πληροφορία γιατι είναι και το Low pass. Όσο για το orthnormals επιβεβαίωσε με αυτό:

http://www.mathworks.com/help/toolbox/wavelet/ref/orthfilt.html


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on July 10, 2012, 17:01:19 pm
ε το εκανα με την wfilter το ιδιο ειναι. Δλδ εσενα βικ σου βγαινουν στην g1 θετικα τα μονα? Δν μπορω να καταλαβω που γινεται το λαθος στους υπολογισμους με τους τυπους. Νομιζω το ειπανε πιο πριν στο τοπικ οτι παιζει να μη θελει το ενα -. Απο την αλλη το εχω τρεξει και με τις δυο περιπτωσεις και οι οπτικες διαφορες ειναι απειροελαχιστες.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Salvation on July 10, 2012, 17:05:33 pm
h0 =

   0.70711
   0.70711

h1 =

  -0.70711
   0.70711

g0 =

   0.70711
   0.70711

g1 =

   0.70711
  -0.70711

Αυτά βγάζω. Τα έβγαλα σύμφωνα με του τύπους του gonzalez στο βιβλίο.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Time_Tripper on July 10, 2012, 17:23:43 pm
χμ τελικα το προβλημα ηταν θεωρητικο παρα υπολογιστικο. το matlab ξεκιναει απο το 1 αλλα το φιλτρο απο το 0 οποτε το 1ο του ματλαμπ πρεπει να το αντιμετωπισουμε ως το 0 και ετσι βγαινουν αναποδα τα προσημα...

εδιτ: παντως οι τυποι του gonzalez ειναι ιδιοι με του ντελο οποτε δν υπαρχει θεμα


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Merlin on July 10, 2012, 18:47:52 pm
τελικά πόσα μηδενικά θέλει να προσθέσουμε στην orthonormalFilterBanks(h) και που??
εμένα χωρις μηδενικά δουλεύει κομπλέ!! αν βάλω όπου και να τα βάλω μου δημιουργεί αντίστοιχες πρώτες μηδενικές στήλες και γραμμές στην τελική ανακατασκευασμένη εικόνα!!

Εγώ το έκανα να κάνει τόσο zero padding , όσο χρειάζεται για το φιλτράρισμα της εκάστοτε στήλης / γραμμής της εικόνας.Δηλαδή την καλώ εκτός από το demo και στις synthesis2d και analysis2d.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Merlin on July 10, 2012, 18:55:28 pm
Έχει σημασία  σε ποια συνάρτηση θα κάνουμε αναδίπλωση συχνότητας(φίλτρο G ή εικόνα W)?Είτε με G(-z),είτε με W(-z),έχω το ίδιο ακριβώς αποτέλεσμα.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Selportion on July 10, 2012, 20:17:36 pm
εγω το έκανα με filter και δεν ειχα τέτοιο πρόβλημα ..
Ναι όντως και σε μένα με filter και έπιασε, με την y0 να έχει σχεδόν όλη την πληροφορία όπως λένε τα παιδιά.

Η εικόνα που θα πάρουμε θα είναι ακριβώς η ίδια με την αρχικη; Δλδ αν παρετε τη διαφορα της farm με τη xhat θα είναι όλο μηδενικα;
Η διαφορά τους φαίνεται να φθάνει σε μια δοκιμή που έκανα μέχρι και 200, αλλά ο μέσος όρος της είναι της τάξης του e-014. Οπτικά βλέπω λίγο χαμηλότερο contrast και μια απαλή υφή παραμορφωσης, αλλά πάλι ίσως απλά κάτι κάνω λάθος. Βασικά έχουν διαφορετικό μέγιστο, έτσι όταν διαιρώ, διαιρώ με διαφορετικό μέγιστο την x και την xhat για να είναι όλες οι τιμές στο [0 1] και να τις δω.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: lynx on July 11, 2012, 02:17:29 am
αφού το Κ είναι το μέγεθος της h0 και αντίστοιχα της g0,

ο τύπος g1(n) = (-1)^n*g0(2K - 1 - n)

δεν σημαίνει ότι η g1 έχει μήκος 2Κ;

όμως, στην εικόνα 7.6 του βιβλίου δεν φαίνεται κάτι τέτοιο..


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on July 11, 2012, 02:38:46 am
αφού το Κ είναι το μέγεθος της h0 και αντίστοιχα της g0,

ο τύπος g1(n) = (-1)^n*g0(2K - 1 - n)

δεν σημαίνει ότι η g1 έχει μήκος 2Κ;

όμως, στην εικόνα 7.6 του βιβλίου δεν φαίνεται κάτι τέτοιο..

https://www.thmmy.gr/smf/index.php?topic=50776.msg887470#msg887470


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Tempus Fugit on July 11, 2012, 02:55:36 am
Η διαφορά τους φαίνεται να φθάνει σε μια δοκιμή που έκανα μέχρι και 200, αλλά ο μέσος όρος της είναι της τάξης του e-014. Οπτικά βλέπω λίγο χαμηλότερο contrast και μια απαλή υφή παραμορφωσης, αλλά πάλι ίσως απλά κάτι κάνω λάθος. Βασικά έχουν διαφορετικό μέγιστο, έτσι όταν διαιρώ, διαιρώ με διαφορετικό μέγιστο την x και την xhat για να είναι όλες οι τιμές στο [0 1] και να τις δω.
Αν τις κανεις uint8 και δεν διαιρεσεις (οταν δεν ειναι double δεν υπαρχει προβλημα να τις εμφανισει σε κλιμακα 0-255 στο imshow) τοτε βγαινουν ιδιες. Σε μορφη double το σφαλμα ειναι θεμα του ελαχιστου σφαλματος του ματλαμπ.


Ναι το y00 έχει την περισσότερη πληροφορία γιατι είναι και το Low pass. Όσο για το orthnormals επιβεβαίωσε με αυτό:

http://www.mathworks.com/help/toolbox/wavelet/ref/orthfilt.html
Αυτά βγάζω. Τα έβγαλα σύμφωνα με του τύπους του gonzalez στο βιβλίο.


Νομίζω απ' ότι είχα δει σε ένα pdf που έλεγε τις λύσεις ότι τα g1 και h1 μπορουν να αντισταφούν από θέμα προσήμου στο
συγκεκριμένο φίλτρο. Αν θυμαμαι καλα τα είχα δει αναποδα απ' ότι τα βγαζα δλδ το ενα να εχει μειον στο πρωτο tap κ το
αλλο στο δευτερο κ αντιστροφα, και δουλευαν κ οι 2 εκδοχες.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: lynx on July 11, 2012, 03:44:34 am
αφού το Κ είναι το μέγεθος της h0 και αντίστοιχα της g0,

ο τύπος g1(n) = (-1)^n*g0(2K - 1 - n)

δεν σημαίνει ότι η g1 έχει μήκος 2Κ;

όμως, στην εικόνα 7.6 του βιβλίου δεν φαίνεται κάτι τέτοιο..

https://www.thmmy.gr/smf/index.php?topic=50776.msg887470#msg887470

ναι το είχα δει αυτό το ποστ σου (πολύ χρήσιμο btw, ευχαριστούμε!)

και πάνω σε αυτό μου ήρθε η απορία..

απλά αναρωτιέμαι δηλ, γιατί το g1 να έχει διπλάσιο μέγεθος από τα υπόλοιπα..


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on July 11, 2012, 12:00:00 pm
αφού το Κ είναι το μέγεθος της h0 και αντίστοιχα της g0,

ο τύπος g1(n) = (-1)^n*g0(2K - 1 - n)

δεν σημαίνει ότι η g1 έχει μήκος 2Κ;

όμως, στην εικόνα 7.6 του βιβλίου δεν φαίνεται κάτι τέτοιο..

https://www.thmmy.gr/smf/index.php?topic=50776.msg887470#msg887470

ναι το είχα δει αυτό το ποστ σου (πολύ χρήσιμο btw, ευχαριστούμε!)

και πάνω σε αυτό μου ήρθε η απορία..

απλά αναρωτιέμαι δηλ, γιατί το g1 να έχει διπλάσιο μέγεθος από τα υπόλοιπα..

Αν χρησιμοποιήσεις μηδενικά , όλα θα έχουν διπλάσιο μέγεθος

Δηλαδή , λες πχ  h0 = [0 0 h], άρα το h0 θα έχει τώρα μέγεθος 2K και τα υπόλοιπα θα έχουν και αυτά διπλάσιο μέγεθος αφού προκύπτουν από το h0. Απ'ότι δοκίμασα βγάζει τα ίδια αποτελέσματα χωρίς μηδενικά.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: lynx on July 11, 2012, 12:35:56 pm


Αν χρησιμοποιήσεις μηδενικά , όλα θα έχουν διπλάσιο μέγεθος

Δηλαδή , λες πχ  h0 = [0 0 h], άρα το h0 θα έχει τώρα μέγεθος 2K και τα υπόλοιπα θα έχουν και αυτά διπλάσιο μέγεθος αφού προκύπτουν από το h0. Απ'ότι δοκίμασα βγάζει τα ίδια αποτελέσματα χωρίς μηδενικά.

α, μάλιστα, αυτό με τα μηδενικά βγαίνει παίρνοντας τον αντίστροφο Μ/Τ z; (για να βγάλεις το h0 από το Η0)


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: Emfanever on July 11, 2012, 13:31:29 pm


Αν χρησιμοποιήσεις μηδενικά , όλα θα έχουν διπλάσιο μέγεθος

Δηλαδή , λες πχ  h0 = [0 0 h], άρα το h0 θα έχει τώρα μέγεθος 2K και τα υπόλοιπα θα έχουν και αυτά διπλάσιο μέγεθος αφού προκύπτουν από το h0. Απ'ότι δοκίμασα βγάζει τα ίδια αποτελέσματα χωρίς μηδενικά.

α, μάλιστα, αυτό με τα μηδενικά βγαίνει παίρνοντας τον αντίστροφο Μ/Τ z; (για να βγάλεις το h0 από το Η0)

Το απλό h είναι η κρουστική απόκριση του H0 και βγαίνει με τη συνάρτηση filter αν βάλεις είσοδο [1 0 0 0 0 ...] ή απευθείας με την impz. Τα μηδενικά εγώ τα έβαλα μόνος μου για να ισχύουν οι τύποι ακριβώς όπως τους έχει , αλλά μάλλον δε χρειάζονται.


Title: Re: [ΨΕΕ] Εργασιες 2012
Post by: lynx on July 11, 2012, 14:02:54 pm
αα, οκ!

ευχαριστώ!