THMMY.gr

Από e-thmmy

=================================================================================


=================================================================================

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

Δευτέρα - Τετάρτη - Παρασκευή 10-12

Ευχαριστώ για την προώθηση!

To avatar ειναι beanchmark?? :P

Όχι. Απλώς είναι μια παλιά ... αμαρτία. Alice Cooper rulez....
Χα! Καλό!

Θα το περάσουμε από έναν συμπιεστή και θα μας δώσει ... Poison!

Η τελευταία τιμή την οποία πέρνουμε για την εκτίμηση της εικόνας λέει οτι είναι η

Nαι έχει γίνει ένα μικρό λάθος, αυτό όμως είναι στο κείμενο και όχι στη σχέση.
Η εκτίμηση γίνεται με βάση τις 5 παραμέτρους α-e
από τα
x[m,n-1]
x[m-1,n]
x[m-1,n-1]
x[m-1,n+1]
x[m-2,n-2]

Το σφάλμα έγινε κατά το copy/paste και διέφυγε. Θα ανεβάσω διορθωτικό κείμενο αύριο το πρωί

Προσέξτε κάτι όμως όλοι. Αυτά που είναi bold
(και θα τα πώ και αύριο στο μάθημα)
Τα σχόλια παίζουν πολύ σημαντικό ρόλο, καθώς δείχνουν κατα πόσο κατανοήσατε τη αντίστοιχη θεωρία.

αυριο θα γίνει το μάθημα τελικά; Ο Καμένος μου έλεγε κάτι περί ημιαργίας..

Για το αν θα γίνει μάθημα αύριο, περιμένω κι εγώ νέα.
Αν δεν υπάρξουν, θα γίνει κανονικά. [ Ο οδηγός σπουδών δεν αναφέρει τίποτε για την 16η Νοέμβρη].

Θα μιλήσω με τον κύριο Στρίντζη και βλέπουμε.

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

Aκούς Καμένε, ακούς? :P :P
Πάει το πινγκ-πονγκ μας..!

Σήμερα, θα γίνει κανονικά μάθημα

Θανάσης Μαδεμλής

 Κανένας χριστιανός και μη που να βρήκε τις εντροπίες της αρχικής και της μετασχηματισμένης εικόνας?
Πόσο τις βρήκατε? Τι τιμές?

Εγώ βρήκα της αρχικής περίπου 7,56 και της μετασχηματισμένης 6,68.
(το περίπου πάει στο οτι δε θυμάμαι και τα 4 δεκαδικα ψηφία).

Κανενας άλλος προσφορά?? :)

7,5683

2,9347

Ουπςς...Εχουμε μια μικρή απόκλιση στη 2η...Γιατι έτσι?? :o

Η πρωτη βγαινει 7,56
Η δευτερη (του ερρορ) βγαινει 4,9...

Μετα το 3ο ερωτημα (και την καλυτερη κωδικοποιηση) οι εντροπιες πεφτουν στο 2,9 και 2,2 για την πρωτη και δευτερη συναρτηση αυτοσυσχετισης αντιστοιχα.

ΥΓ. Καντε ban τον user mademlis απ το τοπικ!! ;D

 Ερώτηση. Στο 2ο ερώτημα για την κωδικοποίηση Huffman χρησιμοποιώντας την εντολή dict στο matlab προκύπτουν αποτελέσματα της μορφής

[ 28]    [1x7  double]
[232]    [1x40 double]

Τι σημαίνουν αυτά τα αποτελέσματα. Πόσα bytes χρειάζονται για το κάθε σύμβολο? Τι σημαίνει double?

Ευχαριστώ!

Αν κατάλαβα καλά σημαίνει οτι πχ στο παράδειγμα που θέτεις για την φωτεινότητα 28 χρειάζονται 7 bits στην κωδικοποίηση.Για την 232 χρειάζονται 40.Το double να σου πώ δεν κατάλαβα τι είναι.

Το έκανε κανείς το 2ο ερώτημα???

Ολα τα νούμερα που βλέπω είναι λιγάκι off. (εκτός από την εντροπία της αρχικής)

Κάντε μια βόλτα απο το http://www.iti.gr/files/optimal%20construction%20of%20reduced.pdf
στον πίνακα της σελίδας 25, για να δείτε τι εντροπίες βγάζουν οι καλές μέθοδοι (μέχρι το 2000). Το λιγότερο 4.3....

Από κει και πέρα, τι άλλο να πω ?

