帮我编写一段代码，使用PCL点云库完成森林点云的单木分割

好的，以下是一段使用PCL点云库完成森林点云单木分割的代码示例：

#include <iostream>
#include <pcl/point_cloud.h>
#include <pcl/octree/octree.h>
#include <pcl/io/pcd_io.h>
#include <pcl/filters/extract_indices.h>

int main (int argc, char** argv)
{
    // 读入森林点云数据
    pcl::PointCloud<pcl::PointXYZ>::Ptr cloud (new pcl::PointCloud<pcl::PointXYZ>);
    if (pcl::io::loadPCDFile<pcl::PointXYZ> ("forest.pcd", *cloud) == -1) 
    {
        std::cout << "Error loading point cloud" << std::endl;
        return (-1);
    }

    // 创建octree并设置八叉树的分辨率为1.0m
    float resolution = 1.0;
    pcl::octree::OctreePointCloudSearch<pcl::PointXYZ> octree (resolution);
    octree.setInputCloud (cloud);
    octree.addPointsFromInputCloud ();

    // 循环每一个叶子节点
    std::vector<int> leafNodeIndices;
    for (size_t i = 0; i < octree.getLeafNodeCount (); ++i)
    {
        leafNodeIndices.clear ();
        octree.getLeafNodeIndices (i, leafNodeIndices);

        // 创建一个指向该叶子节点对应点云的指针
        pcl::PointCloud<pcl::PointXYZ>::Ptr leafNodeCloud (new pcl::PointCloud<pcl::PointXYZ>);
        leafNodeCloud->width = leafNodeIndices.size ();
        leafNodeCloud->height = 1;
        leafNodeCloud->points.resize (leafNodeCloud->width * leafNodeCloud->height);

        for (size_t j = 0; j < leafNodeIndices.size (); ++j)
        {
            leafNodeCloud->points[j].x = cloud->points[leafNodeIndices[j]].x;
            leafNodeCloud->points[j].y = cloud->points[leafNodeIndices[j]].y;
            leafNodeCloud->points[j].z = cloud->points[leafNodeIndices[j]].z;
        }

        // 存储该叶子节点对应的点云
        pcl::io::savePCDFileASCII (std::to_string (i) + ".pcd", *leafNodeCloud);