I have a text file with time "T2-distribution" and Amplitude at each depth.X-axis are my values adjacent to 'T2' row only. I want area under curve between T2 and amplitude at the respective 'DEPTH' (ignore values adjacent to 'DEPTH'- only in row 4)?

Question

KRISHNA NISHCAL BHARATULA 2017-3-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/330343-i-have-a-text-file-with-time-t2-distribution-and-amplitude-at-each-depth-x-axis-are-my-values-adja

评论： KRISHNA NISHCAL BHARATULA 2017-3-17

Hole: 123A
Top: 78.67
T2  3.34   4.26   8.16   9.84  10.01
DEPTH   0      1      2    3     4
78.67   3.38  9.28  7.86  0.00  0.00
80.52   4.68  9.27  0.00  0.00  1.24
91.50   6.66  2.34  1.23  0.00  0.00
95.67   7.88  6.88  5.43  0.00  2.34
100.01  5.44  6.54  8.23  5.01  0.00

2 个评论
显示无隐藏无

KSSV 2017-3-17

Isn't the same? https://in.mathworks.com/matlabcentral/answers/327418-i-have-a-following-text-file-delimited-with-tabs-that-contains-both-strings-and-decimal-numbers

KRISHNA NISHCAL BHARATULA 2017-3-17

yes it is. But I need to find the area under each of those curves which gives 'porosity of rock' at a particular depth. I appreciate your quick response, thank you very much.

:)

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

KSSV 2017-3-17

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/330343-i-have-a-text-file-with-time-t2-distribution-and-amplitude-at-each-depth-x-axis-are-my-values-adja#answer_259119

在 MATLAB Online 中打开

fid = fopen('your txt file') ;
S = textscan(fid,'%s','delimiter','\n') ;
fclose(fid) ;
S = S{1} ;
%%Get positions
idx1 = find(not(cellfun('isempty',strfind(S, 'T2')))); % T2 position 
idx2 = find(not(cellfun('isempty',strfind(S, 'DEPTH')))); % T2 position 
idx3 = idx2:length(S) ;
%%Get the required 
x = cell2mat(cellfun(@str2num,strsplit(S{idx1}),'un',0)) ;
depth = cell2mat(cellfun(@str2num,strsplit(S{idx2}),'un',0)) ;
data = cell2mat(cellfun(@str2num,S(idx3),'un',0)) ;
y = data(:,1) ;
data = data(:,2:end) ;
%%each area 
N = length(depth) ;
Area = zeros(N,1) ;
for i = 1:N
    Area(i) = trapz(x,data(i,:)) ;
end