Я пытаюсь преобразовать приведенный ниже код Matlab в C++, используя codegen. Однако он не работает при сборке, и я получаю ошибку:
"Если не указано" rows ", первый вход должен быть вектором. Если вектор является переменным размером, либо первое измерение, либо второе должно иметь фиксированную длину 1. Вход [] не поддерживается Используйте вход 1-на-0 или 0-на-1 (например, нули (1,0) или нули (0,1)), чтобы представить пустой набор. "
Затем он указывает на [id, m, n] = unique (id); являясь виновником. Почему он не строится и как лучше всего это исправить?
function [L,num,sz] = label(I,n) %#codegen
% Check input arguments
error(nargchk(1,2,nargin));
if nargin==1, n=8; end
assert(ndims(I)==2,'The input I must be a 2-D array')
sizI = size(I);
id = reshape(1:prod(sizI),sizI);
sz = ones(sizI);
% Indexes of the adjacent pixels
vec = @(x) x(:);
if n==4 % 4-connected neighborhood
idx1 = [vec(id(:,1:end-1)); vec(id(1:end-1,:))];
idx2 = [vec(id(:,2:end)); vec(id(2:end,:))];
elseif n==8 % 8-connected neighborhood
idx1 = [vec(id(:,1:end-1)); vec(id(1:end-1,:))];
idx2 = [vec(id(:,2:end)); vec(id(2:end,:))];
idx1 = [idx1; vec(id(1:end-1,1:end-1)); vec(id(2:end,1:end-1))];
idx2 = [idx2; vec(id(2:end,2:end)); vec(id(1:end-1,2:end))];
else
error('The second input argument must be either 4 or 8.')
end
% Create the groups and merge them (Union/Find Algorithm)
for k = 1:length(idx1)
root1 = idx1(k);
root2 = idx2(k);
while root1~=id(root1)
id(root1) = id(id(root1));
root1 = id(root1);
end
while root2~=id(root2)
id(root2) = id(id(root2));
root2 = id(root2);
end
if root1==root2, continue, end
% (The two pixels belong to the same group)
N1 = sz(root1); % size of the group belonging to root1
N2 = sz(root2); % size of the group belonging to root2
if I(root1)==I(root2) % then merge the two groups
if N1 < N2
id(root1) = root2;
sz(root2) = N1+N2;
else
id(root2) = root1;
sz(root1) = N1+N2;
end
end
end
while 1
id0 = id;
id = id(id);
if isequal(id0,id), break, end
end
sz = sz(id);
% Label matrix
isNaNI = isnan(I);
id(isNaNI) = NaN;
[id,m,n] = unique(id);
I = 1:length(id);
L = reshape(I(n),sizI);
L(isNaNI) = 0;
if nargout>1, num = nnz(~isnan(id)); end
Удалите анонимную функцию, сохраненную в переменной vec
и сделайте vec
подфункцией:
function y = vec(x)
coder.inline('always');
y = x(:);
Без опции 'rows'
вход в unique
функцию всегда интерпретируется как вектор, и выход всегда является вектором. Так, например, что-то вроде id = unique(id)
будет иметь эффект id = id(:)
если все элементы id
матрицы были уникальными. Нет никакого вреда в том, чтобы вводить вектор. Поэтому измените строку
[id,m,n] = unique(id);
в
[id,m,n] = unique(id(:));
Просто FYI, если вы используете MATLAB R2013b или новее, вы можете заменить error(nargchk(1,2,nargin))
с narginchk(1,2)
.
Как говорится в сообщении об ошибке, для unique
codegen требуется, чтобы вход был вектором, если не передано "rows".
Если вы посмотрите на отчет (нажмите ссылку "Открыть отчет") и наведите указатель мыши на идентификатор, вы, вероятно, увидите, что его размер не является ни 1-by-N
ни N-by-1
. Требование к unique
можно увидеть, если вы ищете unique
здесь:
Вы можете сделать одно из нескольких вещей:
Сделайте id
вектором и рассматривайте его как вектор для вычисления. Вместо декларации:
id = reshape(1:prod(sizI),sizI);
вы можете использовать:
id = 1:numel(I)
Тогда id
будет вектором строки.
Вы также можете сохранить код как есть и сделать что-то вроде:
[idtemp,m,n] = unique(id(:));
id = reshape(idtemp,size(id));
Очевидно, что это приведет к созданию копии, idtemp
, но это может привести к меньшему изменению кода.