img1 = imread("MRI.jpg");

img2 = imread("CT.jpg");

% Convert the images to grayscale

gray1 = rgb2gray(img1);

gray2 = rgb2gray(img2);

% Compute the correlation coefficient between the two images

corrCoef = corr2(gray1, gray2);

% Blend the two images using the correlation coefficient

alpha = corrCoef; % Weight for the first image

beta = 1 - corrCoef; % Weight for the second image

blendedImg = alpha*img1 + beta*img2;

% Display the results

figure;

subplot(2,2,1); imshow(img1); title('MRI IMAGE');

subplot(2,2,2); imshow(img2); title('CT IMAGE');

subplot(2,2,3); imshow(blendedImg); title('FUSED IMAGE');