%% 时间:2015年5月8日 %% 修改人:961256834(qq) clc; clear all; close all; %% 计时开始 tic; %% 输入图像 ori_image = imread('girl.bmp');   %% 判断该图像是彩色图像还是灰度图像,如果是彩色图像把它转换为灰度图像 [x,y,z] = size(ori_image);            if(z~=1)     ori_image = rgb2gray(ori_image);    end   %% 计算像素的自相关矩阵 fx = [-2 -1 0 1 2];            % x方向梯度算子 Ix = filter2(fx,ori_image);    % x方向滤波 fy = [-2;-1;0;1;2];            % y方向梯度算子) Iy = filter2(fy,ori_image);    % y方向滤波 Ix2 = Ix.^2; Iy2 = Iy.^2; Ixy = Ix.*Iy; clear Ix; clear Iy; %% 产生7*7的高斯窗函数,sigma=2,对图像进行高斯滤波 h = fspecial('gaussian',[7 7],2);       Ix2 = filter2(h,Ix2); Iy2 = filter2(h,Iy2); Ixy = filter2(h,Ixy); %% 计算像素的兴趣值 height = size(ori_image,1);              %计算ori_im矩阵的行数 width = size(ori_image,2);               %计算ori_im矩阵的列数 result = zeros(height,width);            % 纪录角点位置,角点处值为1 R = zeros(height,width); Rmax = 0;                                % 图像中最大的R值 for i = 1:height     for j = 1:width         M = [Ix2(i,j) Ixy(i,j);Ixy(i,j) Iy2(i,j)];                      R(i,j) = det(M)-0.06*(trace(M))^2;                             if R(i,j) > Rmax             Rmax = R(i,j);         end;     end; end; %% 进行非极大抑制,窗口大小3*3 cnt = 0; for i = 2:height-1     for j = 2:width-1                if R(i,j)>0.002*Rmax && R(i,j)>R(i-1,j-1) && R(i,j)>R(i-1,j) && ...            R(i,j)>R(i-1,j+1) && R(i,j)>R(i,j-1) && R(i,j)>R(i,j+1) && ...            R(i,j)>R(i+1,j-1) && R(i,j)>R(i+1,j) && R(i,j) > R(i+1,j+1)             result(i,j) = 1;             cnt = cnt+1;         end;     end; end; [posc,posr] = find(result == 1); cnt;  numbercnt = cnt; disp(['特征点的个数为',num2str(numbercnt),'个。']);                                    %% 结束计时 toc; t = toc; disp(['本程序的运行时间为',num2str(t),'秒。']) %% 先给出原来的图像 figure(1);imshow(ori_image);title('原图像');xlabel('图像的列数');ylabel('图像的行数');axis on; %% 描绘特征点分布 figure(2);plot(posr,posc,'B*');view(0,-90);title('特征点的分布');xlabel('图像的列数');ylabel('图像的行数');axis on; %% 图像中显示特征点 figure(3);imshow(ori_image);hold on;plot(posr,posc,'R*');title('图像中显示特征点');xlabel('图像的列数');ylabel('图像的行数');axis on;

