Three-dimensional Point Spread Function Engineering for High Axial Localization Precision