I have a some x (longitude mid point), y (latitude mid point), z (value) data of 0.5degree size grid cell for a region stored in Xarray variable. I have been trying to overlay it on the topographic grid (load_earth_relief and then grdgradient) by using pygmt.grdsample using same spacing (“0.5m/0.5m”) and cubic-spline interpolation for both. Topography grid is well plotted but the above overlying data is interpolated from the mid-point of the edge cells of my region. So this is leaving a strip kind of outer edge non-computed. Any solution for this?