deok5.m

close all
clear

load sampleMovie

NUM_DEPTH_THRESHOLD = 10;
NUM_NOISE_THRESHOLD = 0.05;
NUM_GRAD_THRESHOLD = 100;

%
colormap(jet);
I = K0(1).cdata(:,:,1);

I2 = I/NUM_DEPTH_THRESHOLD;
I3 = I2*NUM_DEPTH_THRESHOLD;

figure, imshow(I);

figure, imshow(I3), colormap('jet');






I4 = zeros(size(I));



for i =1:max(max(I2)) 
	temp = I2 == i;
	
	[B,L,N,A] = bwboundaries(temp,'noholes');
	

	for j=1:length(B)
		
		a = cell2mat(B(j));
		
		[m n] = size(a);
		
		for k = 1:m
			I4(a(k,1),a(k,2)) = i;
		end
		
		
	end
	
end

I5 = I4 * NUM_DEPTH_THRESHOLD;
figure, image(I5),colormap('jet');

%Appy Sobel to extract edge
hy = fspecial('sobel');
hx = hy';
Iy = imfilter(double(I), hy, 'replicate');
Ix = imfilter(double(I), hx, 'replicate');
gradmag = sqrt(Ix.^2 + Iy.^2);
figure, imshow(gradmag,[]), title('Gradient magnitude (gradmag)')


%Noise Identification 
sum = 0;
[m n] = size(I);
total = m *n;
for i = 0:255
	
	sum = sum + length(find(I == i)); 
	
	percentSum = sum/total;
	
	if percentSum > NUM_NOISE_THRESHOLD 
		noiseIndex = i;
		break;
	end
	
	
end

Inoise = (gradmag > 50) & (I < noiseIndex);
SE = strel('disk',1,0);
Inoise = imdilate(Inoise,SE);
Inoise = imfill(Inoise,'holes');
figure, imshow(Inoise);


[L num] = bwlabel(Inoise);

[B,L,N,A] = bwboundaries(Inoise,'noholes');
	


[row col] = find(Inoise);

Inoise2 = double(I);

for i=1:length(row)
	Inoise2(row(i),col(i)) = NaN;
end
figure, imshow(Inoise2);

Inoise3 = Inoise2;

while(find(isnan(Inoise3))) 
	Inoise3 = nlfilter(Inoise3, [3,3], @neighError);
	figure, imshow(Inoise3);
end



% Inoise3 = double(Inoise2);
% Inoise4 = inpaint_nans(Inoise3,6);
% figure, imshow(Inoise4);
% Inoise5 = Inoise4/256;
% imshow(Inoise5)

% figure, imshow(Inoise);
% 
% 
% 
% Inoise2 = immultiply(I,~Inoise) + uint8(Inoise*255);
% figure, imshow(Inoise2);

I = Inoise3;


%Appy Sobel to extract edge
hy = fspecial('sobel');
hx = hy';
Iy = imfilter(double(I), hy, 'replicate');
Ix = imfilter(double(I), hx, 'replicate');
gradmag = sqrt(Ix.^2 + Iy.^2);
figure, imshow(gradmag,[]), title('Gradient magnitude (gradmag)')


%Noise Identification 
sum = 0;
[m n] = size(I);
total = m *n;
for i = 0:255
	
	sum = sum + length(find(I == i)); 
	
	percentSum = sum/total;
	
	if percentSum > NUM_NOISE_THRESHOLD 
		noiseIndex = i;
		break;
	end
	
	
end

Inoise = (gradmag > 50) & (I < noiseIndex);
SE = strel('disk',1,0);
Inoise = imdilate(Inoise,SE);
Inoise = imfill(Inoise,'holes');
figure, imshow(Inoise);


[L num] = bwlabel(Inoise);

[B,L,N,A] = bwboundaries(Inoise,'noholes');
	


[row col] = find(Inoise);

Inoise2 = double(I);

for i=1:length(row)
	Inoise2(row(i),col(i)) = NaN;
end
figure, imshow(Inoise2);

Inoise3 = Inoise2;

while(find(isnan(Inoise3))) 
	Inoise3 = nlfilter(Inoise3, [3,3], @neighError);
	figure, imshow(Inoise3);
end



% Inoise3 = double(Inoise2);
% Inoise4 = inpaint_nans(Inoise3,6);
% figure, imshow(Inoise4);
% Inoise5 = Inoise4/256;
% imshow(Inoise5)

% figure, imshow(Inoise);
% 
% 
% 
% Inoise2 = immultiply(I,~Inoise) + uint8(Inoise*255);
% figure, imshow(Inoise2);