我一直在尝试将这个等式实现到matlab中：
http://www.mathworks.com/matlabcentral/answers/uploaded_files/46747/Capture.PNG
Omega是图像u1的某个已知区域。
和p1,i是：
http://www.mathworks.com/matlabcentral/answers/uploaded_files/46748/Capture1.PNG
注意：u1是图像（按列排序），应该在此图像的每个像素中计算p1,i。
现在，我以这种方式计算了p1,i：
[f1,xi1] = ksdensity(u1(:),1:255);
p1_u1 = reshape(f1(floor(u1(:))+1),M,N);
现在我的问题是计算前一个等式。
我尝试过使用for循环，但它会永远消失......
还有其他建议吗？也许有一种方法可以使用ksdensity并更改积分内的值？
谢谢！