To perform computer-aided diagnosis of the thyroid nodules on ultrasound images, the location and boundary of nodules should be clearly defined. However, the identification of thyroid nodule boundary is a difficult issue due to the biological characteristics of the nodules, the physics and quality of ultrasound imaging, and the subjective factors and operating conditions of the operator. In this study, we propose a novel and semiautomatic method for detecting the boundary of thyroid nodule based on the Variance-Reduction (V-R) statistics without image preprocessing. The region of interest (ROI) is first automatically generated according to the initial inputs of the nodule's major and minor axes. The boundary candidate pixel points are then extracted by using the V-R statistics from the grayscale values of all pixel points in the ROI. Three filtering methods are further applied to eliminate the outlier pixel points to ensure that the remaining candidate pixel points are located on the nodule boundary. Finally, the remaining pixel points are smoothened and linked together to form the final boundary. The proposed method is validated with ultrasound images of 538 thyroid nodules, with manual delineation by experienced radiologist as gold standard. The effectiveness is evaluated and compared with previous publications using boundary error metrics and overlapping area metrics with the same data set. The results show that the normalized average mean boundary error is 1.02%, the true positive overlapping area ratio achieves 93.66% and false positive overlapping area ratio is limited to 7.68%. In conclusion, our proposed method is reliable and effective in detecting thyroid nodule boundary on ultrasound images.