等值线实验报告.docVIP

一、实验内容 采用网格序列法实现等值线生成。 基本要求： 网格分辨率100x100 网格结点数据为随机生成的[0,1]数据 等值为0.2、0.8 将等值线用分别用不同颜色标示 二、主要技术及实现 1. 画网格 采用for循环，分别画出网格中的横线和竖线。 2. 给网格顶点赋值 使用一个二维数组，对应到网格各顶点，使用 (rand()%11)/10.00为网格顶点赋值。 3. 生成等值线 逐个网格处理，将网格顶点值分为一正三负，三正一负，两正两负三种情况，其中两正两负时需要考虑同侧不同号时的二义性处理。 等值点坐标 利用线性插值法求出。 二义性处理 依据线性插值法判断中心点的正负，画出的等值线应保证中心点与同区域的顶点正负值相同。 三、实验结果 四、实验分析和总结 通过这次实验，学会了如何通过VC++进行画线，在实现等值线生成算法的过程中，加深了对算法的理解。因为在写代码的时候是一边调试一边写，所以很容易判断每一种情况处理的是否正确，为调试工作省去了不少麻烦。但是这样也反映出了边想边做的弊端，就是不能对代码有一个全局的掌控，所以以后会尽量先理清思路，对算法有一个整体上的设计之后再动手。 五、核心源代码 //代码要有详细的注释 //以等值线0.2为例 …… #define TABSIZE 100 //网格数量 #define CELLLEN 10 //网格长度 #define TABSTART 0 //表格起始坐标 #define EV1 0.2 //等值线的值 #define EV2 0.8 //等值线的值 …… // CMFCApplication1View 绘制 void CMFCApplication1View::OnDraw(CDC* /*pDC*/) { CMFCApplication1Doc* pDoc = GetDocument(); ASSERT_VALID(pDoc); if (!pDoc) return; // TODO: 在此处为本机数据添加绘制代码 CClientDC dc(this); CPen pen(PS_SOLID,0.5,RGB(0,0,0)); CPen* oldpen=dc.SelectObject(pen); //dc.MoveTo(20,20); // dc.LineTo(20,100); int i,j; double x0,x1,x2; //画网格 int m=TABSTART,n=TABSTART+(TABSIZE-1)*CELLLEN; for( i=m;i=n;i=i+CELLLEN) { dc.MoveTo(i,m); dc.LineTo(i,n); dc.MoveTo(m,i); dc.LineTo(n,i); } double point[TABSIZE][TABSIZE]; //网格顶点的值 int mark[TABSIZE][TABSIZE]; //网格顶点的正负值 int pointx[TABSIZE][TABSIZE]; //网格中横线上的等值点 int pointy[TABSIZE][TABSIZE]; //网格中竖线上的等值点 //为网格顶点赋值 for( i=0;iTABSIZE;i++) { for(j=0;jTABSIZE;j++) { point[i][j]=(rand()%11)/10.00; //TRACE(...%d, point[i][j], ); } } //判断网格顶点的正负值，正为1，负为0 for( i=0;iTABSIZE;i++) { for(j=0;jTABSIZE;j++) { if(point[i][j]EV1) mark[i][j]=0; else mark[i][j]=1; } } CClientDC dc1(this); CClientDC dc2(this); CPen pen1(PS_SOLID,0.5,RGB(255,0,0)); //等值线0.2使用红色 CPen* oldpen1=dc1.SelectObject(pen1); int a,b; double c,d; //计算网格线上的等值点坐标 for(a=0;aTABSIZE;a++) { for(b=0;bTABSIZE;b++) { if(mark[a][b]+mark[a][b+1]==1) { c=TABSTART+b*CELLLEN; d=TABSTART+(b+1)*CELLLEN; pointx