Let \(A=\begin{bmatrix} A_{ij} \end{bmatrix}\) be an \(n\times n\) operator matrix, where each \(A_{ij}\) is a bounded linear operator on a complex Hilbert space. Among other numerical radius bounds, we show that \(w(A)\le w({\hat{A}}),\) where \({\hat{A}}=\begin{bmatrix} {\hat{a}}_{ij} \end{bmatrix}\) is an \(n\times n\) complex matrix, with \(\begin{aligned} {\hat{a}}_{ij}= {\left\{ \begin{array}{ll} w(A_{ii}) &{}\text {when }i=j, \\ \left\| | A_{ij}|+ | A_{ji}^*| \right\| ^{1/2} \left\| | A_{ji}|+ | A_{ij}^*| \right\| ^{1/2} &{}\text {when }i<j, \\ 0 &{}\hbox {when}\ i>j . \end{array}\right. } \end{aligned}\) This is a considerable improvement of the existing bound \(w(A)\le w({\tilde{A}}),\) where \({\tilde{A}}=\begin{bmatrix} {\tilde{a}}_{ij} \end{bmatrix}\) is an \(n\times n\) complex matrix, with \(\begin{aligned} {\tilde{a}}_{ij}= {\left\{ \begin{array}{ll} w(A_{ii}) &{}\hbox {when}\ i=j, \\ \Vert A_{ij}\Vert &{}\hbox {when}\ i\ne j. \end{array}\right. } \end{aligned}\) Further, applying the bounds, we develop the numerical radius bounds for the product of two operators and the commutator of operators. Also, we develop an upper bound for the spectral radius of the sum of the product of n pairs of operators, which improves the existing bound.