DelaunayTri::DelaunayTri(const std::vector<cv::Point2i> &XYpoints,
const int &width, const int &height, const int &size)
{
this->width = width;
this->height = height;
img = cv::Mat::zeros(cv::Size(width, height), CV_8UC1);
int i;
points.resize(size);
for(i = 0; i < size; ++i)
{
points[i] = std::make_pair(PointCGAL(XYpoints[i].x, XYpoints[i].y), i + 1);
}
dt.insert(points.begin(),points.end());
int idxSize = 0;
size_ = dt.number_of_faces();
indices.resize(size_*3, -1);
centers.resize(size_, cv::Point2d(0, 0));
neighbors.resize(size_*3, 0.0);
radii.resize(size_, 0.0);
for(Delaunay::Finite_faces_iterator it = dt.finite_faces_begin(); it != dt.finite_faces_end(); it++) {
cv::Point2d t[3];
int idx[3];
Delaunay::Face_handle face = it;
int offset = idxSize * 3;
for(i = 0; i < 3; i++) {
t[i] = cv::Point2d(dt.triangle(face)[i].x(), dt.triangle(face)[i].y());
*(idx + i) = face->vertex(i)->info();
indices [offset + i] = face->vertex(i)->info();
}
this->insertMap(idx, idxSize);
cv::Point2d circumCenter = this->GetCircumcenter(t[0], t[1], t[2]);
centers[idxSize] = circumCenter;
radii[idxSize] = distPoint2Point<double>(circumCenter.x, circumCenter.y,
t[1].x, t[1].y);
idxSize++;
}
//this->draw_subdiv();
//cv::imwrite("delaunay.png", img);
this->computeNeighbors();
}
void mud::convertToSpringMass() {
std::vector<Point>::iterator begin = pts.begin();
std::vector<Point>::iterator end= pts.end();
dt.insert(begin,end);
qDebug() <<" n of v" << dt.number_of_vertices();
Edge_iterator eit;
MudParticle tmp_mp1(&mudSp), tmp_mp2(&mudSp);
unsigned int tmp_pmp1, tmp_pmp2;
int nEdge=0;
for (eit = dt.edges_begin();eit != dt.edges_end();eit++) {
Delaunay::Face_handle f = (*eit).first;
int i = (*eit).second;
Delaunay::Vertex_handle vh1 = f->vertex(dt.cw(i));
Delaunay::Vertex_handle vh2 = f->vertex(dt.ccw(i));
float tmpX = (*vh1).point().x();
float tmpY = (*vh1).point().y();
tmp_mp1.pos.setX(tmpX);
tmp_mp1.pos.setY(tmpY);
tmp_mp1.origPos= tmp_mp1.pos;
if (tmpX > -0.9 && tmpY > -0.9 && tmpX < 0.9 && tmpY < 0.9 )
tmp_mp1.type = 0;
else
tmp_mp1.type = 1;
tmpX = (*vh2).point().x();
tmpY = (*vh2).point().y();
tmp_mp2.pos.setX(tmpX);
tmp_mp2.pos.setY(tmpY);
tmp_mp2.origPos= tmp_mp2.pos;
if (tmpX > -0.9 && tmpY > -0.9 && tmpX < 0.9 && tmpY < 0.9 )
tmp_mp2.type = 0;
else
tmp_mp2.type = 1;
//Il faudrait mieux améliorer cette merde.
tmp_pmp1 = utils::findOrAdd(mudPa,&tmp_mp1);
tmp_pmp2 = utils::findOrAdd(mudPa,&tmp_mp2);
mudSp.push_back(new MudSpring (&mudPa, tmp_pmp1, tmp_pmp2));
mudPa[tmp_pmp1]->addSpring(mudSp.size() -1);
mudPa[tmp_pmp2]->addSpring(mudSp.size() -1);
nEdge++;
}
qDebug()<<"mudPa : "<<mudPa.size();
qDebug()<<"nEdge : "<<nEdge<< " et "<<mudSp.size();
}