Generally yes. In your case (relatively small number of points), I would use greenspline since you can use your centroid (x,y,z) values directly as constraints (i.e. no blockm* need) and pick whatever output -R -I you want. You can even give it a mask grid to limit extrapolations outside your bounding polygon, or deal with that later via grdmask and grdmath.