Background: Using the chiral (Kyushu) g-matrix folding model with the densities calculated with Gongny-HFB (GHFB) with the angular momentum projection (AMP), we determined the neutron skin r 208 skin = 0.25 fm from the central values of reaction cross sections σR(EXP) of p+ 208 Pb scattering in Ein = 40 − 81 MeV. As for 48 Ca, the high-resolution E1 polarizability experiment (E1pE) yields r 48 skin (E1pE) = 0.14 − 0.20 fm. Meanwhile, σR(EXP) are available as a function of incident energy Ein for p+ 48 Ca scattering. Aim: Our aim is to determine r 48 skin from the central values of σR(EXP) for p+ 48 Ca scattering by using the Kyushu g-matrix folding model. Results: As for 48 Ca, we first determine neutron radius rn(E1pE) = 3.56 fm from the central value 0.17 fm of r 48 skin (E1pE) and the proton radius rp(EXP) = 3.385 fm of electron scattering, and evaluate matter radius rm(E1pE) = 3.485 fm from the rn(E1pE) and the rp(EXP) of electron scattering. The Kyushu g-matrix folding model with the GHFB+AMP densities reproduces σR(EXP) in 23 ≤ Ein ≤ 25.3 MeV in one-σ level, but slightly overestimates the central values of σR(EXP) there. In 23 ≤ Ein ≤ 25.3 MeV, the small deviation allows us to scale the GHFB+AMP densities to the central values of rp(EXP) and rn(E1pE). The reaction cross sections σR(E1pE) obtained with the scaled densities almost reproduce the central values of σR(EXP) when Ein = 23−25.3 MeV, so that the σR(AMP) obtained with the GHFB+AMP densities and the σR(E1pE) are in one σ of σR(EXP) there. In Ein = 23 − 25.3 MeV, we determine the rm(EXP) from the central values of σR(EXP) and take the average for the rm(EXP). The averaged value is rm(EXP) = 3.471 fm. Eventually, we obtain r 48 skin (EXP) = 0.146 fm from rm(EXP) = 3.471 fm and rp(EXP) = 3.385 fm.