GIS拓扑关系编程报告

编程报告
编程名称：拓扑生成
班级：测绘学院
学号：
姓名：
一．题目：根据结点_弧表，和弧_结点表，求出多边形_弧表,和弧面表。

如图所示：
二.程序设计
1.流程图
2.拓扑生成算法
算法步骤：
（1）得到第一条弧段A，并设置为当前弧段；
（2）判断PL(A)和PR(A)是否为空。

如果都非空，转到第一步，当所有弧段处理完毕后，算法结束；
（3）如果左多边形为空，则创建一个新的多边形P，多边形的第一条弧段为当前弧段，并设置PL(A)=P，设置搜寻起始节点为N0=Ns(A)，搜寻当前节点为NC=NE(A)。

如果右多边形为空，则创建一个新的多边形P，多边形的第一条弧段为当前弧段，并设置PR(A)=P，设置搜寻起始节点N0=NE(A)，搜寻当前节点NC=NS(A)。

（4）判断N0和NC是否相等，如果是，则多边形所有弧段都已经找到，转到第一步。

（5）检查与当前节点相连接的、已经排列好的弧段序列，将当前弧段的下一条弧段A'作为多边形的第二条弧段。

（6）如果NC=NS(A ')，设置PL(A ')=P，令NC=NE(A')；如果NC= NE(A')，设置PR(A ')=P，令NC=NS(A')，转到第四步。

三．主要代码
1．文件读取
void CTuopuDlg::GetherData()//读取文件
{
UpdateData();
int i;
ifstream if1(m_sPathName1);
point=new Point[4];
for(i=0;i<4;i++)
{
if1>>point[i].p>>point[i].po1>>point[i].po2>>point[i].po3;
}
ifstream if2(m_sPathName2);
arc=new arArc[6];
for(i=0;i<6;i++)
{
if2>>arc[i].a>>arc[i].Ns>>arc[i].Ne;
arc[i].Pl=arc[i].Pr=0;
}
}
2.拓扑生成
void CTuopuDlg::OnButtonRun()
{
// TODO: Add your control notification handler code here
UpdateData();
GetherData();
int i,j,Sc,S;
CString datas1[1],datas2[1];
CString str;
char No[20];
char Nc[20];
stp=new structp[6];
int num=0;
for(S=1;S<=6;S++)
{
Sc=S;
if(arc[Sc-1].Pl!=0&&arc[Sc-1].Pr!=0) continue;
if(arc[Sc-1].Pl==0)
{
num++;
stp[S-1].arP.RemoveAll();
stp[S-1].arP.Add(Sc);
arc[Sc-1].Pl=num;
strcpy(No,arc[Sc-1].Ns);
strcpy(Nc,arc[Sc-1].Ne);
step1:if(!strcmp(No,Nc)) continue;
else
{
i=0;
while(i<4)
{
if(!strcmp(Nc,point[i].p))
{
if(Sc==point[i].po3) Sc=point[i].po1;
else if (Sc==point[i].po1) Sc=point[i].po2;
else Sc=point[i].po3;
i=4;
}
i++;
}
stp[S-1].arP.Add(Sc);
i=0;
if(!strcmp(Nc,arc[Sc-1].Ns))
{
arc[Sc-1].Pl=num;
strcpy(Nc,arc[Sc-1].Ne);
i=1;
} if(i==0&&(!strcmp(Nc,arc[Sc-1].Ne))) {
arc[Sc-1].Pr=num;
strcpy(Nc,arc[Sc-1].Ns);
}
goto step1;
}
}
else
{
if(arc[Sc-1].Pr==0)
{
num++;
stp[S-1].arP.RemoveAll();
stp[S-1].arP.Add(Sc);
arc[Sc-1].Pr=num;
strcpy(No,arc[Sc-1].Ne);
strcpy(Nc,arc[Sc-1].Ns);
step2:if(!strcmp(No,Nc))
{
continue;
}
else
{
i=0;
while(i<4)
{
if(!strcmp(Nc,point[i].p))
{
if(Sc==point[i].po3) Sc=point[i].po1;
else if (Sc==point[i].po1) Sc=point[i].po2;
else Sc=point[i].po3;
i=4;
}
i++;
}
stp[S-1].arP.Add(Sc);
i=0;
if(!strcmp(Nc,arc[Sc-1].Ns))
{
arc[Sc-1].Pl=num;
strcpy(Nc,arc[Sc-1].Ne);
i=1;
} if(i==0&&(!strcmp(Nc,arc[Sc-1].Ne)))
{
arc[Sc-1].Pr=num;
strcpy(Nc,arc[Sc-1].Ns);
}
goto step2;
}
}
}
}
stp[0].num=num;
MessageBox("运算完毕，请输出文件！");
}
3．文件输出
void CTuopuDlg::OnPutout()
{
// TODO: Add your control notification handler code here
int i,j;
CString savePath;
char szFilters[]="Data Files (*.txt)|*.txt|All Files(*.*)|*.*||";
CFileDialog dlgFileSave(FALSE,0,0, //参数给为false即为save as OFN_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;i<stp[0].num;i++)
{
int size=stp[i].arP.GetSize();
outfile<<" "<<i+1<<" "<<size;
for(j=0;j<size;j++) {outfile<<" "<<stp[i].arP.GetAt(j);}
outfile<<endl;
}
outfile<<endl;
outfile<<"----------弧_多边形拓扑--------"<<endl;
outfile<<"弧段左右"<<endl;
for(i=0;i<6;i++)
{outfile<<i+1<<" "<<arc[i].Pl<<" "<<arc[i].Pr<<endl;} }
UpdateData(false);
}
四．结果
1，界面
2.文件输出。