sift算法的MATLAB程序完整版.doc

  1. 1、本文档共2页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多

%[image,descriptors,locs]=sift(imageFile)

%

%ThisfunctionreadsanimageandreturnsitsSIFTkeypoints.

%Inputparameters:

%imageFile:thefilenamefortheimage.

%

%Returned:

%image:theimagearrayindoubleformat

%descriptors:aK-by-128matrix,whereeachrowgivesaninvariant

%descriptorforoneoftheKkeypoints.Thedescriptorisavector

%of128valuesnormalizedtounitlength.

%locs:K-by-4matrix,inwhicheachrowhasthe4valuesfora

%keypointlocation(row,column,scale,orientation).The

%orientationisintherange[-PI,PI]radians.

%

%Credits:ThanksforinitialversionofthisprogramtoD.Alvaroand

%J.J.Guerrero,UniversidaddeZaragoza(modifiedbyD.Lowe)

function[image,descriptors,locs]=sift(imageFile)

%Loadimage

image=imread(imageFile);

%IfyouhavetheImageProcessingToolbox,youcanuncommentthefollowing

%linestoallowinputofcolorimages,whichwillbeconvertedtograyscale.

%ifisrgb(image)

%image=rgb2gray(image);

%end

[rows,cols]=size(image);

%ConvertintoPGMimagefile,readablebykeypointsexecutable

f=fopen(tmp.pgm,w);

iff==-1

error(Couldnotcreatefiletmp.pgm.);

end

fprintf(f,P5\n%d\n%d\n255\n,cols,rows);

fwrite(f,image,uint8);

fclose(f);

%Callkeypointsexecutable

ifisunix

command=!./sift;

else

command=!siftWin32;

end

command=[commandtmp.pgmtmp.key];

eval(command);

%Opentmp.keyandcheckitsheader

g=fopen(tmp.key,r);

ifg==-1

error(Couldnotopenfiletmp.key.);

end

[header,count]=fscanf(g,%d%d,[12]);

ifcount~=2

error(Invalidkeypointfilebeginning.);

end

num=header(1);

len=header(2);

iflen~=128

error(Keypointdescriptorlengthinvalid(shouldbe128).);

end

%Createsthetwooutputmatrices(useknownsizeforefficiency)

locs=double(zeros(num,4));

descriptors=double(zeros(num,128));

%Parsetmp.key

fori

文档评论(0)

王水玉 + 关注
实名认证
内容提供者

前面

1亿VIP精品文档

相关文档