设为首页收藏本站转
发帖领奖家园QQ群申请版主版主评价切换到宽版

 找回密码
 立即注册

QQ登录

只需一步，快速开始

扫一扫，访问微社区

晓东CAD家园-论坛»论坛 › ╃CAD编程开发 › ObjectARX › 实体集的最小矩形框。在钣金展开异形件，备料时，也许多 ...
返回列表 发新帖
查看: 1672|回复: 4

[分享] 实体集的最小矩形框。在钣金展开异形件，备料时，也许多得上

[复制链接]
w15994259099

已领礼包: 271个

财富等级: 日进斗金
发表于 2016-12-15 15:19:29 | 显示全部楼层 |阅读模式
【广而告之】积分规则、D豆获取攻略       窝窝空间装扮攻略

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要 登录 才可以下载或查看，没有账号？立即注册

×
  1. /*在网上看到:Lee Mac LISP_Programming写的相关LISP程序，较好，学习其算法思路。
  2. 但在多次均匀旋转角度在[0，180]范围内查询最小面积时，也许会正好错过实体上的轮廓特殊点线，导致计算结果总与实际轮廓线，总差那么一丁点。
  3. 我的程序，加进了对轮廓线上的特殊点角度的计算处理，解决了上述问题!!!
  4. */
  5. void draw_Min_EntSet_Box()
  6. {
  7.     AcDbObjectIdArray objIdArr;
  8.     ads_name ss={0L,0L};
  9.   //选择实体
  10.   if(acedSSGet(NULL,NULL,NULL,NULL,ss)!=RTNORM)
  11.   {
  12.     acutPrintf("\n用户取消!");
  13.     return;
  14.   }
  15.     long len=0;
  16.   acedSSLength(ss,&len);
  17.   for(int i=0;i<len;i++)
  18.   {
  19.     ads_name ent={0L,0L};
  20.     acedSSName(ss,i,ent);
  21.     AcDbObjectId objId;
  22.     acdbGetObjectId(objId,ent);
  23.         objIdArr.append(objId);
  24.   }
  25.   acedSSFree(ss);
  26.    //获得实体组端点线组成的特殊角度:从小到大
  27.     CArray<double,double> angArr;
  28.     getEntsPtAngle(objIdArr,angArr);
  29.   //
  30.     //初始状态:获得实体组的范围
  31.   AcDbExtents ext0;
  32.     Get_AllEnt_Extents(objIdArr,ext0);
  33.     AcGePoint3d centerPt=ext0.minPoint()+(ext0.maxPoint()-ext0.minPoint())/2.0;
  34.   //最小值
  35.   double minAngle=0.0;
  36.     AcDbExtents minExt=ext0;
  37.     //
  38.     int n=2000;
  39.     //acutPrintf("\n进入角度循环!");
  40.   for(int j=1;j<n;j++)
  41.   {
  42.        double curAngle=j*PI/n;
  43.      //acutPrintf("\n当前旋转角:%f",curAngle);
  44.        AcDbExtents curExt;
  45.        rotateEntSet(objIdArr,centerPt,curAngle);
  46.        Get_AllEnt_Extents(objIdArr,curExt);
  47.        rotateEntSet(objIdArr,centerPt,-curAngle);
  48.      double curArea,minArea;
  49.        curArea=(curExt.maxPoint()[0]-curExt.minPoint()[0])*(curExt.maxPoint()[1]-curExt.minPoint()[1]);
  50.        minArea=(minExt.maxPoint()[0]-minExt.minPoint()[0])*(minExt.maxPoint()[1]-minExt.minPoint()[1]);
  51.      //acutPrintf("\n当前面积:%f;最小面积:%f",curArea,minArea);
  52.      if(curArea<minArea)
  53.      {
  54.           minExt.set(curExt.minPoint(),curExt.maxPoint());;
  55.           minAngle=curAngle;
  56.      }
  57.        //
  58.      if(j==n-1)
  59.      {
  60.        acutPrintf("\n暴力浏览:%d 次完毕!",n);
  61.      }
  62.   }
  63.     //
  64.   //acutPrintf("\n特殊点的角度数:%d",angArr.GetSize());
  65.     //特殊角度
  66.   for(int k=0;k<angArr.GetSize();k++)
  67.   {
  68.        double curAngle=angArr[k];
  69.        AcDbExtents curExt;
  70.        rotateEntSet(objIdArr,centerPt,-curAngle);
  71.        Get_AllEnt_Extents(objIdArr,curExt);
  72.        rotateEntSet(objIdArr,centerPt,curAngle);
  73.      double curArea,minArea;
  74.        curArea=(curExt.maxPoint()[0]-curExt.minPoint()[0])*(curExt.maxPoint()[1]-curExt.minPoint()[1]);
  75.        minArea=(minExt.maxPoint()[0]-minExt.minPoint()[0])*(minExt.maxPoint()[1]-minExt.minPoint()[1]);
  76.      //acutPrintf("\n当前面积:%f;最小面积:%f",curArea,minArea);
  77.      if(curArea-minArea<0.0)
  78.      {
  79.           minExt.set(curExt.minPoint(),curExt.maxPoint());;
  80.           minAngle=-curAngle;
  81.       acutPrintf("\n特殊点角度的最小值交换:成功!");
  82.      }
  83.      else
  84.      {
  85.       ;//acutPrintf("\n特殊点的最小值:不交换!");
  86.      }

  88.      if(k==angArr.GetSize()-1)
  89.      {
  90.        acutPrintf("\n特殊点角度浏览:%d 次完毕!",angArr.GetSize());
  91.      }
  92.   }
  93.     //特殊角度
  94.     //
  95.   //建立矩形
  96.     //AcDbExtents minExt;
  97.     AcDbPolyline* pPLine=new AcDbPolyline();
  98.     pPLine->addVertexAt(0,AcGePoint2d(minExt.minPoint()[0],minExt.minPoint()[1]));
  99.     pPLine->addVertexAt(1,AcGePoint2d(minExt.maxPoint()[0],minExt.minPoint()[1]));
  100.     pPLine->addVertexAt(2,AcGePoint2d(minExt.maxPoint()[0],minExt.maxPoint()[1]));
  101.     pPLine->addVertexAt(3,AcGePoint2d(minExt.minPoint()[0],minExt.maxPoint()[1]));
  102.   pPLine->setClosed(Adesk::kTrue);
  103.     //加入
  104.   AcDbBlockTable* pBT=NULL;
  105.   AcDbBlockTableRecord* pBTRec=NULL;
  106.   acdbHostApplicationServices()->workingDatabase()->getBlockTable(pBT,AcDb::kForRead);
  107.   pBT->getAt(ACDB_MODEL_SPACE,pBTRec,AcDb::kForWrite);
  108.   pBT->close();
  109.   pBTRec->appendAcDbEntity(pPLine);
  110.   //转换
  111.   //acutPrintf("\n旋转角:%f",minAngle*180/PI);
  112.     AcGeMatrix3d xform;
  113.   xform.setToRotation(-minAngle,AcGeVector3d(0,0,1),centerPt);
  114.     pPLine->transformBy(xform);
  115.   //
  116.   pPLine->close();
  117.   pBTRec->close();
  118.     //
  119. }


