How to create a surface?

Hi all,

 

I have a simple question. I tried to find out few materials over internet. But I could not gather much information on this. How can I build/construct a surface based on a set of 3D points and edges connecting these points? I tried the example shown in this link here:

 

http://opencascade.blogspot.com/2010/03/surface-modeling-part6.html

 

This is what I have done so far, without positive results. Experienced people kindly give some information on this and let me know what is it that I am missing.

 

//constants for algorithm

    const Standard_Integer aNbIter = 5; //number of algorithm iterations

    const Standard_Integer aNbPnts = 10; //sample points per each constraint

    const Standard_Integer aDeg = 3; //requested surface degree ? const Standard_Integer aMaxDeg = theMaxDeg;

    const Standard_Integer aMaxSeg = 10000;

    const Standard_Real aTol3d = 1.e-04;

    const Standard_Real aTol2d = 1.e-05;

    const Standard_Real anAngTol = 1.e-02; //angular

    const Standard_Real aCurvTol = 1.e-01; //curvature

 

    gp_Pnt p1(0, 0, 0);

    gp_Pnt p2(0, 10, 10);

    gp_Pnt p3(10, 10, 5);

    gp_Pnt p4(10, 0, -5);

 

    Handle(Geom_TrimmedCurve) aSegment1 = GC_MakeSegment(p1, p2);

    Handle(Geom_TrimmedCurve) aSegment2 = GC_MakeSegment(p2, p3);

    Handle(Geom_TrimmedCurve) aSegment3 = GC_MakeSegment(p3, p4);

    Handle(Geom_TrimmedCurve) aSegment4 = GC_MakeSegment(p4, p1);

 

    TopoDS_Edge anEdge1 = BRepBuilderAPI_MakeEdge(aSegment1);

    TopoDS_Edge anEdge2 = BRepBuilderAPI_MakeEdge(aSegment2);

    TopoDS_Edge anEdge3 = BRepBuilderAPI_MakeEdge(aSegment3);

    TopoDS_Edge anEdge4 = BRepBuilderAPI_MakeEdge(aSegment4);

 

    TopoDS_Face myFace;

 

    Handle(GeomAdaptor_HCurve) aCurve1 = new GeomAdaptor_HCurve(

                GeomAdaptor_HCurve(aSegment1));

    Handle(GeomAdaptor_HCurve) aCurve2 = new GeomAdaptor_HCurve(

                GeomAdaptor_HCurve(aSegment2));

    Handle(GeomAdaptor_HCurve) aCurve3 = new GeomAdaptor_HCurve(

                GeomAdaptor_HCurve(aSegment3));

    Handle(GeomAdaptor_HCurve) aCurve4 = new GeomAdaptor_HCurve(

                GeomAdaptor_HCurve(aSegment4));

 

    TColStd_ListOfTransient theBoundaries;

    theBoundaries.Append(aCurve1);

    theBoundaries.Append(aCurve2);

    theBoundaries.Append(aCurve3);

    theBoundaries.Append(aCurve4);

 

    Handle(Geom_Surface) aRes;

 

    GeomPlate_BuildPlateSurface aPlateBuilder;

 

    TColStd_ListIteratorOfListOfTransient anIt (theBoundaries);

    if (anIt.More())

    {

        int i = 1;

        for (; anIt.More(); anIt.Next(), i++)

        {

            const Handle(Standard_Transient)& aCur = anIt.Value();

            if (aCur.IsNull())

            {

                std::cout << "IS NULL" << std::endl;

                Standard_ConstructionError::Raise ("ACISAlgo::MakeSurface()");

                return;

            }

            else if (aCur->IsKind (STANDARD_TYPE (Adaptor3d_HCurveOnSurface)))

            {

                //G1 constraint

                const Handle(Adaptor3d_HCurveOnSurface)& aHCOS = Handle(Adaptor3d_HCurveOnSurface)::DownCast (aCur);

                Handle (GeomPlate_CurveConstraint) aConst =

                        new GeomPlate_CurveConstraint (aHCOS, 1 /*GeomAbs_G1*/, aNbPnts, aTol3d, anAngTol, aCurvTol);

                aPlateBuilder.Add (aConst);

            }

            else if (aCur->IsKind (STANDARD_TYPE (GeomAdaptor_HCurve)))

            {

                //G0 constraint

                const Handle(GeomAdaptor_HCurve)& aHC = Handle(GeomAdaptor_HCurve)::DownCast (aCur);

                Handle (GeomPlate_CurveConstraint) aConst =

                        new GeomPlate_CurveConstraint (aHC, 0 /*GeomAbs_G0*/, aNbPnts, aTol3d);

                aPlateBuilder.Add (aConst);

            }

            else

            {

                Standard_TypeMismatch::Raise ("ACISAlgo::MakeSurface()");

            }

        }

    }

    else

    {

        Standard_ConstructionError::Raise ("ACISAlgo::MakeSurface()");

    }

    //construct

    aPlateBuilder.Perform();

    if (!aPlateBuilder.IsDone())

    {

        // return aRes;

        std::cout << "Failed building plate" << std::endl;

        return;

    }

 

    const Handle(GeomPlate_Surface)& aPlate = aPlateBuilder.Surface(); //approximation (see BRepFill_Filling - when no initial surface was given)

    TopoDS_Shape resShape = BRepBuilderAPI_MakeFace(aPlateBuilder.SurfInit(), 0.00001);

    BRepTools::Write(resShape, "makeFace1.brep");

 

I am not yet sure whether I am on the right path or not.

 

Thanks in advance

Rakesh Patil