Matlab粗糙度曲线画中线,[转载]（zz）Tamura纹理特征的matlab实现（一）---粗糙度...

2023-10-23 19:50
文章标签 matlab 实现 特征 转载 曲线 纹理 粗糙度 zz 中线 tamura

本文主要是介绍Matlab粗糙度曲线画中线,[转载]（zz）Tamura纹理特征的matlab实现（一）---粗糙度...，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

搜了好久也没找到Tamura纹理特征的实现代码，于是自己动手丰衣足食，找出了最原始的Tamura的论文《Textural

Features Corresponding to Visual

Perception》来研读。今天下午实现了粗糙度(coarseness)部分，与大家共享，有错误的地方，还望大家指正。后续的对比度(contrast)、方向度(directionality)、线性度(linelikeness)、规则度(regularity)、粗略度(roughness)会陆续实现并共享。

Tamura纹理特征要比灰度共生矩阵得到的纹理特征更直观，在视觉效果上更有优势，由于比较懒，一直避重就轻，到今天终于要直面问题了，代码如下：

coarseness.m

%graypic为待处理的灰度图片，2^kmax为最大窗口

function Fcrs=coarseness(graypic,kmax)

%获取图片大小

[h,w]=size(graypic);

%平均灰度值矩阵A

A=zeros(h,w,2^kmax);

%计算有效可计算范围内每个点的2^k邻域内的平均灰度值

for i=2^(kmax-1)+1:h-2^(kmax-1)

for j=2^(kmax-1)+1:w-2^(kmax-1)

for k=1:kmax

A(i,j,k)=mean2(graypic(i-2^(k-1):i+2^(k-1)-1,j-2^(k-1):j+2^(k-1)-1));

end

end

end

%对每个像素点，计算在水平和垂直方向上不重叠窗口之间的Ak差

for i=1+2^(kmax-1):h-2^(kmax-1)

for j=1+2^(kmax-1):w-2^(kmax-1)

for k=1:kmax

Eh(i,j,k)=abs(A(i+2^(k-1),j,k)-A(i-2^(k-1),j));

Ev(i,j,k)=abs(A(i,j+2^(k-1),k)-A(i,j-2^(k-1)));

end

end

end

%对每个像素点计算使E达到最大值的k

for i=2^(kmax-1)+1:h-2^(kmax-1)

for j=2^(kmax-1)+1:w-2^(kmax-1)

[maxEh,p]=max(Eh(i,j,:));

[maxEv,q]=max(Ev(i,j,:));

if

maxEh>maxEv

maxkk=p;

else

maxkk=q;

end

Sbest(i,j)=2^maxkk;

%每个像素点的最优窗口大小为2^maxkk

end

end

%所有Sbest的均值作为整幅图片的粗糙度

Fcrs=mean2(Sbest);

调用：

image=rgb2gray(imread('example.jpg'));

f=coarseness(image,5)

Tamura纹理特征的matlab实现(二)---对比度

a4c26d1e5885305701be709a3d33442f.png (2010-01-13 21:40:27)

TE<

a4c26d1e5885305701be709a3d33442f.png转载TE<

var $tag='tamura,纹理特,对比度,it'; var

$tag_code='c57dde215fc1a8b3216796890dd4f569'; var

$r_quote_bligid='59ead5d90100gx4v'; var $worldcup='0'; var

$worldcupball='0'; 标签：

tamura

纹理特

对比度

it

分类： 学术

Tamura纹理特征之对比度(contrast)，公式如下：

a4c26d1e5885305701be709a3d33442f.png

n一般取1/4。

contrast.m

%graypic为待处理的灰度图片

function Fcon=contrast(graypic)

%二维向量一维化

x=graypic(:);

%四阶矩

M4=mean((x-mean(x))^4);

%方差

delta2=var(x,1);

%峰度

alfa4=M4/(delta2^2);

%标准差

delta=std(x,1);

%对比度

Fcon=delta/(alfa4^(1/4));

调用：

image=rgb2gray(imread('example.jpg'));

f=contrast(image)

Tamura纹理特征的matlab实现(三)---方向度

(2010-01-14 20:52:56)

TE<

a4c26d1e5885305701be709a3d33442f.png转载TE<

var $tag='tamura,纹理,方向度,it'; var

$tag_code='4bad3c6c0fc81846213a62dc64c1a996'; var

$r_quote_bligid='59ead5d90100gxkc'; var $worldcup='0'; var

$worldcupball='0'; 标签：

tamura

纹理

方向度

it

分类： 学术

directionality.m

function