Πιθανά σφάλματα:
Ιάσωνα, Αδαμάντιε, βεβαιωθείτε ότι υπολογίζεται πιθανότητες και για τα στοιχεία που έχουν αρνητικό σφάλμα!
Εβη (όπως θα δεις και στο email):
Μήπως υπολογίζεις το σφάλμα στην εικόνα που προκύπτει μετά την εκτίμηση ? (αν και πάλι υπάρχει πρόβλημα)....

Αυτά για αρχή....

Συγκεκριμένα το arithmetic???????????

Δεν σε πιάνω....

Καλά άστο το arithmetic .. ακόμα το ψάχνουμε.
Αλλά ρε συ το 4.9 καλό δεν είναι για το πρώτο ερώτημα???

όχι.....
Its too good to be true.....

^banghead^ ^banghead^

Αυτο με το αρνητικο σφαλμα μαλλον το παραλλειψαμε!!

Ίδρωσα (και σε έπρηξα και σένα κακόμοιρε δάσκαλε) μέχρι να βγάλω εντροπία 4,9 και λές ότι είναι too good to be true????Damn!!  ^banghead^

Ήταν ανάγκη να υπάρχει και αρνητικό σφάλμα??Why God?Why?

Εμάς μας ενδιαφέρει να λύσουμε την άσκηση έτσι ???? την επανάκτηση της εικόνας την έχουμε 'κάπου' γραμμένη....

το 5,1361 φίλε μαδεμλή πώς ακούγεται????

Έχει καταλάβει κανείς τίποτα απο το φοβερό παράδειγμα του Matlab για την αριθμητική κωδικοποίηση?? ^que^

με την εντολη arithenco θα βγαλεις την αριθμητικη κωδικοποιηση.
Στην ουσια θα ειναι ενας πινακας γραμμη με τα bits του τελικου αριθμου που υποτιθεται υπολογιζει η αριθμητικη κωδ.

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

Εξαρτάται.....
Για ποια δεδομένα μιλάμε ?

Αυτό ρωτούσαμε και εμείς παραπάνω αλλά έλα που δεν μας κατάλαβε ο chief.

Εμένα παντώς για κάποιο λόγο δεν τα δέχονταν αυτά τα ορίσματα....

Τα ορίσματα της arithenco ειναι seq και counts(όπως τουλ τα έχει το παράδειγμα).Ας υποθέσουμε οτι το counts αντιπροσωπεύει τις πιθανότητες συμβόλων.Στο seq όμως τι βάζουμε??

Αν μιλάμε για την εντολή του Matlab, τότε τα πράγματα είναι απλά
Θέλει 2 πίνακες σαν ορίσματα, τον Α και τον Β

στον Α υπάρχει η συμβολοσειρά προς κωδικοποίηση, σαν στήλη και το Β της πιθανότητες του κάθε συμβόλου.

Αν τα σύμβολα είναι s1, s2, s3, s4....s122...
τότε το Β περιέχει
Β = [ p(s1) p(s2)  .... p(s122) ...];

αν η συμβολοσειρά προς κωδικοποιίηση είναι η
s11,s15,s11,s12,s113, .....
Τότε ο Α περιέχει
Α = [11 15 11 12 113 ....]

Προσοχή στα σύμβολα μηδενικής πιθανότητας. Πρέπει να αφαιρεθούν!
και περιέχει

Α πρεπει αφου κανουμε την εικονα μας συμβολοσειρα να ανακαταταξουμε και τις πιθανοτητες αντιστοιχα??
Ελεος!! :P

Τελικα τιποτα δεν εχουμε σωστο και νομιζαμε οτι μας εμενε μονο η αναφορα!

ΥΓ.Αν αυτη ειναι η ευκολη 1η εργασια ειμαι περιεργος να δω την 2η! :P

Μη συζητάς.....
Να σου πω,έτσι όπως το βλέπω εγώ, η 2η θα είναι κωδικοποίηση βίντεο και η 3η κωδικοποιήση βιντεο με ήχο!!!Και νόμιζα ότι τα πολυμέσα έχουν δύσκολη εργασία....Συναρπαστικό εξάμηνο,δεν μπορώ να πω! ^kremala^

Ιάσωνα , περίπου
Βασικά, θα πρέπει να γίνει μια κατάλληλη αντιστοίχιση, ώστε να μην υπάρχουν προβλήματα με σύμβολα μηδενικής πιθανότητας

