如何在没有 Matlab 内置函数的情况下编写 3D 直方图代码？

彩色三元组直方图

请记住，绘制整个 RGB 三元组直方图是一项艰巨的任务，因为您可以拥有 256 × 256 × 256 = 16,777,216 种可能的独特颜色组合。我认为一个稍微更易于管理的任务是计算图像中唯一 RGB 值的直方图（因为其余部分无论如何都会为零）。这仍然是一项相当大的任务，但如果图像相当小，则可能是合理的。下面我认为一个不错的分箱替代方法是将图像调整为较小数量的像素。这可以通过使用该imresize函数来完成。我相信这会降低图像的保真度，并且几乎充当可以“有点”模拟分箱行为的舍入函数。在本例中，我将矩阵字符串数组转换为连接通道L1im、L2im和L3im的图像。下面是一个演示，我使用内置于 MATLAB 中的名为saturn.png. A Resize_Factorof 1 将导致最多的 bin 数量，并且 bin 的数量会随着Resize_Factor增加而减少。请记住，如果使用Resize_Factor.

Resize_Factor = 200;
RGB_Image = imread("saturn.png");
[Image_Height,Image_Width,Number_Of_Colour_Channels] = size(RGB_Image);
Number_Of_Pixels = Image_Height*Image_Width;

RGB_Image = imresize(RGB_Image,[Image_Height/Resize_Factor Image_Width/Resize_Factor]);
L1im = RGB_Image(:,:,1);
L2im = RGB_Image(:,:,2);
L3im = RGB_Image(:,:,3);

L1im_String = string(L1im);
L2im_String = string(L2im);
L3im_String = string(L3im);

RGB_Triplets = L1im_String + "," + L2im_String + "," + L3im_String;
Unique_RGB_Triplets = unique(RGB_Triplets);

for Colour_Index = 1: length(Unique_RGB_Triplets)
    RGB_Colour = Unique_RGB_Triplets(Colour_Index);
    Unique_RGB_Triplets(Colour_Index,2) = nnz(RGB_Colour == RGB_Triplets); 
end

Counts = str2double(Unique_RGB_Triplets(:,2));
Scaling_Factor = Number_Of_Pixels/sum(Counts);
Counts = Counts.*Scaling_Factor;

if sum(Counts) == Number_Of_Pixels
    disp("Sum of histogram is equal to the number of pixels");
end

bar(Counts);
title("RGB Triplet Histogram");
xlabel("RGB Triplets"); ylabel("Counts");

Current_Axis = gca;
Scale = (1:length(Unique_RGB_Triplets));
set(Current_Axis,'xtick',Scale,'xticklabel',Unique_RGB_Triplets);
Angle = 90;
xtickangle(Current_Axis,Angle);