D8算法代码.doc

下载文档 降价啦

145
0
约3.21千字
约 6页
2018-03-13 发布于河南
举报
版权申诉
保障服务

D8算法代码.doc

1、本文档共6页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

D8算法代码

#includeiostream #include fstream #include sstream #include string #includectime #includemath.h using namespace std; int main() { int src[100][100]={0},Vector[100][100]={0},Result[100][100]={0}; int row=100,col=100; coutendl; system(pause); /* for(int r=0;rrow;r++) { for(int c=0;ccol;c++) coutsrc[r][c] ; coutendl; }*/ double Sqrt2; float S=0,N=0,E=0,SE=0,NE=0,NW=0,W=0,SW=0; Sqrt2=sqrt(2); time_t begin,end; begin=clock(); //your code for(int i=0;irow;i++) { for(int j=0;jcol;j++) { S= (i!=(row-1)) ? (src[i][j]-src[i+1][j]):-1; SE= ( i!=(row-1) j!=(col-1) ) ? (src[i][j]-src[i+1][j+1])/sqrt(2):-1; N= (i!=0) ? (src[i][j]-src[i-1][j]):-1; E=(j!=(col-1) ) ? (src[i][j]-src[i][j+1]):-1; NE=( i!=0 j!=(col-1) ) ? (src[i][j]-src[i-1][j+1])/sqrt(2):-1; NW=(i!=0 j!=0) ? (src[i][j]-src[i-1][j-1])/sqrt(2):-1; W= (j!=0) ? (src[i][j]-src[i][j-1]):-1; SW=( i!=(row-1) j!=0 ) ? (src[i][j]-src[i+1][j-1])/sqrt(2):-1; //坡降计算 float M=0; M=(MS)?M:S; M=(MSE)?M:SE; M=(MN)?M:N; M=(ME)?M:E; M=(MNE)?M:NE; M=(MNW)?M:NW; M=(MW)?M:W; M=(MSW)?M:SW; //取最大的坡降 if(M0) { if(M==S) { Vector[i][j]=4; } else if(M==SE) { Vector[i][j]=2; } else if(M==N) { Vector[i][j]=64; } else if(M==E) { Vector[i][j]=1; } else if(M==NE) { Vector[i][j]=128; } else if(M==NW) { Vector[i][j]=32; } else if(M==W) { Vector[i][j]=16; } else if(M==SW) { Vector[i][j]=8; } else { Vector[i][j]=0; } } } } int i1=0,j1=0; /*for(i=0;irow;i++) { for(int j=0;jcol;j++) coutVector[i][j] ; coutendl; }*/ for( i=0;irow;i++) { for(int j=0;jcol;j++) { i1=i; j1=j; while(1)//方向为0时退出 { int x=Vector[i1][j1];//取出方向值 if(Vector[i1][j1]==0 || i1=row ||j1=col) break; switch(x) { // temp=0; Vector[i1][j1]|=0x100000;//设置该点已被遍历 case 1: { Result[i1][j1+1]++; j1=j1+1; break; } case 2: { // if(