added final folder with the ultimate code

Author: synchro--

.gitignore

@@ -1,5 +1,4 @@
 *~
-*.mat
 # Packages #
 ############
 # it's better to unpack these files and commit the raw source

Code/compute_oriented_hist

@@ -0,0 +1,32 @@
+function [CBnplus1, dam] = grow_regions_inside_Q2(CB, Q)
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%   [CBnplus1, dam] = grow_regions_inside_Q(CB, Q)
+%
+%   Grows the connected regions of CB that lie inside Q
+%
+%   Inputs:     CB -  Matrix with the labelled regions
+%               Q  -  Binary matrix with the connected region Q
+%
+%   Outputs:   CBnplus1 - Matrix with the grown labelled regions
+%              dam      - Binary matrix identifying the dam points
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+nhood = strel([1 1 1; 1 1 1; 1 1 1]);
+dam   = zeros(size(CB));
+
+CBnplus1 = CB; % para no tocar CB 
+while(~isequal(Q, ((CBnplus1 > 0) & Q) | dam)), % la union de los puntos que pertenecen a un CB o que estan en una presa tienen que ser igual a q. continua hasta que esté lleno, hasta que se hayan llenado del todo. or dam a lo mejor se puede tocar
+    newq  = Q & (CBnplus1 > 0); % 
+    nlist = get_number_list(CBnplus1(newq)); % nqint again
+
+    tmp1 = zeros(size(CB)); % posiciones de dam
+    tmp2 = CBnplus1; % nuevas posiciones de CB. lo inicializamos asi
+    for k=nlist; % para los cb que esten en q
+        imdk = imdilate(CBnplus1==k, nhood) & Q; % no puede salir de q cuando se dilata. 
+        tmp1 = tmp1 + imdk;  % le suma estos puntos
+        tmp2(imdk & (CBnplus1 == 0)) = k; % ahora pertenecen al CB k. 
+    end;
+
+    dam = tmp1 > 1; % se han solapado dos dilataciones
+    CBnplus1 = tmp2; 
+end;

Code/HOG/compute_histogram_bins.m Code/final/compute_histogram_bins.m

@@ -16,7 +16,7 @@
 % gmax=sgolayfilt(gmax,3,7);
 
 %threshold filtering
-thresh=5;
+thresh=8;
 r45(r45<thresh)=0;
 r125(r125<thresh)=0;
 gmax(gmax<thresh)=0;
@@ -40,17 +40,18 @@
 tmp=max(grad45,grad125);
 final_grad=max(gmax,tmp);
 figure(999)
-imshow(uint8(final_grad));
+imshow(uint8(final_grad)),title('Oriented Gradient');
 
 %%%%final post-processing %%%%%
 sgo_grad=sgolayfilt(final_grad,2,7);
-median_grad=medfilt2(final_grad,[4 4]);
+median_grad=medfilt2(final_grad,[3 3]);
+median_grad(median_grad<thresh)=0;
 
 
 figure(1)
-imshow(uint8(sgo_grad));
-title('median');
+imshow(uint8(median_grad));
+title('Median filter');
 
 figure(2)
-imshow(uint8(median_grad));
-title('sgolay filter');
+imshow(uint8(sgo_grad));
+title('Sgolay filter');