def funcGauss2D(varSizeX, varSizeY, varPosX, varPosY, varSd):
varSizeX = int(varSizeX)<EOL>varSizeY = int(varSizeY)<EOL>aryY, aryX = sp.mgrid[<NUM_LIT:0>:varSizeX,<EOL><NUM_LIT:0>:varSizeY]<EOL>aryGauss = (<EOL>(<EOL>np.power((aryX - varPosX), <NUM_LIT>) +<EOL>np.power((aryY - varPosY), <NUM_LIT>)<EOL>) /<EOL>(<NUM_LIT> * np.power(varSd, <NUM_LIT>))<EOL>)<EOL>aryGauss = np.exp(-aryGauss)<EOL>return aryGauss<EOL>
Create 2D Gaussian kernel. Source:
    http://mathworld.wolfram.com/GaussianFunction.html