The ICESat-2 Land and Vegetation Height (ATL08) product is a new control point dataset for large-scale topographic mapping and geodetic surveying. However, its elevation accuracy is typically affected by multiple factors. The study aims to propose a new approach to classify ATL08 terrain estimates into different accuracy levels and extract reliable ground control points (GCPs) from ICESat-2 ATL08. Specifically, the methodology is divided into three stages. First, the ATL08 terrain estimates are matched with the raw ATL03 photon cloud data, and the ATL08 terrain estimates are used to fit a continuous terrain curve. Then, using the fitted continuous terrain curve and raw ATL03 photon cloud data, a pseudo-waveform is generated for grading the ATL08 terrain estimates. Finally, all the ATL08 terrain estimates are graded based on the peak characteristics of the generated pseudo-waveform. To validate the feasibility of the proposed method, four study areas from the National Ecological Observatory Network (NEON), characterized by various terrain features and forest types were selected. High-accuracy airborne lidar data were used to evaluate the accuracy of graded ICESat-2 terrain estimates. The results demonstrate that the method effectively classified all ATL08 terrain estimates into different accuracy levels and successfully extracted high-accuracy GCPs. The root mean square errors (RMSEs) of the first accuracy level in the four selected study areas were 0.99 m, 0.51 m, 1.88 m, and 0.65 m, representing accuracy improvement of 51.7%, 58.2%, 83.1%, and 68.8%, respectively, compared to the original ATL08 terrain estimates before classifying. Additionally, a comparison with the conventional threshold-based GCP extraction method demonstrated the superior performance of our proposed approach. This study introduces a new approach to extract high-quality elevation control points from ICESat-2 ATL08 data, particularly in forested areas.