It mostly depends of what points should move (only neighbours or all of them), but after you decide it you can calculate the closest point in the line for each point and then move it lineary in that direction. Here is an example moving all points to a line defined by three points
clc; % Clear the command window.
close all; % Close all figures (except those of imtool.)
clear; % Erase all existing variables. Or clearvars if you want.
workspace; % Make sure the workspace panel is showing.
countmax=5; % units are in micrometres
view([-37.5 30])
X = zeros(1,countmax);
Y = zeros(1,countmax);
Z = zeros(1,countmax);
n=0;
m=0;
y=1;
d = 10;
%d = 15/1000; %distance between each molecule => 15/1000 is 15 micrometre
s = 6; % size of particle
for i = 2:countmax
for j = 2:countmax
for k = 2:countmax
X(i) = X(i-1) + d; %%1 = distance between each molecules (add an equation for it)
Y(j) = Y(j-1) + d;
Z(k) = Z(k-1) + d;
end
end
end
q= 1;
for n = 1:countmax
for m = 1:countmax
b = repelem(X(n),countmax);
c = repelem(Y(m),countmax);
a = [b;c;Z]; %% Check Coordinates
scatter3(a(1,:),a(2,:),a(3,:),s,'Filled','r') %% 4 value from left changes the size of each point
%xlim([d (countmax-1)*d]) % Axis limitsg
%ylim([d (countmax-1)*d])
%zlim([d (countmax-1)*d])
xlabel('My x label') % Axis Labels
ylabel('My y label')
zlabel('My z label')
hold on
rotate3d on
end
end
grid off;
Image = getframe();
K = sum(sum(rgb2gray(Image.cdata)==255));
percentageofwhite(y) = K/numel(rgb2gray(Image.cdata))*100;
q = percentageofwhite(y);
%delete(findall(gcf,'type','text'))
txt=text(0.0,0.95,sprintf('Water Volume = %0.3f%%',q(1,1)),'Units','normalized');
close all; % Close all figures (except those of imtool.)
clear; % Erase all existing variables. Or clearvars if you want.
workspace; % Make sure the workspace panel is showing.
countmax=5; % units are in micrometres
view([-37.5 30])
X = zeros(1,countmax);
Y = zeros(1,countmax);
Z = zeros(1,countmax);
n=0;
m=0;
y=1;
d = 10;
%d = 15/1000; %distance between each molecule => 15/1000 is 15 micrometre
s = 6; % size of particle
for i = 2:countmax
for j = 2:countmax
for k = 2:countmax
X(i) = X(i-1) + d; %%1 = distance between each molecules (add an equation for it)
Y(j) = Y(j-1) + d;
Z(k) = Z(k-1) + d;
end
end
end
q= 1;
[XX,YY,ZZ] = meshgrid(X,Y,Z);
AllDimensions = [XX(:),YY(:),ZZ(:)];
LinePoints = [30,10,20;
30,20,20;
30,30,20;
30,40,20];
% Find which points should be directed to where
LinePointMap = zeros(length(AllDimensions),1);
for idx=1:length(AllDimensions)
[~,idxMin] = min( sum(abs(AllDimensions(idx,:)-LinePoints).^2,2) );
LinePointMap(idx) = idxMin;
end
Differences = LinePoints(LinePointMap,:)-AllDimensions;
MaxStep = 50;
for idxAnimation = 1:MaxStep
scatter3(AllDimensions(:,1)+idxAnimation/MaxStep*Differences(:,1),...
AllDimensions(:,2)+idxAnimation/MaxStep*Differences(:,2),...
AllDimensions(:,3)+idxAnimation/MaxStep*Differences(:,3),...
s,'Filled','r') %% 4 value from left changes the size of each point
xlim([d (countmax-1)*d]) % Axis limitsg
ylim([d (countmax-1)*d])
zlim([d (countmax-1)*d])
xlabel('My x label') % Axis Labels
ylabel('My y label')
zlabel('My z label')
hold off
rotate3d on
grid off;
pause(0.1)
end
Image = getframe();
K = sum(sum(rgb2gray(Image.cdata)==255));
percentageofwhite(y) = K/numel(rgb2gray(Image.cdata))*100;
q = percentageofwhite(y);
%delete(findall(gcf,'type','text'))
txt=text(0.0,0.95,sprintf('Water Volume = %0.3f%%',q(1,1)),'Units','normalized');
