读取带有NA的textscan文件作为空matlab

Question

我有一个包含 n 行和 5 列的文本文件 D 我想阅读。每行将代表一个向量。但是，如果值等于 NA，我想将其留空：

D:

122 12 10 NA NA
110 10 30 45 87 
110 12 10 NA NA
..

我正在这样做：

bid = fopen(D,'r');
bfile = textscan(bid,'%d %d %d %s %s','TreatAsEmpty',{'NA'})
var=122;
IndVar=bfile {:,1}==var;
rest=bfile(:,[2:end])
vecVar=rest(IndVar)

我想从没有 NA 的行中获取向量： vecVar = [12 10]； [10 30 45 87]； [12 10]；

我试试这个：

rest= rest(~isnan(rest));

我得到错误：

'cell' 类型的输入参数的未定义函数 'isnan'。

我已经玩了一个多小时了仍然无法弄清楚.. 任何帮助将不胜感激！

非常感谢！

Answer 1

如错误消息中所述，isnan 仅适用于数值数组。要对数值数组的元胞数组进行操作，请使用cellfun，它将对数组中的每个元胞应用一个命令。顺便说一句，NaN 内部实际上是一个数字（单数或双数），因此您将希望使用数字版本的 NaN 而不是文本字符串 'NA' 为您的空值，@987654324 @ 是一个内置命令，因此您需要更改该变量名称。

您需要稍微澄清一下，但据我了解，给定包含矩阵的输入文件：

122 12 10 NA NA
110 10 30 45 87 
110 12 10 NA NA

你期望：

vecVar ={...
    [12 10];
    [10 30 45 87 ];
    [12 10]
    };

在我的第一个答案中，我错过了您正在阅读整数的事实。这将给出上面的结果，我认为这就是你所追求的：

bid = fopen(d,'r');
% integers don't have NaN, so read in as float instead.
% also, collect everything into a single matrix because it makes it.
% easier to debug
bfile = textscan(bid,'%f %f %f %f, %f', 'TreatAsEmpty', 'NA', 'CollectOutput', true);
fclose(D);

% convert from a cell to a matrix to make i
bfile = bfile{1};
var=122;
IndVar = bfile(:,1) == var;
rest = bfile(:, 2:end);

% Turn the matrix back into a cell with each row a single cell
[nrows, ncols] = size(rest);
rest = mat2cell(rest, ones(nrows, 1), ncols);

% get rid of the NaN values
rest = cellfun(@(x) x(~isnan(x)), rest, 'UniformOutput', false);

这将生成一个元胞数组，其中每个元素都包含一个可以使用desired_vector = rest{index_of_your_choice};检索的向量

Answer 2

以下是在元胞数组val 中放入每行数值的代码：

    %'Number of columns'
    NC = 5;

    %'Read the strings'
    f = fopen('test.txt', 'r');
    raw= textscan(f, repmat('%s',1,NC));
    fclose(f);
    if numel(raw) ~= NC
            error('Something is wrong with the file format');
    else
            NR = numel(raw{1});
    end;

    %'Convert data'
    for c = 1:NC
            for r = 1:NR
                    raw{c}{r} = sscanf(raw{c}{r}, '%d');
            end;
    end;

    %'Concatenate rows into val'
    val = cell(1,NR);
    buf = cell(1,NC);
    for r = 1:NR
            for c = 1:NC
                    buf{c} = raw{c}{r};
            end;
            val{r} = [buf{:}];
    end;

您不能使用普通矩阵来存储值，因为存在具有不同数量值的行。这意味着行向量只能存储到元胞数组中。