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 and、or以及and-or语法总结
Apr 14 Python
浅析Python多线程下的变量问题
Apr 28 Python
基于python实现的抓取腾讯视频所有电影的爬虫
Apr 22 Python
动感网页相册 python编写简单文件夹内图片浏览工具
Aug 17 Python
Python编程之微信推送模板消息功能示例
Aug 21 Python
Python自定义简单图轴简单实例
Jan 08 Python
pandas数据框,统计某列数据对应的个数方法
Apr 11 Python
python清除字符串前后空格函数的方法
Oct 21 Python
pygame游戏之旅 游戏中添加显示文字
Nov 20 Python
GDAL 矢量属性数据修改方式(python)
Mar 10 Python
Django 拼接两个queryset 或是两个不可以相加的对象实例
Mar 28 Python
详解python程序中的多任务
Sep 16 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对接微信公众平台消息接口开发流程教程
2014/03/25 PHP
PHP SOCKET编程详解
2015/05/22 PHP
PHP排序算法之基数排序(Radix Sort)实例详解
2018/04/21 PHP
jscript之List Excel Color Values
2007/06/13 Javascript
用正则表达式 动态创建/增加css style script 兼容IE firefox
2009/03/10 Javascript
autoIMG 基于jquery的图片自适应插件代码
2011/03/12 Javascript
js实现动态添加、删除行、onkeyup表格求和示例
2013/08/18 Javascript
不同Jquery版本引发的问题解决
2013/10/14 Javascript
利用PM2部署node.js项目的方法教程
2017/05/10 Javascript
Nodejs中Express 常用中间件 body-parser 实现解析
2017/05/22 NodeJs
web前端开发中常见的多列布局解决方案整理(一定要看)
2017/10/15 Javascript
基于vue cli 通过命令行传参实现多环境配置
2018/07/12 Javascript
jQuery实现的页面弹幕效果【测试可用】
2018/08/17 jQuery
微信小程序城市选择及搜索功能的方法
2019/03/22 Javascript
JS中的防抖与节流及作用详解
2019/04/01 Javascript
黑科技 Python脚本帮你找出微信上删除你好友的人
2016/01/07 Python
python3.6 实现AES加密的示例(pyCryptodome)
2018/01/10 Python
python使用tkinter实现简单计算器
2018/01/30 Python
python实现逐个读取txt字符并修改
2018/12/24 Python
基于 Django 的手机管理系统实现过程详解
2019/08/16 Python
使用PyCharm进行远程开发和调试的实现
2019/11/04 Python
Python解析多帧dicom数据详解
2020/01/13 Python
python输出第n个默尼森数的实现示例
2020/03/08 Python
使用Python中tkinter库简单gui界面制作及打包成exe的操作方法(二)
2020/10/12 Python
纯HTML5+CSS3制作生日蛋糕(代码易懂)
2016/11/16 HTML / CSS
英国香水店:The Perfume Shop
2017/03/27 全球购物
美国农场鲜花速递:The Bouqs
2018/07/13 全球购物
哈曼俄罗斯官方网上商店:Harman.club
2020/07/24 全球购物
小学毕业感言150字
2014/02/05 职场文书
小学中秋节活动方案
2014/02/06 职场文书
作风建设演讲稿
2014/05/23 职场文书
软件工程毕业生自荐信
2014/07/04 职场文书
大学生预备党员自我评价
2015/03/04 职场文书
创业计划之特色精品店
2019/08/12 职场文书
2019大学生暑期实习心得总结
2019/08/21 职场文书
MySQL令人大跌眼镜的隐式转换
2021/08/23 MySQL