clc,clear all;clf
NX = 2;
NY = 2;
NZ = 2;
fx = 0;
fy = 0;
fz = 0;
fname = 'Cube_2x2x2_fx-0_fy-0_fz-0_1';
%fnsize = ([num2str(NX),'x',num2str(NY),'x',num2str(NZ)]);
fnrandp = (['fx-',num2str(fx),'_','fy-',num2str(fy),'_','fz-',num2str(fz),'.emf']);
%fname = (['Cube_',fnsize,'_',fnrandp,'_1']);

[x,y,z]=meshgrid(1:NX,1:NY,1:NZ);

r1=fx*rand(NY,NX,NZ);
r2=fy*rand(NY,NX,NZ);
r3=fz*rand(NY,NX,NZ);

xr=x+r1%-.8;
yr=y+r2%-.8;
zr=z+r3%-.8;

figure('Units', 'pixels', ...
    'Position', [100 10 2000 500]);
hold on;

for i = 1:length(z(1,:,:))
subplot(1,3,1); plot3(xr(:,:,i),yr(:,:,i),zr(:,:,i),'ko','MarkerSize',6,...
'MarkerEdgeColor',[0 0 0],'MarkerFaceColor',[0 0 0]);view([12 55]);
hold on; title('Isometric','FontSize',20,'FontName','Helvetica','FontWeight','bold');
set(gca,'Box','on','LineWidth',2,'FontSize',16,'FontName','Helvetica',...
    'XTick',0:1:NX,'YTick',0:1:NY,'ZTick',[0 NZ])
xlabel('X-Axis','FontSize',16,'FontName','Helvetica','FontWeight','bold');
ylabel('Y-Axis','FontSize',16,'FontName','Helvetica','FontWeight','bold');
zlabel('Z-Axis','FontSize',16,'FontName','Helvetica','FontWeight','bold');
xlim([0 NX]);ylim([0 NY]);zlim([0 NZ]);
%
subplot(1,3,2); plot3(xr(:,:,i),yr(:,:,i),zr(:,:,i),'ko','MarkerSize',6,...
'MarkerEdgeColor',[0 0 0],'MarkerFaceColor',[0 0 0]);view([90 0]);
hold on;title('View Down X-Axis','FontSize',20,'FontName','Helvetica','FontWeight','bold');
set(gca,'Box','on','LineWidth',2,'FontSize',16,'FontName','Helvetica',...
    'XTick',0:1:NX,'YTick',0:1:NY,'ZTick',0:1:NZ)
xlabel('X-Axis','FontSize',16,'FontName','Helvetica','FontWeight','bold');
ylabel('Y-Axis','FontSize',16,'FontName','Helvetica','FontWeight','bold');
zlabel('Z-Axis','FontSize',16,'FontName','Helvetica','FontWeight','bold');
xlim([0 NX]);ylim([0 NY]);zlim([0 NZ]);
%
subplot(1,3,3); plot3(xr(:,:,i),yr(:,:,i),zr(:,:,i),'ko','MarkerSize',6,...
'MarkerEdgeColor',[0 0 0],'MarkerFaceColor',[0 0 0]);view([90 90]);
hold on;title('View Down Z-Axis','FontSize',20,'FontName','Helvetica','FontWeight','bold');
set(gca,'Box','on','LineWidth',2,'FontSize',16,'FontName','Helvetica',...
    'XTick',0:1:NX,'YTick',0:1:NY,'ZTick',0:1:NZ)
xlabel('X-Axis','FontSize',16,'FontName','Helvetica','FontWeight','bold');
ylabel('Y-Axis','FontSize',16,'FontName','Helvetica','FontWeight','bold');
zlabel('Z-Axis','FontSize',16,'FontName','Helvetica','FontWeight','bold');
xlim([0 NX]);ylim([0 NY]);zlim([0 NZ]);
end

%cd('C:\Users\Colin\Desktop\Point Clouds\Simple Cubic');
%print -dmeta '2x2x2_fx-0_fy-0_fz-0.emf'

%view([90 90]); 



m = 1;
for i = 1:length(zr(:,1,1))
    for j = 1:length(yr(1,:,1))
        for k = 1:length(xr(1,1,:))
            
            p(m,:) = [m xr(i,j,k) yr(i,j,k) zr(i,j,k)];
            m = m+1;
        end
    end
end

cd('C:\Users\Colin\Desktop\Point Clouds\Simple Cubic')

%dlmwrite(fname,p,'delimiter','\t','newline','pc');

%axis([0,N+1,0,N+1])