% make aprox 500x500 pixels

% % perfect lattices
% image = zeros(10);
% image(10,10) = 1;
% image = repmat(image, 50,50);
% size(image)
% 
% imwrite(image, '/home/wrightc/matlab/dimages/perfect1010' , 'tiff');
% 
% image = zeros(10,20);
% image(10,20) = 1;
% image = repmat(image, 50,25);
% size(image)
% 
% imwrite(image, '/home/wrightc/matlab/dimages/perfect1020' , 'tiff');
% 
% image = zeros(20,10);
% image(20,10) = 1;
% image = repmat(image, 25,50);
% size(image)
% 
% imwrite(image, '/home/wrightc/matlab/dimages/perfect2010' , 'tiff');
% 
% image = zeros(20);
% image(20,20) = 1;
% image = repmat(image, 25,25);
% size(image)
% 
% imwrite(image, '/home/wrightc/matlab/dimages/perfect2020' , 'tiff');
% 
% image = zeros(20,40);
% image(20,40) = 1;
% image = repmat(image, 25,13);
% size(image)
% 
% imwrite(image, '/home/wrightc/matlab/dimages/perfect2040' , 'tiff');
% 
% image = zeros(40,20);
% image(40,20) = 1;
% image = repmat(image, 13,25);
% size(image)
% 
% imwrite(image, '/home/wrightc/matlab/dimages/perfect2040' , 'tiff');

% % staggered perfect lattices
% image = zeros(20);
% image(20,20) = 1;
% image(10,10) = 1;
% image = repmat(image, 25,25);
% size(image)
% 
% imwrite(image, '/home/wrightc/matlab/dimages/stagg1010' , 'tiff');
% 
% image = zeros(20,10);
% image(20,10) = 1;
% image(10,5) = 1;
% image = repmat(image, 25,50);
% size(image)
% 
% imwrite(image, '/home/wrightc/matlab/dimages/stagg1005' , 'tiff');
% 
% image = zeros(20,10)';
% image(10,20) = 1;
% image(5,10) = 1;
% image = repmat(image, 50,25);
% size(image)
% 
% imwrite(image, '/home/wrightc/matlab/dimages/stagg0510' , 'tiff');
% 
% image = zeros(20,40);
% image(20,40) = 1;
% image(10,20) = 1;
% image = repmat(image, 25,13);
% size(image)
% 
% imwrite(image, '/home/wrightc/matlab/dimages/stagg1020' , 'tiff');
% 
% image = zeros(20,40)';
% image(40,20) = 1;
% image(20,10) = 1;
% image = repmat(image, 13,25);
% size(image)
% 
% imwrite(image, '/home/wrightc/matlab/dimages/stagg2010' , 'tiff');

% add in some error
image = zeros(40,20);
image = repmat(image, 13,25);
size(image)
for i = 0:13
    h = 40*i;
    for j = 0:25
        r=3;
        p=.4;
        if(rand<p)
            i1 = h + 30 + floor(r*randn());
            j1 = 20*j + 15 + floor(r*randn());
            image(i1,j1) = 1;
        end
        if(rand<p)
            i1 =abs( h + 10 + floor(r*randn()));
            j1 = abs(20*j + 5 + floor(r*randn()));
            if(j1==0)
                j1=1;
            end
            image(i1,j1) = 1;
        end
    end
end

imwrite(image, '/home/wrightc/matlab/dimages/error2010_3.4' , 'tiff');