This paper employs the reproducing kernel (RK) approximation for evaluation of field theory-based incompatibility tensor in a polycrystalline plasticity simulation. The modulation patterns, which is interpreted as mimicking geometrical-type dislocation substructures, are obtained based on the proposed method. Comparisons are made using FEM and RK based approximation methods among different support sizes and other evaluation conditions of the strain gradients. It is demonstrated that the evolution of the modulation patterns needs to be accurately calculated at each time step to yield a correct physical interpretation. The effect of the higher order strain derivative processing zone on the predicted modulation patterns is also discussed.