The numerical generating method for global homeomorphic grid from multi-dimensional towel-shaped reality region to square computational region is investigated. The numerical iterative method for multi-dimensional equal-parameter transformation is given. Its error analysis and the sufficient condition of numerical grid without overlapping is pressented. Finally, an example about two-dimensional non-overlapping numerical grid generation for Qiantang Estuary is shown.