harris算法 matlab,Harris算法的Matlab代码相关推荐

  1. 测角误差估计算法matlab,Harris角点检测 及 Matlab实验

    目录 1 基础知识 1.1 图像变化的类型 1.2 提取点特征的作用 1.3 什么是好的角点检测算法? 2 Harris 角点检测 2.1 Harris角点检测基本思想 2.2 Harris角点检测: ...

  2. 基于MATLAB图像harris角点提取算法仿真

    figure; ori_im2 = imread('ucsb1.jpg');     % 读取图像 ori_im2=rgb2gray(ori_im2);  fx = [6 0 -6;8 0 -8;6 ...

  3. k均值算法原理详细讲解以及matlab代码实现

    有研究生物电信号处理和机器学习的欢迎加我qq429500506共同交流学习进步. 最近更新文章的频率太低了,主要原因是不想为了发文章而发文章,想潜心研究,写好文章,顺便想说一句开源万岁,最近一个月虽然 ...

  4. matlab中存档算法代码,MATLAB 智能算法超级学习手册中程序代码

    [实例简介] MATLAB 智能算法超级学习手册中程序代码 [实例截图] [核心代码] dc90ef43-7920-434e-bdb8-0636c31c0b44 └── MATLAB 智能算法超级学习 ...

  5. 快速行进算法(fast_marching_kroon)的matlab代码

    快速行进算法(fast_marching_kroon)的matlab代码 快速行进算法用于求解程函方程得到走时场,源代码转载于link 文章目录 快速行进算法(fast_marching_kroon) ...

  6. 货物配送问题的matlab,免疫算法求解配送中心选址问题matlab代码

    免疫算法求解配送中心选址问题matlab代码 所属分类:matlab例程 开发工具:matlab 文件大小:31KB 下载次数:3 上传日期:2020-12-01 16:11:00 上 传 者:代码大 ...

  7. 【智能优化算法-白鲸算法】基于白鲸优化算法求解多目标优化问题附matlab代码

    1 内容介绍 在本文中,一种新的基于群体的元启发式算法灵感来自白鲸的行为鲸鱼,称为白鲸优化(BWO),是为了解决优化问题而提出的.三在 BWO 中建立了探索.开发和鲸落的阶段,对应于成对游泳.猎物和鲸 ...

  8. 智能优化算法:海鸥算法原理及Matlab代码

    智能优化算法:海鸥算法原理及Matlab代码 文章导航 1. 算法原理 2. SOA算法流程 2.1 迁徙行为(exploration ability) 2.2 攻击行为(exploitation a ...

  9. matlab 浮雕算法,photoshop图像滤镜——浮雕算法(含matlab代码)

    一.photoshop图像滤镜--浮雕算法 浮雕的算法是对图像的每一个点进行卷积处理,采用的矩阵如下 [100:00 0:00-1]: 假设原图像为X,处理后的图像为Y ,也就是说,对于坐标为(i,j ...

  10. 【路径优化】基于帝企鹅算法求解TSP问题(Matlab代码实现)

    目录 1 帝企鹅算法 2 旅行商问题(TSP) 3 运行结果 4 参考文献 5 Matlab代码实现 1 帝企鹅算法 帝企鹅优化算法(emperor penguin optimizer,EPO)是Ga ...

最新文章

  1. HTML, CSS, Javascript, jQuery之间的关系
  2. std string与线程安全_这才是现代C++单例模式简单又安全的实现
  3. 码神之路博客项目构建记录
  4. Onew——全球领先的区块链数字资产服务商
  5. 【杂文】Do A Slash
  6. 完成端口(CompletionPort)详解 - 手把手教你玩转网络编程系列之三 1-转
  7. app登录界面背景 css_Google flutter这么火?撸一个APP登录界面(上)
  8. igs无法分配驱动器映射表_明白步进驱动器原理,就知道为什么步进电机要加驱动器了...
  9. 零基础学php rar,php实现rar文件的读取和解压
  10. 金蝶专业版过账提示运行时错误5_金蝶KIS专业版常见故障及解决办法
  11. [转贴]请珍惜你身边爱你的人
  12. 华为 DAYU 详细介绍
  13. Unity3D: 做自己的卡牌游戏 (一)
  14. Python 的 sum():Pythonic 的求和方法
  15. python解决哲学家就餐问题(and型信号量)
  16. Maven Archetype脚手架制作
  17. 电器企业网站建设-电器网站设计建设方案
  18. 后台前端解决方案vue-element-admin的安装及使用
  19. 【原创】性能测试面试题及答案 - 2020年1月19日最新整理
  20. 如何把未压缩的.avi文件批量地转为.yuv文件(yuv420)?

热门文章

  1. 以物言志——溟㠭(展)艺术
  2. java localdate_Java 时间类-Calendar、Date、LocalDate/LocalTime
  3. 论文解读:Sadeepcry:使用自我注意和自动编码器网络的蛋白质结晶倾向预测的深度学习框架
  4. HTTP响应报文体和长连接详解
  5. 如何用python抓取qq音乐_手把手教你使用Python抓取QQ音乐数据(第二弹)
  6. 个人支付方案(免签约)-支付宝当面付 附带源码
  7. 吸附钴Co的螯合树脂
  8. Android Studio 4.2的SDK安装问题
  9. 小米翻译机闪亮登场,它称得上是便携式“翻译官”吗?
  10. c语言 细胞自动机,Leetcode 289. 生命游戏 (细胞自动机)