THMMY.gr

Για όσους έχουν αρχίσει να ασχολούνται με την 1η εργασία..

Στην εκφώνηση δίνεται η αρχική εικόνα gray scale x στο σχήμα 1.
Πρέπει να δουλέψουμε με αυτήν ή κάθε φορά θα παράγουμε δική μας τυχαία εικόνα στο Matlab?
Αν πρέπει να δουλέψουμε με αυτήν με ποιον τρόπο θα την πάρουμε σε μορφή πίνακα στο Matlab ?

Έχει βγάλει κανείς άκρη; Την εικόνα θα τη σχεδιάσουμε με βάση τον τύπο x(n1,n2) και τα α,ω,b είπε θα τα βρούμε από το διάγραμμα.

Update:

Έχει καταλάβει κανείς τι πρέπει να κάνουμε στο 1ο ερώτημα ? :P

Όχι ακριβώς :P
Μόνο παίζω με τις συχνότητες για να πετύχω το σχημα  ::)

Μήπως να του πούμε να δώσει τις τιμές να τελειώνουμε; Τζάμπα παιδευόμαστε. Η τουλάχιστον να δώσει τιμες στα διαγράμματα μήπως από εκεί βγει κάτι

Παιδια για να μπορω να κανω την εργασια τι γνωσεις πρεπει να εχω; και περιπου μεχρι που πρεπει να διαβασω στις σημειωσεις; Με καλυπτουν η πρεπει να ξερω και κατι εξτρα; (δεν πατησα σε κανενα μαθημα)

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

Έχει καταλάβει κάνεις πως ακριβώς βρίσκουμε τις τιμές απ τα διαγράμματα τελικα????

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

Είπε ότι θα δώσει παράταση έως την Κυριακή..
Μια μέρα δηλαδή αν δεν κάνω λάθος..

Νομίζω ότι οι αποστάσεις των κορυφών από το κέντρο του γραφήματος (π.χ. στο heatmap που δίνει)
σου δίνουν τις τιμές για τα u1,u2 που θα χρησιμοποιήσεις στους τύπους για τα ω. Τώρα τα a1,a2 και b
φαντάζομαι τα βάζεις στην τύχη μέχρι να πετύχεις περίπου ίδια εικόνα και μέτρο για τον  μετασχηματισμό Fourier της εικόνας. Αν θέλει ας επιβεβαιώσει και κάποιος άλλος.

Όποιος θέλει να κάνει το Heatmap  δεν ειναι απλα Heatmap αλλα colormap και γίνετε με :

Ωραίος  ;) Θυμάται κανείς τι παίζει με τον πίνακα H για να έχει η τελική εικόνα y πραγματικές τιμές ?

Στο s(n1, n2) = sin(ω1*n2)
sin(ω2*n2) είναι και τα 2 n2?

Το πρώτο είναι ω1 απλά έχει ξεχάσει να το διορθώσει..

Μήπως εννοείς n1?

Ναι προφανώς :P Μπερδεύτηκα !

Παιδιά κανείς που κατάλαβε ποιος θα είναι ο Η και τί θα κάνουμε για να βρούμε τη y?

Επειδή βλέπω δεν υπάρχει κινητικότητα στο θέμα, εγώ την τελείωσα (τώρα για την απόλυτη ορθότητά της δε ξέρω). Αν κάποιος χρειάζεται βοήθεια ας μου στείλει pm :)

το σχημα 1 το κανουμε export στην matlab και απο κει βρισκουμε X?

Δεν κάνεις κάποιο export από την εκφώνηση. Βρίσκεις τις κατάλληλες μεταβλητές u1,u2 από
το μέτρο του DFT που σου δίνει (ή από το heatmap) και προσπαθείς να βγάλεις όσο γίνεται μια εικόνα με το ίδιο μέτρο DFT, "πειράζοντας" και τα a1, a2, b.

Το μέτρο του DFT που το δίνει; Από το heatmap πως το καταλαβαίνουμε;
Τις μεταβλητές u1,u2 τις δίνει. Οι ω11,ω12,ω21,ω22 πως βγαίνουν;

απ οτι εχω καταλάβει βάζεις στο ω11 ας πουμε οπου υ1 ,συμφωνα με τον τυπο που δίνεται, την αποστάση απ το κεντρο μιας κορυφης ως προς τον αξονα υ1, αντιστοιχα μετα για τον αλλον αξονα για το ω12. Έπειτα κάνεις το ίδιο και για μια απ τις εσωτερικές κορυφες

"Να προσδιορίσετε ποιο ή ποια στοιχεία του πίνακα X αντιστοιχούν σε ένα διδιάστατο πραγμα-
τικό ημίτονο της μορφής" σ αυτο τι κανατε ?

+1

Προθεσμια υποβολης ειναι 5/4 βραδυ υποθετω ετσι;

ενας που τ κανε?με max κ ind2sub το κανατε?

Ναι βράδυ, αλλα κυριακή δεν ειχε πει ; αφου μας εδωσε ακομα μια ολοκληρη μερα ::)

Μέχρι βράδυ Κυριακής είναι νομίζω!

Στο 2ο ερώτημα της εργασίας λέει στο γ:"να επιδεικνύει ότι όντως παράγεται το σωστό γράφημα". Αυτό πως ακριβώς θα το κάνουμε?

