[Παράλληλα και Διανεμημένα Συστ.] Εργασια 3

[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 με προσοχή στα ονόματα των αρχείων!

[height=400]

  Όταν ανεβάζουμε φωτογραφίες στις Ανακοινώσεις και Έκτακτα νέα, βάζουμε τη μεγαλύτερη πλευρά 400 (width=400 ή height=400 ). π.χ. [img height=400 (κλείνει η αγκύλη) 

[ReMi0s]

.To Z(i) είναι για κάθε pixel καθώς καθορίζει ''ποσοστό ομοιότητας'' με τα άλλα pixel.
@Andromedas:Τι εννοείς σε αυτό το σημείο,τι πολλαπλασιάζεις;

Επίσης να κάνω μια ερώτηση στο παράδειγμα του matlab ξέρετε τι είναι τα filtSigma και patchSigma.Εγώ καταλαβαίνω ότι το filtSigma είναι ο όρος σ που έχει στους τύπους στο pdf αλλά το patchSigma σε τι αντιστοιχεί,ή μήπως τα έχω μπερδέψει;

δεν είμαι σίγουρος αλλα αυτά τα δύο μαζί πρέπει να κάνουν το σ^2 . Γιατί για διαφορετικό patchSize βγάζεις διαφορετικά αποτελέσματα, οπότε νομίζω το patchSigma είναι πολλαπλασιαστής του filtSigma για να κρατάει αναλογία. Ολα αυτά με επιφύλαξη

[Andromedas]

2)@Andromedas:Τι εννοείς σε αυτό το σημείο,τι πολλαπλασιάζεις;
1) patchSigma ;

Λοιπόν μια γρήγορη απάντηση.

Code:

1)
..
  % create 3D cube
  B = patchCube(I, patchSize); %ανοίγει τον πίνακα ανάλογα με το patchSize πρακτίκα εάν έχεις μια I 3x3 και %patch size 3x3 κάνεις ένα padding γύρω από κάθε pixel άρα κάθε z της εικόνας 3x3x9 είναι το i στοιχείο με το %padding του γύρω γύρω. 
 [m, n, d] = size( B );
  B = reshape(B, [ m*n d ] );
% κάνει τον m x n x z πίνακα σε  m*n x z δηλαδή μια γραμμή είναι το i pixel μαζί με το padd του 


2)
% φτίαχνει την gaussian σε 3x3 (patchSize)
  % gaussian patch
  H = fspecial('gaussian',patchSize, patchSigma);
  H = H(:) ./ max(H(:));