[Fdir,sita]=directionality(graypic)

%sita为各像素点的角度矩阵，在线性度中会用到，所以这里作为结果返回

[h w]=size(graypic);

%两个方向的卷积矩阵

GradientH=[-1 0 1

-1

0 1

-1

0 1];

GradientV=[ 1 1 1

0

0 0

-1

-1 -1];

%卷积，取有效结果矩阵

MHconv=conv2(graypic,GradientH);

MH=MHconv(3:h,3:w);

MVconv=conv2(graypic,GradientV);

MV=MVconv(3:h,3:w)

%向量模

MG=(abs(MH)+abs(MV))./2;

%有效矩阵大小

validH=h-2;

validW=w-2

%各像素点的方向

for i=1:validH

for j=1:validW

sita(i,j)=atan(MV(i,j)/MH(i,j))+(pi/2);

end

end

n=16;

t=12;

Nsita=zeros(1,n);

%构造方向的统计直方图

for i=1:validH

for j=1:validW

for k=1:n

if

sita(i,j)>=(2*(k-1)*pi/2/n)

&&

sita(i,j)

&& MG(i,j)>=t

Nsita(k)=Nsita(k)+1;

end

end

end

end

for k=1:n

HD(k)=Nsita(k)/sum(Nsita(:));

end

%假设每幅图片只有一个方向峰值，为计算方便简化了原著

[maxvalue,FIp]=max(HD);

Fdir=0;

for k=1:n

Fdir=Fdir+(k-FIp)^2*HD(k);%公式与原著有改动

end

调用：

image=rgb2gray(imread('example.jpg'));

[Fdir,sita]=directionality(image)

Tamura纹理特征的matlab实现(四)---线性度

(2010-01-14 21:01:04)

TE<

a4c26d1e5885305701be709a3d33442f.png转载TE<

var $tag='tamura,纹理,线性度,it'; var

$tag_code='e3de83ab0ce44a1975fe0fc93d4a6717'; var

$r_quote_bligid='59ead5d90100gxke'; var $worldcup='0'; var

$worldcupball='0'; 标签：

tamura

纹理

线性度

it

分类： 学术

linelikeness.m

function

Flin=linelikeness(graypic,sita,d)

%d为共生矩阵计算时的像素间隔距离

n=16;

[h,w]=size(graypic);

%构造方向共生矩阵

PDd1=zeros(n,n);

PDd2=zeros(n,n);

PDd3=zeros(n,n);

PDd4=zeros(n,n);

PDd5=zeros(n,n);

PDd6=zeros(n,n);

PDd7=zeros(n,n);

PDd8=zeros(n,n);

for i=d+1:h-d-2

for j=d+1:w-d-2

for m1=1:n

for

m2=1:n

%下方向

if

