THMMY.gr

Μαθήματα Κύκλου Ηλεκτρονικής & Υπολογιστών => Ψηφιακή Επεξεργασία Εικόνας => Topic started by: Ροζ συννεφάκι on April 13, 2017, 10:34:03 am

Title: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Ροζ συννεφάκι on April 13, 2017, 10:34:03 am

Η εκφώνηση της Εργασίας #2 έχει αναρτηθεί στην ενότητα "Υλικό μαθήματος" στο ethmmy.
Προθεσμία υποβολής: 27/04/2017

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Indy on April 13, 2017, 18:24:59 pm

Τί πρέπει να διαβάσουμε (και από πού) για την εργασία 2; Δεν έχω ιδέα για ιστογράμματα κλπ

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: gdims on April 13, 2017, 22:06:30 pm

Quote from: Indy on April 13, 2017, 18:24:59 pm

Τί πρέπει να διαβάσουμε (και από πού) για την εργασία 2; Δεν έχω ιδέα για ιστογράμματα κλπ

Η εργασία δεν φαίνεται να ζητάει πολλά πράγματα από θεωρία. Πρέπει να καταλαβαίνεις τι σημαίνει το ιστόγραμμα και τι είναι ένας σημειακός μετασχηματισμός. Ρίξε μια ματιά στις διαφάνειες του Ντελόπουλου για όλα αυτά, και κάνε αναζητήσεις στο Google όπως "Histogram of an image".
Θα βγάλεις γρήγορα άκρη.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 14, 2017, 01:18:29 am

Quote from: gdims on April 13, 2017, 22:06:30 pm

Ασχολήθηκε κανείς με το ερώτημα 2.2;
Πώς μπορούμε να περάσουμε μαζί με την συνάρτηση πχ @normpdf και ορίσματα όπως μεση τιμή και τυπική απόκλιση;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: gdims on April 14, 2017, 14:08:39 pm

Παιδιά στο 2.1, μπορείτε να δείξετε κανένα ιστόγραμμα της εικόνας εξόδου;

Εγώ δοκίμασα να παράξω ιστόγραμμα:

h = normpdf(v, 0.5, 0.4) / sum(normpdf(v, 0.5, 0.4));

Δηλαδή κανονική κατανομή με μέση τιμή 0.5 και τυπική απόκλιση 0.4, και πήρα το ακόλουθο figure:

