SIFT源码python实现

未使用包，python源码实现。 SIFT 尺度不变特征转换(Scale-invariant feature transform或SIFT)是一种电脑视觉的算法用来侦测与描述影像中的局部性特征，它在空间尺度中寻找极值点，并提取出其位置、尺度、旋转不变量，此算法由 David Lowe在1999年所发表，2004年完善总结。 其应用范围包含物体辨识、机器人地图感知与导航、影像缝合、3D模型建立、手势辨识、影像追踪和动作比对。 此算法有其专利，专利拥有者为英属哥伦比亚大学。 局部影像特征的描述与侦测可以帮助辨识物体，SIFT 特征是基于物体上的一些局部外观的兴趣点而与影像的大小和旋转无关。对于光线、噪声、些微视角改变的容忍度也相当高。基于这些特性，它们是高度显著而且相对容易撷取，在母数庞大的特征数据库中，很容易辨识物体而且鲜有误认。使用 SIFT特征描述对于部分物体遮蔽的侦测率也相当高，甚至只需要3个以上的SIFT物体特征就足以计算出位置与方位。在现今的电脑硬件速度下和小型的特征数据库条件下，辨识速度可接近即时运算。SIFT特征的信息量大，适合在海量数据库中快速准确匹配。 SIFT算法的特点有： 1. SIFT特征是图像的局部特征，其对旋转、尺度缩放、亮度变化保持不变性，对视角变化、仿射变换、噪声也保持一定程度的稳定性； 2. 独特性（Distinctiveness）好，信息量丰富，适用于在海量特征数据库中进行快速、准确的匹配； 3. 多量性，即使少数的几个物体也可以产生大量的SIFT特征向量； 4. 高速性，经优化的SIFT匹配算法甚至可以达到实时的要求； 5. 可扩展性，可以很方便的与其他形式的特征向量进行联合。 SIFT算法可以解决的问题： 目标的自身状态、场景所处的环境和成像器材的成像特性等因素影响图像配准/目标识别跟踪的性能。而SIFT算法在一定程度上可解决： 1. 目标的旋转、缩放、平移（RST） 2. 图像仿射/投影变换（视点viewpoint） 3. 光照影响（illumination） 4. 目标遮挡（occlusion） 5. 杂物场景（clutter） 6. 噪声 SIFT算法的实质是在不同的尺度空间上查找关键点(特征点)，并计算出关键点的方向。SIFT所查找到的关键点是一些十分突出，不会因光照，仿射变换和噪音等因素而变化的点，如角点、边缘点、暗区的亮点及亮区的暗点等。 Lowe将SIFT算法分解为如下四步： 1. 尺度空间极值检测：搜索所有尺度上的图像位置。通过高斯微分函数来识别潜在的对于尺度和旋转不变的兴趣点。 2. 关键点定位：在每个候选的位置上，通过一个拟合精细的模型来确定位置和尺度。关键点的选择依据于它们的稳定程度。 3. 方向确定：基于图像局部的梯度方向，分配给每个关键点位置一个或多个方向。所有后面的对图像数据的操作都相对于关键点的方向、尺度和位置进行变换，从而提供对于这些变换的不变性。 4. 关键点描述：在每个关键点周围的邻域内，在选定的尺度上测量图像局部的梯度。这些梯度被变换成一种表示，这种表示允许比较大的局部形状的变形和光照变化。 本文沿着Lowe的步骤，参考Rob Hess及Andrea Vedaldi源码，详解SIFT算法的实现过程。 未使用包，python源码实现。

import numpy
import Image


class contrast（object）:

	def lowcontrast（self arr）:

	    H W = arr.shape
	    re = []

	    for i in range（1H-1）:
                for j in range（1W-1）:
                    temp = arr[i-1:i+2j-1:j+2]
		    winStd = temp.std（）
		    winMean = temp.mean（）
		    if winStd == 0:
			continue				
		    elif abs（（arr[i][j] - winMean）/winStd） > 0.3:
					re.append（（ij））

            return re


# con  = contrast（）
# im = Image.open（‘f:\\5500\\tt\\s.jpg‘）.convert（‘L‘）
# ima = numpy.array（im）
# re = con.lowcontrast（ima）
# print len（re）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    .......       567  2015-05-27 20:56  SIFT-master\contrast.py

    .......      3020  2015-05-27 20:56  SIFT-master\descriptor.py

    .......      1537  2015-05-27 20:56  SIFT-master\factor.py

    .......      4156  2015-05-27 20:56  SIFT-master\findextrema.py

    .......      2644  2015-05-27 20:56  SIFT-master\harris.py

    .......      1810  2015-05-27 20:56  SIFT-master\hessian.py

    .......      2442  2015-05-27 20:56  SIFT-master\process.py

    .......       267  2015-05-27 20:56  SIFT-master\README.txt

     文件       4281  2018-11-16 13:43  SIFT-master\sift.py

     文件       2952  2018-11-16 13:43  SIFT-master\__pycache__\findextrema.cpython-37.pyc

     文件       2307  2018-11-16 13:43  SIFT-master\__pycache__\harris.cpython-37.pyc

     文件       1946  2018-11-16 13:43  SIFT-master\__pycache__\hessian.cpython-37.pyc

     文件       2820  2018-11-16 13:43  SIFT-master\__pycache__\process.cpython-37.pyc

     目录          0  2018-11-16 13:43  SIFT-master\__pycache__

     目录          0  2018-11-16 13:43  SIFT-master

----------- ---------  ---------- -----  ----

                30749                    15