Digital photoelasticity is used for the experimental evaluation of the strain intensity factor for a rigid line inclusion embedded in an elastic matrix. The multi-parameter stress field equations are derived using the Airy's stress function approach. The ten-step phase shifting technique (PST) is used to obtain the whole field isochromatic and isoclinic data, and the adaptive quality-guided phase unwrapping (AQGPU) algorithm is used for unwrapping the phase map. An over-deterministic linear least squares approach is used to solve the multi-parameter stress field equations in an optimization framework to obtain the inclusion tip parameters and the actual inclusion tip location. The relations for strain intensity factor and stress intensity factor (SIF) with the parameters is presented. Photoelastic experiments are performed on two inclusion configurations where inclusion is kept along and inclined to the direction of loading. The obtained strain intensity factor is compared with the analytical solution, and a good match is observed. © 2019 Elsevier Ltd