2d gaussian function

395 ビュー (過去 30 日間)
lital
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 件のコメント
Andicha Zain
Andicha Zain 2020 年 10 月 17 日
not project but it's my topic thesis.
Walter Roberson
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
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 件のコメント
BBB
BBB 2019 年 11 月 11 日
Hi Walter,
It's perfect and working! Thanks :)
Andicha Zain
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 ?
please help me.

サインインしてコメントする。

その他の回答 (5 件)

Image Analyst
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
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
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
Mahdi 2014 年 7 月 22 日
I faced the same problem, just so others know you can use fspecial('gaussian', hsize, sigma) intrinsic function.
  1 件のコメント
Image Analyst
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
Logan 2011 年 11 月 25 日
where do you get the gaussC function?
  1 件のコメント
Walter Roberson
Walter Roberson 2011 年 11 月 26 日
The code for gaussC is shown in the original Question by lital.

サインインしてコメントする。


Karbala'a Unvi. Science
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 件のコメント
Image Analyst
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
rusgu fcf 2016 年 5 月 6 日
what does the variable 'center' signify in the code of gaussian2d by lital ?
  3 件のコメント
Rym Benchaabane
Rym Benchaabane 2020 年 3 月 29 日
is center a real number ie 0? or a two element array ie [0,0]?
Thank you
Image Analyst
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 Center および File ExchangeImages についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by