How to implement the 'selected' feature of a mesh healer in code.

[JihoKang24]

hello devs

I found your mesh inspector mesh healer

when I use this healing options In this case
after 2 holes is removed when this selected repair

how can I do this work in cpp code?

always appreciate for your helps

[Grantim]

Hello!

With this settings it only use:

MeshLib/source/MRMesh/MRMeshBuilder.h

Lines 91 to 97 in a895000

	/// the function finds groups of mesh vertices located closer to each other than \param closeDist, and unites such vertices in one;
	/// then the mesh is rebuilt from the remaining triangles
	/// \param optionalVertOldToNew is the mapping of vertices: before -> after
	/// \param uniteOnlyBd if true then only boundary vertices can be united, all internal vertices (even close ones) will remain
	/// \return the number of vertices united, 0 means no change in the mesh
	MRMESH_API int uniteCloseVertices( Mesh & mesh, float closeDist, bool uniteOnlyBd = true,
	VertMap * optionalVertOldToNew = nullptr );

and

MeshLib/source/MRMesh/MRFillHoleNicely.h

Lines 52 to 58 in a895000

	/// fills a hole in mesh specified by one of its edge,
	/// optionally subdivides new patch on smaller triangles,
	/// optionally make smooth connection with existing triangles outside the hole
	/// \return triangles of the patch
	MRMESH_API FaceBitSet fillHoleNicely( Mesh & mesh,
	EdgeId holeEdge, ///< left of this edge must not have a face and it will be filled
	const FillHoleNicelySettings & settings );

(for small holes)

But in your case there is no small holes so only uniteCloseVertices is called

[JihoKang24]

hello grantim
this is mesh files in this code result
case_result.zip

this is our code to remove collide face when add two closed mesh

    // remove previous duplicated vertex to work like load mesh
    bool VxMeshLibSolidFix(Mesh& mesh)
    {
	mesh.pack();
	auto pos = mesh.points;

	vector<stVertexKey<Vector3f>> vtrVertices = vector<stVertexKey<Vector3f>>(pos.size());
	BitSetParallelFor(mesh.topology.getValidVerts(), [&](VertId v)
	{
		vtrVertices[(int)v].v3Position = pos[v];
	});

	unordered_map<stVertexKey<Vector3f>, std::vector<VertId>, stHashVert<Vector3f>> umapHash;
	unordered_map<VertId, VertId> umapVert;
	for (VertId v : mesh.topology.getValidVerts())
	{
		auto& stVert = vtrVertices[v];
		umapHash[stVert].push_back(v);
	}

	unordered_map<int, bool> umapRemove;
	for (auto& kvp : umapHash)
	{
		std::sort(kvp.second.begin(), kvp.second.end());
		for (int i = 1; i < kvp.second.size(); i++)
		{
			umapRemove[kvp.second[i]] = true;
		}
	}

	vector<int> vtrVertOffset = vector<int>(mesh.topology.numValidVerts());
	int offset = 0;
	VertCoords verts;
	for (int i = 0; i < mesh.topology.numValidVerts(); i++)
	{
		VertId v = (VertId)i;
		vtrVertOffset[i] = offset;
		if (umapRemove.contains(i))
		{
			offset++;
			continue;
		}

		verts.push_back(mesh.points[v]);
	}

	auto remapVert = [&](int v)
	{
		if (umapRemove.contains(v))
		{
			int rep = umapHash[vtrVertices[v]][0];
			return VertId(rep - vtrVertOffset[rep]);
		}

		return VertId(v - vtrVertOffset[v]);
	};

	auto tris = mesh.topology.getAllTriVerts();
	Triangulation indexList;
	for (auto tri : tris)
	{
		ThreeVertIds tvId;
		for (int i = 0; i < 3; i++)
		{
			tvId[i] = remapVert((int)tri[i]);
		}

		bool bValid = true;
		for (int i = 0; i < 3 && bValid; i++)
		{
			for (int j = 0; j < 3 && bValid; j++)
			{
				if (i == j)
					continue;

				if (tvId[i] == tvId[j])
				{
					bValid = false;
				}
			}
		}

		if (bValid == true)
			indexList.push_back(tvId);
	}

	Mesh tempMesh = Mesh::fromTrianglesDuplicatingNonManifoldVertices(std::move(verts), indexList);
	mesh = tempMesh;
	return true;
    }

    bool VxMeshLibRemoveOveralppingArea(SVXModuleParameters* pstModuleParameter)
    {
	CVXVObjectPrimitive* pCMeshIn = (CVXVObjectPrimitive*)pstModuleParameter->GetVXObject(vxrObjectTypePRIMITIVE, true, 0);
	
	MR::Mesh mesh;
	VxConvertVxToMRMesh(pCMeshIn, mesh);
	FindOverlappingSettings findSettings;
	VxMeshLibSolidFix(mesh);
	MeshSave::toAnySupportedFormat(mesh, "C:\\Users\\asdf\\Desktop\\originMesh.stl");
	findSettings.maxDistSq *= mesh.computeBoundingBox().size().lengthSq();
	FaceBitSet fsOverlap = *MR::findOverlappingTris(mesh, findSettings);
	fsOverlap = MeshComponents::getLargeByAreaComponents(MeshPart(mesh, &fsOverlap), 0.1f, {});

	EdgeLoops holes;
	mesh.deleteFaces(fsOverlap);

	int cnt = 0;
	FillHoleNicelySettings fillHoleParams;
	fillHoleParams.triangulateOnly = true;
	do
	{
		holes = findRightBoundary(mesh.topology);
		cnt = 0;
		for (auto hole : holes)
		{
			for (auto e : hole)
			{
				double dSumPerimeter = mesh.holePerimiter(e);
				if (dSumPerimeter < 4.34)
				{
					fillHoleParams.triangulateParams.metric = getPlaneFillMetric(mesh, e);
					MR::fillHoleNicely(mesh, e, fillHoleParams);
					cnt++;
				}
			}
		}
	} while (cnt > 0);

	MeshBuilder::uniteCloseVertices(mesh, 0.0);
	MeshSave::toAnySupportedFormat(mesh, "C:\\Users\\asdf\\Desktop\\fixSmallHoles.stl");
	VxMeshLibSolidFix(mesh);
	MeshSave::toAnySupportedFormat(mesh, "C:\\Users\\asdf\\Desktop\\afterSolid.stl");

	holes = findRightBoundary(mesh.topology);
	auto getClosestEdge = [&](EdgeLoops eLoops) -> int
	{
		double dMin = 1e30;
		int iMin = 0;
		for (int i = 1; i < eLoops.size(); i++)
		{
			for (auto& e : eLoops[i])
			{
				double dDist = (mesh.orgPnt(eLoops[0][0]) - mesh.orgPnt(e)).lengthSq();
				if (dDist < dMin)
				{
					dMin = dDist;
					iMin = i;
				}
			}
		}

		return iMin;
	};

	std::cout << holes.size() << "\n";
	while (holes.size() > 1)
	{
		int iCloseset = getClosestEdge(holes);
		buildCylinderBetweenTwoHoles(mesh, holes[0][0], holes[iCloseset][0], { .metric = getEdgeLengthStitchMetric(mesh) });

		holes.erase(holes.begin() + iCloseset);
		holes.erase(holes.begin());
	}

	MeshSave::toAnySupportedFormat(mesh, "C:\\Users\\asdf\\Desktop\\stitchhole.stl");
	VxConvertMRMeshToVx(mesh, pCMeshIn);

	return true;
    }

