D8算法代码.docVIP

#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(