[求助]：求线相交问题

shengl

怎么用flag=mdlIntersect_allBetweenElms(pInsO,pInsT,g_num,OneEdp,TwoEdp,NULL,0.00001);求复杂链相交时,得出的交点远远超过我在屏幕上所见的交点;另外如果交点为线上的点时,交点个数为偶数.真让人郁闷!

gp341

这个函数我用过 没有问题 你的限差太小了把
因为咱们的地形图一般都是米做主单位 所以给1就已经很小了

MicroStation迷

把你的参数定义以及构造OneEdp,TwoEdp这两个原素描述符的程序段也发上来让大家帮你分析分析。

shengl

POS_S *qPos,*pPos;
        LOCATE_S *qLoc,*pLoc;
        int flag;
        Dpoint3d pInsO[101],pInsT[101],PL[101];
        int numO,numT,numVT,numVO,i,j,k,n,ba,ii;
        Dpoint3d *ptO,*ptT;
        MSElementDescr *OneEdp,*TwoEdp;
        int Value,Tflag,Oflag;
       
        g_pLocate=malloc(sizeof(LOCATE_S));
        g_pLocate->next=NULL;
        pLoc=g_pLocate;
       
        while(g_pPos->next!=NULL)
        {
                pPos=g_pPos->next;
                qPos=pPos->next;
                      
                       if(mdlElmdscr_read(&OneEdp,pPos->filePos,MASTERFILE,FALSE,NULL)==0L)
                       {
                               mdlDialog_openInfoBox("程序运行遇到错误");
                               return -1;
                       }
                while(qPos!=NULL)
                {       
                        if(mdlElmdscr_read(&TwoEdp,qPos->filePos,MASTERFILE,FALSE,NULL)==0L)
                               {
                                       mdlDialog_openInfoBox("程序运行遇到错误");
                                       return -1;
                               }
                               flag=mdlIntersect_allBetweenElms(pInsO,pInsT,g_num,OneEdp,TwoEdp,NULL,1);

MicroStation迷

进入调试方式查看一下pPos->filePos、qPos->filePos的值是否是您所期望的？
另外，要适时地调用mdlElmdscr_freeAll ()去释放OneEdp和TwoEdp。

shengl

这些问题能过解释果交点为线上的点时,交点个数为偶数的问题吗?

MicroStation迷

交点个数为偶数？
是否求交的两线有一段重叠？
好像有重重叠时会返回两个交点(我有些记不清了)。