The water index (WI) is designed to highlight inland water bodies in remotely sensed imagery. The application of WI for water body mapping is mainly based on the thresholding method. However, there are three primary difficulties with this method: (1) inefficient identification of mixed water pixels; (2) confusion of water bodies with background noise; and (3) variation in the threshold values according to the location and time of image acquisitions. Considering that mixed water pixels usually appear in narrow rivers or shallow water at the edge of lakes or wide rivers, an automated method is proposed for extracting rivers and lakes by combining the WI with digital image processing techniques to address the above issues. The data sources are the Landsat (TM) (Thematic Mapper) and ETM+ (Enhanced Thematic Mapper Plus) images for three representative areas in China. The results were compared with those from existing thresholding methods. The robustness of the new method in combination with different WIs is also assessed. Several metrics, which include the Kappa coefficient, omission and commission errors, edge position accuracy and completeness, were calculated to assess the method's performance. The new method generally outperformed the thresholding methods, although the degree of improvement varied among WIs. The advantages and limitations of the proposed method are also discussed.