how to apply the arithmetic mean filter to a medical image to improve it? I have wrote this code but it did not work, there is an error which I could not figure it out.

%arithmatic mean filter
im=imread('chest.tif');%loading image
figure,imshow(im);
title('original');
[row col]=size(im);%storing size of image
for i=1:1:row-2;%sweeping through rows
    for j=1:1:col-2;%sweeping through columns
        for u=1:1:2;%sweeping through window
            for v=1:1:2;
                %taking arithmetic mean
                im(i,j)=im(i,j)+im(i+u,j+v);
                im(i,j)=im(i,j)/9;
            end
        end       
    end
end
%showing final image
figure,imshow(im);
title('arthmean');

[madhan ravi]

I formatted your code today , next time please select the whole code and press the code button so that it's easy to read.

[Jan]

@Lama: If you mention "there is an error which I could not figure it out" please share any details with the readers: Do you get an error message or do the results differ from your expectations?

[Lama Awawdeh]

ok. thanks for your comment

[Lama Awawdeh]

@jan there is an error message

[VISHNU VARDHAN]

there is an error in the above code.

[VISHNU VARDHAN]

%arithmatic mean filter

clc;

clear all;

close all;

x=imread('1.jpg');

x1=imresize(x,[256,256]);

figure;imshow(x1)

title('original image')%loading image.

x2=rgb2gray(x1);

figure;imshow(x2);title('gray image');

[m,n]=size(x2);%storing size of image.

for i=1:m-3

for j=1:n-3

a=x(i:i+3,j:j+3);

v(i,j)=sum(sum(a));

end

end

y=mat2gray(v);

figure;imshow(y);

title('average image');

this is useful for average of the any hd image.if u want to increase size,change pixcels in imresize command...................

[Jahid Hasan]

Is this loop similar to finding mean as well, how about median, what way I can write a Mexican function? Any suggestions

Answer 1

[KALYAN ACHARJYA]

Avoid multiple for loops, use inbuilt imfilter function for masking operation.

The concept is same in all images. For image enhancement you can perform numerous operation depends on input image.

%arithmetic mean filter
im=rgb2gray(imread('chest.tif'));
h=fspecial('average',3);
filter_image=imfilter(im,h);

[Lama Awawdeh]

according to the task, it is not allow for us,student, to use built in function. that is why i used for loop.

[KALYAN ACHARJYA]

No Error in your code

%arithmetic mean filter
im=rgb2gray(imread('test1.jpg'));
[row col]=size(im);%storing size of image
for i=1:1:row-2;%sweeping through rows
    for j=1:1:col-2;%sweeping through columns
        for u=1:1:2;%sweeping through window
            for v=1:1:2;
                %taking arithmetic mean
                im(i,j)=im(i,j)+im(i+u,j+v);
                im(i,j)=im(i,j)/9;
            end
        end       
    end
end
%showing final image
figure,imshow(im);
title('arthmean');

[KALYAN ACHARJYA]

Replace this

im=imread('chest.tif')

with

im=rgb2gray(imread('chest.tif'))

[Lama Awawdeh]

then whay this message is shown ?

and attach is the image that i must use

Error using rgb2gray>parse_inputs (line 77)
MAP must be a m x 3 array.

Error in rgb2gray (line 52)
[X, threeD] = parse_inputs(X);

[KALYAN ACHARJYA]

The image is not attched properly, attach the image again.

[Lama Awawdeh]

it is a very noisy image of the chest.

[KALYAN ACHARJYA]

tif image format is not supported to attach here, use other third party website, like googe drive, dropbox etc

[Lama Awawdeh]

https://drive.google.com/file/d/1A_gkswca76QgzHMfyMTNfwk0D-hxMoXP/view?usp=sharing

[KALYAN ACHARJYA]

No error shown here

clc;
clear all;
close all;
im=imread('chest.tif');
[row col]=size(im);%storing size of image
for i=1:1:row-2;%sweeping through rows
    for j=1:1:col-2;%sweeping through columns
        for u=1:1:2;%sweeping through window
            for v=1:1:2;
                %taking arithmetic mean
                im(i,j)=im(i,j)+im(i+u,j+v);
                im(i,j)=im(i,j)/9;
            end
        end       
    end
end
%showing final image
figure,imshow(im);
title('arthmean');

[KALYAN ACHARJYA]

[Lama Awawdeh]

thank you so much for your help. i appriciate your time that spent in helping me.

[KALYAN ACHARJYA]

Welcome !

[Lama Awawdeh]

my pleasure, done.

Answer 2

[Jan]

for u=1:1:2;%sweeping through window
    for v=1:1:2;
        %taking arithmetic mean
        im(i,j)=im(i,j)+im(i+u,j+v);
        im(i,j)=im(i,j)/9;
    end
end

This is a 2x2 window starting right and top of the current pixel. You divide the result repeatedly by 9. I assume you mean:

for i = 2:row-1
    for j = 2:col - 1
        a = 0;
        for u = -1:1
            for v = -1:1
                a = a + im(i + u, j + v);
            end
        end
        im(i, j) = a / 9;
    end
end

conv2 and filter2 will me more efficient.