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 10 Python
Python操作MySQL简单实现方法
Jan 26 Python
Python发送email的3种方法
Apr 28 Python
python开发之函数定义实例分析
Nov 12 Python
Python实现包含min函数的栈
Apr 29 Python
使用Python的turtle模块画图的方法
Nov 15 Python
Python 调用PIL库失败的解决方法
Jan 08 Python
Python实现的大数据分析操作系统日志功能示例
Feb 11 Python
python地震数据可视化详解
Jun 18 Python
Python中random模块常用方法的使用教程
Oct 04 Python
利用python调用摄像头的实例分析
Jun 07 Python
python单向链表实例详解
May 25 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函数
2010/02/16 PHP
整理的9个实用的PHP库简介和下载
2010/11/09 PHP
通过Javascript将数据导出到外部Excel文档的函数代码
2012/06/15 Javascript
基于pthread_create,readlink,getpid等函数的学习与总结
2013/07/17 Javascript
jQuery获得IE版本不准确webbrowser的解决方法
2014/02/23 Javascript
jQuery通过Ajax返回JSON数据
2015/04/28 Javascript
JavaScript实现垂直滚动条效果
2017/01/18 Javascript
javascript ES6 新增了let命令使用介绍
2017/07/07 Javascript
JS交互点击WKWebView中的图片实现预览效果
2018/01/05 Javascript
express如何使用session与cookie的方法
2018/01/30 Javascript
JavaScript中利用Array filter() 方法压缩稀疏数组
2018/02/24 Javascript
vue实现裁切图片同时实现放大、缩小、旋转功能
2018/03/02 Javascript
解决vue js IOS H5focus无法自动弹出键盘的问题
2018/08/30 Javascript
浅谈KOA2 Restful方式路由初探
2019/03/14 Javascript
原生JS实现随机点名项目的实例代码
2019/04/30 Javascript
JS如何实现网站中PC端和手机端自动识别并跳转对应的代码
2020/01/08 Javascript
JS实现可视化音频效果的实例代码
2020/01/16 Javascript
Python continue语句用法实例
2014/03/11 Python
Python的迭代器和生成器使用实例
2015/01/14 Python
基于Django的ModelForm组件(详解)
2017/12/07 Python
Python实现基于KNN算法的笔迹识别功能详解
2018/07/09 Python
解决Pandas的DataFrame输出截断和省略的问题
2019/02/08 Python
django2.0扩展用户字段示例
2019/02/13 Python
python实现对列表中的元素进行倒序打印
2019/11/23 Python
如何使用Python发送HTML格式的邮件
2020/02/11 Python
python生成13位或16位时间戳以及反向解析时间戳的实例
2020/03/03 Python
Python搭建Keras CNN模型破解网站验证码的实现
2020/04/07 Python
Python实现AI换脸功能
2020/04/10 Python
Python中的全局变量如何理解
2020/06/04 Python
详解python安装matplotlib库三种失败情况
2020/07/28 Python
PHP两种查询函数array/row的区别
2013/06/03 面试题
*p++ 自增p 还是p所指向的变量
2016/07/16 面试题
药学专业大学生个人的自我评价
2013/11/04 职场文书
主题实践活动总结
2014/05/08 职场文书
学校端午节活动方案
2014/08/23 职场文书
redis实现排行榜功能
2021/05/24 Redis