GIS拓扑生成程序.doc

拓扑生成程序 罗宇 使用VisualC+6.0新建一个AppWizard MFC.exe项目，目名为tuopu。选择创建一个对话框类型的应用程序。第一，项界面设计。打开对话框资源IDD_TUOPU_DIALOG,界面设计如图1：添加完控件后，打开ClassWizard,针对IDD_MY_DIALOG所指向的类CMYDlg,添加DDX成员变量；CStringm_sPathName1;CStringm_sPathName2;CStringm_editshow;系统自动生成如下代码：void CTuopuDlg:DoDataExchange(CDataExchange* pDX)CDialog:DoDataExchange(pDX);/AFX_DATA_MAP(CTuopuDlg)DDX_Text(pDX, IDC_EDIT1, m_sPathName1);DDX_Text(pDX, IDC_EDIT2, m_sPathName2);DDX_Text(pDX, IDC_EDIT3, m_editshow);/AFX_DATA_MAP第二， 代码编写。在tuopuDlg.h定义如下结构体：typedef struct arArc/弧段的结构体int a; char Ns20;char Ne20;int Pl;int Pr;arArc;typedef struct Point/点的结构体 char p20;int po1,po2,po3;Point;typedef struct structp/多边形的结构体 int num;CArray arP;structp;arArc *arc;Point *point;structp *stp;利用类向导为控件添加功能函数：系统自动生成的消息路由为：BEGIN_MESSAGE_MAP(CMyDlg, CDialog)/AFX_MSG_MAP(CMyDlg)OON_BN_CLICKED(IDC_BUTTON_OPEN1, OnButtonOpen1)ON_BN_CLICKED(IDC_BUTTON_OPEN2, OnButtonOpen2)ON_BN_CLICKED(IDC_BUTTON_RUN, OnButtonRun)ON_COMMAND(ID_EXIT, OnExit)ON_COMMAND(IDM_USE, OnUse)ON_BN_CLICKED(IDC_PUTOUT, OnPutout)/AFX_MSG_MAPEND_MESSAGE_MAP()在tuopuDlg.cpp添加各个功能函数代码如下：void CTuopuDlg:OnButtonOpen1() /打开文件1：结点-弧拓扑关系/ TODO: Add your control notification handler code hereCFileDialog fdlg(true,NULL,NULL,OFN_OVERWRITEPROMPT,拓扑关系文件|*.dat;*.txt|,NULL);if(fdlg.DoModal()=IDOK)m_sPathName1=fdlg.GetPathName(); if(m_sPathName1= ) return;UpdateData(false);void CTuopuDlg:OnButtonOpen2() /打开文件2：弧-结点拓扑关系/ TODO: Add your control notification handler code hereCFileDialog fdlg(true,NULL,NULL,OFN_OVERWRITEPROMPT,拓扑关系文件|*.dat;*.txt|,NULL);if(fdlg.DoModal()=IDOK)m_sPathName2=fdlg.GetPathName(); if(m_sPathName2= ) return; (CButton*)GetDlgItem(IDC_BUTTON_RUN)-EnableWindow(true); UpdateData(false);void CTuopuDlg:GetherData()/读取文件UpdateData();int i;ifstream if1(m_sPathName1);point=new Point4;for(i=0;ipointi.ppointi.po1pointi.po2pointi.po3;ifstream if2(m_sPathName2);arc=new arArc6;for(i=0;iarci.aarci.Nsarci.Ne;arci.Pl=arci.Pr=0;void CTuopuDlg:OnButtonRun() /拓扑关系生成 UpdateData(); GetherData();/ TODO: Add your control notification handler code hereint i,j,Sc,S;CString datas11,datas21;CString str;char No20;char Nc20;stp=new structp6;int num=0;for(S=1;S=6;S+) Sc=S;if(arcSc-1.Pl!=0&arcSc-1.Pr!=0) continue;if(arcSc-1.Pl=0) num+; stpS-1.arP.RemoveAll(); stpS-1.arP.Add(Sc); arcSc-1.Pl=num; strcpy(No,arcSc-1.Ns); strcpy(Nc,arcSc-1.Ne); step1:if(!strcmp(No,Nc) continue; else i=0; while(i4) if(!strcmp(Nc,pointi.p) if(Sc=pointi.po3) Sc=pointi.po1; else if (Sc=pointi.po1) Sc=pointi.po2; else Sc=pointi.po3; i=4; i+; stpS-1.arP.Add(Sc); i=0; if(!strcmp(Nc,arcSc-1.Ns) arcSc-1.Pl=num; strcpy(Nc,arcSc-1.Ne); i=1; if(i=0&(!strcmp(Nc,arcSc-1.Ne) arcSc-1.Pr=num; strcpy(Nc,arcSc-1.Ns); goto step1; elseif(arcSc-1.Pr=0)num+; stpS-1.arP.RemoveAll(); stpS-1.arP.Add(Sc); arcSc-1.Pr=num; strcpy(No,arcSc-1.Ne); strcpy(Nc,arcSc-1.Ns); step2:if(!strcmp(No,Nc) continue; else i=0; while(i4) if(!strcmp(Nc,pointi.p) if(Sc=pointi.po3) Sc=pointi.po1; else if (Sc=pointi.po1) Sc=pointi.po2; else Sc=pointi.po3; i=4; i+; stpS-1.arP.Add(Sc); i=0; if(!strcmp(Nc,arcSc-1.Ns) arcSc-1.Pl=num; strcpy(Nc,arcSc-1.Ne); i=1; if(i=0&(!strcmp(Nc,arcSc-1.Ne) arcSc-1.Pr=num; strcpy(Nc,arcSc-1.Ns); goto step2; stp0.num=num;datas10=-面_弧拓扑-rn;datas10+=面 弧数 弧号 rn;for(i=0;inum;i+)int size=stpi.arP.GetSize();str.Format(%d %d ,i+1,size);datas10+=str;for(j=0;jsize;j+)str.Format(%d ,stpi.arP.GetAt(j);datas10+=str;datas10+=rn;datas10+=rn;m_editshow=datas10;datas20=-弧_多边形拓扑-rn;datas20+=弧段 左 右 rn;for(i=0;iEnableWindow(true);UpdateData(false);void CTuopuDlg:OnPutout() /将结果输出txt文件/ TODO: Add your control notification handler code hereint i,j; CString savePath;char szFilters=Data Files (*.txt)|*.txt|All Files(*.*)|*.*|;CFileDialog dlgFileSave(FALSE,0,0, /参数给为false即为save asOFN_FILEMUSTEXIST | OFN_PATHMUSTEXIST, (CString)szFilters, this);if(IDOK=dlgFileSave.DoModal()savePath = dlgFileSave.GetPathName()+.txt;/路径加文件名再加后缀UpdateData(FALSE);fstream outfile;outfile.open(savePath,ios:out|ios:app); outfile-面_弧拓扑-endl; outfile面 弧数 弧号 endl;for(i=0;istp0.num;i+)int size=stpi.arP.GetSize();outfile i+1 size;for(j=0;jsize;j+) outfile stpi.arP.GetAt(j);outfil