Παιδιά, μην τρελλαθούμε κιόλας....
Κολλάτε σε 2-3 προγραμματιστικά κολπάκια και πελαγώνετε....
Δεν φαντάζομαι τι θα γινόταν αν το γράφατε σε C.....

Θεωρώ ότι μάλλον έχετε κουραστεί, οπότε καλύτερα να κοιμηθήτε και τα λέμε αύριο στις 4.... (περίπου)

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

[Προλαβαίνω:

Αυτό σημαίνει είτε ότι:
ο κόπος που θα ρίξεις στο εξάμηνο θα είναι πολύ μεγαλύτερος από τον κόπο μιας απλής εξεταστικής
είτε ότι
Ο κόπος που θα ρίξεις στο εξάμηνο, θα σε απαλλάξει από πολύ κόπο στην περίοδο της εξεταστικής.

Εδώ συνήθως συμβαίνει το 2ο.....
]

Ναι δεν λεω οτι εχουμε να λυσουμε τις τρελες εξισωσεις ή το τρελο προβλημα.
Απλα δεν υπαρχει προφανης τροπος (ως συνηθως βεβαια) να δουμε αν αυτο που βγαζουμε ειναι σωστο και η αληθεια ειναι οτι η εργασια κρυβει παγιδες...δεν εννοω στις εκφωνησεις αλλα οπως γινεται εδω με τους περισσοτερους που κατι δεν εχουμε καταλαβει καλα ή δεν υπολογισαμε.

Οσο για το εξαμηνο δεν εχω παρα να συμφωνησω...τα ιδια τραβηξαμε και περυσι στο 7ο αλλα εγω προσωπικα περασα 7 μαθηματα (αριθμος ρεκορ για μια εξεταστικη) !

Θανάση,έχω ερώτηση που σπάει κόκκαλα!!(ετοιμάσου να σου πάρω τη θέση στο ΙΠΤΗΛ...)

κοίτα,κάνω την αριθμητική και το θέμα είναι ότι μου βγάζει ένα λόγο συμπίεσης περίπου 1,8:1 ενώ ο Huffman πέζει πολύ πιο πάνω( θυμάμαι 17,κάτι:1).Βρήκα ότι η αριθμητική είναι πολύ αποτελεσματική κυρίως σε περιπτώσεις που έχω μικρό αλφάβητο συμβόλων(εδώ προφανώς δεν έχω) αλλά τόσο μεγάλη διαφορά στη συμπίεση μου έκανε εντύπωση.

Οπότε,επί του πρακτέου,ή είναι σωστό κ πάμε για το τρίτο ή έχει γίνει μανούρα με τους πίνακες της αριθμητικής και θα αρχίσω να δέρνω τον υπολογιστή..(προφανώς κ φταίει αυτός κ όχι εγώ,έτσι!??)

Κ ένα τελευταίο,στο τρίτο,το σύστημα στο Matlab το λύνω ή δίνω στο πρόγραμμα απλά εισόδους τις βέλτιστες τιμές των συντελεστών και καθαρίζω (το κατά πόσο καθαρίζω,άλλο θέμα)?Να το κάνω να δουλεύει στη γενικότητα?Προς το παρόν έχω κάνει 5-6 συναρτήσεις,να τις ενοποιήσω όλες σε μία κ να τρέχει ο κώδικας από ένα .m file??

Τι να σε κάνω,έμπλεξες με ηλεκτρολόγους και οι ερωτήσεις θα σκάνε από όλες τις μπάντες...

Thx

Παιδιά το παλεύει κανείς άλλος σε java ή είμαι ο μόνος και να το γυρίσω καλύτερα??
Το πρώτο μέρος όλα καλά πήγαν αλλά μετά για huffman και αριθμητική λίγο που είδα δε βρήκα καμιά σοβαρή υλοποίηση...

LOL nub! :S

ρε συ black ,παίξε Matlab γιατί αλλιώς θα μπλέξεις..γιατί να μη χρησιμοποιήσεις τις embedded κ να τελειώνει η ιστορία???

Mxiro,άντε αγόρι μου κάνε κάνα headhot να γουστάρουμε όλοι μαζί σου..;)

Πάμε ένα ένα και σιγά σιγά.

Τα κύρια μέρη, θα πάνε προς τον αγαπητό Παναγιώτη.....

