Hi,
I have a rectangular cuboid in which I have a cloud of points, each having a mass assigned to it. Input file ('file_%d.txt') provides the following data regarding these points:
column 1: point type (1, 2, 3 …)
column 2: x coordinate
column 3: y coordinate
column 4: z coordinate
column 5: mass of the point
There are over 10k points listed in a file.
What I want to do is to go along z direction which is 100cm long, cut it into 1cm thick slices (along z, so 100 slices in total, in the code they are named regions) and for each slice calculate the total mass of the points which are in that specific slice. Also, I have 3 input files that contain these points and I need to obtain average mass for each slice over these 3 files. Then I would like to export that average mass over each slice as a function of the distance from 0 to 100cm (which corresponds to slice numbers) to either a .txt or .xlsx file.
Last time I did some programming in Matlab was ages ago – the code below is needed for data processing. I got stuck with the for & if loops and would be super grateful for your help guys.
mini=1; inc=1; maxi=3;
regions = zeros(1,100)
dr=1
c=0;
for i=1:length(mini:inc:maxi)
c=c+1;
str = sprintf('file_%d.txt',mini+(i-1)*inc);
temp = importdata(str,' ',9);
data = sortrows(temp.data)
var_bottom=0
var_top=1
fork=1:length(data)
var_bottom=var_bottom+1
var_top=var_top+1
if data(j,4)=<dr*var_top && if data(j,4)>dr*var_bottom
total(k+1)=sum(data(j,4))
???????? – not sure if the for loop above makes sense and, on top of that, my use of && was deemed “illegal” by Matlab. Also, not sure how I should code the rest of what I need
end
fclose(fileID);
c
regions = regions/c;
regions = regions/dr
x = 1:dr:length(regions)
plot(mass,regions, '-.')
for k = 1 : length(mini:inc:maxi)
filename = sprintf('position_vs_mass_%d.xlsx', mini+(i-1)*inc);
xlswrite(filename, [mass(:),regions(:)]);
end