(https://preview.ibb.co/dnz5t5/hist.png) (https://ibb.co/gWk0RQ)

Οι μπάρες είναι το ιστόγραμμα της εικόνας εξόδου, η καμπύλη είναι η ιδανική περίπτωση. Με προβληματίζει λίγο το τμήμα στα δεξιά αν και νομίζω ότι είναι λόγω greedy προσέγγισης.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 14, 2017, 16:44:10 pm

To διάγραμμα που ζήτησες .
(http://i67.tinypic.com/2u8cow8.jpg)

και για μ=0.5 και σ=0.1
(http://i67.tinypic.com/2crr05x.jpg)

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: gdims on April 14, 2017, 17:02:50 pm

Ευχαριστώ όμως δεν μπορώ να βγάλω ξεκάθαρο συμπέρασμα...
Είναι λίγο χοντρές οι γραμμές και δεν βλέπω λεπτομέρειες.
Πόσο μεγάλο L έβαλες;
Σου κάνει τέλεια προσέγγιση;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 14, 2017, 17:05:24 pm

Σορρυ πριν έκανα λάθος, είχα την εικόνα για σ=0.1 , με μπέρδεψες γιατί πουθενά δεν ζητάει σ=0.4
Τώρα αυτό ειναι για 100 διαστήματα. Αμα θέλεις να ανεβάσω και για 20 που έχεις εσύ
http://imgur.com/eBR7FuK

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Goemon on April 14, 2017, 17:08:53 pm

Quote from: Indy on April 13, 2017, 18:24:59 pm

Τί πρέπει να διαβάσουμε (και από πού) για την εργασία 2; Δεν έχω ιδέα για ιστογράμματα κλπ

Το κεφάλαιο 3 του Gonzales είναι επίσης είναι λίγο πολύ αυτά που έλεγε ο Ντελόπουλος. Εγώ αυτό διαβάζω τώρα

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: gdims on April 14, 2017, 17:11:53 pm

Ναι, εγώ απλά δοκιμές κάνω, για να δω πόσο καλά προσεγγίζεται το ιστόγραμμα.
Και αναφέρομαι μόνο στο 2.1, δηλαδή οι εντολές:

L = 100;
v = linspace(0, 1, L);
h = normpdf(v, 0.5, 0.1) / sum(normpdf(v, 0.5, 0.1));
y = histtransform(x,h,v);

σου έδωσαν το ιστόγραμμα που μου έδειξες;
Σίγουρα χρησιμοποιείς greedy αλγόριθμο;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 14, 2017, 17:16:36 pm

Quote from: gdims on April 14, 2017, 17:11:53 pm

Σίγουρος δεν είμαι, αλλα το έκανα με την λογική που γράφει στο pdf.
~~Βάζοντας αυτές τις εντολές δεν παίρνω σωστό αποτέλεσμα~~. Εγω τα γραφίματα που σου έδειξα πριν τα έκανα μέσω του ερωτήματος 2.2 χρησιμοποιώντας και την συνάρτηση histtransform
Βασικά άκυρο, το ίδιο αποτέλεσμα παίρνω. Κάτι είχα γράψει λάθος πρίν

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: gdims on April 14, 2017, 17:20:33 pm

Δεν έχω καταλάβει ακριβώς το σκοπό του 2.2 ακόμα, αλλά μάλλον στοχεύει στο να βελτιώσει τη προσέγγιση. Θα ήθελα να δω το αποτέλεσμα που παίρνεις με τις εντολές που σου έστειλα, ή με τις εντολές που γράφει η εκφώνηση στο τέλος του 2.1.

EDIT: α, αφού παίρνεις το ίδιο μάλλον κάτι κάνω λάθος....

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 14, 2017, 17:25:18 pm

Quote from: gdims on April 14, 2017, 17:20:33 pm

Μείωσα και το Linewidth για λεπτομέρια που θές

Παιδιά με διόρθωσε ο gdims , τα προηγούμενα διαγράμματα που ανέβασα δεν είναι ''σωστά'' , γιατί δεν ακολουθούσαν τον greedy αλγόριθμο. Δεν είχα ακολουθήσει αυτό που γράφει στο pdf
'Ο έλεγχος γίνεται αφού μετασχηματιστούν όλα τα pixels μίας συγκεκριμένης τιμής φωτεινότητας.'
για αυτό τα προηγούμενα διαγράμματα συγκλίνουν τέλεια στις εκάστοτε κατανομές.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Ροζ συννεφάκι on April 15, 2017, 21:58:19 pm

Για τις δυο πρώτες περιπτώσεις (ομοιόμορφες κατανομές) του 2.3 τα αποτελέσματα είναι ίδια;;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 15, 2017, 22:24:36 pm

Quote from: Ροζ συννεφάκι on April 15, 2017, 21:58:19 pm

Για τις δυο πρώτες περιπτώσεις (ομοιόμορφες κατανομές) του 2.3 τα αποτελέσματα είναι ίδια;;

Για ίδιο διάστημα ναι.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Ροζ συννεφάκι on April 16, 2017, 20:34:29 pm

Quote from: ReMi0s on April 15, 2017, 22:24:36 pm

Για ίδιο διάστημα ναι.

Τι εννοείς για ίδιο διάστημα; Αφού το διάστημα της κάθε κατανομής το δίνει και το διάστημα των τιμών της φωτεινότητας στο ιστόγραμμα είναι [0,1]. Έτσι δεν είναι;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 16, 2017, 21:24:31 pm

Quote from: Ροζ συννεφάκι on April 16, 2017, 20:34:29 pm

Εννοώ για ίδιο αριθμό υποδιαστημάτων. Και ναι θεωρώ πως έχεις δίκιο για το [0,1].
Anyway βγαίνει ίδιο αποτέλεσμα λόγο της κανονικοποίησης του h

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: lomptzoit on April 17, 2017, 23:53:47 pm

Για το 2.1 ο τυπος λεει να παρουμε οσα ειναι 0 να τα κανουμε v(1) μετα αμα ισχυει η συνθηκη τα 1 v(1) κοκ.
Πως θα το κανουμε αυτο εφοσον ο πινακας Χ που μας δινει εχει μονο τιμες στο διαστημα [0,1]?

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 18, 2017, 03:18:14 am

Quote from: lomptzoit on April 17, 2017, 23:53:47 pm

το έχεις καταλάβει λάθος.
Παίρνεις όσες τιμές είναι ίδιες και τις μετασχηματίζεις σε v(1), ξεκινόντας απο τις μικρές τιμές φωτεινότητας. Μετά ελέγχεις αν το h(1) έχει γεμίσει, αν δεν έχει γεμίσει συνεχίζεις με τις επόμενες τιμές φωτεινότητας και τις μετασχηματίζεις σε v(1) αλλιώς συνεχίζεις με το v(2) και ούτο καθεξής.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: lomptzoit on April 18, 2017, 23:48:01 pm

Quote from: ReMi0s on April 18, 2017, 03:18:14 am

Τι εννοεις ξεκινωντας απο τις μικρες τιμες φωτεινοτητας?

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 19, 2017, 02:53:42 am

Quote from: lomptzoit on April 18, 2017, 23:48:01 pm

Τι εννοεις ξεκινωντας απο τις μικρες τιμες φωτεινοτητας?

ξεκινώντας απο το 0 και καταλήγωντας στο 1.... τί να εννοώ;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: isitsou on April 20, 2017, 23:14:30 pm

Μπορεί κάποιος να μου εξηγήσει τι παίζει με το function pointer στο ερώτημα 2.2? Το γκουγκλαρα αλλά δεν έβγαλα και πολύ άκρη.Βασικά αυτό που με μπέρδεψε τελείως στην εκφώνηση είναι η φράση "επειδή στη γενική περίπτωση η f είναι μη μηδενική εκτός του διαστήματος [h(1); h(end)]".
Αν έχω καταλάβει σωστά στο ερώτημα 2.2 πρέπει να βρεις το πλήθος των πιξελ με φωτεινότητα v στο διάστημα [ d_k,d_k+1 ] και να το διαιρέσεις με τον συνολικό αριθμό των πιξελ της εικόνας, ώστε να βρεις την πυκνότητα της πιθανότητας στο εν λόγο διάστημα. right??? :-\

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 21, 2017, 00:53:46 am

Quote from: isitsou on April 20, 2017, 23:14:30 pm

function pointer εννοεί ότι μπορείς να καλέσεις την συνάρτηση pdf2hist ως εξης : pdf2hist(d,@myfunction) (h @myfunction για αρχή μπορεί να είναι @normpdf για τις δοκιμές, ωστόσο έτσι δεν μπορείς να βάλεις ορίσματα οπότε θέλει να κάνεις δική σου συνάρτηση που να παίρνει μέση τιμή και απόκλιση).
Απο εκεί και πέρα ολοκληρώνεις το κάθε διάστημα για να πάρεις την πιθανότητα. Ο τύπος λέει P(x1<x<x2)=ολοκλήρωμα της συνάρτησης απο x1 εώς x2. Για την ολοκλήρωση μπορείς να χρησιμοποιήσεις την συνάρτηση integral(f,d1,d2)

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: isitsou on April 21, 2017, 16:13:23 pm

Quote from: ReMi0s on April 21, 2017, 00:53:46 am

Όταν λες " ωστόσο έτσι δεν μπορείς να βάλεις ορίσματα " καταλαβαίνω ότι η κατανομή που θα έχουμε δεν θα είναι πάντα κανονική(ώστε να χρησιμοποιούμε συνέχεια την normpdf). Άρα, η συνάρτηση f θα πρέπει να επιλέγει αναλόγως σε ομοιόμορφη ή κανονική κατανομή και να χρησιμοποιεί την εκάστοτε συνάρτηση πυκνότητας πιθανότητας. Και μετά πάμε και κάνουμε ολοκλήρωση( integral(...) ) στο διάστημα [d_k,d_k+1] με ολοκληρωτέα συνάρτηση αυτή που δίνεται απο την f.Σωστός ή έχω καταλάβει αρλούμπες??? :D

edit: Ή μήπως πρέπει να κάνουμε διαφορετικές f για την κάθε κατανομή? δηλαδή πχ f_norm για κανονική κατανομή και f_uni για ομοιόμορφη κατανομή.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: abadasa13 on April 21, 2017, 19:04:33 pm

Quote from: isitsou on April 21, 2017, 16:13:23 pm

Εγώ πάντως αυτό που θα κάνω -δεν είμαι σίγουρος αν είναι βέβαια σωστό- είναι να φτιάξω μια συνάρτηση η οποία θα παίρνει τα απαραίτητα ορίσματα + ένα ακόμα που θα επιλέγει
αν θα χρησιμοποιείς normpdf ή unifpdf. Και μετά θα υπολογίζω το ολοκλήρωμα με integral σε αυτή τη συνάρτηση. Δεν είμαι σίγουρος για το πόσο σωστό ή υλοποιήσιμο είναι
αυτό που περιγράφω.

Αυτό που δεν κατάλαβα εγώ είναι το τι ρόλο παίζουν τα d. Εννοώ για ποιο λόγο να μην είναι π.χ. 100 ή 256 ή 1000 διαστήματα ίσου μεγέθους και στην κανονική κατανομή;
Πρακτικά, αν είναι απλά μεγάλο το μέγεθος του d όπως και να τα πάρεις τα διαστήματα, το ίδιο αποτέλεσμα θα πάρεις. Όχι;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Ragnar Lothbrok on April 21, 2017, 20:40:24 pm

Quote from: abadasa13 on April 21, 2017, 19:04:33 pm

Αυτό που δεν κατάλαβα εγώ είναι το τι ρόλο παίζουν τα d. Εννοώ για ποιο λόγο να μην είναι π.χ. 100 ή 256 ή 1000 διαστήματα ίσου μεγέθους και στην κανονική κατανομή;
Πρακτικά, αν είναι απλά μεγάλο το μέγεθος του d όπως και να τα πάρεις τα διαστήματα, το ίδιο αποτέλεσμα θα πάρεις. Όχι;

Όχι. Όσο περισσότερα τα διαστήματα, τόσο χειρότερο το τελικό αποτέλεσμα. Μην ξεχνάς ο αλγοριθμος είναι greedy. Για αυτό πρέπει να βρεις το ιδανικό αριθμό διαστήματος. Το τι εννοώ με τον όρο greedy, το αφήνω να στο εξηγήσει ο remios που εχει ταλέντο

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 21, 2017, 21:21:50 pm

Quote from: isitsou on April 21, 2017, 16:13:23 pm

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

function f=my_norm(x,m,s)
f=normpdf(x,m,s)
end

Η ίσως κατευθείαν
f=normpdf(x,0.5,0.1)

και αντίστοιχα για την ομοιόμορφη f=unifpdf(x,0,1) ή (x,0,2).

Ο άλλος τρόπος είναι με symbolic functions.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: abadasa13 on April 21, 2017, 23:23:08 pm

@Ragnar
Α το ξέχασα εντελώς ότι χρησιμοποιούμε μπακαλαλγόριθμο! :D
Επειδή δεν έγραψα τον κώδικα ακόμα, το σκεφτόμουν θεωρητικά.
Ναι οκ τώρα βγάζει νόημα η διαδικασία thanks!

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Kthulu on April 21, 2017, 23:31:28 pm

Quote from: ReMi0s on April 19, 2017, 02:53:42 am

ξεκινώντας απο το 0 και καταλήγωντας στο 1.... τί να εννοώ;

Θα παίρνουμε δλδ διαδοχικά τις τιμές 0/255, 1/255, 2/255 ... κοκ ;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Ragnar Lothbrok on April 22, 2017, 12:51:54 pm

Quote from: Kthulu on April 21, 2017, 23:31:28 pm

Θα παίρνουμε δλδ διαδοχικά τις τιμές 0/255, 1/255, 2/255 ... κοκ ;

Εξαρτάται τι εννοείς :P

Πρέπει να βρεις ενα τροπο να ξεκινας απο τα πιξελ με τη χαμηλοτερη φωτεινοτητα ( 0 ) και να φτανεις μεχρι τα πιξελ με την μεγιστη φωτεινοτητα ( 1 ) . Προφανως αν σαρωσεις την εικονα σου κελι κελι, δεν θα το πετυχεις αυτο, γιατι καθε πιξελ εχει διαφορετικη φωτεινοτητα και δεν ειναι σε διατεταγμενη αυξουσα σειρα.

Εγω αυτο που εκανα ηταν να φτιαξω εναν temporary πινακα, να τα βαλω ολα τα πιξελ σε αυξουσα σειρα, να τρεξω τον αλγοριθμο και μετα να τα επαναφερω στην αρχικη τους μορφη, με τις φωτεινοτητες πια αλλαγμενες.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: xristos_st on April 22, 2017, 17:20:03 pm

Στο ερωτημα 2.2 δεν εχω καταλαβει τι ειναι ακριβως το d και πως το χρησιμοποιουμε στο 2.3. Μπορει καποιος να εξηγησει?

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 22, 2017, 18:46:27 pm

Quote from: xristos_st on April 22, 2017, 17:20:03 pm

Το d είναι ένας πίνακας που περιέχει πχ τις τιμές [0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1]
Αυτός ο πίνακας καθορίζει τα διαστήματα σου. Δηλαδή απο αυτόν τον πίνακα έχεις το διάστημα [0, 0.1] ,[0.1, 0,2], [0.2, 0.3] .....
τα διαστήματα σου χρειάζονται για την ολοκλήρωση.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Indy on April 22, 2017, 19:01:18 pm

Το 1 το έκανα με τις κλίσεις των ευθειών...
Για το 1 β) όμως τί παραμέτρους πρέπει να περάσουμε; Δεν έχω πολυκαταλάβει... καμιά βοήθεια; Δοκίμασα 0.5, 0, 0.5, 1 αλλά μου βγήκε όλη μαύρη. "Ασπρόμαυρη" τί εννοεί; μόνο δύο τιμές;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 22, 2017, 19:42:38 pm

Quote from: Indy on April 22, 2017, 19:01:18 pm

αυτές οι τιμές που έβαλες είναι σωστές. Αν σου βγαίνει όλη μαύρη τότε κάτι έχεις κάνει λάθος.
Ασπρόμαυρη σημαίνει ασπρόμαυρη.... ο πίνακας σου θα έχει μόνο τιμές 0 και 1

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: isitsou on April 22, 2017, 19:56:58 pm

Quote from: ReMi0s on April 21, 2017, 21:21:50 pm

Θενξ ε λοτ για την απάντηση!!! ;)

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Indy on April 22, 2017, 21:46:59 pm

Quote from: ReMi0s on April 22, 2017, 19:42:38 pm

Δεν είχα μπορέσει να παρεβρεθώ στις αντίστοιχες διαλέξεις και δεν είμαι σίγουρος αν κατάλαβα σωστά από τις σημειώσεις και απ' ό,τι βρήκα online για σημειακό μετασχηματισμό, αλλά το έκανα έτσι:

Code:

function Y = pointtransform(X, x1, y1, x2, y2)
    [n, m] = size(X);
    Y = zeros(n, m);
    for i = 1:n
        for j = 1:m
            if X(i, j) < x1
                Y(i, j) = (y1 / x1) * X(i, j);
            elseif X(i, j) < x2
                Y(i, j) = ((y2 - y1) / (x2 - x1)) * X(i, j);
            else
                Y(i, j) = ((1 - y2) / (1 - x2)) * X(i, j);
            end
        end
    end
end

Ήθελε κάτι άλλο να φανταστώ;

Δοκίμασα και αυτό:

Code:

function Y = pointtransform(X, x1, y1, x2, y2)
    [n, m] = size(X);
    Y = zeros(n, m);
    for i = 1:n
        for j = 1:m
            if X(i, j) < x1
                Y(i, j) = (y1 / x1) * X(i, j);
            elseif X(i, j) < x2
                Y(i, j) = ((y2 - y1) / (x2 - x1)) * (X(i, j) - x1) + y1;
            else
                Y(i, j) = ((1 - y2) / (1 - x2)) * (X(i, j) - x2) + y2;
            end
        end
    end
end

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

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 23, 2017, 16:02:36 pm

Δεν κάθομαι να δω τους κώδικες που ανέβασες, ωστόσο η απάντηση στο πρόβλημα σου είναι η εξίσωση της ευθείας απο το γυμνάσιο.
y=λ(x-x1)+y1
λ=(y2-y1)/(x2-x1)

αυτό το προσαρμόζεις για τρεις περιπτώσεις. Για (0,0) με (x1,y1) , (x1,y1) με (x2,y2) και (x2,y2) με (1,1)

ο δεύτερος κώδικας είναι σωστός

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: thanospr on April 24, 2017, 11:41:10 am

Στο τελευταίο ερώτημα πως επιλέγουμε τα διαστήματα d?

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: isitsou on April 24, 2017, 17:36:34 pm

Quote from: thanospr on April 24, 2017, 11:41:10 am

Στο τελευταίο ερώτημα πως επιλέγουμε τα διαστήματα d?

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 24, 2017, 19:26:32 pm

θεωρώ πως το διάστημα θα είναι το [0,1] ωστόσο αλλάζει ο αριθμός των υποδιαστημάτων που θα επιλέξεις. Για τον αριθμό των υποδιαστημάτων λογικά κάνεις δοκιμές , ώστε να προσεγγίζει όσο καλύτερα γίνεται την εκάστοτε κατανομή. Τι εννοώ, πχ αν επιλέξεις 100 υποδιαστήματα τότε λόγο του greedy αλγορίθμου θα δείς ότι αποκλίνει πολύ απο την κατανομή. Ωστόσο αν επιλέξεις μικρότερο αριθμό, πχ 10 υποδιαστήματα τότε θα την προσεγγίζει καλύτερα.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: isitsou on April 24, 2017, 21:36:29 pm

Quote from: ReMi0s on April 24, 2017, 19:26:32 pm

Δηλαδή εσύ στον αλγόριθμό σου μεταβάλλεις μόνο το πλήθος των διαστημάτων ή και το μήκος του κάθε διαστήματος ώστε να "τρέχει" και για μεγάλο πλήθος υποδιαστημάτων? Το ρωτάω γιατί στην εκφώνηση λέει να "προδιαγράψουμε" τον τρόπο με τον οποίο θα επιλέγονται τα διαστήματα.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 25, 2017, 14:45:43 pm

Quote from: isitsou on April 24, 2017, 21:36:29 pm

Εγώ τουλάχιστον παίρνω ίσα υποδιαστήματα, οπότε αν αλλάξεις τον αριθμό αλλάζει και το μήκος τους... Δεν καταλαβαίνω ακριβώς την ερώτηση σου. Πάντως δεν έκανα αλγόριθμο να βρίσκει βέλτιστο διάστημα...
μια απλή εντολή d=linspace(0,1,10) ή d=[0:0.1:1] αρκεί...

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: talkshowhost on April 25, 2017, 18:41:24 pm

Στο 2.1 ξεκινάω από τα pixel της X με τη χαμηλότερη τιμή φωτεινότητας και την αλλάζω σε v(1). Αν ο αριθμός pixel που έχουν αντιστοιχηθεί στην v(1) είναι μικρότερος από τον επιθυμητό χρησιμοποιώ και τα pixel της X με τις επόμενες τιμες φωτεινοτητας. Αλλιώς συνεχίζω στις μεγαλυτέρες τιμές φωτεινότητας και τις αλλάζω σε v(2) κτλ
Σωστά?
'Ομως κάτι πάει λάθος στον κώδικα που έχω γράψει και δεν συνεχίζει στις v(2),v(3), οπότε προκύπτει ιστόγραμμα μίας τιμής..

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: abadasa13 on April 25, 2017, 18:43:58 pm

Quote from: isitsou on April 24, 2017, 21:36:29 pm

Κι εγώ στην αρχή αυτό κατάλαβα που λέει ο isitsou παραπάνω. Ότι δηλαδή αν για 1000 ίσα υποδιαστήματα δε βγαίνει σωστό αποτέλεσμα πάρε το d1-d2 < d2-d3 < d3-d4
μήπως και καταφέρουν και "μοιραστούν" οι τιμές καλύτερα. Δεν το δοκίμασα βέβαια οπότε δεν ξέρω αν λειτουργεί. Λέω να το αφήσω όπως έχει δηλαδή ίσα (μεταξύ τους)
διαστήματα. Εκτός αν αυτό ζητάει... Για άλλη μια φορά ψάχνουμε τα ζητούμενα...

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: abadasa13 on April 25, 2017, 18:47:34 pm

Quote from: talkshowhost on April 25, 2017, 18:41:24 pm

Σωστά! Υπενθυμίζω παίρνεις ΟΛΕΣ τις τιμές με φωτεινότητα π.χ. 0.1 και τις πετάς στο αντίστοιχο v(x). Και μετά κάνεις έλεγχο αν γέμισε το v(x).
Αν δε συνεχίζει στα επόμενα v πιθανότατα έχει ξεχάσει να χρησιμοποιήσεις μεταβλητή που κάνει αυτή τη δουλεία.
Δηλαδή να ελέγχει αν γέμισε το v(x) και αν έχει γεμίσει να γίνεται x = x+1.
Χωρίς να δω τον κώδικά σου δεν μου έρχεται κάποια άλλη ιδέα. :S

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: talkshowhost on April 25, 2017, 19:27:48 pm

Quote from: abadasa13 on April 25, 2017, 18:47:34 pm

Στη συγκριση ποσοστο pixel v >= επιθυμητο ποσοστο pixel v υπήρχε πρόβλημα, επειδή οι αριθμοί είναι πραγματικοί. Τελικά συγκρίνω αριθμούς pixel και μπαίνει κανονικά στην if για να αυξηθεί το x.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Kthulu on April 25, 2017, 20:58:59 pm

Στην histtransform, οσοι το κανατε με ταξινομημένο πίνακα, πως καταφέρατε να επιστρέψετε τον πίνακα με τα στοιχεία στις αρχικές του θέσεις;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Thtas on April 25, 2017, 22:10:47 pm

Quote from: Kthulu on April 25, 2017, 20:58:59 pm

υπαρχει ποιο απλος τροπος για να το κανεις
χρησιμοποιεις την συναρτηση unique(X) που σου δινει ολες τις διαφορετικες φωτεινοτητες ταξινομημενες και μετα εφαρμοζεις στην διαδικασια που λεει στο pdf μια και τωρα γνωριζεις καθε φορα την αμεσως επομενη φωτεινοτητα καθε φορα

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: talkshowhost on April 25, 2017, 22:14:05 pm

για το πρωτο ερωτημα του 2, αν χρησιμοποιησεις την sort του ματλαμπ, μπορεί να επιστρέψει vector που σε καθε θεση εχει τη θεση στο vector που ταξινομεις του αντιστοιχου στοιχειου του ταξινομημενου vector,
δηλαδη [xsorted,I]=sort(x) και x(Ι(i))=xsorted(i)
Δεν ξέρω τι άλλο μπορεί να χρειάζεται ο κώδικάς σου για να εχεις τον πινακα που θες..

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Kthulu on April 25, 2017, 22:19:54 pm

Quote from: talkshowhost on April 25, 2017, 22:14:05 pm

Ειδού το αποτέλεσμα..

Quote from: Thtas on April 25, 2017, 22:10:47 pm

Χμμ καλο ακουγεται. θα το προσπαθησω

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: pilatos on April 25, 2017, 22:54:06 pm

Quote from: ReMi0s on April 21, 2017, 21:21:50 pm

function f=my_norm(x,m,s)
f=normpdf(x,m,s)
end

Η ίσως κατευθείαν
f=normpdf(x,0.5,0.1)

και αντίστοιχα για την ομοιόμορφη f=unifpdf(x,0,1) ή (x,0,2).

To x που μπαινει τι ειναι?

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Kthulu on April 25, 2017, 22:58:07 pm

Quote from: pilatos on April 25, 2017, 22:54:06 pm

To x που μπαινει τι ειναι?

Ο πινακας της εικονας

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: pilatos on April 25, 2017, 23:10:04 pm

E και τοτε η εξοδος f θα ειναι ο πινακας αυτος μονο που αντι για τιμες φωτεινοτητας πλεον θα εχει την πυκνοτητα πιθανοτητας της καθε τιμης φωτεινοτητας. Τι θα ολοκληρωσεις μετα για να σου δωσει την πιθανοτητα της καθε τιμης φωτεινοτητας?

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 26, 2017, 10:01:10 am

Quote from: pilatos on April 25, 2017, 23:10:04 pm

Ολοκληρώνεις την f για να πάρεις την πιθανότητα ή αλλιώς ποσοστό των pixel στο τάδε διάστημα... δηλαδή το h σου.
Έπειτα με την histtransform παίρνεις τις τιμές φωτεινότητας.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: unfo on April 26, 2017, 11:58:38 am

Μπορει καποιος να καταλαβει που εχω λαθος εδω ή γενικα αν ειναι ολο λαθος ?

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: pilatos on April 26, 2017, 12:44:56 pm

Quote from: unfo on April 26, 2017, 11:58:38 am

Μπορει καποιος να καταλαβει που εχω λαθος εδω ή γενικα αν ειναι ολο λαθος ?

Δεν το εκανα με sorted οποτε μου ειναι λιγο δυσκολο να καταλαβω.
Απλα αυτο που παρατηρησα ειναι πως ο ελεγχος count/MN <h(k) πρεπει να γινεται αφου καταχωρουνται οι τιμες μιας σταθμης και οχι πριν συμφωνα παντα με την εκφωνηση.

@Remios Aυτο που δεν μπορω να καταλαβω ειναι οτι η εξοδος της f=normpdf(X,a,b) ειναι πινακας. Πως θα κανω μετα h=integral(f,d1,d2) αν η f δεν ειναι συναρτηση. Καπου το χανω.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: gdims on April 26, 2017, 13:34:56 pm

Παιδιά δεν υπάρχει κανένας λόγος και εγώ δε βλέπω κανένα νόημα στο να περάσεις την εικόνα x μέσα σε κάποια πυκνότητα πιθανότητας f. Η f είναι μια κλασική συνάρτηση, μιας μεταβλητής που ολοκληρώνεται ΚΑΤΑ ΔΙΑΣΤΗΜΑΤΑ....

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: sabpap on April 26, 2017, 14:54:38 pm

Quote from: unfo on April 26, 2017, 11:58:38 am

Μπορει καποιος να καταλαβει που εχω λαθος εδω ή γενικα αν ειναι ολο λαθος ?

Δεν είναι όλο λάθος,κανα-δυο σημεία μόνο!Τσέκαρε το συνημμένο αρχείο.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Indy on April 26, 2017, 15:12:52 pm

Στο 2.2 τί ακριβώς πρέπει να είναι το h που επιστρέφεται; 1xn vector με τιμές τις πιθανότητες για κάθε διάστημα; Ή 2xn vector με πρώτη γραμμή τις τιμές (τα d?) και δεύτερη τις πιθανότητες για κάθε διάστημα;
Επίσης, στην εκφώνηση λέει ότι πρέπει να περιγράψουμε τη μέθοδο αριθμητικής ολοκλήρωσης που υλοποιούμε, ενώ απ' τα λίγα που έχω καταλάβει (και δει εδώ μέσα), αρκεί(?) το h(i) = integral(f, d(i), d(i + 1)); οπότε την ολοκλήρωση την παίρνουμε έτοιμη απ το matlab. Ή θέλει πολλά περισσότερα; Ομολογώ πως μάλλον δεν το χω καταλάβει καλά το ερώτημα.

Επίσης, για το 2.3 πώς περνάμε στην pdf2hist την εικόνα X σαν όρισμα για την f? Πες ότι κάνουμε την

Code:

function Y = uni1(X)
    Y = unipdf(X, 0, 1);
end

Μετά καλώντας την pdf2hist ως Y = pdf2hist(d, @uni1) πώς περνιέται η εικόνα X σαν όρισμα;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: unfo on April 26, 2017, 18:42:50 pm

Quote from: sabpap on April 26, 2017, 14:54:38 pm

Δεν είναι όλο λάθος,κανα-δυο σημεία μόνο!Τσέκαρε το συνημέννο αρχείο.

Σε ευχαριστω πολυ... :D

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: abadasa13 on April 26, 2017, 19:04:13 pm

Quote from: Indy on April 26, 2017, 15:12:52 pm

Code:

function Y = uni1(X)
    Y = unipdf(X, 0, 1);
end

Μετά καλώντας την pdf2hist ως Y = pdf2hist(d, @uni1) πώς περνιέται η εικόνα X σαν όρισμα;

Για το πρώτο πράγμα που ρωτάς, το h πρέπει να είναι ενα 1xn vector με τιμές του ιστογράμματος για κάθε διάστημα. Το d περιγράφει απλά τα διαστήματα όπου αντιστοιχούν οι τιμές του h.
Αυτό για τη μέθοδο αριθμητικής ολοκλήρωσης ούτε εγώ το πολυέπιασα. Αυτό που λες κάνω πάντως. Υποθέτω ότι απλά θέλει να πεις με λόγια ότι "α αν ολοκληρώσω την f από εδώ ως εκεί βγαίνει το ιστόγραμμα". Κάτι παραπάνω δεν ξέρω - δεν κατάλαβα.
Για το 2.3, χρησιμοποιείς την pdf2hist για να δημιουργήσεις το ιστόγραμμα από τη ζητούμενη κατανομή. Και μετά η έξοδος h αυτής της συνάρτησης γίνεται είσοδος στην histtransform όπου εκεί βάζεις την εικόνα X.
Ελπίζω να βοήθησα!

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Indy on April 26, 2017, 21:32:12 pm

Quote from: abadasa13 on April 26, 2017, 19:04:13 pm

Όμως η pdf2hist θέλει ένα function pointer για να κληθεί, το f δηλαδή... και το f αν είναι πχ το @uni1 που όρισα παραπάνω, καλεί πχ τη unipdf... Η unipdf όμως παίρνει σαν πρώτο όρισμα το X... Ποιό X θα είναι αυτό (και πώς θα περνιέται στη unipdf, άρα πιο πριν στη uni1); Μήπως απλά η uni1 να loadάρει manually την εικόνα (lena.bmp) με το όνομά της;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: abadasa13 on April 26, 2017, 22:45:06 pm

Quote from: Indy on April 26, 2017, 21:32:12 pm

Y = unifpdf (x, 0, 1)

Όταν γράφεις αυτό το Y είναι η τιμή που έχει η ομοιόμορφη κατανομή από 0-1 στο x. Άρα το x είναι απλά ο άξονας που έχεις κάτω. Δηλαδή τα διαστήματα που ορίζουν τα d.
Και μετά κάνεις τη διαδικασία που σου περιέγραψα παραπάνω.

Από PDF πας σε Histogram με pdf2hist και μετά μπαίνει στο "κόλπο" η εικόνα σου. Αφού έχεις το επιθυμητό ιστόγραμμα και θες να πειράξεις την εικόνα για να προσεγγίζει αυτό.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Indy on April 26, 2017, 23:59:44 pm

Quote from: abadasa13 on April 26, 2017, 22:45:06 pm

αα όντως αυτό είναι το x της unipdf, σωστά σωστά... Τhanks!

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: unfo on April 27, 2017, 14:25:46 pm

Για την integral και τη unifpdf μου βγαζει το εξης προβλημα στο matlab :

Undefined function or method 'integral' for input arguments of type 'double'.
Undefined function or method 'unipdf' for input arguments of type 'double'.

Ξερει κανεις τι ειναι ?

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Gaara on April 27, 2017, 14:41:14 pm

~~Για την κατωφλιωση στο 0.5 τι κανω?~~
Σωστά κατάλαβα?

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 27, 2017, 16:05:29 pm

Quote from: unfo on April 27, 2017, 14:25:46 pm

Γιατί δεν υπάρχει η κατανομή unipdf... Είναι unifpdf .

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: unfo on April 27, 2017, 16:22:37 pm

Quote from: ReMi0s on April 27, 2017, 16:05:29 pm

Γιατί δεν υπάρχει η κατανομή unipdf... Είναι unifpdf .

Απο γρηγοραδα το εγραψα χωρις f... Ξερω οτι ειναι unifpdf .
Απ οτι καταλαβα παντως η απαντηση στο ερωτημα μου ειναι οτι λειπει καποιο toolbox !

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Indy on April 27, 2017, 17:42:21 pm

Quote from: abadasa13 on April 26, 2017, 22:45:06 pm

Και πώς περνάμε τα d στη unifpdf τότε; Αν κάνουμε Y = pdf2hist(d, @uni1) πρέπει κάπως να περαστεί και το d στη uni1 (που καλεί τη unifpdf)... Κάτι του στυλ Y = pdf2hist(d, @uni1(d)) ? Πώς θα γίνει κάτι τέτοιο;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: driverbulba on April 27, 2017, 18:25:33 pm

Ψάξτε για anonymous functions στο matlab

Έτσι γίνεται το πέρασμα παραμέτρων σε function handle

https://www.mathworks.com/help/matlab/matlab_prog/anonymous-functions.html (https://www.mathworks.com/help/matlab/matlab_prog/anonymous-functions.html)

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Indy on April 27, 2017, 18:28:56 pm

Quote from: driverbulba on April 27, 2017, 18:25:33 pm

thanks

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Kthulu on April 27, 2017, 19:09:55 pm

Quote from: abadasa13 on April 26, 2017, 19:04:13 pm

Και τι v θα παρουμε;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: abadasa13 on April 27, 2017, 19:58:53 pm

Quote from: Kthulu on April 27, 2017, 19:09:55 pm

Και τι v θα παρουμε;

Εγώ ολοκληρώνω από d(i) εως d(i+1). Τα v που λες τα χρησιμοποίησα στα 2.3 για να αντιστοιχίσω το μέσο του διαστήματος d(i) εως d(i+1) στην "μπάρα" του ιστογράμματος.

Δεν ξέρω αν είναι σωστό. Ας επιβεβαιώσει και κάποιος άλλος.

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 27, 2017, 20:32:38 pm

Quote from: abadasa13 on April 27, 2017, 19:58:53 pm

ρε παιδιά μη κολλάτε..
έχεις δημιουργήσει μια συνάρτηση έστω
function f=myuni1(x)
f=unifpdf(x,0,1)
end

έπειτα αντιστοιχείς το f στην συνάρτηση αυτή με f=@myuni1
και καλείς την pdf2hist(d,f).
δεν χρειάζεται πουθενα το x.
μετα απλά ολοκληρώνεις την f για όλα τα διαστήματα του d με την integral(f,d(i),d(i+1))

εσύ όταν ολοκληρώνεις μια συνάρτηση χρειάζεσαι την τιμή του x ή χρειάζεσαι τον τύπο της;

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Kthulu on April 27, 2017, 20:45:07 pm

Οκ βασικά το γράφει στην εκφώνηση, τα V πρέπει να είναι ο μέσος όρος του αντίστοιχου διαστήματος.

Μα είναι απαραίτητο να φτιαξουμε εξτρα συναρτήσεις; δεν αρκεί να ορίσουμε π.χ. f=@(x)unifpdf(x,0,1) ?

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: ReMi0s on April 27, 2017, 20:57:18 pm

Quote from: Kthulu on April 27, 2017, 20:45:07 pm

αμα λειτουργεί έτσι ακόμα καλύτερα :)

Title: Re: [Ψηφιακή Επεξεργασία Εικόνας] 2η Εργασία 2017
Post by: Kthulu on April 27, 2017, 21:09:20 pm

Quote from: ReMi0s on April 27, 2017, 20:57:18 pm

αμα λειτουργεί έτσι ακόμα καλύτερα :)

Λειτουργει ;)