Question

图像直方图归一化，下面的公式如何用matlab编辑

分子为第n个子带中感知系数的个数，分母为第n个子带中系数的总个数。分子分母我均已求得，大小均为<512*768 double>的矩阵，求得的直方图归一化，其纵坐标代表各子带感知系数所占的比例。谢谢回答。
最后可以得到一幅具有L个柄（L为总共选取的子带个数）的归一化直方图，其纵坐标表示各自带中感知系数所占的比例。

您好， 上面的图片即为详细内容。意思就是任意一个子带，记为n子带(n属于L，总共选取L个子带)。实际上这里的N值您可以忽略。分子分母的矩阵大小均为512*768,其中分子该选取子带中是经过处理的系数矩阵，其中的一部分数字已经处理为0，感知系数的个数即为非零个数，分母为该选取子带中未经处理的系数矩阵。





Answer 1

你的问题描述不够清楚。

我是这样理解你的问题的：那个512*768矩阵的第n行是你说的第n个子带中的768个系数的系数值，第n行中有Ct(n)个感知系数，这Ct(n)个感知系数的值相同，假设这个值是t(n)，你现在想统计每行的P(n)

记你说的512*768的矩阵为A，其他变量名如我上面所说，假设你说的n是矩阵的第n行(如果是第n列，将程序中的A转置一下即可)，
编程如下：

for i=1:512
Temp=tabulate(A(i,:));
%tabulate的用法见链接：
%http://www.mathworks.cn/cn/help/stats/tabulate.html?searchHighlight=tabulate
S(i)=Temp(Temp(:,1)==t(n),3) %这个语句得到t(n)在A的第n行出现的频率
end

程序结束，S中储存了各子带系数中感知系数所占的比例。

如对程序不明或我对你的问题理解有误，欢迎追问补充

追问

我在问题那一栏补充了一下解释。
不是您说的512*768矩阵的第n行为第n个子带中的768个系数的系数值。
我的理解是，这个问题是对任一选取子带，分别求得分子分母系数矩阵中的非零总数，相除就得到相应子带的感知系数所占的比例。昨天已计算出各个系数均在（0,1）之间。这是不是就是论文所说的归一化？那个图是怎么用代码实现的？
谢谢您的耐心与细心，期待着您的回复！

追答

你是说，总共有L个子带，每个子带的分母系数对应512*768矩阵(记为A)中的非零元素，分母系数也对应一个512*768矩阵(记为B)中的非零个数，你想对L个子带，求（A中非零个数/B中非零个数），求出来有L个。另外A、B矩阵中元素值在(0,1)之间？

那你可以这样：
把所有子带的A矩阵存储在A(:,:,L)这个三维数组中，A(:,:,n)代表第n个子带的A
同样，B作类似操作，存储在B(:,:,L)中
程序如下：

for n=1:L
temp=[];
temp=reshape(A(:,:,n),1); %将A拉成一个列向量，方便统计
temp=tabulate(temp); %tabulate用法见我第一次的回答
S(n)=temp(temp(:,1)==0,3); %提取A中0元素的频率
S(n)=1-S(n); %得到A中非零元素频率

temp=[];
temp=reshape(B(:,:,n),1);
temp=tabulate(temp);
S(n)=S(n)/(1-temp(temp(:,1)==0,3)); %将A中非零元素频率与B中非零元素频率相比，即感知系数个数与总系数个数比值
end

得到一个L维行向量S，分别对应L个子带的感知系数所占比例。

当然程序还可以有其他设计思路，比如将A,B矩阵先上取整，这样非零元素就全部变成1，然后直接统计1的频率。

追问

谢谢您的回答！我上面说错了，我说求出来的各个子带的感知系数对应比例全在（0,1）之间。不过应该没太多差别，我是用find函数求的。应该效果都是一样的。谢谢您！