开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > OpenGL > 查看源码
CPoint3DCloudDoc.cpp
资源名称：CPoint3DCloud.rar　[点击查看]
上传用户：yanweivga
上传日期：2008-02-19
资源大小：1056k
文件大小：11k
源码类别：

OpenGL

开发平台：

Visual C++

CPoint3DCloudDoc.cpp：源码内容
  1. // CPoint3DCloudDoc.cpp : implementation of the CCPoint3DCloudDoc class
  2. //
  4. #include "stdafx.h"
  5. #include "CPoint3DCloud.h"
  7. #include "CPoint3DCloudDoc.h"
  8. #include "CPoint3DCloudView.h"
  11. #ifdef _DEBUG
  12. #define new DEBUG_NEW
  13. #undef THIS_FILE
  14. static char THIS_FILE[] = __FILE__;
  15. #endif
  17. /////////////////////////////////////////////////////////////////////////////
  18. // CCPoint3DCloudDoc
  20. IMPLEMENT_DYNCREATE(CCPoint3DCloudDoc, CDocument)
  22. BEGIN_MESSAGE_MAP(CCPoint3DCloudDoc, CDocument)
  23. //{{AFX_MSG_MAP(CCPoint3DCloudDoc)
  24. ON_COMMAND(ID_FILE_OPEN, OnFileOpen)
  25. //}}AFX_MSG_MAP
  26. ON_COMMAND(ID_FILE_SEND_MAIL, OnFileSendMail)
  27. ON_UPDATE_COMMAND_UI(ID_FILE_SEND_MAIL, OnUpdateFileSendMail)
  28. END_MESSAGE_MAP()
  30. /////////////////////////////////////////////////////////////////////////////
  31. // CCPoint3DCloudDoc construction/destruction
  33. CCPoint3DCloudDoc::CCPoint3DCloudDoc()
  34. {
  35. // TODO: add one-time construction code here
  37. }
  39. CCPoint3DCloudDoc::~CCPoint3DCloudDoc()
  40. {
  41. }
  43. BOOL CCPoint3DCloudDoc::OnNewDocument()
  44. {
  45. if (!CDocument::OnNewDocument())
  46. return FALSE;
  48. // TODO: add reinitialization code here
  49. // (SDI documents will reuse this document)
  51. return TRUE;
  52. }
  56. /////////////////////////////////////////////////////////////////////////////
  57. // CCPoint3DCloudDoc serialization
  59. void CCPoint3DCloudDoc::Serialize(CArchive& ar)
  60. {
  61. if (ar.IsStoring())
  62. {
  63. // TODO: add storing code here
  64. }
  65. else
  66. {
  67. // TODO: add loading code here
  68. }
  69. }
  71. /////////////////////////////////////////////////////////////////////////////
  72. // CCPoint3DCloudDoc diagnostics
  74. #ifdef _DEBUG
  75. void CCPoint3DCloudDoc::AssertValid() const
  76. {
  77. CDocument::AssertValid();
  78. }
  80. void CCPoint3DCloudDoc::Dump(CDumpContext& dc) const
  81. {
  82. CDocument::Dump(dc);
  83. }
  84. #endif //_DEBUG
  86. /////////////////////////////////////////////////////////////////////////////
  87. // CCPoint3DCloudDoc commands
  89. void CCPoint3DCloudDoc::OnFileOpen() 
  90. {
  91. // TODO: Add your command handler code here
  92. CString m_strPathname;
  93. CString m_strFilename;
  94. CString m_strFileTitle;
  95. CString m_strFileext;//文件的扩展名
  97. LPTSTR lpszFilter = "ASCII Data Files(*.asc , *ASC)|*.asc|ASCII Data Files(*.txt , *TXT)|*.txt|All Files(*.*)|*.*||";
  98. CFileDialog p_dlg(true,".txt",NULL,OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,lpszFilter,NULL);
  99. if (p_dlg.DoModal() == IDOK)
  100. {
  101. m_strPathname = p_dlg.GetPathName();
  102. m_strFilename = p_dlg.GetFileName();
  103. m_strFileTitle = p_dlg.GetFileTitle();
  104. m_strFileext = p_dlg.GetFileExt();//文件扩展名对话框
  106. }
  107. else 
  108. return ;
  109. //m_strFileext.MakeUpper();//扩展名转化成大写字母
  111. m_strFileext.MakeLower();//2005.7.21//
  112. if (m_strFileext!="txt" && m_strFileext!="asc")
  113. {
  114. return ;
  115. }
  116. FILE *p_txt= fopen(m_strPathname,"r");
  117. if (p_txt == NULL)
  118. {
  119.    fclose(p_txt);
  120.    return ;
  121. }
  122. //////////////////////////////////////////////////////////////////////////
  123. //检索文件开头
  124. //   int nSize = sizeof('n');
  125.    TCHAR tHead ;
  126.    int nBegin = 0;
  127.    tHead = fgetc(p_txt);//   tHead = fgetc(p_txt);
  128.    switch(tHead)
  129.    {
  130. case EOF://
  131.    nBegin = 0;
  132.    break;
  133. case '!'://Catia
  134.    OnFileOpenCATIA(p_txt);
  135.    return;
  136.           // nBegin = 50;
  137.    break;
  138. case '#': //Geomagic    
  139. nBegin = 30;
  140.    break;
  141. // case 'n':
  142. // nBegin = 0;
  143. // break;
  144. default:
  145. nBegin = 0;
  146.      break;
  147.    }
  148.    fseek(p_txt , nBegin , SEEK_SET);
  149. //检索文件开头
  150. //////////////////////////////////////////////////////////////////////////
  151. //检索文件结尾 是否有字符串标志
  152. fseek(p_txt ,-14L , SEEK_END);
  153. TCHAR  tchar[14]=_T("");
  154. fgets(tchar , 13 , p_txt);
  155. tchar[13] = '';
  156. CString str = tchar;
  157. int nEndCloud = str.Find("cloud");
  158. int nEnd = -1;
  159. if (nEndCloud == -1)
  160. nEnd = EOF;
  161. else if (nEndCloud != 6)
  162. {
  163.    nEndCloud = nEndCloud - 20;
  164.    nEnd = fseek(p_txt ,nEndCloud , SEEK_END);
  165. }
  166. //检索文件结尾 是否有字符串标志
  167. //////////////////////////////////////////////////////////////////////////
  168. // 判定数据点格式是三个参数 还是六个参数
  169. fseek(p_txt , nBegin , SEEK_SET);
  170. TCHAR tchar1[128] = _T("");
  171. fgets(tchar1 , 127 , p_txt);
  172.     str = tchar1;
  173. int nLine = str.Find('n');
  174. int nNull = str.Find(' ') , nPos = 0;
  175. int sum = 2;
  176. while(nNull < nLine)
  177. {
  178. nPos = str.Find(' ' , nNull + 1);
  179.         if (nPos != -1)
  180. {
  181. nNull = nPos ;
  182. ++sum;
  183. }
  184. else
  185. break;
  186. }
  187. // 判定数据点格式是三个参数 还是六个参数
  188. //////////////////////////////////////////////////////////////////////////
  189. POINT3D prMax , prMin;
  190. prMax.x = -10000000.0;
  191. prMax.y = -10000000.0;
  192. prMax.z = -10000000.0;
  193. prMin.x = 10000000.0;
  194. prMin.y = 10000000.0;
  195. prMin.z = 10000000.0;
  197. //////////////////////////////////////////////////////////////////////////
  198.     fseek(p_txt , nBegin , SEEK_SET);
  199. AfxGetApp()->BeginWaitCursor(); // vector <POINT3D> m_prToothPointsVertexList;//2005.7.28添加//定义一个字符向量
  200. m_prToothPointsVertexList.clear();// POINT3D f_xyz;
  202. POINT7D f_xyz;
  203. POINT3D f_xyzw;
  204. switch(sum) 
  205. {
  206. case 3://
  207. while (fscanf(p_txt,"%f %f %f",&f_xyzw.x,&f_xyzw.y,&f_xyzw.z) != nEnd)
  208. {
  209.            f_xyz = f_xyzw;
  210.    prMax.x=max(prMax.x,f_xyz.x);
  211.    prMax.y=max(prMax.y,f_xyz.y);
  212.    prMax.z=max(prMax.z,f_xyz.z);
  213.    prMin.x=min(prMin.x,f_xyz.x);
  214.    prMin.y=min(prMin.y,f_xyz.y);
  215.    prMin.z=min(prMin.z,f_xyz.z);    
  216.    m_prToothPointsVertexList.push_back(f_xyz);    
  217. }
  218. break;
  219. default://
  220. while (fscanf(p_txt,"%f %f %f %f %f %f",&f_xyz.x,&f_xyz.y,&f_xyz.z,&f_xyz.nx,&f_xyz.ny,&f_xyz.nz) != nEnd)
  221. {
  222.    prMax.x=max(prMax.x,f_xyz.x);
  223.    prMax.y=max(prMax.y,f_xyz.y);
  224.    prMax.z=max(prMax.z,f_xyz.z);
  225.    prMin.x=min(prMin.x,f_xyz.x);
  226.    prMin.y=min(prMin.y,f_xyz.y);
  227.    prMin.z=min(prMin.z,f_xyz.z);    
  228.    m_prToothPointsVertexList.push_back(f_xyz);    
  229. }
  230. break;
  231. }
  233. fclose(p_txt);// vector <POINT3D>::iterator pItVector1,pItVector2;//为字符数组定义游标iterator
  235. //////////////////////////////////////////////////////////////////////////
  236. vector <POINT7D>::iterator pItVector1,pItVector2;
  237. pItVector1= m_prToothPointsVertexList.end()-1;//除了最后一个元素
  238. pItVector2= pItVector1-1;
  239. if (Distance(*pItVector1,*pItVector2)<0.003)
  240. {
  241. m_prToothPointsVertexList.erase(pItVector1);
  242. return ;
  243. }
  244. //////////////////////////////////////////////////////////////////////////
  246. prBoxCenter = (prMax+prMin) / 2;
  247. prBoxSize = (prMax-prMin) / 2;
  248. CCPoint3DCloudView *pView =(CCPoint3DCloudView*) GetView(RUNTIME_CLASS(CCPoint3DCloudView));
  249. m_nDocOpenGL = &(pView->m_nViewOpenGL);// double fmax = max(prBoxSize.x , prBoxSize.y);// fmax = max(prBoxSize.z , fmax);
  250. m_nDocOpenGL->SetBox(tagCVector(prBoxCenter.x , prBoxCenter.y ,prBoxCenter.z) , tagCVector(prBoxSize.x , prBoxSize.y , prBoxSize.z));
  251. AfxGetApp()->EndWaitCursor();
  252.     UpdateAllViews(NULL, 0, 0);
  253. }
  255. void CCPoint3DCloudDoc::OnDraw()
  256. {
  257. if (m_prToothPointsVertexList.empty())
  258. return;
  259. glColor3f(1.0 , 1.0  , 1.0 );// vector<POINT3D>::iterator pItPoint3d;
  260. vector<POINT7D>::iterator pItPoint3d;
  261. glBegin(GL_POINTS);
  262. for (pItPoint3d = m_prToothPointsVertexList.begin() ; 
  263. pItPoint3d != m_prToothPointsVertexList.end() ; ++pItPoint3d)
  264. if (pItPoint3d->w > 0)
  265.       glVertex3f(pItPoint3d->x , pItPoint3d->y ,pItPoint3d->z);
  266.     glEnd();
  267. }
  269. CView* CCPoint3DCloudDoc::GetView(CRuntimeClass *pViewClass)
  270. {
  271. POSITION pos = GetFirstViewPosition();
  272. while (pos != NULL ) 
  273. {
  274. CView* pView = GetNextView(pos);
  275. if (pView->IsKindOf(pViewClass))
  276. return pView;
  277. }
  278.     return NULL;
  279. }
  281. BOOL CCPoint3DCloudDoc::OnFileOpenCATIA(FILE *f_txt)
  282. {
  283. AfxGetApp()->BeginWaitCursor(); // vector <POINT3D> m_prToothPointsVertexList;//2005.7.28添加//定义一个字符向量
  284. m_prToothPointsVertexList.clear();
  285. fseek(f_txt , 0L , SEEK_SET);
  286. //////////////////////////////////////////////////////////////////////////
  287. //catia head
  288. TCHAR tchar[128] = _T("");
  289. CString str ;
  290. int i = 0;
  291. for (i = 0 ; i < 5 ; ++i)
  292. {
  293. memset(tchar , 128 , 0);
  294. fgets(tchar , 127 , f_txt);
  295. tchar[127] = '';
  296. }
  297. /////////////////////////////////////////////////////////////////////////
  298. //获取文件数据格式
  299. str = tchar;
  300.     int nFirstChar = str.Find(''') + 1;
  301.     int nEndChar = str.ReverseFind(''');
  302. int nChar = nEndChar - nFirstChar ;
  303. CString strchar = str.Mid(nFirstChar , nChar);
  304. nChar = (nChar + 1) / 5;
  305. CString *strXYZ = new CString[nChar];
  306. CString *strFormat = new CString[nChar];
  307. for ( i = 0 ; i < nChar ; ++i)
  308. {
  309. strXYZ[i] = strchar.Mid(i * 5 , 1);
  310. strFormat[i] = strchar.Mid(i * 5 + 3 , 1);
  311. }
  312. //获取文件数据格式
  313. //////////////////////////////////////////////////////////////////////////
  314. //点个数
  315. for (i = 0 ; i < 2 ; ++i)
  316. {
  317. memset(tchar , 128 , 0);
  318. fgets(tchar , 127 , f_txt);
  319. tchar[127] = '';
  320. } 
  321. str = tchar;  
  322.     nFirstChar = str.Find('=') + 1;
  323.   nEndChar = str.GetLength() - nFirstChar;
  324.     strchar = str.Mid(nFirstChar , nEndChar);
  325. int nPointSize = atoi(strchar);
  326. m_prToothPointsVertexList.resize(nPointSize);
  327. //
  328. memset(tchar , 128 , 0);
  329. fgets(tchar , 127 , f_txt);
  330. tchar[127] = '';
  332. //catia head
  333. //////////////////////////////////////////////////////////////////////////
  334. int *nXYZ = new int[nChar + 1];
  335.     int j = 0;
  337. POINT3D prMax , prMin;
  338. prMax.x = -10000000.0;
  339. prMax.y = -10000000.0;
  340. prMax.z = -10000000.0;
  341. prMin.x = 10000000.0;
  342. prMin.y = 10000000.0;
  343. prMin.z = 10000000.0;
  345. vector<POINT7D>::iterator pItPoint7D = m_prToothPointsVertexList.begin();
  346.     for (j = 0 ; j < nPointSize ; ++j , ++pItPoint7D)
  347. {
  348. memset(tchar , 128 , 0);
  349. fgets(tchar , 127 , f_txt);
  350. tchar[127] = '';
  351. str = tchar;
  352. for ( i = 0 ; i < nChar ; ++i)
  353. nXYZ[i] = str.Find(strXYZ[i]);
  354. nXYZ[i] = str.GetLength();
  355. for (i = 0 ; i < 3 ; ++i)
  356. {
  357. strchar = str.Mid(nXYZ[i] +1 , nXYZ[i + 1] - nXYZ[i] - 1);
  358. switch(i)
  359. {
  360. case 0:   
  361. pItPoint7D->x = atof(strchar);
  362. prMax.x=max(prMax.x,pItPoint7D->x); 
  363. prMin.x=min(prMin.x,pItPoint7D->x); 
  364. break;
  365. case 1:  
  366. pItPoint7D->y = atof(strchar);
  367. prMax.y=max(prMax.y,pItPoint7D->y);
  368. prMin.y=min(prMin.y,pItPoint7D->y);
  369. break;
  370. case 2:
  371. pItPoint7D->z = atof(strchar);   
  372. prMax.z=max(prMax.z,pItPoint7D->z); 
  373. prMin.z=min(prMin.z,pItPoint7D->z);
  374. break;
  375. }
  376. }
  377. }
  378. //////////////////////////////////////////////////////////////////////////
  380. delete [] strXYZ;
  381. delete []strFormat;
  382. delete []nXYZ;
  384. strXYZ = NULL;
  385. strFormat = NULL;
  386. nXYZ = NULL;
  387.     fclose(f_txt);
  388. prBoxCenter = (prMax+prMin) / 2;
  389. prBoxSize = (prMax-prMin) / 2;
  390. CCPoint3DCloudView *pView =(CCPoint3DCloudView*) GetView(RUNTIME_CLASS(CCPoint3DCloudView));
  391. m_nDocOpenGL = &(pView->m_nViewOpenGL);// double fmax = max(prBoxSize.x , prBoxSize.y);// fmax = max(prBoxSize.z , fmax);
  392. m_nDocOpenGL->SetBox(tagCVector(prBoxCenter.x , prBoxCenter.y ,prBoxCenter.z) , tagCVector(prBoxSize.x , prBoxSize.y , prBoxSize.z));
  393. AfxGetApp()->EndWaitCursor();
  394.     UpdateAllViews(NULL, 0, 0);
  395. return TRUE;
  396. }