Customize Harmonic Potential Fields via Hybrid Optimization over Homotopic Paths