如何在MATLAB中忽略黑色边框

时间:2014-08-29 07:56:55

标签: image matlab image-processing

我在2个非常相似的图像(通过代码RANSAC)之间做了一些比较,之后我将其中一个图像旋转到第一个图像的角度。 问题是在某些图像中你有一个黑色边框,扭曲了主持人和旋转 我是如何让主持人只在没有黑色边框的图像上(忽略它)

function [ output_args ] = ransac( )
%frames( filename)
global numFrames 
a=sprintf('Ransac begin');
disp(a);
for i=1:numFrames 
file_name = sprintf('frames/%0.3i.jpg', i);
file_name2 = sprintf('frames/%0.3i.jpg', i+1);
%file_name = sprintf('frames/008.jpg', i);
%file_name2 = sprintf('frames/049.jpg', i+1);


I1=im2double(imread(file_name2));
I2=im2double(imread(file_name));

% Get the Key Points
Options.upright=true;
Options.tresh=0.0001;
Ipts1=OpenSurf(I1,Options);
Ipts2=OpenSurf(I2,Options);

% Put the landmark descriptors in a matrix
D1 = reshape([Ipts1.descriptor],64,[]);
D2 = reshape([Ipts2.descriptor],64,[]);

% Find the best matches
err=zeros(1,length(Ipts1));
cor1=1:length(Ipts1);
cor2=zeros(1,length(Ipts1));
for i=1:length(Ipts1),
    distance=sum((D2-repmat(D1(:,i),[1 length(Ipts2)])).^2,1);
    [err(i),cor2(i)]=min(distance);
end

% Sort matches on vector distance
[err, ind]=sort(err);
cor1=cor1(ind);
cor2=cor2(ind);

% Make vectors with the coordinates of the best matches
Pos1=[[Ipts1(cor1).y]',[Ipts1(cor1).x]'];
Pos2=[[Ipts2(cor2).y]',[Ipts2(cor2).x]'];
Pos1=Pos1(1:30,:);
Pos2=Pos2(1:30,:);

% Show both images
I = zeros([size(I1,1) size(I1,2)*2 size(I1,3)]);
I(:,1:size(I1,2),:)=I1; I(:,size(I1,2)+1:size(I1,2)+size(I2,2),:)=I2;

% Calculate affine matrix
Pos1(:,3)=1; Pos2(:,3)=1;
M=Pos1'/Pos2';

% Add subfunctions to Matlab Search path
functionname='OpenSurf.m';
functiondir=which(functionname);
functiondir=functiondir(1:end-length(functionname));
addpath([functiondir '/WarpFunctions'])

% Warp the image
I1_warped=affine_warp(I1,M,'bicubic');

% Show the result
%figure,
subplot(1,3,1), imshow(I1);title('Figure 1');
subplot(1,3,2), imshow(I2);title('Figure 2');
subplot(1,3,3), imshow(I1_warped);title('Warped Figure 1');
imwrite(I1_warped,file_name2);
if (mod(i,20)==0 ) 
    disp(sprintf('he make a %d',i)); 
end
end
sprintf('finish');
aaa();
end 

2 个答案:

答案 0 :(得分:1)

  1. 计算黑栏 - this question。只需删除列if 没有任何非零像素。
  2. 用行做同样的技巧。
  3. 完成剩下的工作。
  4. 如果黑色边框不是真的黑色,则会失败 - 例如,深灰色。在这种情况下,应用阈值检测黑色。

    如果此图片的主要部分中有任何黑色列,那么它也将是一个糟糕的解决方案。在这种情况下,您应该检查黑色列的位置,并仅删除相对靠近边界的列。

    这是删除史诗

    内没有黑色行的对称黑匣子的简单版本
    I1=im2double(imread('dart.jpg'));
    
    sizeI = size(I1);
    zeros = floor((sizeI(2) -  min(sum(any(I1))))/2);
    I2 = I1(:, zeros : sizeI(2)-zeros, :);
    nonZero = sum(any(I1,2));
    
    
    sizeI2 = size(I2);
    zerosRows = floor((sizeI(1) -  min(sum(any(I2, 2))))/2);
    I3 = I2(zerosRows : sizeI2(1)-zerosRows, :, :);
    
    subplot(1,3,1), imshow(I1);title('Figure 1');
    subplot(1,3,2), imshow(I2);title('Figure 2');
    subplot(1,3,3), imshow(I3);title('Figure 3');
    

    适用于" good"输入:

    enter image description here

    应用于内部黑线图像 - 效果不佳。 enter image description here

    如果你需要精确检测,这是一个计划:

    1. 使用颜色和方法检查从最左侧到第一列的列 删除所有黑色。
    2. 使用颜色检查从最右侧向后到最后一列的列并删除黑色
    3. 对行做同样的事。
    4. 我不会为这个提供代码,因为它只是OP自己可以做的一些矩阵操作。

答案 1 :(得分:0)

但请你回答不太好 因为你还有一点黑色边框...有人可以修复代码吗? 最后不会是黑色边框