“…The workspace boundary for the needle is a 1D curve in each xy-plane means the z h -plane for z ¼ h, h ∈ 0, h max ½ with a clinical value h ¼ 200 mm. The S is the plan of the needle insertion trajectory, θ k , φ k ðÞ are the rotation angles with respect to y and z axes θ min ≤ θ k ≤ θ max , ϕ min ≤ φ k ≤ φ max [45].…”