CloudCompare forum

Code: Select all

#include "stdafx.h"
#include <iostream>
#include <fstream>
#include <string>
#include <algorithm>
#include <RegistrationTools.h>
#include <SimpleCloud.h>
#include <SimpleMesh.h>

int getLinesCount(std::string filePath) {
	std::ifstream inFile(filePath);
	return std::count(std::istreambuf_iterator<char>(inFile), std::istreambuf_iterator<char>(), '\n');
}

std::vector<double> getPointFromLine(std::string line) {
	std::vector<double> output(3);
	size_t pos = 0, k = 0;
	std::string temp, delimeter = ",";
	while ((pos = line.find(delimeter)) != std::string::npos) {
		temp = line.substr(0, pos);
		output[k] = std::atof(temp.c_str());
		line.erase(0, pos + delimeter.length());
		k += 1;
	}
	output[k] = std::atof(line.c_str());
	return output;
}

CCLib::SimpleCloud* getCloudFromFile(std::string filePath) {
	int linesCount = getLinesCount(filePath);
	CCLib::SimpleCloud *sc = new CCLib::SimpleCloud();
	sc->reserve(linesCount);
	std::ifstream read_file(filePath);
	std::string line;
	size_t j = 0;
	while (std::getline(read_file, line)) {
		std::vector<double> v = getPointFromLine(line);
		CCVector3 *vector = new CCVector3(v[0], v[1], v[2]);
		sc->addPoint(*vector);
	}
	return sc;
}

int main()
{
	CCLib::SimpleCloud *tgt = getCloudFromFile("C:\\pathTo\\file1.csv");
	CCLib::SimpleMesh *tgt_mesh = new CCLib::SimpleMesh(tgt);
	CCLib::SimpleCloud *src = getCloudFromFile("C:\\pathTo\\file2.csv");
	CCLib::ICPRegistrationTools::Parameters *parameters = new CCLib::ICPRegistrationTools::Parameters();
	CCLib::ICPRegistrationTools::ScaledTransformation *scaled_transofrmation = new CCLib::ICPRegistrationTools::ScaledTransformation();
	double RMS = 0.00; unsigned int point_count = 0;
	std::cout << "Start registration" << std::endl;
	CCLib::ICPRegistrationTools::Register(tgt, tgt_mesh, src, *parameters, *scaled_transofrmation, RMS, point_count);
	return 0;
}

Code: Select all

	CCLib::SimpleCloud *tgt = getCloudFromFile("C:\\pathTo\\file1.csv");
	CCLib::SimpleMesh *tgt_mesh = new CCLib::SimpleMesh(tgt, true);
	tgt_mesh->resize(getLinesCount("C:\\pathTo\\file1.csv"));
	std::cout << tgt_mesh->capacity() << " " << tgt_mesh->size() << std::endl;
	std::ofstream out("C:\\somefile.txt");
	for (int i = 0; i < tgt_mesh->size(); i += 1) {
		CCVector3 *A = new CCVector3();
		CCVector3 *B = new CCVector3();
		CCVector3 *C = new CCVector3();
		tgt_mesh->getTriangleVertices(i, *A, *B, *C);
		out << std::to_string(A->x) << " " << std::to_string(A->y) << " " << std::to_string(A->z) << "\t";
		out << std::to_string(B->x) << " " << std::to_string(B->y) << " " << std::to_string(B->z) << "\t";
		out << std::to_string(C->x) << " " << std::to_string(C->y) << " " << std::to_string(C->z) << std::endl;1
		delete A,B,C;
	}
	out.close();
	

Code: Select all

int main()
{
	std::cout << " Start ipc_useCC " << std::endl;
	CCLib::SimpleCloud *tgt = getCloudFromFile("C:\\pathTo\\file1.csv");
	std::cout << "Target point cloud ready" << std::endl;
	CCLib::SimpleCloud *src = getCloudFromFile("C:\\pathTo\\file2.csv");
	std::cout << "Source point cloud ready" << std::endl;
	CCLib::ICPRegistrationTools::ScaledTransformation *scaled_transofrmation = new CCLib::ICPRegistrationTools::ScaledTransformation();
	double RMS = 0.00; unsigned int pointcount = 0;
	CCLib::ICPRegistrationTools::Parameters *params = new CCLib::ICPRegistrationTools::Parameters();
	std::cout << "...Start registration..." << std::endl;
	std::chrono::high_resolution_clock::time_point t1 = std::chrono::high_resolution_clock::now();
	tgt->enableScalarField();
	src->enableScalarField();
	CCLib::ICPRegistrationTools::Register(tgt, 0, src, *params, *scaled_transofrmation, RMS, pointcount);
	std::chrono::high_resolution_clock::time_point t2 = std::chrono::high_resolution_clock::now();
	auto duration = std::chrono::duration_cast<std::chrono::seconds>(t2 - t1).count();
	std::cout << "Registration complete. Elapsed time: " << duration2 << std::endl;
	return 0;
}