Ξεκινάμε:
1) Πρώτα από όλα, από τους καινούριους στο ΙΠΤΗΛ, κάποιοι και κάποιες πέρυσι ήταν στη θέση σας, σωστά ?
2) Δεν πρέπει να αποθαρρύνουμε κανέναν να ασχοληθεί με τη γλώσσα που ξέρει και γουστάρει.
3) Αν κατάφερες να βγάλεις μέσο μήκος λέξης με την αριθμητική κωδικοποίηση, ή την κωδικοποίηση Huffman μικρότερο απο την εντροπία, τότε, ή θα παρουμε το Νομπελ ή θα βρούμε ένα bug ΤΕΡΑΣΤΙΟ.....

Γιώργο,
 Αν θέλεις, δώσε μου μερικές ώρες για να σου πω περι Java.

Απλα μπορεις αντι να τα υπολογισεις στο χερι να βαλεις το συστημα στο ματλαμπ και να σου κανει τη δουλεια.
Ειναι απλα ενας πινακας 5χ5 .

Ναι, απλά ένας πίνακας 5x5 που λύνεται σε λίγα msec από το Matlab  ;)

Τι θα γίνει με το "Too good AND true" της εντροπίας 4.95 της μετασχηματισμένης εικόνας του 1ου ερωτήματος;
Θανάση μας είχες πει ότι θα είναι κοντά στο 5. Το 4.95 είναι ΠΟΛΥ κοντά στο 5...  :(

Οι φωτεινότητες των pixels είναι από 0 έως 255, έτσι δεν είναι; Ρωτάω γιατί δεν έχω αντιμετωπίσει ακόμα το πρόβλημα που μας έλεγες για το 256 σε int8... Ίσως επειδή δεν έκανα ακόμα Huffman και Αριθμητική...  ::)

Τι 5 ;!;!

Εγώ βγάζω εντροπία 7,5683 για την αρχική εικόνα, 6,8149 για την γραμμ. εκτίμηση και 6,8522
, 6,8207 μέσο μήκος λέξης για τις δυο κωδικοποιήσεις... βαθιά νυχτωμένος;;;

Κάποιος άλλος έβγαλε συμπίεση με τον Huffman 17:1 !!! ...το jpeg βγάζει εύκολα 20:1, αξίζει εντέλει να το χρησιμοποιούμε;

Και μια ερώτηση που δεν είχα σκεφτεί πιο πριν:
Μόλις παρατήρησα ότι κατά τον υπολογισμό της γραμμικής εκτίμησης της αρχικής εικόνας τα αποτελέσματα βγαίνουν από -22 έως 258... τι ακριβώς τιμή θα βάζουμε για τα pixel που η τιμή τους ξεφεύγει από το (0-255);
Όσα έχουν τιμή <0 τα βάζουμε 0 και όσαν έχουν >255 τα βάζουμε 255;;;
Ή μήπως κάνω κάποιο λάθος;

Προφανώς επειδή οι φωτεινότητες βγαίνουν από -22 έως 258 βγαίνει λάθος και η εντροπία...

οι φωτεινότητες των pixels είναι από -255 μέχρι 255 (δεν εμφανίζονται όλες) γιατί στέλνεις τ δφάλμα,όχι την αρχική εικόνα

Hint: αποθήκευσε τον error-πίνακά σου ως double και απλά πόλωσε τον κατά 1. Η συνέχεια με την εντροπία είναι απλά ...fantastico.

Θανάση, το αρχειάκι που έχεις ανεβάσει στο eTHMMY myrearrange.m έχει ένα λάθος( ; ) προς το τέλος. Γράφεις " for i=1:sx " αλλά, τι είναι το sx ; Εν τέλει, στο τελευταίο for τι έπρεπε να κάνει;

Το sx πρέπει να γίνει len

Την είπαμε τη σωστή εντροπία, είναι κοντα στο 5.6 για τα δεδομενά a-e


Με Huffman 17:1 ?
Αυτό σημαίνει ή ότι πάμε για Νόμπελ, ή μιλάμε για ΤΕΡΑΣΤΑ πατάτα.
Πόσοι ποντάρετε στο πρώτο ?

Το λάθος που κάνεις Δημήτρη, είναι στο εξής: Από τη στιγμή που ξέρω τον τρόπο πρόβλεψης, μεταδίδω το σφάλμα εκτίμησης. Έτσι, εντροπία υπολογίζουμε στο σφάλμα πρόβλεψης. Για αυτό "τα μαύρα μεσάνυχτα"

Έξυπνο, αλλά δεν δουλεύει πάντα :-)


