关于osgdelaunay的疑问

answer1103

osg::Vec2Array *LinearConstraint::getRoadTexcoords(const osg::Vec3Array *points)  {
    // need to create a vec2 array from the coordinates that fits the road
    osg::Vec3Array::const_iterator tritr;//常迭代器，这样就不能修改points里的内容了
    osg::ref_ptr<osg::Vec2Array> tcoords= new osg::Vec2Array ;

    for (tritr=points->begin(); tritr!=points->end();tritr++) {
        osg::Vec2 tci(-1.,-1.);
        int ib=0;
        // osg::Vec3Array *varr=dynamic_cast<osg::Vec3Array*>(getVertexArray());
        bool ptfound=false;
        for (osg::Vec3Array::iterator vit=_edgecoords->begin(); vit!= _edgecoords->end() && !ptfound; vit++) {
            if ((*vit)==(*tritr)) {
                tci=_tcoords->at(ib);
                ptfound=true;
            }
            ib++;
        }///找到特定顶点的纹理坐标

        if (!ptfound) { // search for surrounding points and interpolate
            ib=0;
            osg::Vec3 pminus=(_edgecoords->back()); // need pminus for interpolation
            int ibm1=_edgecoords->size()-1;
            for (osg::Vec3Array::iterator vit=_edgecoords->begin(); vit!= _edgecoords->end() /*&& !ptfound*/; vit++) {
                osg::Vec3 pplus=(*vit)-(*tritr);
                osg::Vec3 dpm=pminus-(*tritr);
                pplus.set (pplus.x(),pplus.y(),0);
                dpm.set (dpm.x(),dpm.y(),0);
                float dprod=pplus*dpm/(pplus.length() * dpm.length());
                if (dprod<-0.9999) { // *tritr lies between....
                    osg::Vec2 tminus=_tcoords->at(ibm1);
                    osg::Vec2 tplus=_tcoords->at(ib);
                    float frac=(dpm.length()/(dpm.length()+pplus.length()));
                    tci=tminus+((tplus-tminus)*frac);
                    ptfound=true;
                }
                ibm1=ib;
                ib++;   
                pminus=(*vit);
            }
        }
        tcoords->push_back(tci);
    }
    // some extra points are not interpolated as they lie between 2 interpolated vertices
    for (tritr=points->begin(); tritr!=points->end();tritr++) {
        int ib=tritr-points->begin();
        osg::Vec2 tci=tcoords->at(ib);
        if (tci.x()<-.99 && tci.y()<-.99) {
            // search through each of the primitivesets
            osg::Vec3Array::const_iterator ptitr;
            //    osg::notify(osg::WARN) << "Not calculated " <<  (*tritr).x() <<"," << (*tritr).y() << std::endl;
            for (ptitr=points->begin(); ptitr!=points->end();ptitr++) {
            }
        }
    }
    return tcoords.release();
}

这个函数里面有些地方不是很好懂，特别是   if (!ptfound) 循环里面的这段代码，如果说它是想得到道路周围的点的纹理坐标，那为什么在里面的
if (dprod<-0.9999) { // *tritr lies between....
                    osg::Vec2 tminus=_tcoords->at(ibm1);
                    osg::Vec2 tplus=_tcoords->at(ib);
                    float frac=(dpm.length()/(dpm.length()+pplus.length()));
                    tci=tminus+((tplus-tminus)*frac);
                    ptfound=true;
                }
里面没有push_back?请研究过这段代码的前辈给我点指导，谢谢了。

array

我不明白这里为什么要有push_back，您想要实现什么功能？

answer1103

2# array
您好，我现在想通过一个中心线来构造一条道路，所以在研究这个例子，我不是很明白这个函数的功能，就是这个函数得到的纹理坐标是哪些点的纹理坐标？还有那个-0.9999是为了得到地形上无限靠近道路的点的纹理坐标吗？如果是的话，那我觉得那个里面就应该有一个push_back把得到的纹理坐标添加进去啊，请您不吝赐教。

array

抱歉我没有研究过相关的代码，在我看来这似乎是为了获取某个计算平均值并把它累加到一些变量中

answer1103

4# array
噢，还是谢谢您，我再去研究下，感觉这个对数学的要求比想象中的还高。