Regional gravity field modeling based on rectangular harmonic analysis