方格網程序及算法

1、西 南 交 通 大 學 線 路 工 程 信 息 技 術 編 程 作 業(yè)部分代碼：/Doc類，主要處理文件讀寫，以及運算； 25 void CDTMDoc:OnFileOpen() /響應函數，響應“打開”這一命令消息CMainFrame *pMainFrm = new CMainFrame;pMainFrm = (CMainFrame *)AfxGetApp()-m_pMainWnd;pMainFrm-GetMenu()-GetSubMenu(0)-EnableMenuItem(ID_FILE_OPEN,MF_BYPOSITION | MF_GRAYED);/以下利用CFileDialog類創(chuàng)

2、建一個打開文件對話框CFileDialog OpenFileDlg(TRUE, _T(DEM Files (*DEM), NULL, OFN_HIDEREADONLY | OFN_OVERWRITEPROMPT |OFN_ALLOWMULTISELECT|OFN_ENABLESIZING,_T(Digital Terrain Data Files (*.txt;*.DTD;*.dat)|*.txt;*.DTD;*.dat)|),NULL);CString path;/DoModal()函數成功則返回IDOK或IDCANCEL,其返回值對應于用來關閉對話框的控件的IDif(OpenFileDlg

3、.DoModal() = IDOK )/點擊對話框確定按鈕得到路徑v_OrigPoint.clear();XYLine.clear();poOnLine.clear();path = OpenFileDlg.GetPathName();/得到打開文件的路徑elseAfxMessageBox(_T(Cancel!);return;ifstream file;setlocale(LC_CTYPE, .936);/處理中文路徑問題file.open(path);/setlocale(LC_CTYPE, );/成對使if (!file)AfxMessageBox(_T(Reading file fai

4、led! Maybe the file does not exist!);return;/Open file failed,return.OrigPoint *temp_OrigPoint = new OrigPoint;AfxMessageBox(_T(讀取DTM數據時間可能較長，n請耐心等待!);OnStatusBar(_T(Reading.);char c;while(TRUE)/ /利用Break跳出filetemp_OrigPoint-Mark1; /第一列12,22for (int i = 0;ic;filetemp_OrigPoint-x;filetemp_OrigPoint-y

5、;filetemp_OrigPoint-z;filetemp_OrigPoint-Mark2; /最后一列if (file)v_OrigPoint.push_back(*temp_OrigPoint);elsefile.close(); break;m_PointNUM=v_OrigPoint.size();if (m_PointNUM = 0)/不支持文件格式AfxMessageBox(_T(不支持文件格式！n打開失敗！);return;FileIsOpen = TRUE;Analysis();pMainFrm-m_wndStatusBar.SetPaneInfo(1,ID_State_Do

6、ing,SBPS_NORMAL,150);/100為寬度可調整pMainFrm-m_wndStatusBar.SetPaneText(1,_T(Waitting Location!),true);/ TODO: Add your command handler code herevoid CDTMDoc:OnStatusBar(CString str)/狀態(tài)欄進度CMainFrame *pMainFrm = new CMainFrame;pMainFrm = (CMainFrame *)AfxGetApp()-m_pMainWnd;CStatusBar* StatusBar = (CStatu

7、sBar*)pMainFrm-GetStatusBar();int index =StatusBar-CommandToIndex(ID_State_Doing);StatusBar-SetPaneInfo(index,ID_State_Doing,SBPS_NORMAL,100);/100為寬度可調整StatusBar-SetPaneText(index,str,true);BOOL CDTMDoc:Analysis()/analysis the regionm_left = v_OrigPoint.at(0).x;m_right = v_OrigPoint.at(0).x;m_bottom

8、 = v_OrigPoint.at(0).y;m_top = v_OrigPoint.at(0).y;m_lowest = v_OrigPoint.at(0).z;m_highest = v_OrigPoint.at(0).z;for (int i = 0; i v_OrigPoint.at(i).x)m_left = v_OrigPoint.at(i).x;if (m_right v_OrigPoint.at(i).y)m_bottom = v_OrigPoint.at(i).y;if (m_top v_OrigPoint.at(i).z)m_lowest = v_OrigPoint.at(

9、i).z;if (m_highest m_pMainWnd;CMDIChildWnd *pChild = (CMDIChildWnd *) pFrame-GetActiveFrame();CDTMView *pView = (CDTMView *) pChild-GetActiveView();pView-Invalidate();/更新viewreturn TRUE;bool comparex(const OrigPoint& lhs, const OrigPoint& rhs)if (lhs.x = rhs.x)return lhs.y rhs.y;elsereturn lhs.x rhs

10、.x;bool comparey(const OrigPoint& lhs, const OrigPoint&rhs)if (lhs.y = rhs.y)return lhs.x rhs.x;elsereturn lhs.y rhs.y;void CDTMDoc:Mesh()M = (m_right- m_left)/DTL+1;/設置網格N = (m_top - m_bottom)/DTL+1;if (M*N = 0)AfxMessageBox(_T(先讀取數據);return;PointinPane *PoinPa;/二維數組PoinPa= new PointinPane*M;for(in

11、t i=0; iM; +i)PoinPai = new PointinPaneN;m_PointNUM = v_OrigPoint.size();int i = 0;int j = 0;for (int k=0;km_PointNUM;k+)i = int(v_OrigPoint.at(k).x-m_left)/DTL);j = int(v_OrigPoint.at(k).y-m_bottom)/DTL);PoinPaij.point.push_back(v_OrigPoint.at(k);Node= new Point3D*M+1;for(int i=0; iM+1; +i)Nodei =

12、new Point3DN+1;for (int i=0;iM+1;i+)for (int j=0;jN+1;j+)Nodeij.x = i*200+m_left;Nodeij.y = j*200+m_bottom;Nodeij.z = 0;for (int i=1;iM;i+)for (int j=1;jN;j+)double x = Nodeij.x;double y = Nodeij.y;double l1,l2,l3,l4;double p1,p2,p3,p4;double z1,z2,z3,z4;PointinPane V;V = PoinPai-1j-1;int pid1 = Fin

13、dNearest(V,x,y);if (pid1 = -1)z1 = 0;p1 = 0;z1 = 0;elsel1 = (V.point.at(pid1).x-x)*(V.point.at(pid1).x-x)+(V.point.at(pid1).y-y)*(V.point.at(pid1).y-y);p1 = 1/l1;z1 = V.point.at(pid1).z;V = PoinPai-1j;int pid2 = FindNearest(V,x,y);if (pid2 = -1)z2 = 0;p2 = 0;z2 = 0;elsel2 = (V.point.at(pid2).x-x)*(V

14、.point.at(pid2).x-x)+(V.point.at(pid2).y-y)*(V.point.at(pid2).y-y);p2 = 1/l2;z2 = V.point.at(pid2).z;V = PoinPaij-1;int pid3 = FindNearest(V,x,y);if (pid3 =-1)z3 = 0;p3 = 0;z3 = 0;elsel3 = (V.point.at(pid3).x-x)*(V.point.at(pid3).x-x)+(V.point.at(pid3).y-y)*(V.point.at(pid3).y-y);p3 = 1/l3;z3 = V.po

15、int.at(pid3).z;V = PoinPaij;int pid4 = FindNearest(V,x,y);if (pid4 = -1)z4 = 0;p4 = 0;z4 = 0;elsel4 = (V.point.at(pid4).x-x)*(V.point.at(pid4).x-x)+(V.point.at(pid4).y-y)*(V.point.at(pid4).y-y);p4 = 1/l4;z4 = V.point.at(pid4).z;if (p1+p2+p3+p4 =0)/相鄰的四個網格中沒有找到點Nodeij.z = -200;elsedouble z = (z1*p1+z

16、2*p2+z3*p3+z4*p4)/(p1+p2+p3+p4);Nodeij.z = z;int CDTMDoc:FindNearest(PointinPane V,double x,double y)double mm=200*200*2;double l=0;int pid=-1;if (V.point.size()=0)return -1;for (int i=0;iV.point.size();i+)l=(V.point.at(i).x-x)*(V.point.at(i).x-x)+(V.point.at(i).y-y)*(V.point.at(i).y-y);if (lm_pMain

17、Wnd;CMDIChildWnd *pChild = (CMDIChildWnd *) pFrame-GetActiveFrame();CDTMView *pView = (CDTMView *) pChild-GetActiveView();/pView-m_locationState = FALSE;/關閉定線功能poOnLine.clear();if (XYLine.size()1)double SBase = 0;for (int i=1;iXYLine.size();i+)double xx = XYLine.at(i).x-XYLine.at(i-1).x;double yy =

18、XYLine.at(i).y-XYLine.at(i-1).y;double x0 = XYLine.at(i-1).x;double y0 = XYLine.at(i-1).y;double ss =sqrt(xx*xx+yy*yy);PointOnLine point;for (int j=0;j*divss;j+)point.x = x0 + (j*div/ss)*xx;point.y = y0 + (j*div/ss)*yy;point.s =SBase+j*div;poOnLine.push_back(point);point.x = XYLine.at(i).x;point.y =

19、 XYLine.at(i).y;point.s = SBase+ss;poOnLine.push_back(point);SBase +=ss;elsereturn;/內插int J;int I;lHeighest = m_lowest-5000;lLowest = m_highest+5000;for (int i = 0;ipoOnLine.size();i+)J = (poOnLine.at(i).x)/DTL;I = (poOnLine.at(i).y)/DTL;if (J1|IM-1|IN-1)poOnLine.at(i).z = -200;continue;double x = p

20、oOnLine.at(i).x-DTL*J;double y = poOnLine.at(i).y-DTL*I;double z =NodeJI.z*(1-x/DTL)*(1-y/DTL)+NodeJ+1I.z*(x/DTL)*(1-y/DTL)+NodeJ+1I+1.z*(x/DTL)*(y/DTL)+NodeJI+1.z*(1-x/DTL)*(y/DTL);poOnLine.at(i).z = z;if (zlHeighest)lHeighest = z;if (zInvalidate();/ TODO: Add your command handler code herevoid CDT

21、MDoc:OnUpdateProfileLongitudinal(CCmdUI *pCmdUI)if (XYLine.size()1)pCmdUI-Enable(TRUE);elsepCmdUI-Enable(FALSE);/ TODO: Add your command update UI handler code here/View類主要負責圖像顯示、響應鼠標消息等void CDTMView:OnDraw(CDC* pDC)CDTMDoc* pDoc = GetDocument();ASSERT_VALID(pDoc);if (!pDoc)return;OldPdc = pDC;GetCl

22、ientRect(clientRect);/得到窗口大小CBrush brush(RGB(0, 0, 0);pDC-SelectObject(brush);windowx =pDoc-m_right-pDoc-m_left;/邏輯坐標寬度范圍windowy =pDoc-m_top-pDoc-m_bottom;if (scal FillSolidRect(trueRect,RGB(0,0,0);if (windowx = 0|windowy = 0)windowx = clientRect.right;windowy = clientRect.bottom;elsesizeTotal.cx =

23、clientRect.Width()*scal;sizeTotal.cy = clientRect.Height()*scal; SetScrollSizes(MM_TEXT, sizeTotal); /調節(jié)滾動條trueScale = max(windowx/(scal*clientRect.Width(),windowy/(scal*clientRect.Height();LOGFONT lf;memset( &lf, 0, sizeof( LOGFONT ) );pDC-GetCurrentFont()-GetLogFont(&lf);/獲取系統當前字體CFont font1;lf.lf

24、CharSet=GB2312_CHARSET;/定義字符lf.lfHeight = MyFontHeight*trueScale*0.8;lf.lfWidth=0;lstrcpy(lf.lfFaceName,Myscript);font1.CreateFontIndirect(&lf);pDC-SelectObject(&font1);pDC-SetBkMode(TRANSPARENT);/設置字體輸出背景為透明pDC-SetTextColor(RGB(0,255,255); /設定文本顏色/設置映射關系CPoint pos=GetScrollPosition();pDC-SetMapMode

25、(MM_ISOTROPIC);pDC-SetWindowExt(windowx,windowy);pDC-SetViewportExt(trueRect.Width()-40,-trueRect.Height()+40);pDC-SetViewportOrg(0-pos.x+20,trueRect.Height()-pos.y-20);CString strelev;/用于高程顯示字符串；offx = int(pDoc-m_left);offy = int(pDoc-m_bottom)+1;for (int i = 0; pDoc-m_PointNUM!=0&im_PointNUM;i+)in

26、t pennum = int(pDoc-v_OrigPoint.at(i).z-pDoc-m_lowest+0.5*pDoc-hh)/(pDoc-hh);/計算畫筆顏色pDC-SelectObject(&colorPenpennum);if (pDoc-v_OrigPoint.at(i).Mark2 = 31)pDC-MoveTo(pDoc-v_OrigPoint.at(i).x-offx,pDoc-v_OrigPoint.at(i).y-offy);strelev.Format(_T(%.2lf),pDoc-v_OrigPoint.at(i).z);if (pDoc-v_OrigPoint.

27、at(i).Mark2 = 32)if (i1&(pDoc-v_OrigPoint.at(i-1).Mark2 =33)pDC-MoveTo(pDoc-v_OrigPoint.at(i).x-offx,pDoc-v_OrigPoint.at(i).y-offy);elsepDC-LineTo(pDoc-v_OrigPoint.at(i).x-offx,pDoc-v_OrigPoint.at(i).y-offy);if (IsShowHeight = TRUE&pDoc-v_OrigPoint.at(i).Mark1 = 12&i%20 = 0)pDC-TextOut(pDoc-v_OrigPo

28、int.at(i).x-offx,pDoc-v_OrigPoint.at(i).y-offy,strelev);if (pDoc-v_OrigPoint.at(i).Mark2 = 33)/pDC-LineTo(pDoc-v_OrigPoint.at(i).x-offx,pDoc-v_OrigPoint.at(i).y-offy); if (IsGrid = TRUE)/是否繪制網格DrawGrid(pDC);CPen redPen;redPen.CreatePen(PS_SOLID,5,RGB(255,0,0);pDC-SelectObject(redPen);CString str;for

29、 (int i=0;iXYLine.size();i+)if (i = 0)pDC-MoveTo(pdocument-XYLine.at(0);pDC-SetTextColor(RGB(255,255,0);pDC-TextOut(pdocument-XYLine.at(0).x,pdocument-XYLine.at(0).y,_T(起點);else if(i=2)str.Format(_T(JD %d),i-1);pDC-LineTo(pdocument-XYLine.at(i);pDC-TextOut(pdocument-XYLine.at(i-1).x,pdocument-XYLine

30、.at(i-1).y,str);elsepDC-LineTo(pdocument-XYLine.at(i);int Index = pdocument-XYLine.size()-1;if (Index0&m_locationState = TRUE)int nOldMode=pDC-GetROP2();pDC-SetROP2(R2_NOT);pDC-MoveTo(pdocument-XYLine.at(Index).x,pdocument-XYLine.at(Index).y);pDC-LineTo(OldPoint.x,OldPoint.y);pDC-SetROP2(nOldMode);o

31、ldscal = scal;if (m_locationState=FALSE&Index0)pDC-TextOut(pdocument-XYLine.at(Index).x,pdocument-XYLine.at(Index).y,_T(終點);/ TODO: add draw code for native data here/ CDTMView printingBOOL CDTMView:OnPreparePrinting(CPrintInfo* pInfo)/ default preparationreturn DoPreparePrinting(pInfo);void CDTMVie

32、w:OnBeginPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)/ TODO: add extra initialization before printingvoid CDTMView:OnEndPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)/ TODO: add cleanup after printing/ CDTMView diagnostics/ CDTMView message handlersBOOL CDTMView:OnMouseWheel(UINT nFlags, short zDelt

33、a, CPoint pt)/ TODO: Add your message handler code here and/or call defaultint Index = pdocument-XYLine.size()-1;if (Index=0)CDC *pDC=GetDC();CPoint pos=GetScrollPosition();pDC-SetMapMode(MM_ISOTROPIC);pDC-SetWindowExt(windowx,windowy);pDC-SetViewportExt(trueRect.Width()-40,-trueRect.Height()+40);pD

34、C-SetViewportOrg(0-pos.x+20,trueRect.Height()-pos.y-20);int nOldMode=pDC-GetROP2();pDC-SetROP2(R2_NOT);pDC-MoveTo(pdocument-XYLine.at(Index).x,pdocument-XYLine.at(Index).y);pDC-LineTo(OldPoint.x,OldPoint.y);pDC-MoveTo(pdocument-XYLine.at(Index).x,pdocument-XYLine.at(Index).y);CPoint point1 = pt;pDC-

35、DPtoLP(&point1);pDC-LineTo(point1.x,point1.y);OldPoint = point1;pDC-SetROP2(nOldMode);oldscal = scal;scal += zDelta/500.00;Invalidate();/return CScrollView:OnMouseWheel(nFlags, zDelta, pt);return TRUE;BOOL CDTMView:OnEraseBkgnd(CDC* pDC)/ TODO: Add your message handler code here and/or call defaultr

36、eturn TRUE;/return CScrollView:OnEraseBkgnd(pDC);void CDTMView:OnVScroll(UINT nSBCode, UINT nPos, CScrollBar* pScrollBar)/ TODO: Add your message handler code here and/or call defaultInvalidate();CScrollView:OnVScroll(nSBCode, nPos, pScrollBar);void CDTMView:OnSize(UINT nType, int cx, int cy)Invalid

37、ate();CScrollView:OnSize(nType, cx, cy);/ TODO: Add your message handler code herevoid CDTMView:OnHScroll(UINT nSBCode, UINT nPos, CScrollBar* pScrollBar)/ TODO: Add your message handler code here and/or call defaultInvalidate();CScrollView:OnHScroll(nSBCode, nPos, pScrollBar);void CDTMView:DrawGrid

38、(CDC *pDC)pDC-SetBkMode(TRANSPARENT);/設置字體輸出背景為透明pDC-SetTextColor(RGB(0,255,255); /設定文本顏色LOGFONT lf;memset( &lf, 0, sizeof( LOGFONT ) );pDC-GetCurrentFont()-GetLogFont(&lf);/獲取系統當前字體CFont font1;lf.lfCharSet=GB2312_CHARSET;/定義字符lf.lfHeight = MyFontHeight*trueScale;lf.lfWidth=0;lstrcpy(lf.lfFaceName,M

39、yscript);font1.CreateFontIndirect(&lf);pDC-SelectObject(&font1);CDTMDoc* pDoc = GetDocument();CPen pen;pen.CreatePen(PS_DOT,0,RGB(0,60,60);pDC-SelectObject(pen);int m = pDoc-M;int n = pDoc-N;double dtl = pDoc-DTL;double x0 = 0;double y0 = 0;double xSpan = pDoc-m_right-pDoc-m_left;double ySpan = pDoc

40、-m_top-pDoc-m_bottom;for (int i=0;iMoveTo(x0,y0+i*dtl);pDC-LineTo(x0+dtl*m,y0+i*dtl);for (int i=0;iMoveTo(x0+i*dtl,y0);pDC-LineTo(x0+i*dtl,y0+dtl*n);for (int i=1;im;i+ )for (int j=1;jNodeij.z !=-200)CString str;str.Format(_T(%.2lf),pDoc-Nodeij.z);pDC-TextOut(xx,yy,str);void CDTMView:OnLinelocation()

41、CMainFrame *pMainFrm = (CMainFrame *)AfxGetApp()-m_pMainWnd;pMainFrm-m_wndStatusBar.SetPaneInfo(1,ID_State_Doing,SBPS_NORMAL,80);/100為寬度可調整CString str;if (m_locationState = FALSE)m_locationState = TRUE;str.Format(_T(Locating.);pMainFrm-m_wndStatusBar.SetPaneText(1,str,true);elsem_locationState = FAL

42、SE;str.Format(_T(Waitting Location!);pMainFrm-m_wndStatusBar.SetPaneText(1,str,true);int Index = pdocument-XYLine.size()-1;if (Index=0)OldPoint = pdocument-XYLine.at(Index);Invalidate();/ TODO: Add your command handler code herevoid CDTMView:OnMButtonDown(UINT nFlags, CPoint point)/ TODO: Add your m

43、essage handler code here and/or call defaultCScrollView:OnMButtonDown(nFlags, point);void CDTMView:OnLButtonDown(UINT nFlags, CPoint point)/ TODO: Add your message handler code here and/or call defaultif (m_locationState = TRUE)CDC *pDC=GetDC();CPoint pointtemp = point;CPoint pos=GetScrollPosition()

44、;pDC-SetMapMode(MM_ISOTROPIC);pDC-SetWindowExt(windowx,windowy);pDC-SetViewportExt(trueRect.Width()-40,-trueRect.Height()+40);pDC-SetViewportOrg(0-pos.x+20,trueRect.Height()-pos.y-20);CPen redPen;redPen.CreatePen(PS_SOLID,5,RGB(255,0,0);pDC-SelectObject(redPen);int m =pdocument-M;int n =pdocument-N;

45、double dtl =pdocument-DTL;pDC-DPtoLP(&pointtemp);if ( pointtemp.x(m-1)*dtl|pointtemp.x(n-1)*dtl|pointtemp.yXYLine.push_back(pointtemp);OldPoint = pointtemp;for (int i=0;iXYLine.size();i+)LOGFONT lf;memset( &lf, 0, sizeof( LOGFONT ) );pDC-GetCurrentFont()-GetLogFont(&lf);/獲取系統當前字體CFont font1;lf.lfCha

46、rSet=GB2312_CHARSET;/定義字符lf.lfHeight = MyFontHeight*trueScale*0.8;lf.lfWidth=0;lstrcpy(lf.lfFaceName,Myscript);font1.CreateFontIndirect(&lf);pDC-SetBkMode(TRANSPARENT);/設置字體輸出背景為透明pDC-SetTextColor(RGB(255,255,0); /設定文本顏色pDC-SelectObject(&font1);CString str;if (i = 0)pDC-MoveTo(pdocument-XYLine.at(0);pDC-TextOut(pdocument-XYLine.at(0).x,pdocument-XYLine.at(0).y,_T(起點);else if(i=2)str.Format(_T(JD %d),i-1);pDC-LineTo(pdocument-X