ansaurus

Question

How to read data in chunks from notepad file in Matlab?

Answer 1

+4 A:

First of all, I suggest you don't use variable names such as TIME1,TIME2 etc, since that gets messy quickly. Instead, you can e.g. use a cell array with five rows (one for each well), and one or two columns. In the sample code below, wellData{2,1} is the time for the second well, wellData{2,2} is the corresponding Oil Rate SC - Yearly.

There might be more elegant ways to do the reading; here's something quick:

%# open the file
fid = fopen('Reportq.rwo');

%# read it into one big array, row by row
fileContents = textscan(fid,'%s','Delimiter','\n');
fileContents = fileContents{1};
fclose(fid); %# don't forget to close the file again

%# find rows containing TABLE NUMBER
wellStarts = strmatch('TABLE NUMBER',fileContents);
nWells = length(wellStarts);

%# loop through the wells and read the numeric data
wellData = cell(nWells,2);
wellStarts = [wellStarts;length(fileContents)];

for w = 1:nWells 
    %# read lines containing numbers
    tmp = fileContents(wellStarts(w)+5:wellStarts(w+1)-1);
    %# convert strings to numbers
    tmp = cellfun(@str2num,tmp,'uniformOutput',false);
    %# catenate array
    tmp = cat(1,tmp{:});
    %# assign output
    wellData(w,:) = mat2cell(tmp,size(tmp,1),[1,1]);
end

Jonas 2010-08-25 01:01:26

Thanks Jonas. This works good.

Harpreet 2010-08-25 18:12:45

Answer 2

A:

I used the readline.m code as a function (available to download for free) to keep my over all code less complex. Here's what I did to achieve my objective:

well=readline('Reportq.rwo',[1:130]);
wellData=cell2mat(cellfun(@str2num,well,'UniformOutput', false));
q_yearly1(:,j)=wellData(1:21,2);
q_yearly2(:,j)=wellData(22:42,2);
q_yearly3(:,j)=wellData(43:63,2);
q_yearly4(:,j)=wellData(64:84,2);
q_yearly5(:,j)=wellData(85:105,2);

It works pretty well, and the code readline.m is very general. Here's the link to download the readline.m code.

Harpreet 2010-08-25 18:12:15

ansaurus

tags:

views:

answers:

How to read data in chunks from notepad file in Matlab?

related questions