Αυτό είναι ένα θέμα. Αλλά είναι καθαρά θέμα σύμβασης.
Δείτε πόσο και πώς επηρρεάζει το τελικό αποτέλεσμα.

Από τη στιγμή που μιλάμε για σφάλμα, δεν είναι σωστό να μιλάμε για φωτεινότητα, αλλά για διαφορά φωτεινοτήτων, ή σφάλμα φωτεινοτήτων κλπ

Δηλαδ για τιμές <0 ότι το pixel έχει φωτεινότητα 0 και για τιμές >255 ότι έχει 255...;;
Ένας νομπελίστας τι θα έκανε;;  :D

Tην εντροπία της μετασχηματισμένης εικόνας τη βρήκαμε 5,6116...

Δάσκαλε, την πιάσαμε την ιδέα της μετασχηματισμένης?

Πολύ καλή ερώτηση....
Σκέψου ότι αν για τον xyz λογο, κατα τη μετάδοση ενός πίνακα σφάλματος, τα τελευταία σφάλματα, δεν φτάσουν, και εσύ πρέπει να δείξεις κάτι  άμεσα, πως θα δείξεις φωτεινότητα -2 ή φωτεινότητα 260 ?
Πρακτικά, τα πάντα είναι θέμα σύμβασης:
Ακόμα κι αυτή η μεταξύ μας επικοινωνία. Έχουμε συμφωνήσει, (σύμβαση) ότι η λέξη ΝΑΙ σημαίνει θετική καταφατική απάντηση, και όχι κάτι άλλο. Έχουμε συμφωνήσει, ότι την ένταση της φωτεινότητας θα την δείχνουμε με έναν ακέραιο [0...255] ή [0...65535]. Γιατί να μην συμφωνήσουμε να πειράζουμε την εκτίμηση σε περίπτωση που δώσει μη έγκυρη τιμή φωτεινότητας, προς την πλησιέστερη έγκυρη τιμή ?

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

Κανονικά, βέβαια, οι τιμές αυτές όμως πρέπει να είναι ελάχιστες. Εφαρμόζοντας τώρα το θέμα, με τα δεδομένα a-e και αυτά που υπολογίζονται με την πρώτη εντροπία, στην πρώτη περίπτωση μιλάμε για 30 τιμές και τη δεύτερη για 300 σε σύνολο 65000pixel, συνεπώς το αποτέλεσμα στην εντροπία να μην είναι θεαματικό.

Αυτό πάντως, το αυθαίρετα, δεν ισχύει σε κάθε πρότυπο αποθήκευσης ?

Το βρήκα το λάθος μου.... θέλει λιγάκι προσοχή με τις μετατροπές uin8 σε double και τελικά οι entropy και H=imhist() που τόσο είχα χαρεί όταν τις βρήκα στο matlab είναι εντελώς άχρηστες στο πρόβλημά μας. Θέλει να γίνει "χειροκίνητα" η δουλειά...

Flawless....

Στην 3η περίπτωση λύνοντας το σύστημα οι συντελεστές βγαίνουν
α1=0.3126 α2=0.5143  α3=-0.0094 α4=0.2418 α5=-0.0692
οπότε και η εντροπία του σφάλματος σε αυτή την περίπτωση βγαίνει 4.99 . (Η εντροπία στην 1η περίπτωση βγαίνει το σωστό 5.611 άρα δεν υπάρχει λάθος στον υπολογισμό της)
Ωστόσο όμως το 4.99 απέχει ... από το 4.3 που δεν ξέρω πώς αλλά μου καρφώθηκε ως το σωστό που πρέπει να βγει!
Κάποιο αξιόπιστο άτομο ας με διαφωτίσει...(hint-hint: mademlis)!

4,99 βγαζεις για την 2η συναρτηση αυτοσυσχετισης? (οι οροι δλδ προεκυψαν απ αυτην?)

αν ναι τοτε εχουμε ελπιδες να το εχουμε σωστο!

Εδιτ: Τωρα βλεπω και τους ορους...ιδιους τους εχουμε!

 ^wav^

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

Παιδια σας βλεπω ολους να εχετε ασχοληθει σοβαρα με την εργασια και νιωθω...καπως που δεν εχω κανει τιποτα!
Ειπα να κατεβασω λοιπον το αρχειο για Huffman και Αριθμητικη μιας και εχω το Matlab 6.5.
Ξερει κανεις πως το εγκαθιστω στο Matlab?

