Gray level residual (GLR), also known as correlation residual field or image difference, reflects the intensity differences between reference and deformed images. In finite element-based global digital image correlation (DIC), GLR field can be readily computed and has been successfully used for crack or damage detection. The computation of pixel-wise GLR field requires dense displacement results, which can be obtained using mesh interpolation in global DIC. To obtain the GLR field in subset-based local DIC, in which results are obtained at discrete calculation points, a simple, efficient and practical computation method for GLR field is presented, which accounts for the influence of illumination variations unavoidably occurred in real DIC tests. The proposed method can be easily implemented on both CPU and GPU to realize real-time visualization. Experiments showed that the GLR field can quantitatively and intuitively evaluate the pixelwise matching quality and provide clearer and finer locations for new components or textures (e.g., cracks) presented in deformed images. By fully leveraging the pixelwise matching quality information provided by the GLR field, the capability and robustness of local DIC is expected to be further enhanced.