I want to write code in matlab to find a characteristic scale of an image patch. I have an image patch X, I need to maximize function F(X,sigma) so that sigma is the characteristic sigma.

For my experiments, F can be the laplacian of the image, sigma is the "power" of a gaussian filter on the image. F(x,sigma) should return for each x (image patch) its response over scales of increasing sigmas (i.e. a discrete value for each X and sigma). My questions are:

1. The laplacian is obtained by convolving laplacian kernel with the image, where does the sigma takes place ? (is it by smoothing the image before (the laplace) with G(sigma)?)
2. How can I define F(X,sigma) to return discrete value (X is nXm), ? i.e. what does this function do exactly.