[Ψηφιακή Επεξεργασία Εικόνας] 1η Εργασία 2017

[Thmmy.gr portal]

  Thmmy.gr portal
   Forum
   Downloads
   Ενεργ. Λογαριασμού
   Επικοινωνία
  
  Χρήσιμα links
   Σελίδα τμήματος
   Βιβλιοθήκη Τμήματος
   Elearning
   Φοιτητικά fora
   Πρόγραμμα Λέσχης
   Πρακτική Άσκηση
   Ηλεκτρονική Εξυπηρέτηση Φοιτητών
   Διανομή Συγγραμμάτων
   Ψηφιακό Καταθετήριο Διπλωματικών
   Πληροφορίες Καθηγητών
   Instagram @thmmy.gr
   mTHMMY
  
  Φοιτητικές Ομάδες
   ACM
   Aristurtle
   ART
   ASAT
   BEAM
   BEST Thessaloniki
   EESTEC LC Thessaloniki
   EΜΒ Auth
   IAESTE Thessaloniki
   IEEE φοιτητικό παράρτημα ΑΠΘ
   SpaceDot
   VROOM
   Panther
  

THMMY.gr Web

Helios_MultiScribbles2SMF Default ThemeSMFone_Blue

[Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads]

Είσαι πρωτοετής;... Καλώς ήρθες! Μπορείς να βρεις πληροφορίες εδώ. Βοήθεια για τους καινούργιους μέσω χάρτη.

Κατεβάστε εδώ το Android Application για εύκολη πρόσβαση στο forum.

Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads με προσοχή στα ονόματα των αρχείων!

[Δείτε περισσότερα εδώ...]

Συμβουλές καλής χρήσης του φόρουμ: Youtube embed code and links, Shoutbox, Notify, ...

Δείτε περισσότερα εδώ...

[Επικο Burger]

Νομίζω κάτι δεν κατάλαβες καλά!

Ο νεκρής ζώνης είναι συμμετρικός ως προς το 0. Στην εκφώνηση θες το [0,w] να πηγαίνει στο i=0. Άρα (υποθέτω) το [-w,0] πηγαίνει στο -1. Αυτό είναι μία απλή floor.