Ρε παιδια μπορει καποιος να εξηγησει ακριβως πως βρισκουμε τα ω11 ω12 ω21 ω22 α1 α2 και β απο το διαγραμμα επειδη δεν εχω καταλαβει τιποτα; 

Δοκίμασε να βάζεις διάφορες τιμές στα ω11 κλπ και να φτιάχνεις το διάγραμμα του σχήματος 2. Θα δεις ότι ανάλογα τις τιμές μετακινούνται οι ώσεις.

το x(n1,n2)=a1sin(w11n1)sin(w12*n2)+a2sin(w21n1)sin(w22n2)+b πρεπει να το γραψω ετσι στο matlab; επειδη το γραφω ετσι βαζω τιμες στα ω11 ω12 ω21 ω22 α1 α2 και β, n1=256 και n2=364 και μετα βαζω τις εντολες που μας εχει στην εργασια και μου βγαζει error γιατι δεν ειναι πινακας το Χ. Πως μπορω να το x(n1,n2)=a1sin(w11n1)sin(w12*n2)+a2sin(w21n1)sin(w22n2)+b να το κανω πινακα;

Εχω σκαλωσει ασχημα και ειμαι εντελως ασχετος. Διαβασα απο τις διαφανειες του και διαβασα και καποια πραγματα απο το ιντερνετ για Digital Image Processing στο Matlab. Καποιος που μπορει να μου δωσει καιποιες κατευθυντηριες γραμμες;

Για να το κάνεις πίνακα γράφεις:

n1=256;
n2=364;
for i=1:n1
    for j=1:n2
        u(i,j)=a1sin(w11n1)sin(w12*n2)+a2sin(w21n1)sin(w22n2)+b
    end
end

Μπορεί κάποιος να εξηγήσει πως βρίσκουμε το μέτρο του dft του sin(ω1ν1)sin(ω2ν2):

Καποιος τις τιμες ω11 ω12 ω21 ω22 α1 α2 και β επειδη εχουν σπασει τα νευρα μου η τουλαχιστον την μεθοδολογια;
Δεν βγαζω με τιποτα το διαγραμμα. Σαν ω11 πιανω την αποσταση της μιας εξωτερικης κορυφης απο το υ1 και το ω12 απο το υ2 και μετα το ω21 ειναι η αποσταση απο τη μιας εσωτερικης κορυφης απο το υ1 και το ω22 απο το υ2; το α1 ειναι το υψος της εξωτερικης κορυφης, το α2 το υψος της εσωτερικης και το β της κεντρικης της μεγαλης;

ΔΕΝ ΒΓΑΖΩ ΑΚΡΗ!

Δεν παίρνεις τις αποστάσεις για τα ω11.. αλλά για τα u11... και μετά τα αντικαθιστάς στους τύπους που συνδέει τα ω με τα u.

στο ζητουμενο 3 κανατε padding στην εικονα χ,χρειαζεται?

Ναι, θα χρειαστεί και στην x και στο h

Μπορεί κάποις να βοηθήσει στο πως βρίσκουμε το dft του sin(ω1n1)*sin(ω2*n2)? Έχω κολλήσει για τα καλά.

Με την εντολή fft

για ω1,ω2 θα βαλεις συγκεκριμενα u1,u2 οτι 8ες απ το range ποθ σου λεει
μετα X=fft2(x),Xc=fftshift(X),mesh(abs(Xc)), παιξε με τιμες υ1,υ2 κ θα καταλαβεις

Στο ζητουμενο 3 εχω καταλαβει οτι πρεπει να μεγαλωσει το μετρο των εσωτερικων κορυφων του αρχικου Χ και στην ουσια να εξαλειφτουν οι εξωτερικες κορυφες(απο το διαγραμμα του DFT του), Με ενα βαθυπερατο φιλτρο Η θα μπορουσε να γινει αυτο. Λογικα με συνελιξη στο χωρο της συχνοτητας Υ=Χ*Η. Απλα δεν ξερω πως να το κανω στο matlab. Τι τιμες να βαλω στο Η και που ετσι ωστε να το πετυχω. Καποιος καποια βοηθεια;

Αν εννοείς με Matlab σου απάντησαν οι συνάδελφοι. Αν εννοείς με το χέρι (άρα το 1ο ζητούμενο)
εγώ πήρα τον ορισμό του 2D-DFT με τα 2 αθροίσματα. Κι εγώ όμως έχω απορία για το μετά.
Μπορούμε να σπάσουμε το διπλό άθροισμα στους επιμέρους όρους και να βγει ότι τελικά ο 2D-DFT ισούται με
το γινόμενο DFT(sin(ω1n1))DFT(sin(ω2n2)) ???

ναι με σταθερα ω1,ω2 γινεται

Το φιλτράρισμα στο πεδίο των συχνοτήτων υλοποιείται με πολλαπλασιασμό στοιχείο προς στοιχείο Υ=ΗX και όχι
με συνέλιξη. Όπως λες πάρε ένα βαθυπερατό φίλτρο που να κόβει τις συχνότητες που πρέπει και να κρατάει τον dc όρο και τις τέσσερις κοντινές συχνότητες. Το πού ακριβώς θα βάλεις τους άσσους έχει να κάνει με το πώς θα υλοποιήσεις το φιλτράρισμα (την freqfilter.m δηλαδή).

για τα a1,a2 χρησιμοποιησατε οτι χ(η1,η2)=[0,1] για να τα βρειτε η
μονο το γραφημα?

Έχει κανείς ιδέα τι  κάνουμε στο 2γ?