Solved

Change loop to matrix operation in matlab (image processing)

Posted on 2008-10-29
1
1,775 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
[X]
Welcome to Experts Exchange

Add your voice to the tech community where 5M+ people just like you are talking about what matters.

  • Help others & share knowledge
  • Earn cash & points
  • Learn & ask questions
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

Independent Software Vendors: We Want Your Opinion

We value your feedback.

Take our survey and automatically be enter to win anyone of the following:
Yeti Cooler, Amazon eGift Card, and Movie eGift Card!

Question has a verified solution.

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

Suggested Solutions

Title # Comments Views Activity
How to get all the API from website? 11 156
Not seen Link button 5 76
Is there any way to convert exponential value to number in sql server 5 72
Is online banking safe? 11 116
The purpose of this article is to demonstrate how we can use conditional statements using Python.
In this post we will learn how to connect and configure Android Device (Smartphone etc.) with Android Studio. After that we will run a simple Hello World Program.
The viewer will learn how to clear a vector as well as how to detect empty vectors in C++.

740 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