109 lines
3.3 KiB
C++
109 lines
3.3 KiB
C++
#include <igl/opengl/glfw/Viewer.h>
|
|
#include <igl/get_seconds.h>
|
|
#include <igl/readDMAT.h>
|
|
#include <igl/writeDMAT.h>
|
|
#include <igl/triangulated_grid.h>
|
|
#include <igl/cubic.h>
|
|
#include <igl/cycodebase/box_cubic.h>
|
|
#include <igl/cycodebase/point_spline_squared_distance.h>
|
|
|
|
int main(int argc, char * argv[])
|
|
{
|
|
IGL_TICTOC_LAMBDA;
|
|
Eigen::MatrixXd P;
|
|
Eigen::MatrixXi C;
|
|
igl::readDMAT(argc>1?argv[1]: TUTORIAL_SHARED_PATH "/libigl-cursive-P.dmat",P);
|
|
igl::readDMAT(argc>1?argv[2]: TUTORIAL_SHARED_PATH "/libigl-cursive-C.dmat",C);
|
|
|
|
// Cheezy conversion of (P,C) to polyline using 50 samples per curve
|
|
const int ns = 50;
|
|
Eigen::MatrixXd V(ns*C.rows(),P.cols());
|
|
Eigen::MatrixXi E((ns-1)*C.rows(),2);
|
|
for(int c = 0;c<C.rows();c++)
|
|
{
|
|
for(int i = 0;i<ns;i++)
|
|
{
|
|
const double t = double(i)/(ns-1);
|
|
// Evaluate cubic Bezier at parameter t
|
|
Eigen::RowVectorXd Vct;
|
|
igl::cubic(P(C.row(c),Eigen::all),t,Vct);
|
|
V.row(c*ns + i) = Vct;
|
|
if(i>0)
|
|
{
|
|
E.row(c*(ns-1) + (i-1))<< c*ns + (i-1), c*ns + i;
|
|
}
|
|
}
|
|
}
|
|
|
|
Eigen::MatrixXd B1,B2;
|
|
igl::cycodebase::box_cubic(P,C,B1,B2);
|
|
Eigen::MatrixXd BV(C.rows()*4,2);
|
|
BV<<
|
|
B1.col(0), B1.col(1),
|
|
B1.col(0), B2.col(1),
|
|
B2.col(0), B2.col(1),
|
|
B2.col(0), B1.col(1);
|
|
Eigen::MatrixXi BE(C.rows()*4,2);
|
|
for(int c = 0;c<C.rows();c++)
|
|
{
|
|
BE.row(c*4 + 0)<< 0*C.rows() + c, 1*C.rows() + c;
|
|
BE.row(c*4 + 1)<< 1*C.rows() + c, 2*C.rows() + c;
|
|
BE.row(c*4 + 2)<< 2*C.rows() + c, 3*C.rows() + c;
|
|
BE.row(c*4 + 3)<< 3*C.rows() + c, 0*C.rows() + c;
|
|
}
|
|
|
|
Eigen::MatrixXd GV;
|
|
Eigen::MatrixXi GF;
|
|
{
|
|
Eigen::RowVector2d min_corner = P.colwise().minCoeff();
|
|
Eigen::RowVector2d max_corner = P.colwise().maxCoeff();
|
|
// samples on x-axis
|
|
const int sx = 512;
|
|
const int sy = sx * (max_corner(1)-min_corner(1)) / (max_corner(0)-min_corner(0));
|
|
Eigen::RowVector2d diag = (max_corner - min_corner).eval();
|
|
const double bbd = diag.norm();
|
|
diag /= diag.norm();
|
|
min_corner -= 0.1*bbd*diag;
|
|
max_corner += 0.1*bbd*diag;
|
|
|
|
igl::triangulated_grid(sx,sy,GV,GF);
|
|
// Scale and translate grid to fit bounding box
|
|
GV.col(0) = (GV.col(0).array() * (max_corner(0)-min_corner(0))) + min_corner(0);
|
|
GV.col(1) = (GV.col(1).array() * (max_corner(1)-min_corner(1))) + min_corner(1);
|
|
}
|
|
|
|
|
|
Eigen::VectorXd sqrD,S;
|
|
Eigen::VectorXi I;
|
|
Eigen::MatrixXd K;
|
|
for(int r = 0;r<10;r++)
|
|
{
|
|
tictoc();
|
|
igl::cycodebase::point_spline_squared_distance(
|
|
GV, P, C, sqrD, I, S, K);
|
|
printf("%d points in %g secs\n",GV.rows(),tictoc());
|
|
}
|
|
Eigen::VectorXd D = sqrD.array().sqrt();
|
|
|
|
igl::opengl::glfw::Viewer vr;
|
|
vr.data().set_mesh(GV,GF);
|
|
vr.data().set_data(D);
|
|
vr.data().show_lines = false;
|
|
vr.append_mesh();
|
|
vr.data().set_mesh(V,E(Eigen::all,{0,1,1}).eval());
|
|
vr.data().set_points(P,Eigen::RowVector3d(1,0.7,0.2));
|
|
Eigen::MatrixXi CE(C.rows()*2,2);
|
|
CE<< C(Eigen::all,{0,1}),
|
|
C(Eigen::all,{2,3});
|
|
vr.data().set_edges(P,CE,Eigen::RowVector3d(1,0.7,0.2));
|
|
vr.data().point_size = 10;
|
|
//vr.data().set_points(BV,Eigen::RowVector3d(0,1,0));
|
|
//vr.data().set_edges(BV,BE,Eigen::RowVector3d(0,1,1));
|
|
vr.append_mesh();
|
|
vr.data().set_mesh(BV,BE(Eigen::all,{0,1,1}).eval());
|
|
vr.data().line_color = Eigen::RowVector4f(1,1,1,1);
|
|
|
|
vr.launch();
|
|
return 0;
|
|
}
|