正文
显示图像的SIFT flow描述子
小程序:扫一扫查出行
【扫一扫了解最新限行尾号】
复制小程序
【扫一扫了解最新限行尾号】
复制小程序
close all;% 模拟figure 5im = zeros(401,401,3);im(:,:,:) = 0;im(2:200, 2:200, 2) = 255;im(202:400, 202:400, 2) = 255;im(2:200, 202:400, :) = 255;im(202:400, 2:200, :) = 255;im=im2double(im);% im = imread('0000.png');figure;imshow(im);cellsize=3;gridspacing=1;winSize = 3;addpath(fullfile(pwd,'mexDenseSIFT'));addpath(fullfile(pwd,'mexDiscreteFlow'));%不加噪声,(201,201)sift分布sift = mexDenseSIFT(im,cellsize,gridspacing);sift = imnormalize(double(sift));x = 1:128;y = squeeze(sift(100, 100, :));% figure,plot(x, y);% xlabel('The dimension of the descriptor in SIFT flow','FontSize',12);% ylabel('Distribution','FontSize',12);% set(gca, 'FontSize', 12);%计算该点sift的均值和方差M1 = mean(y)V1 = var(double(y))% %加噪声,(201,201)sift分布% im_noise =imnoise(im, 'gaussian', 0.02);% sift_noise = mexDenseSIFT(im_noise,cellsize,gridspacing);% x = 1:128;% y = squeeze(sift_noise(201, 201, :));% figure,plot(x, y);% figure, imshow(im_noise);% %计算该点sift的均值和方差% M2 = mean(y)% V2 = var(double(y))%不加噪声,(100,100)sift分布y1 = squeeze(sift(201, 201, :));figure,plot( x, y1, 'r', x, y, 'b', x, y, '--g', 'LineWidth', 2);xlabel('Dimension Degree of SIFT Flow Descriptor','FontSize',20);ylabel('Distribution','FontSize',20);set(gca, 'FontSize', 20);axis([0, 130, -0.1, 1.1]);hleg1 = legend('B', 'A', 'C');%计算该点sift的均值和方差M3 = mean(y)V3 = var(double(y))% %%加噪声,(100,100)sift分布% y = squeeze(sift_noise(100, 100, :));% figure,plot(x, y);% %计算该点sift的均值和方差% M4 = mean(y)% V4 = var(double(y))