When solving acoustic radiation problem in the routine boundary point method (BPM), the inverse of the particular solution matrix of surface normal velocities is needed and the matrix must be full rank and reversible. However, the type of the particular solution sources is single and the sources are not linear independent, so the locations of the particular solution sources inside the vibrating body must be selected carefully. But, in the routine BPM, the locations are determined by an experiential formula, so the method may be invalid for a complicated vibrating body. In this paper, the construction method of the particular solution sources is improved first, and the singular value decomposition (SVD) technique and the Moore-Penrose pseudoinverse are adopted to realize the inversion. As a consequence, the particular solution matrix of surface normal velocities can be non-full rank and the locations of the particular sources can be determined easily. On the basis of the improved BPM, the spherical wave sources of different orders are proposed to be the particular solution sources. Here, all particular solution sources are located on only one point inside the vibrating body, so the problem of the locations of the particular solution sources is thoroughly solved, and such particular solution sources that are useless for the calculation results are discarded. The theoretical model is established at first, and then the proposed method is used to realize the nearfield acoustic holography (NAH). Subsequently, an experiment is investigated to validate the feasibility and correctness of the proposed method and its application to acoustic holography.