Homoclinic points of 2D and 4D maps via the parametrization method