( https://en.wikipedia.org/wiki/Quantization_(signal_processing)#Mid-riser_and_mid-tread_uniform_quantizers )

Εγώ έτσι το κάνω τουλάχιστον, τώρα εγώ κατάλαβα λάθος ή εσύ, ποιος ξέρει;

Και γω αυτο καταλαβα

[xristos_st]

Τι εννοει με το πλεγμα δειγματοληψιας?Τον καθε πινακα R G B?
Δηλαδη τι πρεπει να παρουσιασουμε στην αναφορα?

[Επικο Burger]

Τι εννοει με το πλεγμα δειγματοληψιας?Τον καθε πινακα R G B?
Δηλαδη τι πρεπει να παρουσιασουμε στην αναφορα?

Και εγω ψιλο σκαλωσα με το τι ηθελε, τελικα εβαλα απλα το bayer pattern σε πινακα, συν την εξισωση περιγραφης που εβγαλα.

[Charles]

Κιεγω αυτο κανω αλλα τα νουμερα βγαινουν οπως ναναι, αλλα αρνητικα αλλα πολυ μεγαλυτερα του ενα. Τα χ1,y1,x2,y2 δεν ειναι  συν/νες απο την αρχικη εικονα Ν0χΜ0 και τα χ,y οι συν/νες του pixel της νεας εικονας που ψαχνεις την τιμη του?

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

[Ροζ συννεφάκι]

Γράφω τους χρόνους μου για reference:

NN Interpolation:  0.1133 sec (matrix operations),  ~4 sec με for loops

Bilinear Interpolation:  0.1070 sec, ~4 sec με for loops

Save PPM: 0.107819 sec για μια εικόνα

Πώς ακριβώς τρέχεις τόσο γρήγορα το bilinear; Εμένα μου έπαιρνε πάρα πολύ ώρα και λέω κάτι θα κάνω λάθος.. Και έβαλα να κάνει τις πράξεις για το interpolation για ένα pixel και θέλει  0.0039s για το ένα χρώμα. Δηλαδή για M,N ίσα με τα M0,N0 θα θέλει γύρω στα 20 λεπτά για όλα τα χρώματα.. Κάνεις κάτι άλλο εκτός από την απλή εφαρμογή της μεθόδου;

[Gaara]

Για την κβαντιση:

Απο float να παω σε x bit int?

[Goemon]

Για την κβαντιση:

Απο float να παω σε x bit int?

Ρωτάς αν αυτό είναι η κβάντιση? λίγο πολύ ναι.

Πχ για x = 0.456 και w = 0.2   ( δλδ στάθμες κβαντισμού = 1/0.2 = 5 δλδ 3bit αρκούν)  τότε q = 3

[isitsou]

Δεν βγάζω άκρη με το τρίτο ΚΑΙ ΤΕΛΕΥΤΑΙΟ μέρος της εργασίας.  Αυτός είναι ο κώδικας μου και μου βγάζει κυριολεκτικά ότι να 'ναι εικόνα. Ανεβάζω και το αρχείο ppm(το οποίο παίρνει τιμές 0-7, όπως νομίζω πως θα πρεπε να παίρνει). Αν μπορεί κάποιος να μου δώσει τα φώτα του γτ δεν βλέπω να παίρνει παράταση και είναι κρίμα μετα απο 100ώρες στην υπόλοιπη εργασία να κολλήσω σε κάτι απλό.     

Code:

function saveasppm(x,filename,K)
[w, h, rgb] = size(x);
 
fid = fopen( filename, 'wt' );
 
fprintf( fid, 'P6 ' );
fprintf( fid, '%d %d %d\n', w, h,K);
for i=1:w
    for j=1:h
        if j~=h %otan graftei olh h prwth seira allakse seira  
        for k=1:3
            fprintf( fid, '%d ',x(i,j,k));  
        end
        else
            fprintf( fid, '%d ',x(i,j,1)); 
            fprintf( fid, '%d ',x(i,j,2));
            fprintf( fid, '%d \n',x(i,j,3));
        end
    end
end
 
fclose( fid ); 
end

PS: Δεν έχω ασχοληθεί ακόμα με τα περιθώρια της εικόνας οπότε μου την βγάζει σε μαύρο πλαίσιο... 
Επισης το ανεβάζω σε txt γιατί δεν με αφήνει σε μορφή ppm.

[ReMi0s]

Δεν βγάζω άκρη με το τρίτο ΚΑΙ ΤΕΛΕΥΤΑΙΟ μέρος της εργασίας.  Αυτός είναι ο κώδικας μου και μου βγάζει κυριολεκτικά ότι να 'ναι εικόνα. Ανεβάζω και το αρχείο ppm(το οποίο παίρνει τιμές 0-7, όπως νομίζω πως θα πρεπε να παίρνει). Αν μπορεί κάποιος να μου δώσει τα φώτα του γτ δεν βλέπω να παίρνει παράταση και είναι κρίμα μετα απο 100ώρες στην υπόλοιπη εργασία να κολλήσω σε κάτι απλό.    

Code:

function saveasppm(x,filename,K)
[w, h, rgb] = size(x);
 
fid = fopen( filename, 'wt' );
 
fprintf( fid, 'P6 ' );
fprintf( fid, '%d %d %d\n', w, h,K);
for i=1:w
    for j=1:h
        if j~=h %otan graftei olh h prwth seira allakse seira  
        for k=1:3
            fprintf( fid, '%d ',x(i,j,k));  
        end
        else
            fprintf( fid, '%d ',x(i,j,1)); 
            fprintf( fid, '%d ',x(i,j,2));
            fprintf( fid, '%d \n',x(i,j,3));
        end
    end
end
 
fclose( fid ); 
end

PS: Δεν έχω ασχοληθεί ακόμα με τα περιθώρια της εικόνας οπότε μου την βγάζει σε μαύρο πλαίσιο...  
Επισης το ανεβάζω σε txt γιατί δεν με αφήνει σε μορφή ppm.

Με μία ματιά είδα διάφορα, αρχικά έχεις ανάποδα το w και το h στο header.
Κατα δεύτερον θέλει fwrite για τα χρώματα h fprintf με %c
Επίσης στα χρώματα δεν θέλει ούτε κενα ούτε νέες γραμμές ΚΑΘΟΛΟΥ
Γενικά δες περίπου στην σελίδα 4 έχουμε γράψει αρκετά για το Save

[panbamid]

Επειδη απο οτι ειδα εχει ρωτηθει και πιο πριν και δεν εχει απαντηθει, μπορει καποιος να μας πει πως εκανε την διγραμμικη παρεμβολη για το πρασινο;

[Επικο Burger]

Με μία ματιά είδα διάφορα, αρχικά έχεις ανάποδα το w και το h στο header.
Κατα δεύτερον θέλει fwrite για τα χρώματα h fprintf με %c
Επίσης στα χρώματα δεν θέλει ούτε κενα ούτε νέες γραμμές ΚΑΘΟΛΟΥ
Γενικά δες περίπου στην σελίδα 4 έχουμε γράψει αρκετά για το Save

Βασικο τα χρωματα να ειναι RGBRGBRGB... σε μια σειρα συνεχομενα χωρις κανενα κενο και νεες γραμμες οπως λεει ο ReMi0s.

Επίσης, νομιζω πιο καλα βγαινει χωρις for, απλα γραφοντας ολα τα χρωματα της εικονας σε ενα vector, δηλαδη [R(1,1) G(1,1) B(1,1) R(1,2) G(1,2) B(1,2) ...], και γραφοντας αυτό στο αρχείο.

[ReMi0s]

Βασικο τα χρωματα να ειναι RGBRGBRGB... σε μια σειρα συνεχομενα χωρις κανενα κενο και νεες γραμμες οπως λεει ο ReMi0s.

Επίσης, νομιζω πιο καλα βγαινει χωρις for, απλα γραφοντας ολα τα χρωματα της εικονας σε ενα vector, δηλαδη [R(1,1) G(1,1) B(1,1) R(1,2) G(1,2) B(1,2) ...], και γραφοντας αυτό στο αρχείο.

Έχεις δίκιο για αυτό με τον vector. Έτσι το έχω κάνει και είναι πολύ γρήγορο.

Επειδη απο οτι ειδα εχει ρωτηθει και πιο πριν και δεν εχει απαντηθει, μπορει καποιος να μας πει πως εκανε την διγραμμικη παρεμβολη για το πρασινο;

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

Στην περίπτωση που είσαι μέσα σε πράσινο και κανεις interpolation για να βρεις την τιμή στον νέο πίνακα.
Έστω f0 αριστερα, f1 πανω, f2 δεξιά kai f3 κάτω. Και εννοείται τα κοντινότερα.... Εδώ γίνεται ο χαμός..

I1=f1*(y2-j)/(y2-y1)+f2*(j-y1)/(y2-y1);
I2=f0*(y3-j)/(y3-y0)+f3*(j-y0)/(y3-y0);
I=I1*(x3-i)/(x3-x1)+I2*(i-x1)/(x3-x1);


Έστω οτι τώρα είσαι μέσα σε κόκκινο η μπλε και ψάχνεις πάλι τα πράσινα.
Όπου bayerI=floor(i) και bayerJ=floor(j)

I1=Green(bayerI-1,bayerJ)*(bayerJ+1-j)+Green(bayerI,bayerJ+1)*(j-bayerJ);
I2=Green(bayerI,bayerJ-1)*(bayerJ+1-j)+Green(bayerI+1,bayerJ)*(j-bayerJ);
I=I1*(bayerI-i)+I2*(i-bayerI+1);

Ξανα λέω πως δεν ξέρω αν είναι θεωρητικά ορθή προσέγγιση.

[Επικο Burger]

Επειδη απο οτι ειδα εχει ρωτηθει και πιο πριν και δεν εχει απαντηθει, μπορει καποιος να μας πει πως εκανε την διγραμμικη παρεμβολη για το πρασινο;

Το θέμα νομιζω δεν ειναι τοσο το πρασινο, όσο οι τιμες των κοκκινων και μπλε, στα πρασινα pixel.

Το βασικό με το πρασινο ειναι να ελέγξεις τι πρασινο έχεις.

Πχ στην απλή υλοποιηση με for:  
1)Αν είσαι σε πρασινο σε γραμμη με κοκκινα-πρασινα, το κοκκινο εινα ο MO του δεξιου και αριστερου πιξελ, ενώ του μπλε είναι ο ΜΟ του πάνω και του κάτω.
2) Αν είσαι σε πρασινο σε γραμμη με μπλε-πρασινα, το κοκκινο ειναι ο ΜΟ του πάνω και του κάτω πιξελ, ενώ του μπλέ ειναι ο ΜΟ του δεξιου και αριστερού.

