matlab直方图归一化_matlab归一化函数normalize

matlab直方图归一化_matlab归一化函数normalize直方图规定化直方图均衡化的优点是能自动增强整个图像的对比度,但它的具体增强效果不易控制,处理的结果总是得到全局的均衡化的直方图.实际工作中,有时需要变换直方图使之成为某个特定的形状,从而有选择地增强某个灰度值范围内的对比度,这时可采用比较灵活的直方图规定化方法.直方图规定化增强处理的步骤如下:令Pr(r)和Pz(z)分别为原始图像和期望图像的灰度概率密度函数。如果对原始图像和期望图像均作直

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE稳定放心使用

直方图规定化

直方图均衡化的优点是能自动增强整个图像的对比度,但它的具体增强效果不易控制,处理的结果总是得到全局的均衡化的直方图.实际工作中,有时需要变换直方图使之成为某个特定的形状,从而有选择地增强某个灰度值范围内的对比度,这时可采用比较灵活的直方图规定化方法.直方图规定化增强处理的步骤如下:

令Pr(r)和Pz(z)分别为原始图像和期望图像的灰度概率密度函数。如果对原始图像和期望图像均作直方图均衡化处理,应有:

matlab直方图归一化_matlab归一化函数normalize

matlab直方图归一化_matlab归一化函数normalize

matlab直方图归一化_matlab归一化函数normalize

由于都是进行均衡化处理,处理后的原图像概率密度函数Ps(S)及理想图像概率密度函数PV(V)是相等的。

于是,我们可以用变换后的原始图像灰度级S代替(2)式中的V。

即:Z = G – 1(S) 

这时的灰度级Z 便是所希望的图像的灰度级。

此外,利用(1)与(3)式还可得到组合变换函数

                 Z = G – 1[T(r)]

对连续图像,重要的是给出逆变换解析式。对离散图像而言,有

matlab直方图归一化_matlab归一化函数normalize

matlab直方图归一化_matlab归一化函数normalize

matlab直方图归一化_matlab归一化函数normalize


clear all;
close all;
clc;

I = imread('2-t.bmp');
L=256;            
NK=zeros(L,1);               %存储原图像直方图数据
Rk_CDF_Normal=zeros(L,1);    %存储直方图规定化后的图像的直方图
Rk_pre=zeros(L,1);           %存储原图像累积直方图数据
[row,col]=size(I);

%%
%获取像素灰度级
for i = 1:row
    for j = 1:col
        num = double(I(i,j))+1;
        NK(num) = NK(num)+1;
    end
end

%计算直方图概率估计
Ps = NK./numel(I);%存储原图像直方图概率数据
%计算累积直方图
for level=1:L
    if level==1
        Rk_pre(level)=Ps(level);
    else
        Rk_pre(level)=Rk_pre(level-1)+Ps(level);
    end
end

%%
%规定化直方图,在这里要得到2-s图像的灰度直方图
H = imread('2-s.jpg');
H = rgb2gray(H);
nk_normal = zeros(L,1);
Rk_normal = zeros(L,1);
[row,col] = size(H);

%获取像素灰度级
for i = 1:row
    for j = 1:col
        num_normal = double(I(i,j))+1;
        nk_normal(num_normal) = nk_normal(num_normal)+1;
    end
end

%计算直方图概率估计
Ps_normal = nk_normal./numel(H);
for level=1:L
    if level==1
        Rk_normal(level)=Ps_normal(level);
    else
        Rk_normal(level)=Rk_normal(level-1)+Ps_normal(level);
    end
end

%计算规定化累积直方图
for level=1:L
    if level==1
        Rk_CDF_Normal(level)=Rk_normal(level);
    else
        Rk_CDF_Normal(level)=Rk_CDF_Normal(level-1)+Rk_normal(level);
    end
end

%%
%计算原图像与目标图像累计直方图数值的差的绝对值
double ScMin=zeros(L,L);
for y=1:L
    for x=1:L
        ScMin(x,y)=abs(Rk_pre(y)'-Rk_CDF_Normal(x)');
    end
end
%建立映射
HisM=zeros(L:1);
for level_x=1:L
    min = 0;
    minV=ScMin(1,level_x);
    for level_y=1:L
        if(minV>ScMin(level_y,level_x))
            minV=ScMin(level_y,level_x);
            min = level_y;
        end
    end
    HisM(level_x)= min;
end

%将原图像的每个像素灰度转换为直方图均衡化后的灰度
[row,col]=size(I);
New=I;
for x = 1:row
    for y = 1:col
        Num = double(I(x,y))+1;
        if Num==L
            New(x,y)=HisM(L);
        end
    end
end
figure;imshow(New),title('规定化后图像');

matlab直方图归一化_matlab归一化函数normalize

matlab直方图归一化_matlab归一化函数normalize

matlab直方图归一化_matlab归一化函数normalize

  

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/180271.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • 扫清盲点,如何正确的从HttpClient 3.x系统升级到HttpClient 4.x

    扫清盲点,如何正确的从HttpClient 3.x系统升级到HttpClient 4.x如果周期比较长的项目,或者这个项目开发人员换过了好几拨人,很有可能出现一些奇怪的问题,比如一个项目中出现了多种Spring注入bean的方式,不同版本的jar冲突等等爬虫项目有的时候更是过犹不及,拿模拟登陆来说,开发人员的迭代,每个人的风格和技术各不相同,模拟登陆的方式也是五花八门,早在之前看到过一个项目的源码,其中使用HttpClient也是各种风格,虽然官方已经强烈建议使用HttpClie…

  • es与数据库的同步方案

    es与数据库的同步方案一、双写模式我们采取MySQL作为主要的数据存储,利用MySQL的事务特性维护数据一致性,使用ElasticSearch进行数据汇集和查询,此时es与数据库的同步方案就尤为重要。保证es与数据库的同步方案:1、首先添加商品入数据库,添加商品成功后,商品入ES,若入ES失败,将失败的商品ID放入redis的缓存队列(或MQ),且失败的商品ID入log文件(若出现redis挂掉,可从日志中取异…

  • 格拉姆矩阵(Gram matrix)详细解读

    格拉姆矩阵(Gram matrix)详细解读目录基础知识-向量的内积Grammatrix介绍Grammatrix的应用-风格迁移一、基础知识-向量的内积1.1向量的内积定义:也叫向量的点乘,对两个向量执行内积运算,就是对这两个向

  • 阿里云服务器使用freessl配置免费证书Nginx

    阿里云服务器使用freessl配置免费证书Nginx环境:阿里云服务器购买的域名服务器:linux+nginxSSL证书:FreeSSL申请的免费证书步骤1、申请ssl证书具体参考二哥的“五分钟搞定HTTPS配置,二哥手把手教”https://blog.csdn.net/qing_gee/article/details/90031376博客,这位大佬写的很详细了2、配置阿里云服务器2.1、上传证书登陆阿里云控制台,搜索“ssl证书应用安全”,上传原有证书,注意一定要将证书转换为pem格式2.2、开启服务器443.

  • Layui分页_pagehelper分页使用

    Layui分页_pagehelper分页使用本文介绍了LayUI分页,LayUI动态分页,LayUIlaypage分页,LayUIlaypage刷新当前页,分享给大家,具体如下:效果图:一、引用js依赖主要是jquery-1.11.3.min.js和layui.all.js,json2.js用来做json对象转换的二、js分页方法封装(分页使用模板laytpl)1、模板渲染/***分页模板的渲染方法*@paramtemp…

    2022年10月28日
  • datatrip 激活码[在线序列号][通俗易懂]

    datatrip 激活码[在线序列号],https://javaforall.cn/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号