floatdata2d Ima,Ima2, Kernel;

ucdata2d ImaUc;

SetDim(Ima, 512, 512);

SetDim(Ima2, 512, 512);

@@

for (i = 200 : 1 : 400, j = 200 : 1 : 400)

{

Ima[i][j] = 128;

}

@@

ImaUc = Ima;

Save(Assign.bmp, gray, ImaUc);

@@

for (i = 0 : 1 : 512, j = 0 : 1 : 512)

if(sqrt((i-255)^2 + (j-255)^2) <30)

{

Ima2[i][j] = 70;

}

else if(sqrt((i-255)^2 + (j-255)^2) <60)

{

Ima2[i][j] = 120;

}

else if(sqrt((i-255)^2 + (j-255)^2) <80)

{

Ima2[i][j] = 170;

}

else if(sqrt((i-255)^2 + (j-255)^2) <100)

{

Ima2[i][j] = 220;

}

else

{

Ima2[i][j] = 255;

}

@@

ImaUc = Ima2;

Save(Assign1.bmp, gray, ImaUc);