% την ανοίγει σε 9x1  και την πολλαπλασιάζει με κάθε i pixel μαζί με το pad του (εδώ εάν κάνεις αθροισμά έχεις δυσδιάστατη συνελιξή )
  % apply gaussian patch on 3D cube
  B = bsxfun( @times, B, H' );

δεν είμαι σίγουρος αλλα αυτά τα δύο μαζί πρέπει να κάνουν το σ^2 . Γιατί για διαφορετικό patchSize βγάζεις διαφορετικά αποτελέσματα, οπότε νομίζω το patchSigma είναι πολλαπλασιαστής του filtSigma για να κρατάει αναλογία. Ολα αυτά με επιφύλαξη

Όχι το σ^2 είναι μόνο του και εφαρμόζεται στο
 D = exp( -D.^2 / filtSigma ); και στην gaussian
Με άλλο patch size αλλάζει η απόδοση καθώς μεγαλώνει ή μικραίνει το παράθυρό που παίρνει τον "μέσο όρο" και αλλάζει και την gaussian.

ελπίζω να βοήθησα

[ReMi0s]

για μεγαλύτερο patchsize πρεπει να παίρνεις καλύτερα αποτελέσματα?
εμένα δουλεύει καλύτερα για 3x3

[Andromedas]

για μεγαλύτερο patchsize πρεπει να παίρνεις καλύτερα αποτελέσματα?
εμένα δουλεύει καλύτερα για 3x3

Δεν είμαι σίγουρος μπορείς να το τεστάρει και να μας πεις 

[ReMi0s]

Δεν είμαι σίγουρος μπορείς να το τεστάρει και να μας πεις  

Στο original φαίνεται να παίρνει καλύτερα αποτελέσματα..
Η ερώτηση μου είναι αν το έχει καταφέρει κανένας.
Επίσης στο original όταν το τρέχεις για διαφορετικά patchSize, πρέπει να αλλάζεις και το patchSigma?
Έχουν σχέση μεταξύ τους;
patchSize = [3 3];
  patchSigma = 5/3;

Δλδ για [5 5] πρέπει να γίνει 5/5 , και [7 7] να γίνει 5/7? Η δεν έχουν σχέση;

επίσης είναι κανείς που έχει αρκετή Ram για να τρέξει το 256*256 χωρίς cuda, για να μας πει χρόνους ? Δεν το επιτρέπει στο διάδη

[kaspas]

Τους πίνακες Β, Η (αυτούς που έχει στον αλγόριθμό του στην matlab αναφέρομαι) που έχει τους πήρατε δεδομένους και υλοποιήσατε όλα τα αποκάτω ή δουλεύεται με τον πίνακα της εικόνας;

[ReMi0s]

Τους πίνακες Β, Η (αυτούς που έχει στον αλγόριθμό του στην matlab αναφέρομαι) που έχει τους πήρατε δεδομένους και υλοποιήσατε όλα τα αποκάτω ή δουλεύεται με τον πίνακα της εικόνας;

Εγώ προσπάθησα να χρησιμοποιήσω τον πίνακα B αλλα κουράστηκα γιατί δεν μου έβγαινε (μεγάλο μέγεθος, διαφορετικά μεγέθη πινάκων μεσα στο cuda kernel). Οπότε το υλοποίησα επεξεργάζοντας τον πίνακα J.

[fanaks]

Τους πίνακες Β, Η (αυτούς που έχει στον αλγόριθμό του στην matlab αναφέρομαι) που έχει τους πήρατε δεδομένους και υλοποιήσατε όλα τα αποκάτω ή δουλεύεται με τον πίνακα της εικόνας;

Εγώ προσωπικά δουλεύω απευθείας πάνω στον πίνακα της εικόνας,αλλά διαβάζοντας και τα παραπάνω αναρωτιέμαι αν ο πίνακας της εικόνας θα πρέπει να είναι 2D ή 3D,λογικά δεν θα πρέπει να είναι 2 διαστάσεων όπως και στο παράδειγμα και να γράψουμε το kernel πρόγραμμα ώστε να χειρίζεται δισδιάστατο πίνακα?

[Aristos]

το matlab χρησιμοποιεί μονοδιάστατους πίνακες για να αποθηκεύσει οποιονδήποτε πίνακα στη μνήμη. παίζεις απλά με τα indexes, όπως και στην δεύτερη εργασία

[lodi]

ξέρει κανένας να πει πώς θα κάνουμε debugging στον Διάδη την cuda στο matlab ?

[ReMi0s]

ξέρει κανένας να πει πώς θα κάνουμε debugging στον Διάδη την cuda στο matlab ?

γενικά το debugging ειναι δύσκολο γιατι δεν σου δίνει κατατοπιστικά errors. ειδικά στον διάδη θα πάρει πολύ ώρα...
το καλύτερο που εχεις να κάνεις είναι να βαζεις σε σχόλια σημεία του κώδικα και σιγά σιγά να βρεις το λάθος.
Όταν παίρνεις Cuda_error_illegal_adress χρειάζεται δυστυχώς να κάνεις restart το matlab.


Το κομμάτι με την shared memory το καταφερε κανεις;
Είναι πολύ μικρή για να χωρέσει ολόκληρο οποιονδήποτε πίνακα εικόνας...

[lodi]

γενικά το debugging ειναι δύσκολο γιατι δεν σου δίνει κατατοπιστικά errors. ειδικά στον διάδη θα πάρει πολύ ώρα...
το καλύτερο που εχεις να κάνεις είναι να βαζεις σε σχόλια σημεία του κώδικα και σιγά σιγά να βρεις το λάθος.
Όταν παίρνεις Cuda_error_illegal_adress χρειάζεται δυστυχώς να κάνεις restart το matlab.


Το κομμάτι με την shared memory το καταφερε κανεις;
Είναι πολύ μικρή για να χωρέσει ολόκληρο οποιονδήποτε πίνακα εικόνας...

μπορείς να χρησιμοποιήσεις την τεχνική blocking για αυτό. Σπας τον πίνακα-εικόνα σε blocks και τα φορτώνεις ένα ένα στην shared memory όπου κάνεις την επεξεργασία

[ReMi0s]

Ψάχνωντας για την γειτονιά κάθε σημείου, όταν αυτό βρίσκεται σε ακριανό σημείο τι πρέπει να κάνουμε για τους γείτονες που δεν υπάρχουν;
Εστω δλδ το Pixel (0,0) και γειτονια σημείων 3x3, μόνο τα (0,1),(1,0) και (1,1) υπάρχουν στην γειτονιά αυτή, για τα υπόλοιπα σημεία μπορούμε να κάνουμε κάτι; Αν τα βάλουμε 0 είναι λάθος γιατί το 0 αντιπροσωπεύει το μαύρο..

[fanaks]

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

[ReMi0s]

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

αυτό που λες απαιτεί να αποθηκεύσω όλες τις γειτονιές σε πίνακα , που είναι χρονοβόρο

Επίσης τι πράξη κάνετε μεταξύ των γειτόνων; Εννοώ πως ερμηνεύσατε τον τύπο του pdf |f(Ni)-f(Nj)|^2