Για κόκκινα και μπλε έχεις κάθε φορα, για κάθε interpolated χρώμα, τα πάνω κατω δεξια και αριστερα.

Έλπιζω να βοηθησα και να μην χαωσα περισσοτερο

[Επικο Burger]

Πώς ακριβώς τρέχεις τόσο γρήγορα το bilinear; Εμένα μου έπαιρνε πάρα πολύ ώρα και λέω κάτι θα κάνω λάθος.. Και έβαλα να κάνει τις πράξεις για το interpolation για ένα pixel και θέλει  0.0039s για το ένα χρώμα. Δηλαδή για M,N ίσα με τα M0,N0 θα θέλει γύρω στα 20 λεπτά για όλα τα χρώματα.. Κάνεις κάτι άλλο εκτός από την απλή εφαρμογή της μεθόδου;

20 λεπτα για το bilinear είναι σίγουρα πολλά

Ο αργός μου χρονος, τα 4 sec, ειναι με υλοποιηση με 2 nested for loops και ελεγχο για καθε pixel.

Ο γρηγορος χρονος ειναι χωρις καθολου for, μόνο με πολλαπλασιασμους πινάκων και μασκες.

[ReMi0s]

Το θέμα νομιζω δεν ειναι τοσο το πρασινο, όσο οι τιμες των κοκκινων και μπλε, στα πρασινα pixel.

