clc
clear

filePath = 'D:\Desktop\ENGGEN 131\MATLAB Project\SuppliedCode\lunar';
extension = 'jpg';
filenames = GenerateImageList(filePath, extension);
imageList = ReadImages(filePath, filenames);

tic


imageMedian = double(RemoveAction(imageList));
imageStacked = double(cat(4,imageList{:}));
distanceArray = zeros(size(imageMedian, 1), size(imageMedian, 2), size(imageStacked, 4));



for i = 1:size(imageStacked,4)
	distanceArray(:,:,i) = (imageMedian(:,:,1) - imageStacked(:,:,1,i)).^2 +...
	(imageMedian(:,:,2) - imageStacked(:,:,2,i)).^2 + ...
	(imageMedian(:,:,3) - imageStacked(:,:,3,i)).^2;
end


[~, maxPostion] = max(distanceArray,[],3);

[rows, cols, colours] = ndgrid(1:size(imageStacked,1), 1:size(imageStacked,2), 1:3);

maxPositionColour = repmat(maxPostion, 1, 1, 3);
colourLinearPosition = sub2ind(size(imageStacked), rows, cols, colours, maxPositionColour);
finalImage = uint8(imageStacked(colourLinearPosition));

toc