Looping through a 3D matrix

95 views (last 30 days)
Robert
Robert on 4 Nov 2015
Commented: TastyPastry on 5 Nov 2015
I have a 3D matrix 35,43,144
I want to run the following code such that it performs the function element by element through the 3rd dimension (using 144 data points) I want this to be done for every pixel giving me 1505 outputs (35x43)
I just can't seem to figure out the correct loop inputs
p = load ('JPLCSRmean.mat')
GRACE = p.datamean
ncol = size(GRACE,2);
for k = 1:144
j = 1:43
i = 1:35
[c,l] = wavedec(GRACE(i,j,k),4,'dmey');
al = appcoef(c,l,'dmey',1);
%[cD1,cD2,cD3,cD4,cD5,cD6] = detcoef(c,l,[1,2,3,4,5,6]);
D = [];
%D(:,6)= wrcoef('a',c,l,'dmey',6);
%D(:,5) = wrcoef('a',c,l,'dmey',5);
D(:,4) = wrcoef('a',c,l,'dmey',4);
D(:,3) = wrcoef('a',c,l,'dmey',3);
D(:,2) = wrcoef('a',c,l,'dmey',2);
D(:,1) = wrcoef('a',c,l,'dmey',1);
filename = sprintf('%s_%04f.txt', 'GRACE_APROXIMATIONS', k)
csvwrite(filename,D)
end

Answers (1)

TastyPastry
TastyPastry on 4 Nov 2015
for i = 1:144
for j = 1:43
for k = 1:35
myVal = myData(i,j,k);
%do something
end
end
end
  2 Comments
TastyPastry
TastyPastry on 5 Nov 2015
You could declare a tracking variable outside the entire loop and increment it each time you print a file.

Sign in to comment.

Categories

Find more on Loops and Conditional Statements in Help Center and File Exchange

Tags

Products

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!