Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- %%
- clear all; close all; clc;
- %H = Height
- %W = Width
- %I = Slice;
- OpenfilePath = 'C:\Users\TBL\Documents\MATLAB\1708_AMORE\01_KB\backhand'
- max_step = 500;
- %%
- for step = 1:max_step
- % 2 -> 0002
- read_file_name = strcat('Image',num2str(step-1,'%04d'),'.jpg');
- img = imread([OpenfilePath,'\' ,read_file_name]);
- %
- img3d(:,:,step) = img;
- fprintf('%i\n', step);
- end
- [H,W,I] = size(img3d);
- %%
- avgTop = zeros(I,W); %top view
- for j = 1:W
- for i = 1:I
- avgTop(i,j) = sum(double(img3d(:,j,i)))/H;
- end
- end
- avgTop = avgTop(:,[65:W-65]);
- %result1 = avgTop([50:W-50],:);
- [H1,W1] = size(avgTop);
- figure(); imshow(avgTop,[]);
- %%
- avgSide = zeros(H,I); %top view
- for j = 1:I
- for i = 1:H
- avgSide(i,j) = sum(double(img3d(i,:,j)))/W;
- end
- end
- avgSide = avgSide([20:H-400],:);
- % result2 = avgSide([50:H-50],:);
- [H2,W2] = size(avgSide);
- figure(); imshow(avgSide,[]);
- %%
- % avgSide2 = zeros(H,I); %top view
- % avgSide2 = img3d(:,50,:);
- % avgSide2 = squeeze(avgSide2);
- %
- % figure();
- % imshow(avgSide2,[]);
- %%
- % for i = 1:H
- % for j = 1:W
- % if result(i,j) < 8
- % result(i,j) = 0;
- % else
- % result(i,j) = 254;
- % end
- % end
- % end
- %
- % figure;
- % imshow(result,[]);
- %
- % totalPixels = H*W;
- % wrinklePixels = 0;
- %
- % for i = 1:H
- % for j = 1:W
- % if result(i,j) ~= 0
- % wrinklePixels = wrinklePixels + 1;
- % end
- % end
- % end
- %
- % wrinklePercentage = wrinklePixels*100/totalPixels
Advertisement
Add Comment
Please, Sign In to add comment