Το βασικό με το πρασινο ειναι να ελέγξεις τι πρασινο έχεις.

Πχ στην απλή υλοποιηση με for:  
1)Αν είσαι σε πρασινο σε γραμμη με κοκκινα-πρασινα, το κοκκινο εινα ο MO του δεξιου και αριστερου πιξελ, ενώ του μπλε είναι ο ΜΟ του πάνω και του κάτω.
2) Αν είσαι σε πρασινο σε γραμμη με μπλε-πρασινα, το κοκκινο ειναι ο ΜΟ του πάνω και του κάτω πιξελ, ενώ του μπλέ ειναι ο ΜΟ του δεξιου και αριστερού.

Για κόκκινα και μπλε έχεις κάθε φορα, για κάθε interpolated χρώμα, τα πάνω κατω δεξια και αριστερα.

Έλπιζω να βοηθησα και να μην χαωσα περισσοτερο

Αυτό που λές ισχύει μόνο για όταν η έξοδος είναι ίδιας ανάλυσης με την είσοδο. Είναι δηλαδή υποπερίπτωση του γενικού αλγορίθμου.
Ναι θα βγάλεις καλή εικόνα, ωστόσο δεν είναι ο θεωρητικά ορθότερος τρόπος, γιατί αμα το σκεφτείς είναι απλά μέσος όρος της γειτονιάς. Ο bilinear έχει αποστάσεις για βάρη.