在Matlab中运行长度编码
我很新的MatLab，我有运行长度编码代码，但它似乎不工作，你能帮我吗？
我有这个input：
ChainCode = 11012321170701000700000700766666666666665555555544443344444333221322222322 我想把它变成RLE输出：
(1,2), (0,1), (1,1), (2,1), (3,1), (2,1), (1,2), (7,1), (0,1), (7,1), (0,1), (1,1), (0,3), (7,1), (0,5), (7,1), (0,2), (7,1), (6,13), (5,8), (4,4), (3,2), (4,5), (3,3), (2,2), (1,1), (3,1), (2,5), (3,1), (2,2)
这是我的代码：
lengthcode = 1; N = 1; for i = 2:length(ChainCode) if x(i)==x(i-1) N = N + 1; valuecode(N) = x(i); lengthcode(N) = lengthcode(N) + 1; else N = 1; lengthcode = 1; end i = i + 1; end
但是这不起作用，我仍然对如何打印输出感到困惑。
我希望你能帮助我。 谢谢。
这是一个没有循环，cellfun或arrayfun的紧凑解决scheme：
chainCode = '11012321170701000700000700766666666666665555555544443344444333221322222322'; numCode = chainCode - '0'; % turn to numerical array J=find(diff([numCode(1)-1, numCode])); relMat=[numCode(J); diff([J, numel(numCode)+1])];
通过坚持原来的实施，以下简单的变化应该工作。
chainCode = '11012321170701000700000700766666666666665555555544443344444333221322222322'; numCode = chainCode - '0'; % turn to numerical array relMat = []; numCode = [numCode nan]; % dummy ending N = 1; for i = 1:length(numCode)-1 if numCode(i)==numCode(i+1) N = N + 1; else valuecode = numCode(i); lengthcode = N; relMat = [relMat; valuecode lengthcode]; N = 1; end end
你可以格式化你喜欢的输出。 例如作为一个序列：
relMatT = relMat'; relSeq = relMatT(:)';
或者将string格式化为build议的输出：
relString = []; for i = 1:length(relMat) relString = [relString, sprintf('(%d, %d), ', relMat(i,1), relMat(i,2))]; end
作为扩展，如果您的源序列中包含字母数字，则应修改上述内容以便比较string而不是数字。
更新 ：要计算原始relMat中唯一代码对的出现次数， relMat尝试查找对并逐行计算零差异。 例如：
relMatUnique = unique(relMat, 'rows'); % find unique pairs nPairs = length(relMatUnique); nOccur = zeros(nPairs, 1); for i = 1:nPairs pairInMat = bsxfun(@minus, relMat, relMatUnique(i,:)); % find pair in relMat nOccur(i) = sum(~sum(pairInMat, 2)); end relMatOccur = [relMatUnique nOccur]; % unique pairs and number of occurrences
% Convert string to numeric array ChainCodeString = '11012321170701000700000700766666666666665555555544443344444333221322222322'; ChainCodeArray = ChainCodeString - '0'; % Initialize CurrentRleValue = ChainCodeArray(1); CurrentRleCount = 1; RleCodeIndex = 1; for i = 2 : length(ChainCodeArray) if ChainCodeArray(i)==ChainCodeArray(i-1) % Increment current run-length count CurrentRleCount = CurrentRleCount + 1; else % Store current run-length valuecode(RleCodeIndex) = CurrentRleValue; lengthcode(RleCodeIndex) = CurrentRleCount; RleCodeIndex = RleCodeIndex + 1; % Initialize next run-length CurrentRleValue = ChainCodeArray(i); CurrentRleCount = 1; end; end;
你可以避免for循环：
chainCode = '11012321170701000700000700766666666666665555555544443344444333221322222322'; numCode = chainCode - '0'; % turn to numerical array % detect edges (changes) edges = arrayfun( @(x,y) x ~= y, ... numCode(1:end-1), ... numCode(2:end)); % get indexes idx = find(edges); % create tuples relMat = cell2mat(arrayfun( ... @(b,e) [ numCode(b) ; e-b+1 ], ... [ 1 (idx + 1) ], ... [ idx length(numCode) ], ... 'UniformOutput', false));
