How to specify boundary condition as DIRICHLET in pcl poisson surface reconstruction?
I'm doing a job of surface reconstruction. But I met an issue :
I want to use the DIRICHLET boundary condition in poisson, but it seems that the poisson of pcl doesn't support specifying boundary condition, it just uses NEUMANN boundary condition always.
So I wonder how to use the DIRICHLET boundary condition in pcl poisson.
BTW：My goal is to calculate the volume of a container. But my pointcloud isn't watertight, so I need the algorithm to 'image' the surface of holes. CloudCompare supports specifying boundary condition, and it works well. But in pcl, the effect of NEUMANN boundary condition is terrible.
The mesh generated by pcl poisson(NEUMANN condition) like below:
The mesh generated by cloudcompare(specified DIRICHLET condition) like below:
The original PoissonRecon code is at this github repository. You can also find there prebuild executables for window command line (--bType to set the conidtion). This is available in the command line executable starting from version 9.0.
[--bType ] This integer specifies the boundary type for the finite elements. Valid values are: 1: Free boundary constraints 2: Dirichlet boundary constraints 3: Neumann boundary constraints The default value for this parameter is 3 (Neumann).
CloudCompare uses version 7.
PCL (1.12.0 at the moment of this post) uses version 4 of PoissonRecon.
Open3D (0.14.1 at the moment of this post) includes a wrapper over version 12, which supports both boundary conditions. It is however hard-coded to use NEUMANN. You should be able to easily change the enum and compile a version of Open3D that uses DIRICHLET condition (never tried this myself).
Alternatively (if you can't use the original console app or recompile Open3D), you can try to work with what you've got.
- You can try to identify the "imaginary faces" based on the area of the faces (smaller density - larger triangle area) and remove them. The original repository offers a SurfaceTrimmer tool (another console project) that does just that (based on density value).
- Close the remaining open mesh using either some hole-closing-method or convex-hull.
- why cmake not compiling for pcl
- Consistent normal calculation of a point cloud
- How to compile Point Cloud Library (PCL) with examples (Ubuntu)?
- send pcl point cloud over tcp in c++
- Why is SSE alignement necessary while doing SIMD instructions?
- Error:Could NOT find PkgConfig (missing: PKG_CONFIG_EXECUTABLE)
- Divide a point cloud into many small surfaces in PCL
- ImportError: DLL load failed with pybind11 and PCL
- After cropping a Point Cloud using pcl::CropBox filter, can't visualize that cropped Point Cloud
- Return quality of fit for PointCloudLibrary RANSAC model fitting
- how to import all the points from a .pcd file into a 2d python array
- Align two Point Clouds with PCL only by translation
- How to handle with 'Policy CMP0074 is not set' warning while using Cmake to build PCL tutorials?
- Given one set of 2D points as groud truth, how to match it with one observed set
- PCLVisualizer flashing when using .spinOnce()
- Occlusion estimation in pointcloud using pcl voxelgridOcclusionEstimation
- Which pcl filter to use to downsample a point cloud
- PCL not finding QHull includes
- How to specify boundary condition as DIRICHLET in pcl poisson surface reconstruction?
- pcl::io::loadPCDFile Could not find file 'gArBaGe'
- PCL Feature matching
- Results of RANSAC plane segmentation depends on the size of point cloud?
- Python - Display 3D Point Cloud
- PCL Own Point Type can't use pcl functions and operators
- Metric to compare two point clouds similarity
- Convex Hull calculation in Point Cloud Library fails in 2 as well as 3 dimensions
- A strange phenomenon when using pcl::ApproximateGrid:
- Use Cmake with PCL and OpenCV
- Can anyone explain the difference between organized and unorganized point cloud?
- Generate image from an unorganized typical scan-pattern lidar's point cloud data