Αν ειναι απλα ενα .m file δεν εχει εγκατασταση...
ειναι απλα μια συναρτηση που την χρησιμοποιεις.
Αρκει να την εχεις στο Current Directory.

Αν δεν την εχεις εκει ή βαλτην σε εναν φακελο με αλλες συναρτησεις...ή βαλτην οπου θες και πηγαινε στο File->Set Path και προσεθεσε το path του φακελου που εχεις βαλει την συναρτηση!

thnx ;D
 done

Σου καρφώθηκε το 4.35 ε ?
Κακώς......
Το 4.35 είνα η εντροπία που δίνει μια μέθοδος από το state-of-the-art, μεταξύ των καλύτερων μεθόδων (S+P)
Το είχα πει κάποια στιγμή, γιατί ακούστηκαν νούμερα όπως 0.29, 1.5 κ.α. :-)

Κατά τα άλλα, αυτην τη στιγμή, δεν βλέπω προφανές λάθος.

Υπόψιν, ότι όλα αυτά χρειάζονται καλό σχολιασμό.
πχ Γιατί το τρίτο μοντέλο εκτίμησης δίνει καλύτερα αποτελέσματα ?
αλλά και σε συνδυασμό με Huffman και αριθμητική.

Απλώς να ενημερώσω σε περίπτωση που υπάρχει και άλλος. Το πρώτο και το τρίτο κομμάτι της άσκησης βγαίνει σχετικά εύκολα σε java. (Μάλιστα είναι ευκολότερο να βγάλεις την εντροπία χειροκίνητα, γιατί σε διαφορετικη περίπτωση πρέπει να μπλέξεις με Java Advanced Imaging(JAI)).
Για Huffman και την αριθμητική υπάρχει τελικά μία λύση που λέγεται Java Native Interface (JNI), δηλαδή να φορτώσεις από  java πρόγραμμα αρχεία πχ από c. Βέβαια σε αυτή τη περίπτωση θα πει κάποιος ας το κανες από την αρχή σε c, πράγμα που ήθελα να αποφύγω!
Θανάση, αν έχεις καμία καλύτερη ιδέα, είναι ευπρόσδεκτη!

Σχετικά με το δεύτερο ερώτημα στο οποίο ζητάει δύο πράγματα:για τον μέγιστο απαιτούμενο αριθμό bytes χρειάζεται να υπολογίσουμε το μέσο μήκος από την κωδικοποιήση Huffman για να το συγκρίνουμε με το μέγιστο αριθμό που θεωρητικά είναι H+1;Επίσης για τον πραγματικό λόγο συμπίεσης θα πρέπει να βρούμε τα συνολικό αριθμό των bits που προκύπτουν από την κωδικοποίηση και να τα διαιρέσουμε με τον αριθμό των bits της αρχικής εικόνας (256*256*8 );
Γενικά τι αποτελέσματα βγάλατε σ'αυτό το ερώτημα;Για τα δεδομένα a,b,c,d,e:
Το μέσο μήκος μου βγαίνει 4,6489 και ο συνολικός αριθμός bits 4999

Και κάτι ακόμα... η συνάρτησης arithenco τι έξοδο μας δίνει;Γιατί από το help δεν πολυκατάλαβα :P

Πάμε από την αρχή.
Aπό τη στιγμή που ζητάει να γίνουν κάποιοι υπολογισμοί σε bytes, καλύτερα όλα να είναι σε bytes. Εξάλλου, αυτή είναι η μονάδα αποθήκευσης.

Τα νούμερα που δίνεις είναι λιγάκι, "ότι να ναι!" και δεν συμβαδίζουν μεταξύ τους

Η arithenco δίνει έξοδο το bitstream προς αποθήκευση.

Για την 1η εκτίμηση με τα δεδομένα a,b,c,d,e:
Εμένα η εντροπία της μετασχηματισμενης εικόνας μου βγαίνει 6,19
και δεν βλεπω καποιο λαθος στον κωδικα μου.Τον κοιταω εδω και 2 μερες και δεν βρισκω καποιο λαθος ετσι ωστε να βγάλω 5,61.
Ασε που αθροίζω και τις πιθανοτητες και μου βγαίνει 1.
1)Σωστο η λάθος?
επισης είδα οτι οι 'διαφορες φωτεινότητας'(σωστα?)ξεκινανε απο το -190 και φτάνουν στο 196.
2)που μπορει να ειναι το λαθος?