钣金

相关帖子

论坛插件加载方法
发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；
如果你在论坛求助问题，并且已经从坛友或者管理的回复中解决了问题，请把帖子标题加上【已解决】；
如何回报帮助你解决问题的坛友，一个好办法就是给对方加【D豆】，加分不会扣除自己的积分，做一个热心并受欢迎的人！
回复

使用道具 举报

w15994259099

已领礼包: 271个

财富等级: 日进斗金
 楼主| 发表于 2016-12-15 15:35:57 | 显示全部楼层
  1. //实体组的旋转变换
  2. bool rotateEntSet(AcDbObjectIdArray& objIdArr,AcGePoint3d centerPt,double angle)
  3. {
  4.   for(int i=0;i<objIdArr.logicalLength();i++)
  5.   {
  6.         AcGeMatrix3d xform;
  7.     xform.setToRotation(angle,AcGeVector3d(0,0,1),centerPt);
  8.     AcDbEntity* pEnt=NULL;
  9.     acdbOpenAcDbEntity(pEnt,objIdArr,AcDb::kForWrite);
  10.     pEnt->transformBy(xform);
  11.     pEnt->close();
  12.   }
  13.   //
  14.   return true;
  15. }

  17. // 20 获得块引用实体的范围
  18. Acad::ErrorStatus Get_Ent_Extent(const AcDbObjectId& idEnt, AcDbExtents& extents)
  19. {
  20.   Acad::ErrorStatus es;
  21.   AcDbEntity *pEnt = NULL;
  22.   if (Acad::eOk == acdbOpenObject(pEnt, idEnt, AcDb::kForRead))
  23.   {
  24.     if(pEnt->isKindOf(AcDbBlockReference::desc()))
  25.     {
  26.    
  27.         AcDbBlockReference *pBlkRef = AcDbBlockReference::cast(pEnt);
  28.             es = pBlkRef->geomExtentsBestFit(extents);
  29.     }

  31.     else
  32.     {
  33.       es = pEnt->getGeomExtents(extents);
  34.     }
  35.     pEnt->close();
  36.   }
  37.   
  38.   return (es);
  39. }

  41. // 21 获得实体组的范围
  42. void Get_AllEnt_Extents(const AcDbObjectIdArray& aridEnt, AcDbExtents& extents)
  43. {
  44.   for (int i = 0; i<aridEnt.length();i++)
  45.   {
  46.     AcDbExtents tem;
  47.     if (Get_Ent_Extent(aridEnt, tem) == Acad::eOk)
  48.     {
  49.       extents.addExt(tem);//最大框
  50.     }
  51.   }
  52. }
  53. //获得实体组端点线组成的特殊角度:从小到大
  54. bool getEntsPtAngle(AcDbObjectIdArray objIdArr,CArray<double,double>& angArr)
  55. {
  56.   //点集--------------------------------------------------------
  57.   CArray<AcGePoint3d,AcGePoint3d> ptArr;
  58.   for(int i=0;i<objIdArr.logicalLength();i++)
  59.   {
  60.        AcDbEntity* pEnt=NULL;
  61.      acdbOpenAcDbEntity(pEnt,objIdArr,AcDb::kForRead);

  63.      //1 AcDbLine
  64.      if(pEnt->isKindOf(AcDbLine::desc()))
  65.      {
  66.        AcDbLine* pLine=AcDbLine::cast(pEnt);
  67.        AcGePoint3d pt1,pt2;
  68.        pt1=pLine->startPoint();
  69.        pt2=pLine->endPoint();

  71.        //数组中还没有点时
  72.            if(ptArr.GetSize()<1)
  73.        {
  74.          if(pt1.distanceTo(pt2)<0.01)
  75.          {
  76.                   ptArr.Add(pt1);
  77.          }
  78.          else
  79.          {
  80.                   ptArr.Add(pt1);
  81.                   ptArr.Add(pt2);
  82.          }
  83.        }
  84.        else
  85.        {
  86.          //pt1
  87.          for(int h=0;h<ptArr.GetSize();h++)
  88.          {
  89.                if(pt1.distanceTo(ptArr[h])<0.01)
  90.            {
  91.              break;
  92.            }

  94.            if(h==ptArr.GetSize()-1)
  95.            {
  96.                        ptArr.Add(pt1);
  97.            }
  98.          }
  99.          //pt2
  100.          if(pt1.distanceTo(pt2)>0.01)
  101.          {
  102.            for(int j=0;j<ptArr.GetSize();j++)
  103.            {
  104.              if(pt2.distanceTo(ptArr[j])<0.01)
  105.              {
  106.                break;
  107.              }
  108.             
  109.              if(j==ptArr.GetSize()-1)
  110.              {
  111.                ptArr.Add(pt2);
  112.              }
  113.            }
  114.          }
  115.          else
  116.          {
  117.            ;
  118.          }
  119.        }
  120.      }
  121.      //1 AcDbLine

  123.        //2 AcDbArc
  124.        else if(pEnt->isKindOf(AcDbArc::desc()))
  125.      {
  126.        AcDbArc* pArc=AcDbArc::cast(pEnt);
  127.        AcGePoint3d pt1,pt2;
  128.        pt1=pArc->center()+pArc->radius()*AcGeVector3d(cos( pArc->startAngle()),sin( pArc->startAngle()),0);  
  129.        pt2=pArc->center()+pArc->radius()*AcGeVector3d(cos( pArc->endAngle()),sin( pArc->endAngle()),0);
  130.   
  131.           //数组中还没有点时
  132.            if(ptArr.GetSize()<1)
  133.        {
  134.          if(pt1.distanceTo(pt2)<0.01)
  135.          {
  136.                   ptArr.Add(pt1);
  137.          }
  138.          else
  139.          {
  140.                   ptArr.Add(pt1);
  141.                   ptArr.Add(pt2);
  142.          }
  143.        }
  144.        else
  145.        {
  146.          //pt1
  147.          for(int h=0;h<ptArr.GetSize();h++)
  148.          {
  149.                if(pt1.distanceTo(ptArr[h])<0.01)
  150.            {
  151.              break;
  152.            }

  154.            if(h==ptArr.GetSize()-1)
  155.            {
  156.                        ptArr.Add(pt1);
  157.            }
  158.          }
  159.          //pt2
  160.          if(pt1.distanceTo(pt2)>0.01)
  161.          {
  162.            for(int j=0;j<ptArr.GetSize();j++)
  163.            {
  164.              if(pt2.distanceTo(ptArr[j])<0.01)
  165.              {
  166.                break;
  167.              }
  168.             
  169.              if(j==ptArr.GetSize()-1)
  170.              {
  171.                ptArr.Add(pt2);
  172.              }
  173.            }
  174.          }
  175.          else
  176.          {
  177.            ;
  178.          }
  179.        }
  180.      }
  181.      //2 AcDbArc

  183.        //3 AcDbPolyline
  184.        else if(pEnt->isKindOf(AcDbPolyline::desc()))
  185.      {
  186.        AcDbPolyline* pPLine=AcDbPolyline::cast(pEnt);
  187.        AcGePoint3d pt;
  188.            for(int i=0;i<pPLine->numVerts();i++)
  189.        {
  190.         pPLine->getPointAt(i,pt);
  191.         if(ptArr.GetSize()<1)
  192.         {
  193.              ptArr.Add(pt);
  194.         }
  195.         else
  196.         {
  197.           //
  198.                    for(int j=0;j<ptArr.GetSize();j++)
  199.            {
  200.              if(pt.distanceTo(ptArr[j])<0.01)
  201.              {
  202.                break;
  203.              }
  204.             
  205.              if(j==ptArr.GetSize()-1)
  206.              {
  207.                ptArr.Add(pt);
  208.              }
  209.            }
  210.           //
  211.         }
  212.        }
  213.            
  214.   
  215.      }
  216.      //3 AcDbPolyline
  217.      //
  218.      pEnt->close();
  219.   }
  220.     //acutPrintf("\n共需计算:%d个点!",ptArr.GetSize());
  221.   //点集--------------------------------------------------------
  222.   if(ptArr.GetSize()<2)
  223.   {
  224.     acutPrintf("\n点数太少!");
  225.     return false;
  226.   }
  227.     //
  228.     //暴力求角度集
  229.   for(i=0;i<ptArr.GetSize()-1;i++)
  230.   {
  231.         for(int j=i+1;j<ptArr.GetSize();j++)
  232.     {
  233.             double dAngle;
  234.       //始向量:+x轴，右手旋转
  235.             dAngle=AcGeVector3d(100,0,0).angleTo(AcGeVector3d(ptArr-ptArr[j]),AcGeVector3d(0,0,1));//[0,360]
  236.             angArr.Add(dAngle);
  237.             //acutPrintf("\n角度:%f",dAngle);
  238.     }
  239.    
  240.   }
  241.   
  242.   //
  243.   return true;
  244. }