appreciate for your helps!!

[Grantim]

Could you please try

 MeshBuilder::uniteCloseVertices(mesh, 0.0);
 buildCylinderBetweenTwoHoles(mesh, { .metric = getMinAreaMetric(mesh) }); // works only if there are exactly two holes left

after this line?

	MeshSave::toAnySupportedFormat(mesh, "C:\\Users\\asdf\\Desktop\\afterSolid.stl");

[JihoKang24]

hello grantim

I've confirmed that selecting "Hole" in the Mesh Healer parameters and clicking "Repair selected" produces different results compared to the code below:
`

do
{
	holes = findRightBoundary(mesh.topology);
	cnt = 0;
	for (auto hole : holes)
	{
		for (auto e : hole)
		{
			double dSumPerimeter = mesh.holePerimiter(e);
			if (dSumPerimeter < 4.34)
			{
				fillHoleParams.triangulateParams.metric = getPlaneFillMetric(mesh, e);
				MR::fillHoleNicely(mesh, e, fillHoleParams);
				cnt++;
			}
		}
	}
} while (cnt > 0);

MeshBuilder::uniteCloseVertices(mesh, 0.0);

`

I'd like to achieve the same border line repair behavior that the Mesh Healer tool does when "Hole" is selected. As shown below, the border line becomes clean and smooth after the repair:

before:

after:

How can I replicate this result via code?

always appreciate for your helps

[Grantim]

Hello!
Appling following code to afterSolid.stl should be equivalent to MeshHealer plus Stitch Holes from MeshInspector

Could you please try

 MeshBuilder::uniteCloseVertices(mesh, 0.0);
 buildCylinderBetweenTwoHoles(mesh, { .metric = getMinAreaMetric(mesh) }); // works only if there are exactly two holes left

after this line?

	MeshSave::toAnySupportedFormat(mesh, "C:\\Users\\asdf\\Desktop\\afterSolid.stl");

[JihoKang24]

I'm so late I solve this problem change hash percision
always appreciate your helps!!!!!