1) Δες ξανά τη σχέση εκτίμησης, γενικότερα,
2) Πρόσεξε τις τιμές που παίρνει η συνάρτηση εκτίμησης, για 4-5 σημεία, και επαλήθευσέ τα με το χέρι. (και φυσικά το σφάλμα εκτίμησης.)
Γενικά δες βήμα βήμα τι σου βγάζει, και πόσο σωστό είναι.
Δόξα το θεώ, και το Matlab Και η MSVC έχουν πολύ καλούς debuggers, για τέτοιες παρακολουθήσεις.

Ρωτήσαμε τον θανάση, και το "Να υπολογίσετε τον μέγιστο απαιτούμενο αριθμό bytes για την μετάδοση της παραπάνω εικόνας"(στο ερώτημα 2) αναφέρεται ΜΟΝΟ στην κωδικοίηση Huffman και όχι στην αριθμιτική κωδικοποίηση...

Aν όμως κάνουμε thresholding στις τιμές του πίνακα της εκτίμησης τότε η κωδικοποίηση δεν θα είναι αντιστρέψιμη δηλαδή δεν θα ισχύει x(m,n)=^x-e για κάθε   m,n.
Έ Θανάση;

deleted

Αυτό είναι το σωστό!


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

Ορίστε? :???: και γιατί αυτό?

Σίγουρα καταλάβατε σωστά? Γιατί, αν βρούμε και τα bytes της αριθμητικής, τι πειράζει? και κυριότερον, πώς θα συγκρίνουμε τις δύο κωδικοποιήσεις? αφού στο κάτω κάτω για να βρούμε τον λόγο συμπίεσης, δεν χρειαζόμαστε και τον απαιτούμενο αριθμό bytes?

Αυτά ακριβώς είναι και η δική μου απορία. Αν δεν βρούμε τα bytes της αριθμητικής τότε τι θα βρούμε με την εφαρμογή της μεθόδου; Γιατί δε μένει και τίποτα άλλο!

Το ΜΕΓΙΣΤΟ απαιτούμενο είπα...Όχι τα bytes της αριθμιτικής γενικά...

Μα και στον Huffman ένα νούμερο δε βγαίνει;; Τι θα πει μέγιστο απαιτούμενο;

deleted

Πρέπει να το εφαρμόσεις για την μετασχηματισμένη, αλλιώς δεν έχει νόημα άλλωστε. Επίσης να προσέξεις με τις αρνητικές τιμές του πίνακα error.

deleted

Δεν έχει νόημα αν διαβάσεις το 3ο ερώτημα ;)


Θα έλεγα πολύ καλύτερα ;)

Worst case analysis..

Μεταξύ ποιων τιμών κυμαίνεται το μέσο μήκος λέξης για Huffman..? Για την αριθμητική είναι fix, οπότε περισσεύει.

Χμ, κάτι σαν να πιάνω από αυτά :-\...
Εννοείς δηλαδή να πάρω την περίπτωση για την οποία χρησιμοποιείται το μέγιστο μήκος λέξης (αυτό που αντιστοιχεί στη λιγότερο πιθανή διαφορά φωτεινότητας) σε κάθε pixel όλης της εικόνας;

Μα η εικόνα μου είναι δεδομένη, οπότε δε χρειάζεται καμιά πρόβλεψη, ξέρω πόσες φορές θα εμφανιστεί κάθε διαφορά, ξέρω και το μήκος λέξης της κάθε διαφοράς, οπότε με το άθροισμα των γινομένων του πλήθους των διαφορών επί το μήκος λέξης της κάθε διαφοράς θα βγει ένα στάνταρ νούμερο!

Ή κάνω λάθος;;;

Όχι, εννοώ να πάρεις το θεωρητικό μέγιστο L για την κωδικοποίηση Huffman... για τη συγκεκριμένη εικόνα.

