In the past decade, the use of three-dimensional forest information from airborne Light Detection and Ranging (LiDAR) has become widespread in forest inventories. Accurate Individual Treetop Detection (ITD) and crown boundary delineation using LiDAR data are critical for obtaining precise inventory metrics. To address this need, we introduced a novel growing tree region (GTR)-driven ITD method that utilizes canopy height models (CHM) derived from very low-resolution airborne LiDAR data. The GTR algorithm consists of three key stages: (i) preserving all height layers through incremental cutting and stacking of CHM; (ii) employing a three-layer concept to identify individual treetops; and (iii) refining the detected treetops using a distance-based filter. Our method was tested in five temperate forests across Central Europe and was compared against the widely-used local maxima (LM) search combined with an optimized variable window filtering (VWF) technique. Our results showed that the GTR method outperformed LM with VWF, particularly in forests with high canopy density. The achieved root mean square accuracies were 74% for the matching rate, 19% for commission errors, and 27% for omission errors. In comparison, the LM with the VWF method resulted in a matching rate of 71%, commission errors of 20%, and omission errors of 31%. To facilitate the application of our algorithm, we developed an R package called TREETOPS, which seamlessly integrates with the lidR package, ensuring compatibility with existing treetop-based segmentation methods. By introducing TREETOPS, we provide the most accurate open-source tool for detecting treetops using low-resolution LiDAR-derived CHM.