(sita(i,j)>=(2*(m1-1)*pi/2/n)

&&

sita(i,j)

&&

(sita(i+d,j)>=(2*(m2-1)*pi/2/n)

&&

sita(i+d,j)

PDd1(m1,m2)=PDd1(m1,m2)+1;

end

%上方向

if

(sita(i,j)>=(2*(m1-1)*pi/2/n)

&&

sita(i,j)

&&

(sita(i-d,j)>=(2*(m2-1)*pi/2/n)

&&

sita(i-d,j)

PDd2(m1,m2)=PDd2(m1,m2)+1;

end

%右方向

if

(sita(i,j)>=(2*(m1-1)*pi/2/n)

&&

sita(i,j)

&&

(sita(i,j+d)>=(2*(m2-1)*pi/2/n)

&&

sita(i,j+d)

PDd3(m1,m2)=PDd3(m1,m2)+1;

end

%左方向

if

(sita(i,j)>=(2*(m1-1)*pi/2/n)

&&

sita(i,j)

&&

(sita(i,j-d)>=(2*(m2-1)*pi/2/n)

&&

sita(i,j-d)

PDd4(m1,m2)=PDd4(m1,m2)+1;

end

%右下方向

if

(sita(i,j)>=(2*(m1-1)*pi/2/n)

&&

sita(i,j)

&&

(sita(i+d,j+d)>=(2*(m2-1)*pi/2/n)

&&

sita(i+d,j+d)

PDd5(m1,m2)=PDd5(m1,m2)+1;

end

%右上方向

if

(sita(i,j)>=(2*(m1-1)*pi/2/n)

&&

sita(i,j)

&&

(sita(i-d,j+d)>=(2*(m2-1)*pi/2/n)

&&

sita(i-d,j+d)

PDd6(m1,m2)=PDd6(m1,m2)+1;

end

%左下方向

if

(sita(i,j)>=(2*(m1-1)*pi/2/n)

&&

sita(i,j)

&&

(sita(i+d,j-d)>=(2*(m2-1)*pi/2/n)

&&

sita(i+d,j-d)

PDd7(m1,m2)=PDd7(m1,m2)+1;

end

%左上方向

if

(sita(i,j)>=(2*(m1-1)*pi/2/n)

&&

sita(i,j)

&&

(sita(i-d,j-d)>=(2*(m2-1)*pi/2/n)

&&

sita(i-d,j-d)

PDd8(m1,m2)=PDd8(m1,m2)+1;

end

end

end

end

end

f=zeros(1,8);

g=zeros(1,8);

for i=1:n

for j=1:n

f(1)=f(1)+PDd1(i,j)*cos((i-j)*2*pi/n);

g(1)=g(1)+PDd1(i,j);

f(2)=f(2)+PDd2(i,j)*cos((i-j)*2*pi/n);

g(2)=g(2)+PDd2(i,j);

f(3)=f(3)+PDd3(i,j)*cos((i-j)*2*pi/n);

g(3)=g(3)+PDd3(i,j);

f(4)=f(4)+PDd4(i,j)*cos((i-j)*2*pi/n);

g(4)=g(4)+PDd4(i,j);

f(5)=f(5)+PDd5(i,j)*cos((i-j)*2*pi/n);

g(5)=g(5)+PDd5(i,j);

f(6)=f(6)+PDd6(i,j)*cos((i-j)*2*pi/n);

g(6)=g(6)+PDd6(i,j);

f(7)=f(7)+PDd7(i,j)*cos((i-j)*2*pi/n);

g(7)=g(7)+PDd7(i,j);

f(8)=f(8)+PDd8(i,j)*cos((i-j)*2*pi/n);

g(8)=g(4)+PDd8(i,j);

end

end

tempM=f./g;

Flin=max(tempM);%取8个方向的线性度最大值作为图片的线性度

调用：

image=rgb2gray(imread('example.jpg'));

Flin=linelikeness(image,sita,4)

%sita为directionality.m返回的结果

Tamura纹理特征的matlab实现(五)---规则度

(2010-01-14 21:04:09)

TE<

a4c26d1e5885305701be709a3d33442f.png转载TE<

var $tag='tamura,纹理,规则度,it'; var

$tag_code='39947dd3859a0df41f3f4de84b8d3520'; var

$r_quote_bligid='59ead5d90100gxkj'; var $worldcup='0'; var

$worldcupball='0'; 标签：

tamura

纹理

规则度

it

分类： 学术

regularity.m

function

Freg=regularity(graypic,windowsize)

%windowsize为计算规则度的子窗口大小

[h,w]=size(graypic);

k=0;

for i=1:windowsize:h-windowsize

for j=1:windowsize:w-windowsize

k=k+1;

crs(k)=coarseness(graypic(i:i+windowsize-1,j:j+windowsize-1),5);

%粗糙度

con(k)=contrast(graypic(i:i+windowsize-1,j:j+windowsize-1));

%对比度

[dire(k),sita]=directionality(graypic(i:i+windowsize-1,j:j+windowsize-1));%方向度

lin=linelikeness(graypic(i:i+windowsize-1,j:j+windowsize-1),sita,4)*10;

%线性度，*10与crs、con、dire同量级化

end

end

%求上述各参数的标准差

Dcrs=std(crs,1);

Dcon=std(con,1);

Ddir=std(dire,1);

Dlin=std(lin,1);

%规则度

Freg=1-(Dcrs+Dcon+Ddir+Dlin)/4/100;

调用：

image=rgb2gray(imread('example.jpg'));

Freg=regularity(image,64)

Tamura纹理特征的matlab实现(六)---粗略度

(2010-01-14 21:11:07)

TE<

a4c26d1e5885305701be709a3d33442f.png转载TE<

var $tag='tamura,纹理,粗略度,it'; var

$tag_code='fc001e7750cac6917c3061ca68bff398'; var

$r_quote_bligid='59ead5d90100gxkl'; var $worldcup='0'; var

$worldcupball='0'; 标签：

tamura

纹理

粗略度

it

分类： 学术

粗略度计算比较简单，是粗糙度和对比度两个属性的合成：

Frgh=Fcrs+Fcon

至此Tamura纹理特征的六个属性终于实现了，工作了一天半，好辛苦。属于经典算法，并没有创新性，所以还要在此基础上做进一步的研究和应用。

这篇关于Matlab粗糙度曲线画中线,[转载]（zz）Tamura纹理特征的matlab实现（一）---粗糙度...的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


http://www.chinasem.cn/article/270100。

相关文章

Idea实现接口的方法上无法添加@Override注解的解决方案

Idea实现接口的方法上无法添加@Override注解的解决方案

《Idea实现接口的方法上无法添加@Override注解的解决方案》文章介绍了在IDEA中实现接口方法时无法添加@Override注解的问题及其解决方法,主要步骤包括更改项目结构中的Languagel... 目录Idea实现接China编程口的方法上无法添加@javascriptOverride注解错误原因解决方
阅读更多...
轻松上手MYSQL之JSON函数实现高效数据查询与操作

轻松上手MYSQL之JSON函数实现高效数据查询与操作

《轻松上手MYSQL之JSON函数实现高效数据查询与操作》：本文主要介绍轻松上手MYSQL之JSON函数实现高效数据查询与操作的相关资料,MySQL提供了多个JSON函数,用于处理和查询JSON数... 目录一、jsON_EXTRACT 提取指定数据二、JSON_UNQUOTE 取消双引号三、JSON_KE
阅读更多...
MySql死锁怎么排查的方法实现

MySql死锁怎么排查的方法实现

《MySql死锁怎么排查的方法实现》本文主要介绍了MySql死锁怎么排查的方法实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧... 目录前言一、死锁排查方法1. 查看死锁日志方法 1：启用死锁日志输出方法 2：检查 mysql 错误
阅读更多...
CSS3中使用flex和grid实现等高元素布局的示例代码

CSS3中使用flex和grid实现等高元素布局的示例代码

《CSS3中使用flex和grid实现等高元素布局的示例代码》：本文主要介绍了使用CSS3中的Flexbox和Grid布局实现等高元素布局的方法，通过简单的两列实现、每行放置3列以及全部代码的展示，展示了这两种布局方式的实现细节和效果，详细内容请阅读本文，希望能对你有所帮助... 过往的实现方法是使用浮动加
阅读更多...
Go Mongox轻松实现MongoDB的时间字段自动填充

Go Mongox轻松实现MongoDB的时间字段自动填充

《GoMongox轻松实现MongoDB的时间字段自动填充》这篇文章主要为大家详细介绍了Go语言如何使用mongox库,在插入和更新数据时自动填充时间字段,从而提升开发效率并减少重复代码,需要的可以... 目录前言时间字段填充规则Mongox 的安装使用 Mongox 进行插入操作使用 Mongox 进行更
阅读更多...
MySQL修改密码的四种实现方式

MySQL修改密码的四种实现方式

《MySQL修改密码的四种实现方式》文章主要介绍了如何使用命令行工具修改MySQL密码,包括使用`setpassword`命令和`mysqladmin`命令,此外,还详细描述了忘记密码时的处理方法,包... 目录mysql修改密码四种方式一、set password命令二、使用mysqladmin三、修改u
阅读更多...
Java实现状态模式的示例代码

Java实现状态模式的示例代码

《Java实现状态模式的示例代码》状态模式是一种行为型设计模式,允许对象根据其内部状态改变行为,本文主要介绍了Java实现状态模式的示例代码,文中通过示例代码介绍的非常详细,需要的朋友们下面随着小编来... 目录一、简介1、定义2、状态模式的结构二、Java实现案例1、电灯开关状态案例2、番茄工作法状态案例
阅读更多...
一文教你使用Python实现本地分页

一文教你使用Python实现本地分页

《一文教你使用Python实现本地分页》这篇文章主要为大家详细介绍了Python如何实现本地分页的算法,主要针对二级数据结构,文中的示例代码简洁易懂,有需要的小伙伴可以了解下... 在项目开发的过程中，遇到分页的第一页就展示大量的数据，导致前端列表加载展示的速度慢，所以需要在本地加入分页处理，把所有数据先放
阅读更多...
SpringMVC前后端传值的几种实现方式

SpringMVC前后端传值的几种实现方式

《SpringMVC前后端传值的几种实现方式》本文主要介绍了SpringMVC前后端传值的方式实现,包括使用HttpServletRequest、HttpSession、Model和ModelAndV... 目录一、从Controller层到JSP界面1、使用HttpServletRequest的方式2、使
阅读更多...
树莓派启动python的实现方法

树莓派启动python的实现方法

《树莓派启动python的实现方法》本文主要介绍了树莓派启动python的实现方法,文中通过图文介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧... 目录一、RASPBerry系统设置二、使用sandroidsh连接上开发板Raspberry Pi三、运
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2