In this paper, we present an accurate and efficient isogeometric topology optimization method that integrates the non-uniform rational B-splines based isogeometric analysis and the parameterized level set method for minimal compliance problems. The same NURBS basis functions are used to parameterize the level set function and evaluate the objective function, and therefore the design variables are associated with the control points. The coefficient matrix that parameterizes the level set function is set up by a collocation method that uses the Greville abscissae. The zero-level set boundary is obtained from the interpolation points corresponding to the vertices of the knot spans. Numerical examples demonstrate the validity and efficiency of the proposed method.