Με την αριθμητική κωδικοποίηση και το μέσο μήκος λέξης της εικόνας τι γίνεται?
Δηλαδή ξέρουμε θεωρητικά ότι η αριθμητική κωδικοποίηση πετυχαίνει μέσο μήκος ίσο με την εντροπία....Αλλά αυτό το μέσο μήκος πρέπει κάπως να το υπολογίσουμε για να εξακριβώσουμε την ισότητα??
Και αν ναι, πως θα γίνει αυτό με το matlab αφού το μόνο αποτέλεσμα που έχουμε από την αριθμητική κωδικοποίηση είναι μία γραμμή με τα bit που αντιστοιχούν διαδοχικά στα σύμβολα...Πως ξέρουμε πόσα bits αντιστοιχούν σε κάθε σύμβολο για να μπορέσουμε να βρούμε το μέσο μήκος L?? :( :(

Αχ.....

Πρόκειται για  την απλή μέθοδο των τριών

Για να κωδικοποιήσεις αυτά τα Ν σύμβολα, χρειάστηκαν Μ bits....
Kατα μέσο για 1 σύμβολο θέλω.... ?

@marauber: Την ανισοτητα που γραψαμε θελει στο Huffman βρε !

@BoB: Και γω το ιδιο πιστευα στην αρχη (οτι στελνουμε την αρχικη εικονα και οτι το 1ο ερωτημα ειναι ανεξ. του 2ου). Λαθος! Στελνεις το σφαλμα. Και υποτιθεται ο δεκτης γνωριζοντας την συναρτηση εκτιμησης και το σφαλμα ξαναφτιαχνει την εικονα. Η συμπιεση τοτε ειναι καλυτερη.

και μενα μου βγαινει 5,6487 :) αλλα με εντροπια 5,6157 και αυτη η διαφορα πρεπει να ειναι γιατι χρησιμοποιησα ακεραιοποιηση και οχι round για τη στρογγυλοποιηση.

Απ'ότι μου είπε ο Θανάσης, ανάλογα ποια μέθοδο χρησιμοποιούμε για στρογγυλοποίηση βγαίνει ελάχιστα διαφορετικό αποτέλεσμα στην εντροπία που συζητάτε..

Χρησιμοποιήσαμε αρχικά την int16 και μετά την round, και νομίζω βρήκαμε το ίδιο αποτέλεσμα, 5,6113

Κάποιοι πήραν την fix (???) και έβγαλαν ελαφρώς διαφορετικό..

Προσοχή ότι ακαιρεοποίηση με το int8 δεν είναι σωστό..

εγω την κανω σε c++ και χρησιμοποιω εκει την (int)(x^[m,n]).
γιατι ειναι λάθος με την int8 του matlab?
ισχυει το ιδιο και για την int της c++?

H int8 του Matlab είναι ισοδύναμη με τo char της C++
Γενικά οι αντιστοιχίες είναι

Matlab --  C
int8 --> char
int16 --> short int
int32 --> long int


Είναι λάθος η χρήση της Int8 γιατί πολύ απλά, δεν επαρκεί για να περιγράψει τα λάθη (χρειάζονται τουλάχιστον 9 Bitγια κώδικα σταθερού μήκους.) Εσύ χρησιμοποιείς την int (Πιιθανότατα 32βιτ) Που είναι υπεραρκετή για την δουλειά αυτή.

:P

Εdit

το παρακάτω μνμ γράστηκε πριν την απάντηση του mademlis, τo στέλνω anyway..

=================================================================

Δες τι εύρος τιμών μπορείς να αποθηκεύσεις σε απλό int..

Απ'ότι θυμάμαι οι τιμές που θέλαμε να αποθηκεύσουμε δεν χωρούσαν όλες σε 8 bits του int8..

Οπότε και πήραμε την int16, ώστε να μπορούμε να χρησιμοποιήσουμε μεγαλύτερους αριθμούς..
(το είπε και ο Θανασης στο μάθημα, ότι int8 λάθος, ιnt 16 ίσως σωστό, round πολύ πιο σωστό)


Απ'ότι βλέπω, στη c++ ο int θέλει 2/4 bytes (ανάλογα με το μηχάνημα), άρα 16/32 bits, οπότε μάλλον είσαι οκ..

thx επιστρέφω με αποτελέσματα για τις αλλες 2 αυτοσυσχετισεις για επαληθευση!

λοιπον:
1)για την 1η αυτοσυσχέτιση:
a=0.95
b=0.95
c=-0.9025
d=e=0
H=5.14798
L_huffman=5.18163

2)για τη 2η:
a=0.31256494410545
b=0.51428707915223
c=-0.00938181267524
d=0.24184708046186
e=-0.06920833725009
H=4.98906
L_huffman=5.02556

anyone με ιδια ή οχι αποτελεσματα?
πως σας ακουγονται?

Σχεδόν τα ίδια με απειροελάχιστες αποκλίσεις στα bold σημεία. Αλλά μιλάμε για πολύ πολύ μικρές αποκλίσεις...