In this paper, we obtained the mean mass density parameter (MMDP) ρm(Z,N) based on the nuclear mass AME2020 and nuclear charge radius CR2013 databases. On this basis, we study the dependence relationship between neutron number N and the two neighboring isotopes MMDP difference δρ(Z,N), and then obtained the empirical formula for δρ(Z,N). By combining this empirical formula with the AME2020 and CR2013 databases, the root-mean-square deviation (RMSD) for 606 nuclear charge radii with Z,N≥ 21 between the calculated and experimental values is 0.0046fm. Considering the odd–even staggering and isospin term correction, the study found that the RMSD is reduced to 0.0040fm. This work also uses Levenberg–Marquart (L-M) neural network approach to study the nuclear charge radius. The RMSD of nuclear charge radius is only 0.0030 fm, where the accuracy is increased by about 25%. Furthermore, the predicted values of nuclear charge radius obtained using this modified empirical formula and Neural Network model are competitive with the results of other models. These results indicate that nuclear charge radius relation and Neural Network model based on δρ(Z,N) proposed in this work have a certain accuracy, simplicity and reliability.