在实际计算时,三种方法计算的金字塔组数noctaves,尺度空间坐标σ,以及每组金字塔内的层数S是一样的。同一时候,如果图像为640*480的标准图像。

金字塔层数:

当中o_min = 0,对于分辨率为640*480的图像N=5。

每组金字塔内图像数:
S=3。即在做极值检測时使用金子塔内中间3张图像。
对于LoG每组金字塔内有S+2张图像(S=-1。0,1,2,3),须要做S+1次高斯模糊操作(后一张图像由前一张做高斯模糊得到)。而DoG每组金字塔有S+3张高斯图像。得到S+2张DoG图像。
尺度空间系数:

当中,S表示每组金字塔内图像层数,n为当前高斯层数。取0-4。DoG须要5个尺度系数得到6张GSS图像。而LoG仅仅须要前4个尺度系数得到5张图像。

LoG

高斯核使用正太分布(高斯函数)计算模糊模版,N维空间正太分布方程为:

 

于是,二维高斯模板上的距离中心点为(x,y)的元素相应高斯计算公式为:

 

规范化的高斯拉普拉斯图像为

 

终于构造LoG金字塔有5层。每层有S+2=5张图像,每层金字塔内每张图像尺度是前一张的k倍。即构成的连续尺度序列:

当中o为当前金字塔层数,n为在当前金字塔层中图像张数。
因为卷积计算性质:

在计算时。通过对前一张图像做尺度系数为的卷积操作。能够降低卷积计算次数。故在金字塔每层内的S+2张图像,须要S+1次卷积操作,每次LoG核的尺度系数为:

 

watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQveGlhb3dlaV9jcXU=/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/SouthEast" alt="" />

图1. LoG金字塔示意图。 左側为图像尺度空间系数,红色矩形框为实际參与比較的空间系数,右側为在上一张图像做LoG卷积操作的尺度系数

DoG

因为LoG与Gauss核具有例如以下关系:

即,

 

因此。LoG算子能够用高斯差分算子DoG(Difference of Guassians)表示。

 

于是通过高斯金字塔每层内相邻两张图像相减能够得到DoG金字塔。对于最后须要S张图像寻找极值点,须要S+2张DoG图像。S+3张高斯图像。详细关系如图2.所看到的。

watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQveGlhb3dlaV9jcXU=/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/SouthEast" alt="" />

图 2. 由高斯金字塔得到DoG金字塔及其相应的尺度空间系数示意图。

LoG & DoG

一个直观的比較结果,使用分别计算LoG和DoG,得到:

watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQveGlhb3dlaV9jcXU=/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/SouthEast" alt="" />

能够看到。LoG比DoG明显须要很多其它的加法运算和乘法运算。尽管DoG须要在每层金字塔多做一次高斯操作(即为了得到S+2张DoG图须要S+3张高斯模糊图),但通过减法代替LoG核的计算过程,显著降低了运算次数,大大节省了运算时间。

