close all;clc,clear;
%% image name and read image
%双焦面c2图像
fileName_c2 = 'G:\大仪器项目\我的工作\data\标定曲线测试\data2\two focal\c2-1\Image1_1.tif';
% 扫描图像
fileName_scan = 'G:\大仪器项目\我的工作\data\标定曲线测试\data2\scan-c2\scan-c2-copy\0_0_0_200_0.tif';
img_c2 = imread(fileName_c2);
img_scan = imread(fileName_scan);
%% find features
img_c2_Points = detectFASTFeatures(img_c2);
img_scan_Points = detectFASTFeatures(img_scan);
%% image show
figure(1);
imshow(img_c2);
title('100 Feature Points');
hold on;
plot(selectStrongest(img_c2_Points,100));
figure(2);
imshow(img_scan);
title('300 Feature Points');
hold on;
plot(selectStrongest(img_scan_Points, 300));
%% Extract feature descriptor
[img_c2_Features, img_c2_Points] = extractFeatures(img_c2, img_c2_Points);
[img_scan_Features, img_scan_Points] = extractFeatures(img_scan, img_scan_Points);
%% find Putative point matches
Point_Pairs = matchFeatures(img_c2_Features, img_scan_Features);
%% display matched features
matchedImg_c2Points = img_c2_Points(Point_Pairs(:, 1), :);
matchedImg_scanPoints = img_scan_Points(Point_Pairs(:, 2), :);
figure(3);
showMatchedFeatures(img_c2, img_scan, matchedImg_c2Points, matchedImg_scanPoints, 'montage');
title('Putatively Matched Points(Including Outliers)');