【计算机视觉】匹配地理标记图像

1、提取特征

先对图像使用SIFT特征提取代码进行处理，并且将特征保存在和图像同名（但文件名是.sift，而不是.jpg）的文件路径下。

download_path = "D:\imagepi"
path = "D:\imagepi"

imlist = sift.get_imlist(download_path)
nbr_images = len(imlist)

featlist = [imname[:-3] + 'sift' for imname in imlist]
for i, imname in enumerate(imlist):
    sift.process_image(imname, featlist[i])

2、使用局部描述子匹配

我们可以对所有的组合图像进行逐个匹配，我们将每对图像间的匹配特征数保存在matchscores数组中。

matchscores = zeros((nbr_images, nbr_images))

for i in range(nbr_images):
    for j in range(i, nbr_images): #only compute upper triangle
        print( 'comparing ', imlist[i], imlist[j])
        l1, d1 = sift.read_features_from_file(featlist[i])
        l2, d2 = sift.read_features_from_file(featlist[j])
        matches = sift.match_twosided(d1, d2)
        nbr_matches = sum(matches>0)
        print(