3D打印技术之切片引擎(4)

3D打印技术之切片引擎(4)

文章转自:http://blog.csdn.net/fourierFeng/article/details/46753651

这一篇文章我讲一下多边打印的问题,多边打印是切片引擎的一项关键的技术。 

双边打印

图1 双边打印

首先,它可以保证打印实体表面免受内部填充的冲击,保证外观的真实度;其次,在上层在相对于下层倾斜较大时,多边打印可以很好的起到支撑的作用,避免上层塌陷。 但是,目前来说,我的多边打印还不够普适,对一些不规范的模型,以及模型中非常尖锐的特征效果并不好,对绝大部分的较为平滑的模型是完全没有问题的。 

下面就简单说一下它的原理:假设边界中的任意相邻的向量AB和BC,这里要找的是点d(角ABC中心线上的一点),看下图:

找出边界环中每一个相邻向量的d点,工作就基本完成了。所以原理非常简单,只不过是很多琐碎的细节需要处理好,比如说向量Bd的方向问题,B点和d点的欧氏距离等等,不能再说了,再说就有误导人之嫌,其实这段时间细想,我的切片引擎的这些想法并无优秀可言,可是还是想把那段时间的工作记录下来,仅此而已,闲话说几句,这段时间在做与测绘相关的算法,本来认为可以写成博文与大家分享的东西,结果硬是变成了核心期刊上的论文了,哎,只是觉得那个东西离论文水平的创新还有不小的距离啊,看来,国内的论文质量……,呵呵。下面呈上代码,其中getInnerBoundary函数有些不够规范,用了goto跳转,而且一些细节也不是最科学,当初主要是赶进度,为了适应公司的建模能力不足。

void getAngularBisector(float3 &bisector,float3 point,float3 leftPoint,float3 rightPoint)
{
    float3 v1,v2;
    float norm1,norm2,angle;
    get_vector_diff(v1,leftPoint,point);
    get_vector_diff(v2,rightPoint,point);
    if(v1[0]*v2[1]==v1[1]*v2[0])
    {
        if(v1[1]==v2[1])
        {
            bisector[0]=0;
            bisector[1]=1;
            bisector[2]=0;
        }
        else if(v1[0]==v2[0])
        {
            bisector[0]=1;
            bisector[1]=0;
            bisector[2]=0;
        }
        else
        {
            bisector[0]=point[0]-1;
            bisector[1]=((v1[0]-v2[0])+(v1[1]-v2[1])*point[1])/(v1[1]-v2[1]);
            bisector[2]=0;
        }

    }
    else
    {
        getNormalizeVector(v1);
        getNormalizeVector(v2);
        if(v1[0]*v2[1]==v1[1]*v2[0])
        {
            if(v1[1]==v2[1])
            {
                bisector[0]=0;
                bisector[1]=1;
                bisector[2]=0;
            }
            else if(v1[0]==v2[0])
            {
                bisector[0]=1;
                bisector[1]=0;
                bisector[2]=0;
            }
            else
            {
                bisector[0]=point[0]-1;
                bisector[1]=((v1[0]-v2[0])+(v1[1]-v2[1])*point[1])/(v1[1]-v2[1]);
                bisector[2]=0;
            }

        }
        else
        {
            get_vector_sum(bisector,v1,v2);
        }
    }
}


void getInnerPoint(float3 &innerPoint,float3 point,float3 leftPoint,float3 rightPoint,float margin)
{
    float3 bisector;
    float distance,rate;
    distance=margin;
    getAngularBisector(bisector,point,leftPoint,rightPoint);    rate=distance/sqrt(pow(bisector[0],2)+pow(bisector[1],2)+pow(bisector[2],2));
    innerPoint[0]=point[0]+bisector[0]*rate;
    innerPoint[1]=point[1]+bisector[1]*rate;
    innerPoint[2]=point[2]+bisector[2]*rate;

    if(get_vector3_det(point,rightPoint,innerPoint)<0)
    {

        innerPoint[0]=point[0]-bisector[0]*rate;
        innerPoint[1]=point[1]-bisector[1]*rate;
        innerPoint[2]=point[2]-bisector[2]*rate;
    }
}



void getInnerPoint(Phasor *innerPhasor,Phasor *phasor,Phasor*leftPhasor,Phasor *rightPhasor,float margin,int fillMaterial)
{
    getInnerPoint(innerPhasor->beginPoint,phasor->beginPoint,leftPhasor->beginPoint,phasor->endPoint,margin);

    getInnerPoint(innerPhasor->endPoint,phasor->endPoint,phasor->beginPoint,rightPhasor->endPoint,margin);

    innerPhasor->material=fillMaterial;
}




void getInnerBoundary(Phasor *&innerPhasors,Phasor *phasors,int phasor_num,vector<vector<int> > &closedSet
    ,float lineHeight,int fillMaterial,bool *innerBoundaryStatus)
{
    innerPhasors=new Phasor[phasor_num];

    int index,indexMain,phasorIndex,phasorIndexLeft,phasorIndexRight;

    float margin=lineHeight;
againScan:
    if(margin<0.1)
    {
        *innerBoundaryStatus=false;
        innerPhasors=NULL;
        return;
    }

    for(indexMain=0;indexMain!=closedSet.size();++indexMain)
    {
        int count=0;
        for(index=0;index!=closedSet[indexMain].size();++index)
        {
            phasorIndex=closedSet[indexMain][index];

            if(index>0)
            {
                phasorIndexLeft=closedSet[indexMain][index-1];
            }
            else
            {
                phasorIndexLeft=closedSet[indexMain][closedSet[indexMain].size()-1];
            }
            phasorIndexRight=closedSet[indexMain][(index+1)%closedSet[indexMain].size()];

            getInnerPoint(innerPhasors+phasorIndex
                ,phasors+phasorIndex
                ,phasors+phasorIndexLeft
                ,phasors+phasorIndexRight
                ,margin,fillMaterial);

            for(int i=0;i!=phasor_num;++i)
            {
                if(i==phasorIndex)
                {
                    continue;
                }

                if(get_vector_distance2(innerPhasors[phasorIndex].beginPoint,phasors[phasorIndex].beginPoint)>
                    get_vector_distance2(innerPhasors[phasorIndex].beginPoint,phasors[i].beginPoint))
                {
                    count++;
                    break;
                }
            }

            if(count>closedSet[indexMain].size()/5+1)
            {
                margin-=0.05;
                goto againScan;
            }
        }
    }
}
展开阅读全文

没有更多推荐了,返回首页