编辑一个Matlab小程序关于提高CT图像清晰度识别肿瘤结构并改变其颜色并截图带注释

以下是一个简单的Matlab小程序，用于提高CT图像的清晰度、识别肿瘤结构并改变其颜色，并截取带有注释的图像。具体步骤如下：

1. 读入CT图像并显示

img = imread('CT_image.jpg');
imshow(img);

2. 对图像进行锐化处理，提高图像清晰度

sharp_img = imsharpen(img);
imshow(sharp_img);

3. 对图像进行阈值分割，将肿瘤结构分离出来

threshold = graythresh(sharp_img);
bw_img = imbinarize(sharp_img, threshold);
imshow(bw_img);

4. 对分割后的二值图像进行形态学处理，填充空洞并去除噪点

se = strel('disk', 5);
bw_img = imclose(bw_img, se);
bw_img = imfill(bw_img, 'holes');
bw_img = bwareaopen(bw_img, 50);
imshow(bw_img);

5. 将肿瘤结构标记出来，并改变其颜色

labeled_img = bwlabel(bw_img);
rgb_img = label2rgb(labeled_img, 'jet', 'k', 'shuffle');
imshow(rgb_img);

6. 在标记后的图像上添加注释，如肿瘤大小、位置等信息

stats = regionprops(labeled_img, 'Area', 'Centroid');
for i = 1:length(stats)
    area = stats(i).Area;
    centroid = stats(i).Centroid;
    text(centroid(1), centroid(2), ['Area: ' num2str(area)], 'Color', 'white', 'FontSize', 12, 'HorizontalAlignment', 'center');
end

7. 将带有注释的图像截取下来，并保存

annotated_img = getframe(gca);
imwrite(annotated_img.cdata, 'annotated_CT_image.jpg');

完整代码如下：

% 读入CT图像并显示
img = imread('CT_image.jpg');
subplot(2,2,1);
imshow(img);
title('Original Image');

% 对图像进行锐化处理，提高图像清晰度
sharp_img = imsharpen(img);
subplot(2,2,2);
imshow(sharp_img);
title('Sharpened Image');

% 对图像进行阈值分割，将肿瘤结构分离出来
threshold = graythresh(sharp_img);
bw_img = imbinarize(sharp_img, threshold);
subplot(2,2,3);
imshow(bw_img);
title('Binary Image');

% 对分割后的二值图像进行形态学处理，填充空洞并去除噪点
se = strel('disk', 5);
bw_img = imclose(bw_img, se);
bw_img = imfill(bw_img, 'holes');
bw_img = bwareaopen(bw_img, 50);
subplot(2,2,4);
imshow(bw_img);
title('Processed Image');

% 将肿瘤结构标记出来，并改变其颜色
labeled_img = bwlabel(bw_img);
rgb_img = label2rgb(labeled_img, 'jet', 'k', 'shuffle');
figure;
imshow(rgb_img);
title('Labeled Image');

% 在标记后的图像上添加注释，如肿瘤大小、位置等信息
stats = regionprops(labeled_img, 'Area', 'Centroid');
for i = 1:length(stats)
    area = stats(i).Area;
    centroid = stats(i).Centroid;
    text(centroid(1), centroid(2), ['Area: ' num2str(area)], 'Color', 'white', 'FontSize', 12, 'HorizontalAlignment', 'center');
end

% 将带有注释的图像截取下来，并保存
annotated_img = getframe(gca);
imwrite(annotated_img.cdata, 'annotated_CT_image.jpg');