function centerGPos = cenP(I )
%cenP finds the core of the PSF
%
[row, col] = find(I > );
minRow = min(row); maxRow = max(row);
minCol = min(col); maxCol = max(col);
pos(:, :, ) = repmat((minRow:maxRow)', [1, maxCol - minCol + 1]);
pos(:, :, ) = repmat((minCol:maxCol), [maxRow - minRow + , ]);
centerGPos = [sum(sum(pos(:, :, ).*double(I(minRow:maxRow, minCol:maxCol))))/sum(sum(double(I(minRow:maxRow, minCol:maxCol)))), ...
sum(sum(pos(:, :, ).*double(I(minRow:maxRow, minCol:maxCol))))/sum(sum(double(I(minRow:maxRow, minCol:maxCol))))];
centerGPos = int16(centerGPos);
end