灰度共生矩阵Matlab中实现

							<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&amp;q=%D7%AA%D4%D8&amp;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

    for j =
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

end



%--------------------------------------------------------------------------

%3.计算四个共生矩阵P,取距离为1,角度分别为0,45,90,135

%--------------------------------------------------------------------------

P = zeros(16,16,4);

for m = 1:16

    for n =
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

end



%%---------------------------------------------------------

% 对共生矩阵归一化

%%---------------------------------------------------------

for n = 1:4

    P(:,:,n) =
P(:,:,n)/sum(sum(P(:,:,n)));

end



%--------------------------------------------------------------------------

%4.对共生矩阵计算能量、熵、惯性矩、相关4个纹理参数

%--------------------------------------------------------------------------

H = zeros(1,4);

I = H;

Ux =
H;     
Uy = H;

deltaX= H;  deltaY = H;

C =H;

for n = 1:4

    E(n) =
sum(sum(P(:,:,n).^2)); %%能量

    for i =
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)^2
P(i,j,n)+I(n);  %%惯性矩

          

           
Ux(n) = iP(i,j,n)+Ux(n); %相关性中μx

           
Uy(n) = j
P(i,j,n)+Uy(n); %相关性中μy

       
end

    end

end

for n = 1:4

    for i =
1:16

       
for j = 1:16

           
deltaX(n) = (i-Ux(n))^2P(i,j,n)+deltaX(n); %相关性中σx

           
deltaY(n) = (j-Uy(n))^2
P(i,j,n)+deltaY(n); %相关性中σy

           
C(n) =
ijP(i,j,n)+C(n);            

       
end

    end

    C(n) =
(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))  % 输出数据;







0



                            <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 &amp;&amp; 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&nbsp;Toolbox&nbsp;大全(四)</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">抢沙发&gt;&gt;</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&nbsp;:</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="#">更多&gt;&gt;</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&amp;src=blogicp&amp;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>昵&nbsp;&nbsp;&nbsp;称:</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>&nbsp;&nbsp;&nbsp;--></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">&lt;&nbsp;前一篇</span><a href="http://blog.sina.com.cn/s/blog_ad4d75200102v9mr.html">[转载]【转载】MATLAB&nbsp;Toolbox&nbsp;大全(四)</a></div>
								  <div class="SG_floatR"><span class="SG_txtb">后一篇&nbsp;&gt;</span><a href="http://blog.sina.com.cn/s/blog_ad4d75200102vb48.html">2014未完成</a></div>
				</div>
	<div class="clearit"></div>
			
</div>
  • 2
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值