The spin-1 Haldane chain is an example of the symmetry-protected-topological (SPT) phase in one dimension. Experimental realization of the spin chain materials usually involves both the uniaxialtype, D(S z ) 2 , and the rhombic-type,, single-ion anisotropies. Here, we provide a precise ground-state phase diagram for spin-1 Haldane chain with these single-ion anisotropies. Using quantum numbers, we find that the Z2 symmetry breaking phase can be characterized by double degeneracy in the entanglement spectrum. Topological quantum phase transitions take place on particular paths in the phase diagram, from the Haldane phase to the Large-Ex, Large-Ey, or Large-D phases. The topological critical points are determined by the level spectroscopy method with a newly developed parity technique in the density matrix renormalization group [Phys. Rev. B 86, 024403 (2012)], and the Haldane-Large-D critical point is obtained with an unprecedented precision, (D/J)c=0.9684713(1). Close to this critical point, a small rhombic single-ion anisotropy |E|/J ≪ 1 can destroy the Haldane phase and bring the system into a y-Néel phase. We propose that the compound [Ni(HF2)(3-Clpy)4]BF4 is a candidate system to search for the y-Néel phase.Introduction. Quantum magnetism of integer-spin chains has been attracting attention for decades. It was stimulated by the Haldane conjecture [1] that the lowest excitation in the antiferromagnetic Heisenberg model are gapped if and only if the spin S is an integer. Experimental evidences for the Haldane gap were discovered in several S=1 quasi-one dimensional (Q1D) materials, e.g., [7, 8], and [Ni(C 2 H 8 N 2 ) 2 NO 2 ]BF 4 (NENB) [9]. Due to the crystal field and the spin-orbit coupling, the microscopic effective Hamiltonian for the Q1D spin chains involves the single-ion anisotropies,