img = imread('CAM1_C001H001S0014001202.png');
edges = edge(img, 'sobel');
[H, T, R] = hough(edges);
P = houghpeaks(H, 1000, 'threshold', ceil(0.1 * max(H(:))));
lines = houghlines(edges, T, R, P, 'FillGap', 5, 'MinLength', 3);
img = imfill(img, 'holes');
img =imclose(img,strel('disk',8));
img =imopen(img,strel('disk',8));
title('After Pre-Processing')
xy = [lines(k).point1; lines(k).point2];
angle = atan2d(diff(xy(:, 2)), diff(xy(:, 1)));
if abs(angle) > 80 && abs(angle) < 100
vertical_lines = [vertical_lines; xy(1, 1)];
vertical_lines = sort(vertical_lines);
for k = 1:length(vertical_lines)
if vertical_lines(k) - last_x >= 25
goodFilteredLines = [goodFilteredLines, vertical_lines(k)];
last_x = vertical_lines(k);
for k = 1:length(goodFilteredLines)
if k == 3 || k == 9 || k == 21
plot([goodFilteredLines(k), goodFilteredLines(k)], [0, 1024], 'LineWidth', 1, 'Color', 'k');
capturedImage = frame.cdata;
imgGray = (rgb2gray(capturedImage));
imwrite(imgGray, 'processedImage.png');