i am using following code to find area of layers. I am trying but masking somewhere some mistake. Any help from community will be highly appreciated.
Attached file are the xyz-coordinates. I measure the surface area of one layer using bellow code. This coding is combining all layers. I want to find area of each layer individually. Regards to all for their cooperation and guidance.
vertices = load(' attached file ')vertices = round(sortrows(vertices,2));Vx = round(vertices(:,1));Vy = round(vertices(:,2));Vz = round(vertices(:,3));G_surface_area =0layer = []surface_area = 0for i = min(Vy):10:max(Vy) % This for loop, i am using to measure by 10 spacing.
minVy = i for j = i:i + 9 % This measure is for layer to find its area.
layer_new = vertices(any(Vy== j,2),:); layer = [layer; {layer_new}]; end layer = cell2mat(layer) Lx = layer(:,1); Ly = layer(:,2); Lz = layer(:,3); surface_area = polyarea(Lx,Lz) G_surface_area = G_surface_area + surface_area end
Best Answer