1 files changed, 47 insertions, 0 deletions

diff --git a/SD-VBS/benchmarks/stitch/src/matlab/main.m b/SD-VBS/benchmarks/stitch/src/matlab/main.m
new file mode 100755
index 0000000..9be8504
--- /dev/null
+++ b/SD-VBS/benchmarks/stitch/src/matlab/main.m
@@ -0,0 +1,47 @@
+Icur_list=readData('capitol');
+Icur=Icur_list{1};
+%Icur=Icur(1:150,:,:);
+[x y v]=harris(Icur);
+interestPntsCur=getANMS(x, y, v, 24);
+Fcur=extractFeatures(Icur, interestPntsCur(:,1), interestPntsCur(:,2));
+
+for id=2:length(Icur_list)
+    Iprev=Icur;
+    Fprev=Fcur;
+    interestPntsPrev=interestPntsCur;
+    Icur=Icur_list{id};
+    %Icur=Icur(1:150,:,:);
+    [x y v]=harris(Icur);
+    %showInterestPoints(I, x, y)
+    interestPntsCur=getANMS(x, y, v, 16);
+    %showInterestPoints(I, interestPnts(:,1), interestPnts(:,2))
+    Fcur=extractFeatures(Icur, interestPntsCur(:,1), interestPntsCur(:,2));
+    matchedPntsIdx=matchFeatures(Fprev, Fcur);
+    matchedPntsPrev=interestPntsPrev(matchedPntsIdx(:,1),:);
+    matchedPntsCur=interestPntsCur(matchedPntsIdx(:,2),:);
+    
+    %subplot(1,2,1);showInterestPoints(Iprev, matchedPntsPrev(:,1), matchedPntsPrev(:,2));axis image
+    %subplot(1,2,2);showInterestPoints(Icur, matchedPntsCur(:,1), matchedPntsCur(:,2));axis image
+    %drawnow
+    %pause;
+    %printImage(['featureMatching' int2str(id-1) '_' int2str(id)]);
+    [retH retPntsIdx]=ransac(matchedPntsPrev, matchedPntsCur, 10000, 100);  
+    subplot(2,2,1);showInterestPoints(Iprev, matchedPntsPrev(retPntsIdx,1), matchedPntsPrev(retPntsIdx,2));axis image
+    subplot(2,2,2);showInterestPoints(Icur, matchedPntsCur(retPntsIdx,1), matchedPntsCur(retPntsIdx,2));axis image
+    T=maketform('projective', retH');
+    Tfn_Mat=T;
+    [nr nc nd]=size(Iprev);
+    Itrans=imtransform(Iprev,T,'XData', [1 nc], 'YData', [1 nr],'FillValues',[0;0;0]);
+    subplot(2,2,3);imshow(Itrans)    
+    subplot(2,2,4);imagesc(rgb2gray(abs(Icur-Itrans)))
+    drawnow
+    pause;
+    printImage(['ransac' int2str(id-1) '_' int2str(id)]);
+end
+
+    %[Itrans2 XData YData]=imtransform(Iprev,T,'FillValues',[0;0;0]);
+    %for the fast iteration
+    %[Itrans2 XData2 YData2]=imtransform(Iprev,T,'XData', [-2*nc 2*nc], 'YData', [-2*nr 2*nr],'FillValues',[0;0;0]);
+
+    %h=imshow(Itrans);
+