Getting two shapes When there are more then One top level shape

Hi,

I am creating a OpenSceneGraph plugin to load IGES files using OpenCascade. I am using OpenCascade to read and triangulate IGES files. My problem is that if a model contains more than one top level images I am getting some shapes twice.

Following is a screenshot of CAD exchanger for this model
[URL=http://imageshack.com/i/f7jhnlp][IMG]http://imageshack.com/a/img547/7215...

and here is how i am getting image looks like in Osg Viewer
[URL=http://imageshack.com/i/j652dbp][IMG]http://imageshack.com/a/img690/600/...

Here is how i am reading this IGES file using XDE

Handle(TDocStd_Document) doc;
XCAFApp_Application::GetApplication()->NewDocument("MDTV-XCAF", doc);

IGESCAFControl_Reader reader;
reader.SetColorMode(true);
reader.SetNameMode(true);
reader.SetLayerMode(true);

//IGESControl_Reader Reader;
reader.ReadFile( (Standard_CString)filePath.c_str() );
// transfer data from reader to doc
if(!reader.Transfer(doc))
{
cout return NULL;
}

// To get a node considered as an Assembly from an XDE structure, you can use the Label of the node.
Handle (XCAFDoc_ShapeTool) assembly;
assembly = XCAFDoc_DocumentTool::ShapeTool(doc->Main());

// To query, edit, or initialize a Document to handle Colors of XCAF
Handle(XCAFDoc_ColorTool) colors;
colors = XCAFDoc_DocumentTool::ColorTool(doc->Main());

this is where i read colors from document

/// Try to Get shape color from top level shape
Quantity_Color color;
if (!colors->GetColor(frshapes.Value(i), color))
{
/// If no color found in top level shape then iterate children
/// Sometimes colors are not attached to top level shapes
/// So we will iterate children of this shape to find a color
if (frshapes.Value(i).HasChild())
{
TDF_ChildIterator it;
for(it.Initialize(frshapes.Value(i)); it.More(); it.Next())
{
if(colors->GetColor(it.Value(), XCAFDoc_ColorGen, color) ||
colors->GetColor(it.Value(), XCAFDoc_ColorSurf, color) ||
colors->GetColor(it.Value(), XCAFDoc_ColorCurv, color) )
{
/// if we find a color save it for this shape and break
shapeColor = osg::Vec3(color.Red(), color.Green(), color.Blue());
break;
}
}
}
}

triangulation and shape healing is as follows

BRepTools::Clean(aShape);

ShapeFix_Shape fixer(aShape);
fixer.Perform();
BRepBuilderAPI_Sewing sew;
sew.Add(fixer.Shape());
sew.Perform();
aShape = sew.SewedShape();

/// call to incremental mesh on this shape
BRepMesh::Mesh(aShape, 1);

Finally iterating over faces and getting vertices and indices

for (TopExp_Explorer ex(aShape, TopAbs_FACE); ex.More(); ex.Next())
{
TopoDS_Face face = TopoDS::Face(ex.Current());
TopLoc_Location loc;
/// triangulate current face
Handle (Poly_Triangulation) triangulation = BRep_Tool::Triangulation(face, loc);

if (!triangulation.IsNull())
{
int noOfNodes = triangulation->NbNodes();
// Store vertices. Build vertex array here
for(unsigned int j = 1; j NbNodes(); j++)
{
// populate vertex list
// Ref: http://www.opencascade.org/org/forum/thread_16694/?forum=3
gp_Pnt pt = (triangulation->Nodes())(j).Transformed(loc.Transformation());
vertexList->push_back(osg::Vec3(pt.X(), pt.Y(), pt.Z()));

colorList->push_back(shapeColor);
}

/// now we need to get face indices for triangles
// get list of triangle first
const Poly_Array1OfTriangle& triangles = triangulation->Triangles();

//No of triangles in this triangulation
noOfTriangles = triangulation->NbTriangles();

Standard_Integer v1, v2, v3;
for (unsigned int j = 1; j {
/// If face direction is reversed then we add verticews in reverse order
/// order of vertices is important for normal calculation later
if (face.Orientation() == TopAbs_REVERSED)
{
triangles(j).Get(v1, v3, v2);
}
else
{
triangles(j).Get(v1, v2, v3);
}
triangleStrip->push_back(index + v1 - 1);
triangleStrip->push_back(index + v2 - 1);
triangleStrip->push_back(index + v3 - 1);
}
index = index + noOfNodes;
}
}
}
Previous posts by users on these forums have helped great so far. I just started with OpenCascade.
Please Help