(转载请注明作者和出处:http://blog.csdn.net/xiaowei_cqu 未经同意请勿用于商业用途)

【计算机视觉】SIFT中LoG和DoG比較的更多相关文章

  1. 【计算机视觉】SIFT中LoG和DoG比较

    <SIFT原理与源码分析>系列文章索引:http://www.cnblogs.com/tianyalu/p/5467813.html 在实际计算时,三种方法计算的金字塔组数noctaves ...

  2. 计算机视觉-sift(1)原理

    1999年由David Lowe首先发表于计算机视觉国际会议(International Conference on Computer Vision,ICCV),2004年再次经David Lowe整 ...

  3. SIFT中的尺度空间和传统图像金字塔

    SIFT中的尺度空间和传统图像金字塔 http://www.zhizhihu.com/html/y2010/2146.html 最近自己混淆了好多概念,一边弄明白的同时,也做了一些记录,分享一下.最近 ...

  4. matlab中log函数与rssi转距离

    我们通常所说的log是指以10为底的对数,而MATLAB中的log却不是这样.Matlab中的log函数在默认情况下是以e为底,即loge,如果需要计算以10为底的对数,那么需要用log10()函数. ...

  5. LoadRunner中log的使用总结

    LoadRunner中log的使用总结 1.log的设置方式. 在 runtime setting中可以设置log的生成方式: 默认的log方式: Enable logging选中,log optio ...

  6. Android中Log机制详解

    Android中Log的输出有如下几种: Log.v(String tag, String msg);        //VERBOSELog.d(String tag, String msg);   ...

  7. 0-Android系统各层中LOG的使用

     Android系统各层中LOG的使用 , ANDROID_LOG_DEFAULT, /* only for SetMinPriority() */ ANDROID_LOG_VERBOSE, ANDR ...

  8. Sift中尺度空间、高斯金字塔、差分金字塔(DOG金字塔)、图像金字塔

    转自:https://blog.csdn.net/dcrmg/article/details/52561656 一. 图像金字塔 图像金字塔是一种以多分辨率来解释图像的结构,通过对原始图像进行多尺度像 ...

  9. 计算机视觉-sift(2)代码理解

    之前结合不同人的资料理解了sift的原理,这里通过opencv中的代码来加深对sift的实现的理解. 使得能够从原理性理解到源码级的理解.不过该博文还是大量基于<赵春江, opencv2.4.9 ...

随机推荐

  1. 什么是JavaScript框架-------share

    摘要:现代网站和web应用程序趋向于依赖客户端的大量的javascript来提供丰富的交互.特别是通过不刷新页面的异步请求来返回数据或从服务器端的脚本(或数据系统)中得到响应.在这篇文章中,你将会了解 ...

  2. 点击增加删除class

    <!DOCTYPE html><html lang="en"><head> <meta charset="UTF-8" ...

  3. 【简●解】[HAOI2007] 理想的正方形

    [简●解][HAOI2007] 理想的正方形 可恶的\(DP\). [题目大意] 有一个\(a*b\)的整数组成的矩阵,现请你从中找出一个\(n*n\)的正方形区域,使得该区域所有数中的最大值和最小值 ...

  4. luogu 1113 杂务--啥?最长路?抱歉,我不会

    P1113 杂务 题目描述 John的农场在给奶牛挤奶前有很多杂务要完成,每一项杂务都需要一定的时间来完成它.比如:他们要将奶牛集合起来,将他们赶进牛棚,为奶牛清洗乳房以及一些其它工作.尽早将所有杂务 ...

  5. 条款14:在资源管理类中心copying行为(Think carefully about copying behavior in resource-manage classes)

    NOTE: 1.复制RAII 对象必须一并赋值它所管理的资源,所以资源的copying行为决定RAII对象的copying行为. 2.普遍而常见的RAII class copying 行为是: 抑制c ...

  6. windows下安装oracle客户端和php扩展

    先来抱怨下 ,按这玩楞费了我大半天的时间,一路的坑! 我的电脑是win7 64位的 第一步  打开php.ini  把 extension=php_oci8_12c.dll extension=php ...

  7. 大数据学习——关于hive中的各种join

    准备数据 2,b 3,c 4,d 7,y 8,u 2,bb 3,cc 7,yy 9,pp 建表: create table a(id int,name string) row format delim ...

  8. Struts2的Action配置的各项默认值

    1 如果没有为action指定class,默认是ActionSupport 2 如果没有为action指定method,默认执行action中的execute()方法 3 如果没有指定result的n ...

  9. Coloring Brackets (区间DP)

    Once Petya read a problem about a bracket sequence. He gave it much thought but didn't find a soluti ...

  10. 『NYIST』第八届河南省ACM竞赛训练赛[正式赛一]-CodeForces 236A,虽然很水,但有一个很简单的函数用起来方便

    A. Boy or Girl time limit per test 1 second memory limit per test 256 megabytes input standard input ...