[每日一码] （4）获得准确的AcDb3dSolid的包围盒

admin

1. void asdkasdktightboundbox()
2. {
3.     ads采用name eName;
4.     ads采用point pt;
7.     if( RTNORM != acedEntSel("\nPlease select a solid ", eName, pt ))
8.         return;
11.     AcDbObjectId id;
12.     acdbGetObjectId( id, eName );
13.     AcDb3dSolid*pSolid = NULL;
16.     acdbOpenObject(pSolid, id, AcDb::kForRead );
17.     if( pSolid == NULL ) {
18.         acutPrintf("\nObject not a solid");
19.         return;
20.     }
23.     double volume;
24.     AcGePoint3d centroid;
25.     double momInertia[3];
26.     double prodInertia[3];
27.     double prinMoments[3];
28.     AcGeVector3d prinAxes[3];
29.     double radiiGyration[3];
32.     AcDbExtents extents;
35.     Acad::ErrorStatus es;
36.     es = pSolid->getMassProp(
37.                 volume,
38.                 centroid,
39.                 momInertia,
40.                 prodInertia,
41.                 prinMoments,
42.                 prinAxes,
43.                 radiiGyration,
44.                 extents);
45.     AcGePoint3d max = extents.maxPoint();
46.     AcGePoint3d min = extents.minPoint();
49.     int xAxis, zAxis;
50.     double bigVol = (max.x - min.x ) * ( max. y - min.y ) * (max.z - min.z);
51.     acutPrintf("\ncalculated original box ");
52.     acutPrintf("\nvolume %e", bigVol );
53.     acutPrintf("\nOriginal box is %e the size of the object", bigVol/volume );
56.     // some work to make x the largest prin axis, and z the smallest.
57.     if(prinMoments[0] > prinMoments[1] )
58.     {
59.         xAxis = 0;
60.         zAxis = 1;
61.     }
62.     else
63.     {
64.         xAxis = 1;
65.         zAxis = 0;
66.     }
70.     if(prinMoments[zAxis] > prinMoments[2] )
71.         zAxis =  2;
72.     else if(prinMoments[xAxis] < prinMoments[2] )
73.         xAxis = 2;
77.     AcGeMatrix3d mat;
78.     AcGeVector3d yAxisVec = prinAxes[zAxis].crossProduct( prinAxes[xAxis]);
79.     mat.setCoordSystem( centroid, prinAxes[xAxis], yAxisVec, prinAxes[zAxis]);
80.     acutPrintf("\nxAxis %e %e %e", prinAxes[xAxis].x, prinAxes[xAxis].y,prinAxes[xAxis].z);
81.     acutPrintf("\nyAxis %e %e %e", yAxisVec.x, yAxisVec.y, yAxisVec.z);
82.     acutPrintf("\nzAxis %e %e %e", prinAxes[zAxis].x, prinAxes[zAxis].y,prinAxes[zAxis].z);
85.     AcDbEntity*pNewEnt = NULL;
88.     es = pSolid->getTransformedCopy( mat.inverse(), pNewEnt );
89.     if( pNewEnt == NULL )
90.     {
91.         acutPrintf("\nCannot transform solid");
92.         acutPrintf("\nError: %s", acadErrorStatusText(es));
93.         pSolid->close();
94.         return;
95.     }
98.     AcDbExtents newExtents;
99.     pNewEnt->getGeomExtents( newExtents );
100.     delete pNewEnt;
101.     //delete pNewEnt;
102.     max = newExtents.maxPoint();
103.     min = newExtents.minPoint();
104.     double deltaX = (max.x - min.x );
105.     double deltaY = ( max. y - min.y );
106.     double deltaZ = (max.z - min.z );
107.     double smallVol = deltaX * deltaY * deltaZ;
110.     acutPrintf("\ncalculated smaller box ");
111.     acutPrintf("\nvolume %e", smallVol );
112.     acutPrintf("\nAligned box is %e the size of the object", smallVol/volume );
113.     acutPrintf("\nOriginal box is %e the size of the aligned box", bigVol/smallVol);
116.     AcDb3dSolid*pNewSolid = new AcDb3dSolid;
117.     pNewSolid->createBox(deltaX, deltaY, deltaZ );
118.     AcGePoint3d newBoxCenter = min + 0.5 * (max - min);
119.     pNewSolid->transformBy( mat * AcGeMatrix3d::translation( newBoxCenter - AcGePoint3d::kOrigin ) );
120.     AcDbBlockTableRecord*pRecord;
121.     acdbOpenObject(pRecord, pSolid->ownerId(), AcDb::kForWrite );
122.     pSolid->close();
124.     if( pRecord == NULL )
125.     {
126.         acutPrintf("\nCannot open BTR");
127.         delete pNewSolid;
128.         return;
129.     }
132.     pRecord->appendAcDbEntity(pNewSolid);
133.     pNewSolid->close();   
134.     pRecord->close();
137. }

复制代码