论坛插件加载方法
发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；
如果你在论坛求助问题，并且已经从坛友或者管理的回复中解决了问题，请把帖子标题加上【已解决】；
如何回报帮助你解决问题的坛友，一个好办法就是给对方加【D豆】，加分不会扣除自己的积分，做一个热心并受欢迎的人！
回复 支持 反对

使用道具 举报

w15994259099

已领礼包: 271个

财富等级: 日进斗金
 楼主| 发表于 2016-12-15 15:39:31 | 显示全部楼层
【广而告之】积分规则、D豆获取攻略       窝窝空间装扮攻略
代码不太简捷，漂亮，也没有优化，请朋友们多多指教!
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；
如果你在论坛求助问题，并且已经从坛友或者管理的回复中解决了问题，请把帖子标题加上【已解决】；
如何回报帮助你解决问题的坛友，一个好办法就是给对方加【D豆】，加分不会扣除自己的积分，做一个热心并受欢迎的人！
回复 支持 反对

使用道具 举报

LoveArx

已领礼包: 13个

财富等级: 恭喜发财
发表于 2016-12-15 20:44:10 | 显示全部楼层
很好的代码，旋转那个里面

acdbOpenAcDbEntity(pEnt,objIdArr<i>,AcDb::kForWrite);
    pEnt->transformBy(xform);

