matlab 計(jì)算灰度圖像的一階矩,二階矩,三階矩實(shí)例
一階矩，定義了每個(gè)顏色分量的平均強(qiáng)度
二階矩，反映待測(cè)區(qū)域顏色方差，即不均勻性
三階矩，定義了顏色分量的偏斜度，即顏色的不對(duì)稱性
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) %求增強(qiáng)對(duì)比度后的圖像一階矩
Jstd = std(std(J)) %求原圖像的二階矩，因?yàn)橐淮蝧td函數(shù)表示按列求標(biāo)準(zhǔn)差，兩次std表示求整個(gè)矩陣的標(biāo)準(zhǔn)差
Kstd = std(std(K)) %求增強(qiáng)對(duì)比度后的圖像二階矩
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) %求增強(qiáng)對(duì)比度后的圖像三階矩
部分函數(shù)說(shuō)明：
mean2(A) : 求矩陣A的均值
std（x,flag,dim）: 求x的標(biāo)準(zhǔn)偏差
std(x,0,1) : 0表示求標(biāo)準(zhǔn)差時(shí)除n-1，1表示按列劃分
std(x,1,2) : 1表示求標(biāo)準(zhǔn)差時(shí)除n，2表示按行劃分
補(bǔ)充知識(shí)：圖像的重心和二階矩
圖像的重心
圖像實(shí)際上就是個(gè)矩陣，每個(gè)位置的元素就是該處的像素。 這里碰到了求圖像重心的問(wèn)題，特此總結(jié)：
計(jì)算公式：
其中(xi,yi)是像素點(diǎn)的坐標(biāo)，pi是該點(diǎn)的像素值。
以下是利用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+'); %十字標(biāo)出重心位置
圖像的二階矩
我們這里只討論二階矩的問(wèn)題。
二階矩最終是形成了一個(gè)二階矩陣，如下：
計(jì)算方法：
其中(r0,c0)是重心坐標(biāo)。
以上這篇matlab 計(jì)算灰度圖像的一階矩,二階矩,三階矩實(shí)例就是小編分享給大家的全部?jī)?nèi)容了，希望能給大家一個(gè)參考，也希望大家多多支持腳本之家。
相關(guān)文章
Django csrf 驗(yàn)證問(wèn)題的實(shí)現(xiàn)
csrf是通過(guò)偽裝來(lái)自受信任用戶的請(qǐng)求來(lái)利用受信任的網(wǎng)站。這篇文章主要介紹了Django csrf 驗(yàn)證問(wèn)題的實(shí)現(xiàn)，小編覺(jué)得挺不錯(cuò)的，現(xiàn)在分享給大家，也給大家做個(gè)參考。一起跟隨小編過(guò)來(lái)看看吧2018-10-10
今天小編就為大家分享一篇python 讀取修改pcap包的例子，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2019-07-07
python制作小說(shuō)爬蟲(chóng)實(shí)錄
本文給大家介紹的是作者所寫(xiě)的第一個(gè)爬蟲(chóng)程序的全過(guò)程，從構(gòu)思到思路到程序的編寫(xiě)，非常的細(xì)致，有需要的小伙伴可以參考下2017-08-08
Python實(shí)現(xiàn)常用文本內(nèi)容提取
在日常工作和學(xué)習(xí)中,我們經(jīng)常需要從PDF、Word文檔中提取文本,本文將介紹如何使用Python編寫(xiě)一個(gè)文本內(nèi)容提取工具,有需要的小伙伴可以參考下2025-03-03
Python中index()函數(shù)與find()函數(shù)的區(qū)別詳解
這篇文章主要介紹了Python中index()函數(shù)與find()函數(shù)的區(qū)別詳解,Python index()方法檢測(cè)字符串中是否包含子字符串 str ,如果指定beg開(kāi)始和end結(jié)束范圍,則檢查是否包含在指定范圍內(nèi),需要的朋友可以參考下2023-08-08
這段時(shí)間老是“不務(wù)正業(yè)”的搞一些東西玩。之前的貪吃蛇，俄羅斯方塊激發(fā)了我研究游戲算法的興趣。經(jīng)過(guò)1個(gè)星期的構(gòu)思，連連看的連接算法終于出爐了。再過(guò)一段時(shí)間就基于這個(gè)算法使用JavaScript推出網(wǎng)頁(yè)版的連連看。下面是說(shuō)明及代碼。2008-11-11
每天上班最痛苦的事情就是早起早起早起！這是大部分上班族的痛苦，但是不上班又是不可能的啦，因?yàn)槎际菫榱烁沐X。本文用Python制作了三款有趣的鬧鐘，感興趣的可以學(xué)習(xí)一下2022-05-05
教你利用Selenium+python自動(dòng)化來(lái)解決pip使用異常
今天帶大家來(lái)學(xué)習(xí)怎么利用Selenium+python自動(dòng)化解決pip使用異常,文中介紹的非常詳細(xì),對(duì)正在學(xué)習(xí)python的小伙伴們有很大的幫助,需要的朋友可以參考下2021-05-05
Python生成隨機(jī)數(shù)的方法有很多，譬如random模塊、NumPy庫(kù)、secrets模塊等。本文為大家整理了五個(gè)生成隨機(jī)數(shù)的方法（應(yīng)該是最全的了吧），希望對(duì)大家有所幫助2023-02-02
python中threading.Semaphore和threading.Lock的具體使用
python中的多線程是一個(gè)非常重要的知識(shí)點(diǎn),本文主要介紹了python中threading.Semaphore和threading.Lock的具體使用,具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下2023-08-08