获取直方图均衡的切割点 - matlab

Question

我正在尝试获取用于创建直方图均衡化的切割点（通过histeq 或以其他方式），即获取用于确定每个 bin 边缘的原始像素值。

为了简化下面的例子，我将使用一维矩阵：

如果原图是：[ 0.2 , 0.25 , 0.34 , 0.4 , 0.4 , 0.4 , 0.6 , 0.6 , 0.6 ]

直方图均衡后的图像（到3个bin）是：[ 0.2 , 0.2 , 0.2 , 0.4 , 0.4 , 0.4 , 0.6 , 0.6 , 0.6 ]

如何获得 3 个 bin 的边缘？

即如何获得以下向量（代表切割点）：[0.37 0.5]？

非常感谢！！！

附：得到向量[0.34 0.4] 或[0.4 0.6]（分别是histeq 之后的特定bin 中包含的最小值和最大值）也很令人满意。

Answer 1

（如果我的问题是正确的）我可以想象做如下的事情：

dVec = diff(eqdata);
fVec = find(dVec);
cuttingPoints = eqdata(fVec) + dVec(fVec)/2;

例子：

eqdata =  [ 0.2 , 0.2 , 0.2 , 0.4 , 0.4 , 0.4 , 0.6 , 0.6 , 0.6 ];

输出：

0.3000    0.5000

编辑：如果数据未排序：

dVec = diff(sort(unique(eqdata)));
fVec = find(dVec);
cuttingPoints = eqdata(fVec) + dVec(fVec)/2;

Answer 2

这是一个适用于向量和矩阵（已排序或未排序）的解决方案。它使用unique 和median 函数。

您将在 bin 值之间获得中间值（此处使用 median）。您还将拥有分配给每个图像元素的 bin 索引。

以下代码

%data
A = [ 0.2 , 0.25 , 0.34 , 0.4 , 0.4 , 0.4 , 0.6 , 0.6 , 0.6 ]; % vector
%A = rand(5,5); % matrix (other possible input)
bin_number = 3;    

%histogram
J = histeq(A,bin_number); %histogram equalization
[Ju, ia, ic] = unique(J(:)); %uniquing

%values of in-between bins (here using median)
cutting_index = median([Ju(1:end-1) Ju(2:end)]');

%get the assignment
assign_index = reshape(ic,size(A));

会产生

cutting_index = [0.2500    0.7500]
assign_index = [1     1     1     2     2     2     3     3     3]'