Improved position estimation for a square photodetector having terminals at its corners is provided by first estimating an event position in a coordinate system aligned with the diagonals of the square, and then performing a coordinate rotation to provide an estimated position in a coordinate system aligned with the sides of the square. These two steps can also be mathematically combined. It is believed that the improved accuracy of this position determination approach is a result of a greater degree of linearity in charge division between terminals along diagonal axes than along orthogonal axes of a square detector. Flood images provided by this method have been compared to flood images provided by the conventional Anger method, and substantially less image distortion is observed with the present method than with the Anger method. This reduced distortion can enable automatic crystal segmentation from raw flood image data, which is particularly beneficial in connection with simplifying positron emission tomography (PET) system calibration.