[BDim,BD]=N_bandtrace(BW)

function [BDim,BD]=N_bandtrace(BW)
%α߽
BD = [];
BWT_tm = BW;
BDim=repmat(logical(uint8(0)),size(BW,1),size(BW,2));

%ǰһθٺб߽
while sum(sum(BWT_tm)) ~= 0
    [BD_tm,BWT_tm] = bandtrace(BWT_tm);
    if length(BD) < length(BD_tm)
        BD = BD_tm;
    end
end
%ɨ
for i = 1 : length(BD)
    BDim(BD(i,1),BD(i,2)) = 1;
end