Construction of a high spatial resolution and high precision marine gravity field in coastal areas is constrained by the low quality and sparse coverage of altimetry data, except for limited shipborne and airborne gravity surveys. To address this problem, a mean sea surface height constraint factor (MSSHCF) method based on the ordinary kriging method and the remove-restore technique is proposed from the perspective of interpolation. In this method, the data is standardized during the interpolation process to reduce the error and mean sea surface as variables related to the marine gravity field are added to the semi-variance function in ordinary kriging to obtain a marine gravity field with a spatial resolution of 1′ × 1′. Validation experiments show that the MSSHCF method more closely agrees with the referenced SS V28, DTU17 global marine gravity models than the ordinary kriging method. Our results were further validated against shipborne data; the accuracy of the MSSHCF method is 0.13 and 0.33 mGal higher than that of the ordinary kriging method in two experimental areas. The effects of ocean depth and offshore distance on the results were also assessed. These results show that the proposed method is more accurate than the ordinary kriging method, when the distance and depth varied. Therefore, our study demonstrates that the MSSHCF method is an innovative and feasible tool for extracting gravity fields along coastal, beach, and island areas.