% for j = 1:4
    load 2stl;
%     a = [num2str(j),'.mat'];
%     load(a);
    gridX = 1:1500;
    gridY = 1:500;
    gridZ = 1:2;

    gridDATA = zeros(500,1500,2);
    globalx(globalx>0.5) = 1;
    globalx(globalx<1) = 0;

    figure; colormap(gray); imagesc(1-globalx);
    caxis([0 1]); axis equal; axis off; drawnow;
    
    for i = gridZ        
        gridDATA(:,:,i) = globalx;
    end
    
    output = 'v_small.stl';
%     output = [num2str(j),'.stl'];
    CONVERT_voxels_to_stl(output,gridDATA,gridY,gridX,gridZ,'binary');
% end


%     load whole;
%     gridX = 1:600;
%     gridY = 1:200;
%     gridZ = 1:20;`    `                                                                       
%     gridDATA = zeros(200,600,20);
%     whole(whole>0.1) = 1;
%     whole(whole<1) = 0;
%     
%     for i = gridZ        
%         gridDATA(:,:,i) = whole;
%     end
%     
%     CONVERT_voxels_to_stl('a.stl',gridDATA,gridY,gridX,gridZ,'binary');


% gridX = 1:600;
% gridY = 1:200;
% gridZ = 1:200;
% load density;
% density(density>0.5) = 1;
% density(density<1) = 0;
% gridDATA = zeros(600,200,200);
% d = 20;
% for i = 1:600
%     for j = 1:200
%         for k = 1:200
%             dis = round(sqrt((j-99.5)^2+(k-99.5)^2));
%             if dis<=100
%                 % area left
%                 if -d/100*j+100+d<k
%                     if d/100*j+100-d>k
%                         gridDATA(i,j,k) = density(i,dis);
%                     end
%                 end
%
%                 % area right
%                 if -d/100*j+100+d>k
%                     if d/100*j+100-d<k
%                         gridDATA(i,j,k) = density(i,dis);
%                     end
%                 end
%
%                 % area up
%                 if -100/d*j+100/d*(100+d)<k
%                     if 100/d*j-100/d*(100-d)<k
%                         gridDATA(i,j,k) = density(i,dis);
%                     end
%                 end
%                 % area down
%                 if -100/d*j+100/d*(100+d)>k
%                     if 100/d*j-100/d*(100-d)>k
%                         gridDATA(i,j,k) = density(i,dis);
%                     end
%                 end
%
% %                 gridDATA(i,j,k) = density(i,dis);
%             end
%         end
%     end
% end
% CONVERT_voxels_to_stl('a.stl',gridDATA,gridX,gridY,gridZ,'binary');