Solved

Change loop to matrix operation in matlab (image processing)

Posted on 2008-10-29
1
1,767 Views
Last Modified: 2016-02-10
Hei!

I have a problem with image processing in matlab for a lab course. I have only basic matlab knowledge and hope to be able to explain my problem with my basic school english...
The code attached is a loop in matlab to calculate a "moving average" on a picture. The first and second loop is for the rows (german: Zeilen) and columns (Spalten) of the picture (P), the third and fourth loop is for the rows and columns of the filter matrix (kern). Then it sums up the values of the filter matrix and divides it by the sum of the loadings. ewr and ewc are needed to calculate the new border of the image. (zero padding or replicate)

My question:
I need a way to do this with fast matrix/vector operations if possible.
I am not allowed to use any imfilter, conv,  conv2, block-processing or sliding-window functions.
What i can imagine is to exchange the 2 inner loops by a matrix-multiplication. I would multiply the filter matrix with the right submatrix of the image. But to replace the outer loops i need to know the position/index to know wich submatrix i have to multiply...
I read something about the not so well known accumarray() function but i don't know how i can use this for 2D operation.
Maybe someone can give me a hint...

Thanks in advance!

Greetings from Vienna, Austria.

beta
for ip = (ewr+1):(rc(1)+(ewr)) %Zeilen

    for np = (ewc+1):(rc(2)+ewc)%Spalten

        for Z = -ewr:ewr

            for S = -ewc:ewc

            PF(ip,np)=PF(ip,np)+P(ip+Z,np+S)*kern(Z+ewr+1,S+ewc+1);

            end

        end

        PF(ip,np)=PF(ip,np)/(sum(sum(kern,1),2));  

    end

end

Open in new window

0
Comment
Question by:strupet
1 Comment
 

Accepted Solution

by:
strupet earned 0 total points
ID: 22838649
i did it with a matrix multiplication (see attached code) but it seems much slower than the previous version with four loops.
for ip = (row_ext+1):(rc(1)+row_ext) %rows

    for np = (col_ext+1):(rc(2)+col_ext)%columns

        PF(ip,np)=sum(sum(kern.*P( (ip-row_ext):(ip+row_ext) , (np-col_ext):(np+col_ext)),1),2)/sum(sum(kern,1),2);

    end

end

Open in new window

0

Featured Post

Is Your Active Directory as Secure as You Think?

More than 75% of all records are compromised because of the loss or theft of a privileged credential. Experts have been exploring Active Directory infrastructure to identify key threats and establish best practices for keeping data safe. Attend this month’s webinar to learn more.

Question has a verified solution.

If you are experiencing a similar issue, please ask a related question

Suggested Solutions

Title # Comments Views Activity
mapAB Challlenge 35 123
oracle query help 18 99
python sqlite question 11 46
Help to convert powershell script into a gui 9 133
When we want to run, execute or repeat a statement multiple times, a loop is necessary. This article covers the two types of loops in Python: the while loop and the for loop.
Whether you’re a college noob or a soon-to-be pro, these tips are sure to help you in your journey to becoming a programming ninja and stand out from the crowd.
This theoretical tutorial explains exceptions, reasons for exceptions, different categories of exception and exception hierarchy.
The goal of the tutorial is to teach the user how to use functions in C++. The video will cover how to define functions, how to call functions and how to create functions prototypes. Microsoft Visual C++ 2010 Express will be used as a text editor an…

947 members asked questions and received personalized solutions in the past 7 days.

Join the community of 500,000 technology professionals and ask your questions.

Join & Ask a Question

Need Help in Real-Time?

Connect with top rated Experts

21 Experts available now in Live!

Get 1:1 Help Now