Here is my first attempt at a 2 - D DOG.

Do[Plot3D[(1/((1*s)*(Sqrt[2*Pi]))) 
            Exp[-((((x^2) + 
                        (y^2)))/((2*(((1*s))^2))))] - (1/((1/
                      s)*(Sqrt[2*Pi]))) 
            Exp[-((((x^2) + 
                        (y^2)))/((2*(((1/s))^2))))], {x, -6, 
          6}, {y, -6, 6}, PlotPoints -> 75, 
        PlotRange -> {- .5,  .6}, 
        ViewPoint -> {1.539,  -2.840,  1.008}], {s, 1/2, 2, 
        1.5}]

\(ret = 
      Import["c:/zbs/html/content/vision/model/mathematic/Jesusface.bmp"]\)

Show[ret]

$ImportFormats

$ExportFormats

Shallow[InputForm[ret]]

retarray = ret[[1, 1] ];

Dimensions[retarray]

Export["c:/zbs/html/content/vision/model/mathematic/ganglion.txt", 
      retarray, "Table"]

MatrixQ[retarray]

(1/((1/2)*(Sqrt[2*Pi]))) 
        Exp[-((((((x - 7))^2) + (((y - 7))^2)))/((2*(((1/2))^2))))] - (1/((2)*(Sqrt[2*Pi]))) 
        Exp[-((((((x - 7))^2) + (((y - 7))^2)))/((2*(((2))^2))))]

s=(1/((1/
                2)*(Sqrt[2*Pi])))

se is the standard deviation for the excitatory portion of the gaussian in the 2D DOG.
si is the standard deviation for the inhibitory protion of the gausian the the 2D DOG

se  =  3

si  =  6

width is the radius of the receptive field in x and y (or i and j) .   It is sensitive to the size of the inhibitory region .   So currently these cells are all on center cells.

width  =  4*si

ganglion = 
      Table[(1/((se)*(Sqrt[2*Pi]))) 
              Exp[-((((((i - width))^2) + (((j - width))^2)))/((2*
                          (se^2))))] - (1/((si)*(Sqrt[2*Pi]))) 
              Exp[-((((((i - width))^2) + (((j - width))^2)))/((2*(((si))^2))))] // N, {i, 2*width + 1}, {j, 
          2*width + 1}]

Dimensions[ganglion]

Export["c:/zbs/html/content/vision/model/mathematic/cell.txt", 
      ganglion, "Table"]

Show[SurfaceGraphics[ganglion], 
      PlotRange -> {- .03,  .07}]

<< LinearAlgebra`MatrixManipulation`

width

subretarray = SubMatrix[retarray, {1, 1}, {2*width + 1, 2*width + 1}]; 
    MatrixForm[subretarray]

Dimensions[subretarray]

output = ganglion . subretarray

Dimensions[output]