By the end of December 2020, an earthquake with Mw about 6.4 hit Sisak-Moslavina County, Croatia. The town of Petrinja was the most affected region with major power outage and many buildings collapsed. The damage also affected neighbor countries such as Bosnia and Herzegovina and Slovenia. As a light of this devastating event, a deformation map due to this earthquake could be generated by using remote sensing imagery from Sentinel-1 SAR data. InSAR could be used as deformation map but still affected with noise factor that could problematize the exact deformation value for further research. Thus in this study, 17 SAR data from Sentinel-1 satellite is used in order to generate the multi-temporal interferometry utilize Stanford Method for Persistent Scatterers (StaMPS). Mean deformation map that has been compensated from error factors such as atmospheric, topographic, temporal, and baseline errors are generated. Okada model then applied to the mean deformation result to generate the modeled earthquake, resulting the deformation is mostly dominated by strike-slip with 3 meter deformation as right lateral strike-slip. The Okada sources are having 11.63 km in length, 2.45 km in width, and 5.46 km in depth with the dip angle are about 84.47° and strike angle are about 142.88° from the north direction. The results from this modeling can be used as learning material to understand the seismic activity in the latest 2020 Petrinja, Croatia Earthquake.