直方图规定化

一、直方图规定化的原理和步骤

---

二、具体实现

clear all;
close all;
clc;
%% 计算给定的直方图累积分布函数
r=127;                              
x=-r:r+1;
sigma=20;
y1=exp(-((x-80).^2)/(2*sigma^2));
y2=exp(-((x+80).^2)/(2*sigma^2));
y=y1+y2;                        %双峰高斯函数，任意函数都可以
y=y/sum(y);                   %归一化，使函数符合概率分布的sum(y)==1这样一个规律
G=zeros(1,256);               %函数的累积直方图
G=cumsum(y);

%%  计算待处理的累积分布函数
img=imread('pout.tif');
[m n]=size(img);
hist=imhist(img);       %imhist能返回每个灰度级的像素值,返回的是列向量
p=hist'/(m*n)