An area-type shear difference method, which has been previously developed by one of the authors, can analyze stresses in a member of structures using photoelastic data (principal stress differences and principal stress directions) obtained from photoelastic experiments. The method uses an approximate equation for a directly integrated stress equation of equilibrium between two points in any direction. In this study, most of the procedures in the stress analysis using the method were automated, and the stress analysis method was applied to photoelastic data obtained from an experiment on a circular disk subjected to diametric compression. The results indicated that the stress components obtained by the method were in good agreement with those obtained theoretically.