建议判断下返回值，避免异常情况，其他地方也是如此，尽量都判断返回值再下一步处理。
AcadErrorStatus es = acdbOpenAcDbEntity(pEnt,objIdArr<i>,AcDb::kForWrite);
if (es == Acad::eOk || es == Acad::eWasOpenForWrite)
   pEnt->transformBy(xform);
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；
如果你在论坛求助问题，并且已经从坛友或者管理的回复中解决了问题，请把帖子标题加上【已解决】；
如何回报帮助你解决问题的坛友，一个好办法就是给对方加【D豆】，加分不会扣除自己的积分，做一个热心并受欢迎的人！
回复 支持 反对

使用道具 举报

黑蜻蜓
发表于 2019-7-7 10:09:37 | 显示全部楼层
【广而告之】积分规则、D豆获取攻略       窝窝空间装扮攻略
谢谢版主分享！
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；
如果你在论坛求助问题，并且已经从坛友或者管理的回复中解决了问题，请把帖子标题加上【已解决】；
如何回报帮助你解决问题的坛友，一个好办法就是给对方加【D豆】，加分不会扣除自己的积分，做一个热心并受欢迎的人！
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|申请友链|Archiver|手机版|小黑屋|辽公网安备|晓东CAD家园 ( 辽ICP备15016793号 )

GMT+8, 2024-11-6 03:44 , Processed in 0.236130 second(s), 39 queries , Gzip On.

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表