图像的结构张量MATLAB代码
来源:互联网 发布:c语言else是什么意思 编辑:程序博客网 时间:2024/06/06 00:52
转自:http://www.cnblogs.com/tiandsp/p/3307378.html
根据结构张量能区分图像的平坦区域、边缘区域与角点区域。
此算法也算是计算机科学最重要的32个算法之一了。链接的文章中此算法名称为Strukturtensor算法,不过我搜索了一下,Strukturtensor这个单词好像是德语,翻译过来就是structure tensor结构张量了。
此处所说的张量不是相对论或黎曼几何里的张量,黎曼几何的张量好多论文都叫张量场了。也不是数学界还没研究明白的对矩阵进行扩展的高阶张量,主要是张量分解。这里的结构张量就是一个矩阵,一个对图像像素进行组织的数据结构而已。
像素组织而成的矩阵如下:
这个公式太常见了,在harris角点检测中就用到了。其中Ix,Iy就是原对原图像在x和y方向求得的偏导。
然后求矩阵E的行列式K和迹H。然后根据K和H的关系就能区分图像的区域模式了。
模式分以下三类:
平坦区域:H=0;
边缘区域:H>0 && K=0;
角点区域:H>0 && K>0;
harris角点检测就用到了第三类判断。
当然,在实际应用的时候H和K的值肯定都不会是理想,所以我用的都是近似判断。
处理结果如下:
原图:
平坦区域:
边缘区域:
角点区域(好像也不全角点,求角点还是harris好了):
结构张量行列式与迹的关系:
其中红框为平坦区域,黄框为边缘区域,铝框为角点区域。
matlab代码如下:
clear all; close all; clc;img=double(imread('lena.jpg'));[m n]=size(img);imshow(img,[])[Ix Iy]=gradient(img);Ix2=Ix.^2;Iy2=Iy.^2;Ixy=Ix.*Iy;k=1;lambda=zeros(m*n,2);for i=1:m for j=1:n st=[Ix2(i,j) Ixy(i,j);Ixy(i,j) Iy2(i,j)]; %结构张量 K=det(st); %求行列式 H=trace(st); %求迹 %所有的判断都是近似的 % if H<50 %认为是平坦区域 % if H>50 && abs(K)<0.01*10^(-9) %认为是边缘区域 if H>50 && abs(K)>0.01*10^(-9) %认为是角点区域 img(i,j)=255; end lambda(k,:)=[K H]; k=k+1; endendfigure;plot(lambda(:,1),lambda(:,2),'.');ylabel('trace');xlabel('det');figure;imshow(img,[])
0 0
- 图像的结构张量MATLAB代码
- MATLAB与图像处理(一):图像中结构张量的求法(Structure Tensor)
- MATLAB与图像处理(一):图像中结构张量的求法(Structure Tensor)
- 结构张量
- matlab 稀疏张量的HOSVD 推荐系统
- matlab练习程序(结构张量structure tensor)
- matlab练习程序(结构张量structure tensor)
- 图像的Mallat算法分解(Matlab代码)
- 【转】图像复原的matlab代码
- 基于PCNN的图像增强matlab代码
- Matlab代码实现图像的傅里叶变换
- 张量的展开与matlab下的工具包操作
- 张量的展开与matlab下的工具包操作
- matlab tensor张量库的安装与测试
- 结构张量扩散
- 结构张量structure tensor
- 部分matlab处理图像代码
- 图像细化matlab代码实现
- string.replaceAll()中的特殊字符($ \)与matcher.appendReplacement
- Linux下的tar压缩解压缩命令详解
- 弹出软键盘挤压上一个activity页面变形
- C.土豪银行
- java工厂模式示例
- 图像的结构张量MATLAB代码
- 基于Java NIO的Socket通信
- Linux系统中NCBI BLAST+本地化教程
- mybatis-范围和生命周期
- [BZOJ1853][SCOI2010][容斥原理]幸运数字
- linux压缩和解压缩命令大全
- 第一个只出现一次的字符
- J2EE-MVC分层及常用词解释
- 字符串处理 explode--使用一个字符串分割另一个字符串