# 2d gaussian function

조회 수: 276 (최근 30일)
lital 2011년 8월 3일
댓글: Walter Roberson 2020년 10월 17일
I need to plot a 2d gaussian function, where x and y corresponds to the image pixels, my code uses a nested for loop which makes my program run extremely slow, is there a way to write this in a more faster way?
(right now it takes about 8-10 sec to run on 1920*1080 size matrix and i need to produce 173,340 2d gaussian functions images which is too much time...)
thanks in advance for any help.
this is my current code:
function mat = gauss2d(mat, sigma, center)
gsize = size(mat);
for r=1:gsize(1)
for c=1:gsize(2)
mat(r,c) = gaussC(r,c, sigma, center);
end
end
function val = gaussC(x, y, sigma, center)
xc = center(1);
yc = center(2);
exponent = ((x-xc).^2 + (y-yc).^2)./(2*sigma);
val = (exp(-exponent));
##### 댓글 수: 10이전 댓글 8개 표시이전 댓글 8개 숨기기
Andicha Zain 2020년 10월 17일
not project but it's my topic thesis.
Walter Roberson 2020년 10월 17일
It is already Saturday in Atlantic Time and everywhere east of that. There is zero chance that anyone is going to write the code for you in the next few hours before you meet with your professor.
However, if you ask specific questions about MATLAB then there is a possibility that someone might answer in the next few hours.

댓글을 달려면 로그인하십시오.

### 채택된 답변

Walter Roberson 2011년 8월 3일
function mat = gauss2d(mat, sigma, center)
gsize = size(mat);
[R,C] = ndgrid(1:gsize(1), 1:gsize(2));
mat = gaussC(R,C, sigma, center);
##### 댓글 수: 5이전 댓글 3개 표시이전 댓글 3개 숨기기
BBB 2019년 11월 11일
Hi Walter,
It's perfect and working! Thanks :)
Andicha Zain 2020년 9월 21일
편집: Walter Roberson 2020년 10월 15일
hello can anyone help me
how to write program The Influence of Coherence Length and Gaussian Beam on the Multibeam Interference by Matlab ?

댓글을 달려면 로그인하십시오.

### 추가 답변 (5개)

Image Analyst 2011년 11월 26일
Is that 8-10 seconds for all 173,340 Gaussians? That's not bad. I have a demo that randomly places Gaussians in a larger image using fspecial() and indexing, not one pixel at a time like you did. It takes about 7.4 seconds on an old computer for 173,340 randomly placed Gaussians in a 1920x1080 image. Here's the code, in case you're interested:
% Demo to randomly place Gaussians in an image.
% By ImageAnalyst
clc; % Clear the command window.
close all; % Close all figures (except those of imtool.)
imtool close all; % Close all imtool figures.
clear; % Erase all existing variables.
workspace; % Make sure the workspace panel is showing.
tic;
% Set up some parameters.
fontSize = 20;
backgroundGrayLevel = 128;
windowSize = 50; % Could be random if you want.
sigma = 10; % Could be random if you want.
numberOfGaussians = 173340;
rows = 1080;
columns = 1920;
% Create one Gaussian.
g = fspecial('gaussian', windowSize, sigma);
grayImage = backgroundGrayLevel * ones(rows, columns);
% Create random signs so that the Gaussians are
% randomly brighter or darker than the background.
s = 2*randi(2, [1 numberOfGaussians])-3;
% Note: g and grayImage are floating point images, not uint8,
% though you could modify the program to have them be uint8 if you wanted.
% Get a list of random locations.
randomRow = randi(rows-windowSize+1, [1 numberOfGaussians]);
randomCol = randi(columns-windowSize+1, [1 numberOfGaussians]);
% Place the Gaussians on the image at those random locations.
for k = 1 : numberOfGaussians
grayImage(randomRow(k):randomRow(k)+windowSize-1, randomCol(k):randomCol(k)+windowSize-1) = ...
grayImage(randomRow(k):randomRow(k)+windowSize-1, randomCol(k):randomCol(k)+windowSize-1) + ...
s(k) * g;
end
toc;
% Display the final image.
imshow(grayImage, []);
caption = sprintf('%d Gaussians, Randomly Placed', numberOfGaussians);
title(caption, 'FontSize', fontSize);
axis on;
% Enlarge figure to full screen.
set(gcf, 'units','normalized','outerposition',[0 0 1 1])
set(gcf,'name','Demo by ImageAnalyst','numbertitle','off')
##### 댓글 수: 2없음 표시없음 숨기기
Muhammed Sameed 2020년 4월 24일
This is a really good piece of code. But how would you do the reverse i.e. from a noisy image with N randomly placed Gaussians, extract the centroid location and sigma of each Gaussian? I am interested in N=2 for the moment, but a generalized function would be helpful.
Walter Roberson 2020년 4월 24일
Hmmm, would that be like finding the values of x(n), y(n), width(n), height(n) such that the sum of the gaussians generated by those parameters is everywhere less than (image + 1/2*EPS(image)), where EPS(image) is 1 for integer-valued images and eps(image) for floating point? Which is to say, that the sum of the gaussians can exceed the value of the image at any given location, but must not exceed it by enough that the total would become the next representable number ? For example on an integer image, if an image location was 42, and the sum of gaussians predicted 42.42 there, then that would be okay because uint8(42.42) would be 42, but predicting 42.52 would not be okay because uint8(42.52) would be 43.

댓글을 달려면 로그인하십시오.

Mahdi 2014년 7월 22일
I faced the same problem, just so others know you can use fspecial('gaussian', hsize, sigma) intrinsic function.
##### 댓글 수: 1이전 댓글 -1개 표시이전 댓글 -1개 숨기기
Image Analyst 2014년 7월 22일
Thanks for calling it out specially - it was kind of buried in my demo code and hard to notice.

댓글을 달려면 로그인하십시오.

Logan 2011년 11월 25일
where do you get the gaussC function?
##### 댓글 수: 1이전 댓글 -1개 표시이전 댓글 -1개 숨기기
Walter Roberson 2011년 11월 26일
The code for gaussC is shown in the original Question by lital.

댓글을 달려면 로그인하십시오.

Karbala'a Unvi. Science 2014년 12월 28일
Dear Sir, I am interested about the code that you wrote about the 2D Gaussian. I have a problem that I want to an image data to be distributed in another image ( image A is the Original, image B is the data one) so that when you see image A you find that there is a noise in it ( where that noise is image B)... I hope that is a good information to help me in building the code... Help will be appreciated.. Thank you in advance
##### 댓글 수: 1이전 댓글 -1개 표시이전 댓글 -1개 숨기기
Image Analyst 2014년 12월 28일
This is not an answer. Who are you talking to? You should have put it as a comment under their answer. Why don't you just try your best and then post your code as a new question? And explain it better there. I don't even understand what you want. I don't know if B is a "data" image or a "noise" image and if you can just add B to A.

댓글을 달려면 로그인하십시오.

rusgu fcf 2016년 5월 6일
what does the variable 'center' signify in the code of gaussian2d by lital ?
##### 댓글 수: 3이전 댓글 1개 표시이전 댓글 1개 숨기기
Rym Benchaabane 2020년 3월 29일
is center a real number ie 0? or a two element array ie [0,0]?
Thank you
Image Analyst 2020년 3월 29일
It's a 2-element vector, as you can see by looking inside his gaussC() function.
I recommend you use my code or Walter's code since lital's had problems.

댓글을 달려면 로그인하십시오.

### 카테고리

Help CenterFile Exchange에서 Images에 대해 자세히 알아보기

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by