img1 = zeros(1000,6000, 'uint8');
img1(:,:) = 255;

t=50;
step=.006;
for i=100:2900
    img1(250:750,i) = (127.5 + cos(2*pi/t * i)*127.5);
    t = t-step;
end

for i=3100:5900
    img1(250:750,i) = (127.5 + cos(2*pi/t * i)*127.5);
    t = t-step;
end

image = [img1(:,1:3000); img1(:,3001:6000)];
imshow(image);
imwrite(image,'spatial_res.bmp');