<h2 id="t_ad4d75200102v9ms" class="titName SG_txta">[转载]Matlab中实现灰度共生矩阵</h2>
<span class="time SG_txtc">(2014-12-08 20:18:02)</span><div class="turnBoxzz"><a href="javascript:;" class="SG_aBtn SG_aBtn_ico SG_turn" action-type="reblog" action-data="{srcBlog:0, blogId:'7c2012a90100y5yo'}"><cite><img class="SG_icon SG_icon111" src="http://simg.sinajs.cn/blog7style/images/common/sg_trans.gif" width="15" height="15" align="absmiddle">转载<em class="arrow">▼</em></cite></a></div> </div>
<div class="articalTag" id="sina_keyword_ad_area">
<table>
<tbody><tr>
<td class="blog_tag">
<script>
var $tag='转载';
var $tag_code='b60b8621e290ec1937052e00bbb0002d';
var $r_quote_bligid='7c2012a90100y5yo';
var $worldcup='0';
var $worldcupball='0';
</script>
<span class="SG_txtb">标签:</span>
<h3><a href="http://search.sina.com.cn/?c=blog&q=%D7%AA%D4%D8&by=tag" target="_blank">转载</a></h3>
</td>
<td class="blog_class">
<span class="SG_txtb">分类:</span>
<a target="_blank" href="http://blog.sina.com.cn/s/articlelist_2907534624_2_1.html">科研</a>
</td>
</tr>
</tbody></table>
</div>
<!-- 正文开始 -->
<div id="sina_keyword_ad_area2" class="articalContent newfont_family">
谢谢分享~~~<div class="blogzz_abstract borderc" style="padding-top:15px;margin:20px 0; border:none; border-top:1px dotted #ccc;"><div class="blogzz_ainfo" style="margin-bottom:12px;"><span style="margin-right:25px;"><strong>原文地址:</strong><a target="_blank" href="http://blog.sina.com.cn/s/blog_7c2012a90100y5yo.html" title="Matlab中实现灰度共生矩阵">Matlab中实现灰度共生矩阵</a></span><span><strong>作者:</strong><a href="http://blog.sina.com.cn/u/2082476713" title="louis" target="_blank">louis</a></span></div><div class="blogzz_acon">%**************************************************************************<br>
%
图像检索——纹理特征
%基于共生矩阵纹理特征提取,d=1,θ=0°,45°,90°,135°共四个矩阵
%所用图像灰度级均为256
%参考《基于颜色空间和纹理特征的图像检索》
%function : T=Texture(Image)
%Image
输入图像数据
%T
: 返回八维纹理特征行向量
%**************************************************************************
% function T = Texture(Image)
Gray = imread(‘d:result5.bmp’);
[M,N,O] = size(Gray);
M = 128;
N = 128;
%--------------------------------------------------------------------------
%1.将各颜色分量转化为灰度
%--------------------------------------------------------------------------
% Gray =
double(0.3Image(:,:,1)+0.59Image(:,:,2)+0.11*Image(:,:,3));
%--------------------------------------------------------------------------
%2.为了减少计算量,对原始图像灰度级压缩,将Gray量化成16级
%--------------------------------------------------------------------------
for i = 1:M
1:N
for n = 1:256/16
if
(n-1)*16<=Gray(i,j)&Gray(i,j)<=(n-1)*16+15
Gray(i,j) = n-1;
end
end
end
%--------------------------------------------------------------------------
%3.计算四个共生矩阵P,取距离为1,角度分别为0,45,90,135
%--------------------------------------------------------------------------
P = zeros(16,16,4);
for m = 1:16
1:16
for i = 1:M
for j = 1:N
if
j<N&Gray(i,j)==m-1&Gray(i,j+1)==n-1
P(m,n,1) = P(m,n,1)+1;
P(n,m,1) = P(m,n,1);
end
if
i>1&j<N&Gray(i,j)==m-1&Gray(i-1,j+1)==n-1
P(m,n,2) = P(m,n,2)+1;
P(n,m,2) = P(m,n,2);
end
if
i<M&Gray(i,j)==m-1&Gray(i+1,j)==n-1
P(m,n,3) = P(m,n,3)+1;
P(n,m,3) = P(m,n,3);
end
if
i<M&j<N&Gray(i,j)==m-1&Gray(i+1,j+1)==n-1
P(m,n,4) = P(m,n,4)+1;
P(n,m,4) = P(m,n,4);
end
end
end
if m==n
P(m,n,:) = P(m,n,:)*2;
end
end
%%---------------------------------------------------------
% 对共生矩阵归一化
%%---------------------------------------------------------
for n = 1:4
P(:,:,n)/sum(sum(P(:,:,n)));
end
%--------------------------------------------------------------------------
%4.对共生矩阵计算能量、熵、惯性矩、相关4个纹理参数
%--------------------------------------------------------------------------
H = zeros(1,4);
I = H;
Ux =
H;
Uy = H;
deltaX= H;
C =H;
for n = 1:4
sum(sum(P(:,:,n).^2)); %%能量
1:16
for j = 1:16
if P(i,j,n)~=0
H(n) = -P(i,j,n)log(P(i,j,n))+H(n); %%熵
end
I(n) = (i-j)^2P(i,j,n)+I(n);
Ux(n) = iP(i,j,n)+Ux(n); %相关性中μx
Uy(n) = jP(i,j,n)+Uy(n); %相关性中μy
end
end
for n = 1:4
1:16
for j = 1:16
deltaX(n) = (i-Ux(n))^2P(i,j,n)+deltaX(n); %相关性中σx
deltaY(n) = (j-Uy(n))^2P(i,j,n)+deltaY(n); %相关性中σy
C(n) =
ijP(i,j,n)+C(n);
end
(C(n)-Ux(n)*Uy(n))/deltaX(n)/deltaY(n);
%相关性
end
%--------------------------------------------------------------------------
%求能量、熵、惯性矩、相关的均值和标准差作为最终8维纹理特征
%--------------------------------------------------------------------------
a1 = mean(E)
b1 = sqrt(cov(E))
a2 = mean(H)
b2 = sqrt(cov(H))
a3 = mean(I)
b3 = sqrt(cov(I))
a4 = mean©
b4 = sqrt(cov©)
sprintf(‘0,45,90,135方向上的能量依次为: %f, %f, %f,
%f’,E(1),E(2),E(3),E(4))
sprintf(‘0,45,90,135方向上的熵依次为: %f, %f, %f,
%f’,H(1),H(2),H(3),H(4))
sprintf(‘0,45,90,135方向上的惯性矩依次为: %f, %f, %f,
%f’,I(1),I(2),I(3),I(4))
sprintf(‘0,45,90,135方向上的相关性依次为: %f, %f, %f,
%f’,C(1),C(2),C(3),C(4))
<div class="upBox upBox_add">
<p class="count" id="goldPan-num">0</p>
<p class="link" id="goldPan-give"><img class="SG_icon SG_icon214" src="http://simg.sinajs.cn/blog7style/images/common/sg_trans.gif" width="20" height="16" title="赠金笔" align="absmiddle">赠金笔</p>
</div>
</div>
<div class="clearit"></div>
</div>
<div class="articalInfo">
<!-- 分享到微博 {$t_blog} -->
<div class="IL">
阅读<span id="r_ad4d75200102v9ms" class="SG_txtb">(93)</span><em class="SG_txtb">┊</em>
<a href="#commonComment">评论</a> <span id="c_ad4d75200102v9ms" class="SG_txtb">(0)</span><em class="SG_txtb">┊</em> <a href="javascript:;" onclick="$articleManage('ad4d75200102v9ms',5);return false;">收藏</a><span id="f_ad4d75200102v9ms" class="SG_txtb">(0)</span>
<em class="SG_txtb">┊</em><a href="#" onclick="scope.articel_quote_alert && scope.articel_quote_alert('7c2012a90100y5yo');return false;">转载原文</a> <span id="fn_[转载]Matlab中实现灰度共生矩阵" class="SG_txtb"></span><em class="SG_txtb">┊</em>
<a onclick="return false;" href="javascript:;"><cite id="d1_digg_ad4d75200102v9ms">喜欢</cite></a><a id="d1_digg_down_ad4d75200102v9ms" href="javascript:;"><b>▼</b></a>
<em class="SG_txtb">┊</em><a href="http://blog.sina.com.cn/main_v5/ria/print.html?blog_id=blog_ad4d75200102v9ms" target="_blank">打印</a><em class="SG_txtb">┊</em><a id="q_ad4d75200102v9ms" onclick="report('ad4d75200102v9ms');return false;" href="#">举报</a>
</div>
<div class="IR">
<table>
<tbody><tr><!--
-->
</tr>
</tbody></table>
</div>
</div>
<div class="clearit"></div>
<div class="blogzz_zzlist borderc" id="blog_quote" style="display:none">加载中,请稍候......</div>
<div class="articalfrontback SG_j_linedot1 clearfix" id="new_nextprev_ad4d75200102v9ms">
<div><span class="SG_txtb">前一篇:</span><a href="http://blog.sina.com.cn/s/blog_ad4d75200102v9mr.html">[转载]【转载】MATLAB Toolbox 大全(四)</a></div>
<div><span class="SG_txtb">后一篇:</span><a href="http://blog.sina.com.cn/s/blog_ad4d75200102vb48.html">2014未完成</a></div>
</div>
<div class="clearit"></div>
<div id="loginFollow"></div>
<div class="allComm">
<div class="allCommTit">
<div class="SG_floatL">
<strong>评论</strong>
<span id="commAd_1" style="display:none;">
<span style="margin-left:15px; width:220px; display:inline-block;"><a target="_blank" href="http://blog.sina.com.cn/lm/8/2009/0325/105340.html">重要提示:警惕虚假中奖信息</a></span>
</span>
</div>
<div class="SG_floatR"><a class="CP_a_fuc" href="#post">[<cite>发评论</cite>]</a></div>
</div>
<ul id="article_comment_list" class="SG_cmp_revert"><li><div class="noCommdate"><span class="SG_txtb">做第一个评论者吧! <img class="SG_icon SG_icon134" src="http://simg.sinajs.cn/blog7style/images/common/sg_trans.gif" width="18" height="18" title="" align="absmiddle"><a href="#post">抢沙发>></a></span></div></li></ul>
<div class="clearit"></div>
<div class="myCommPages SG_j_linedot1" style="display: none;">
<div class="SG_page" id="commentPaging" style="display:none;">
<ul class="SG_pages">
</ul>
</div>
<div class="clearit"></div>
</div>
<a name="post"></a>
<div class="writeComm">
<div class="allCommTit">
<div class="SG_floatL">
<strong>发评论</strong>
<span></span>
</div>
<div class="SG_floatR"></div>
</div>
<div class="wrCommTit">
<div class="SG_floatL" id="commentNick" style="">用户5148481065 :</div>
</div>
<div class="formTextarea">
<div style="float:left;" id="commonComment">
<iframe id="postCommentIframe" frameborder="0" style="border:1px solid #C7C7C7;
height:158px;width:448px;maring-top:1px;background-color:white;" src="http://blog.sina.com.cn/main_v5/ria/blank2.html"></iframe>
<textarea id="commentArea" tabindex="1" style="display:none;"></textarea>
</div>
<div id="mobileComment" style="float:left;display:none;">
<textarea id="mbCommentTa" style="width:438px;height:150px;border:1px solid #C7C7C7;line-height:18px;padding:5px;"></textarea>
</div>
<div class="faceblk" id="faceWrap">
<div id="smilesSortShow" class="faceline1">
<div class="facestyle" id="recomm_1542359428673"><a href="#" key="302"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/302-25.gif" alt="小新小浪" title="小新小浪"></a><a href="#" key="308"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/308-25.gif" alt="炮炮兵" title="炮炮兵"></a><a href="#" key="315"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/315-25.gif" alt="张富贵" title="张富贵"></a><a href="#" key="316"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/316-25.gif" alt="旺狗" title="旺狗"></a><a href="#" key="331"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/331-25.gif" alt="悠嘻猴" title="悠嘻猴"></a><a href="#" key="351"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/351-25.gif" alt="酷巴熊" title="酷巴熊"></a></div><span class="SG_more"><a href="#">更多>></a></span><div class="clearit"></div></div>
<ul id="smilesRecommended" class="faceline01"><li><a href="#"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/E___0321EN00SIGT.gif" alt="就不买你" title="就不买你" height="50" width="50"></a></li><li><a href="#"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/E___0320EN00SIGT.gif" alt="股市" title="股市" height="50" width="50"></a></li><li><a href="#"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/E___0319EN00SIGT.gif" alt="发霉" title="发霉" height="50" width="50"></a></li><li><a href="#"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/E___0318EN00SIGT.gif" alt="陈水边" title="陈水边" height="50" width="50"></a></li><li><a href="#"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/E___0317EN00SIGT.gif" alt="裁员" title="裁员" height="50" width="50"></a></li><li><a href="#"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/E___0316EN00SIGT.gif" alt="音乐" title="音乐" height="50" width="50"></a></li><li><a href="#"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/E___0315EN00SIGT.gif" alt="贴你" title="贴你" height="50" width="50"></a></li><li><a href="#"><img src="http://www.sinaimg.cn/uc/myshow/blog/misc/gif/E___0314EN00SIGT.gif" alt="抢车位" title="抢车位" height="50" width="50"></a></li></ul>
</div>
<div class="clearit"></div>
</div>
<div class="formLogin">
<div class="SG_floatL">
<p id="commentlogin" style="display:none;"><span>登录名:</span><input type="text" style="width: 115px;" id="login_name" tabindex="2"> <span>密码:</span><input type="password" style="width: 115px;" id="login_pass" tabindex="3"> <a href="https://login.sina.com.cn/getpass.html" target="_blank">找回密码</a> <a href="https://login.sina.com.cn/signup/signup.php?entry=blog&src=blogicp&srcuid=2907534624" target="_blank">注册</a> <input type="checkbox" id="login_remember"><label for="login_remember" style="display:inline-block;" title="建议在网吧/公用电脑上取消该选项">记住登录状态</label></p><p id="commentloginM" style="display:none;"><span>昵 称:</span><input type="text" style="width: 115px;" id="comment_anonyous" value="新浪网友" tabindex="2" disabled=""></p><p id="quote_comment_p"><!--<input type="checkbox" id="bb"> <label for="bb"><img height="18" align="absmiddle" width="18" title="" src="http://simg.sinajs.cn/blog7style/images/common/sg_trans.gif" class="SG_icon SG_icon110">分享到微博 <img height="15" align="absmiddle" width="15" title="新" src="http://simg.sinajs.cn/blog7style/images/common/sg_trans.gif" class="SG_icon SG_icon11"></label> --></p>
<p id="geetest-box"></p>
</div>
<span style="display: none; color: rgb(153, 153, 153); margin-left: 10px;" id="login_remember_caution"></span>
<!--<div class="SG_floatR" id="anonymity_cont"><input type="checkbox" id="anonymity"/><label for="anonymity">匿名评论</label></div>-->
</div>
<div class="formBtn">
<a href="javascript:;" onclick="return false;" class="SG_aBtn" tabindex="5"><cite id="postcommentid">发评论</cite></a>
<p class="SG_txtc">以上网友发言只代表其个人观点,不代表新浪网的观点或立场。</p>
</div>
</div>
</div>
<div class="clearit"></div>
<div class="articalfrontback articalfrontback2 clearfix">
<div class="SG_floatL"><span class="SG_txtb">< 前一篇</span><a href="http://blog.sina.com.cn/s/blog_ad4d75200102v9mr.html">[转载]【转载】MATLAB Toolbox 大全(四)</a></div>
<div class="SG_floatR"><span class="SG_txtb">后一篇 ></span><a href="http://blog.sina.com.cn/s/blog_ad4d75200102vb48.html">2014未完成</a></div>
</div>
<div class="clearit"></div>
</div>