matlab 计算灰度图像的一阶矩,二阶矩,三阶矩实例


Posted in Python onApril 22, 2020

一阶矩,定义了每个颜色分量的平均强度

matlab 计算灰度图像的一阶矩,二阶矩,三阶矩实例

二阶矩,反映待测区域颜色方差,即不均匀性

matlab 计算灰度图像的一阶矩,二阶矩,三阶矩实例

三阶矩,定义了颜色分量的偏斜度,即颜色的不对称性

matlab 计算灰度图像的一阶矩,二阶矩,三阶矩实例

close all;clear all;clc;
J = imread('lena.jpg');
K = imadjust(J,[70/255 160/255],[]);
figure;
subplot(121),imshow(J);
subplot(122),imshow(K);
[m,n] = size(J);
mm = round(m/2);
mn = round(n/2);
[p,q] = size(K);
pp = round(p/2);
qq = round(q/2);
J = double(J);
K = double(K);
colorsum = 0.0;
Javg = mean2(J) %求原图像一阶矩
Kavg = mean2(K) %求增强对比度后的图像一阶矩
Jstd = std(std(J)) %求原图像的二阶矩,因为一次std函数表示按列求标准差,两次std表示求整个矩阵的标准差
Kstd = std(std(K)) %求增强对比度后的图像二阶矩
for i=1:mm
 for j=1:mn
  colorsum = colorsum+(J(i,j)-Javg)^3;
 end
end
Jske = (colorsum/(mm*mn))^(1/3) %求原图像的三阶矩
colorsum = 0.0;
for i=1:pp
 for j=1:qq
  colorsum = colorsum + (J(i,j)-Kavg)^3;
 end
end
Kske = (colorsum/(pp*qq))^(1/3) %求增强对比度后的图像三阶矩

部分函数说明:

mean2(A) : 求矩阵A的均值

std(x,flag,dim): 求x的标准偏差

std(x,0,1) : 0表示求标准差时除n-1,1表示按列划分

std(x,1,2) : 1表示求标准差时除n,2表示按行划分

补充知识:图像的重心和二阶矩

图像的重心

图像实际上就是个矩阵,每个位置的元素就是该处的像素。 这里碰到了求图像重心的问题,特此总结:

计算公式:

matlab 计算灰度图像的一阶矩,二阶矩,三阶矩实例

其中(xi,yi)是像素点的坐标,pi是该点的像素值。

以下是利用matlab求图像重心

I = imread(‘1.jpg');
I = rgb2gray(I);
imshow(I);
I = double(I);
[rows,cols] = size(I);
x = ones(rows,1)*[1:cols];
y = [1:rows]'*ones(1,cols);
area = sum(sum(I));
meanx = sum(sum(I.*x))/area;
meany = sum(sum(I.*y))/area;
hold on;
plot(meanx,meany,'r+'); %十字标出重心位置

图像的二阶矩

我们这里只讨论二阶矩的问题。

二阶矩最终是形成了一个二阶矩阵,如下:

matlab 计算灰度图像的一阶矩,二阶矩,三阶矩实例

计算方法:

matlab 计算灰度图像的一阶矩,二阶矩,三阶矩实例

其中(r0,c0)是重心坐标。

以上这篇matlab 计算灰度图像的一阶矩,二阶矩,三阶矩实例就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
详细介绍Python的鸭子类型
Sep 12 Python
python操作mysql数据库
Mar 05 Python
关于pip的安装,更新,卸载模块以及使用方法(详解)
May 19 Python
使用python 和 lint 删除项目无用资源的方法
Dec 20 Python
Python绘制七段数码管实例代码
Dec 20 Python
Python实现识别手写数字 Python图片读入与处理
Mar 23 Python
python数字图像处理之骨架提取与分水岭算法
Apr 27 Python
python里运用私有属性和方法总结
Jul 08 Python
Django使用中间件解决前后端同源策略问题
Sep 02 Python
Tensorflow的常用矩阵生成方式
Jan 04 Python
如何使用python socket模块实现简单的文件下载
Sep 04 Python
Python 数据科学 Matplotlib图库详解
Jul 07 Python
python根据完整路径获得盘名/路径名/文件名/文件扩展名的方法
Apr 22 #Python
matlab中二维插值函数interp2的使用详解
Apr 22 #Python
python 一维二维插值实例
Apr 22 #Python
Numpy一维线性插值函数的用法
Apr 22 #Python
python数据处理——对pandas进行数据变频或插值实例
Apr 22 #Python
利用4行Python代码监测每一行程序的运行时间和空间消耗
Apr 22 #Python
Jupyter Notebook折叠输出的内容实例
Apr 22 #Python
You might like
PHP 获取MySQL数据库里所有表的实现代码
2011/07/13 PHP
Tips 带三角可关闭的文字提示
2010/10/06 Javascript
Jquery实现简单的动画效果代码
2012/03/18 Javascript
jquery 元素控制(追加元素/追加内容)介绍及应用
2013/04/21 Javascript
jquery+javascript编写国籍控件
2015/02/12 Javascript
JS实现的左侧竖向滑动菜单效果代码
2015/10/19 Javascript
基于jQuery实现网页打印功能
2015/12/01 Javascript
JavaScript 对象字面量讲解
2016/06/06 Javascript
不间断循环滚动效果的实例代码(必看篇)
2016/10/08 Javascript
js 输入框 正则表达式(菜鸟必看教程)
2017/02/19 Javascript
使用bootstraptable插件实现表格记录的查询、分页、排序操作
2017/08/06 Javascript
微信小程序swiper组件用法实例分析【附源码下载】
2017/12/07 Javascript
Vue组件中slot的用法
2018/01/30 Javascript
react 创建单例组件的方法
2018/04/26 Javascript
Vuex的API文档说明详解
2020/02/05 Javascript
vue用ant design中table表格,点击某行时触发的事件操作
2020/10/28 Javascript
[51:11]2014 DOTA2国际邀请赛中国区预选赛5.21 LGD-CDEC VS DT
2014/05/22 DOTA
[00:47]DOTA2荣耀之路6:天火,天火!
2018/05/30 DOTA
[52:44]VGJ.T vs infamous Supermajor小组赛D组败者组第一轮 BO3 第一场 6.3
2018/06/04 DOTA
[46:58]完美世界DOTA2联赛PWL S3 Forest vs LBZS 第一场 12.17
2020/12/19 DOTA
Python获取邮件地址的方法
2015/07/10 Python
深入理解 Python 中的多线程 新手必看
2016/11/20 Python
python学习之面向对象【入门初级篇】
2017/01/21 Python
python使用多进程的实例详解
2018/09/19 Python
python networkx 包绘制复杂网络关系图的实现
2019/07/10 Python
解决python replace函数替换无效问题
2020/01/18 Python
python 线性回归分析模型检验标准--拟合优度详解
2020/02/24 Python
Python flask框架如何显示图像到web页面
2020/06/03 Python
纯CSS3实现Material Design效果
2017/03/09 HTML / CSS
亚马逊中国官方网站:amazon.cn
2017/05/25 全球购物
英国领先的在线高尔夫商店:Gamola Golf
2019/11/16 全球购物
C++如何引用一个已经定义过的全局变量
2014/08/25 面试题
中文专业学生自我评价范文
2014/02/06 职场文书
异地年检委托书范本
2014/09/24 职场文书
死亡赔偿协议书
2015/01/28 职场文书
java开发双人五子棋游戏
2022/05/06 Java/Android