Filtering Point Clouds with a Passthrough Filter

This tutorial demonstrates how to execute passthrough filtering.

Note

This tutorial is applicable for execution for both within inside and outside a Docker* image. It assumes that the pcl-oneapi-tutorials Deb package is installed, and the user has copied the tutorial directory from /opt/intel/pcl/oneapi/tutorials/ to a user-writable directory.

  1. Prepare the environment:

    cd <path-to-oneapi-tutorials>/passthrough
    
    Copy to clipboard
  2. oneapi_passthrough.cpp should be in the directory with following content:

     1#include <pcl/oneapi/filters/passthrough.h>
     2#include <pcl/io/pcd_io.h>
     3#include <pcl/point_types.h>
     4#include <pcl/point_cloud.h>
     5
     6
     7int main (int argc, char** argv)
     8{
     9
    10  std::cout << "Running on device: " << dpct::get_default_queue().get_device().get_info<sycl::info::device::name>() << "\n";
    11
    12  // Read Point Cloud
    13  pcl::PointCloud<pcl::PointXYZ>::Ptr cloud_( new pcl::PointCloud<pcl::PointXYZ>() );
    14  int result = pcl::io::loadPCDFile("using_kinfu_large_scale_output.pcd", *cloud_);
    15  if (result != 0)
    16  {
    17    pcl::console::print_info ("Load pcd file failed.\n");
    18    return result;
    19  }
    20
    21  // Prepare Point Cloud Memory (output)
    22  pcl::PointCloud<pcl::PointXYZ>::Ptr cloud_filtered_( new pcl::PointCloud<pcl::PointXYZ>() );
    23
    24  // GPU calculate
    25  pcl::oneapi::PassThrough<pcl::PointXYZ> ps;
    26  ps.setInputCloud(cloud_);
    27  ps.setFilterFieldName ("z");
    28  ps.setFilterLimits (0.0, 1.0);
    29  ps.filter(*cloud_filtered_);
    30
    31  // print log
    32  std::cout << "[oneapi passthrough] PointCloud before filtering: " << cloud_->size() << std::endl;
    33  std::cout << "[oneapi passthrough] PointCloud after filtering: " << cloud_filtered_->size() << std::endl;
    34}
    
    Copy to clipboard
  3. Source the Intel® oneAPI Base Toolkit environment:

    source /opt/intel/oneapi/setvars.sh
    
    Copy to clipboard
  4. (Optional) Setup proxy setting to download test data:

    export http_proxy="http://<http_proxy>:port"
    export https_proxy="http://<https_proxy>:port"
    
    Copy to clipboard
  5. Build the code:

    mkdir build && cd build
    cmake ../
    make -j
    
    Copy to clipboard
  6. Run the binary:

    ./oneapi_passthrough
    
    Copy to clipboard
  7. Expected results example:

    [oneapi passthrough] PointCloud before filtering: 993419
    [oneapi passthrough] PointCloud after filtering: 328598
    
    Copy to clipboard

Code Explanation

Load the example PCD into a PointCloud<PointXYZ>.

  // Read Point Cloud
  pcl::PointCloud<pcl::PointXYZ>::Ptr cloud_( new pcl::PointCloud<pcl::PointXYZ>() );
  int result = pcl::io::loadPCDFile("using_kinfu_large_scale_output.pcd", *cloud_);
  if (result != 0)
  {
    pcl::console::print_info ("Load pcd file failed.\n");
    return result;
  }
Copy to clipboard

Prepare output buffer for filtered result.

  // Prepare Point Cloud Memory (output)
  pcl::PointCloud<pcl::PointXYZ>::Ptr cloud_filtered_( new pcl::PointCloud<pcl::PointXYZ>() );
Copy to clipboard

Starts to compute the result.

  // GPU calculate
  pcl::oneapi::PassThrough<pcl::PointXYZ> ps;
  ps.setInputCloud(cloud_);
  ps.setFilterFieldName ("z");
  ps.setFilterLimits (0.0, 1.0);
  ps.filter(*cloud_filtered_);
Copy to clipboard

Result (output log):

  std::cout << "[oneapi passthrough] PointCloud before filtering: " << cloud_->size() << std::endl;
  std::cout << "[oneapi passthrough] PointCloud after filtering: " << cloud_filtered_->size() << std::endl;
Copy to clipboard