A novel way of determining the Hamiltonian of the interacting boson model (IBM) is proposed. Based on the fact that the potential energy surface of the mean-field model, e.g., the Skyrme model, can be simulated by that of the IBM, parameters of the IBM Hamiltonian are obtained. By this method, the multifermion dynamics of surface deformation can be mapped, in a good approximation, onto a boson system. The validity of this process is examined for Sm and Ba isotopes, and an application is presented to an unexplored territory of the nuclear chart, namely, the right lower corner of 208Pb.