Advertisement

L3-004. 肿瘤诊断

阅读量:

在诊断肿瘤疾病时,计算肿瘤体积是很重要的一环。给定病灶扫描切片中标注出的疑似肿瘤区域,请你计算肿瘤的体积。

输入格式:

输入第一行给出4个正整数:M、N、L、T,其中M和N是每张切片的尺寸(即每张切片是一个M×N的像素矩阵。最大分辨率是1286×128);L(<=60)是切片的张数;T是一个整数阈值(若疑似肿瘤的连通体体积小于T,则该小块忽略不计)。

最后给出L张切片。每张用一个由0和1组成的M×N的矩阵表示,其中1表示疑似肿瘤的像素,0表示正常像素。由于切片厚度可以认为是一个常数,于是我们只要数连通体中1的个数就可以得到体积了。麻烦的是,可能存在多个肿瘤,这时我们只统计那些体积不小于T的。两个像素被认为是“连通的”,如果它们有一个共同的切面,如下图所示,所有6个红色的像素都与蓝色的像素连通。

Figure 1

输出格式:

在一行中输出肿瘤的总体积。

输入样例:

复制代码

输出样例:

26

复制代码
 #include<queue>

    
 #include<iostream>
    
 #include<string.h>
    
 using namespace std;
    
 int d[6][3]={{1,0,0},{-1,0,0},{0,1,0},{0,-1,0},{0,0,1},{0,0,-1}};  //上下左右前后六个方向的搜索
    
 int map[65][130][1300];   //三维数组储存数据
    
 int l,n,m,t;
    
 int ans=0;    //计算肿瘤体积
    
 struct node
    
 {
    
     int x,y,z;   //定义三维结构体
    
 };
    
  
    
 int judge(int x,int y,int z)   //判断坐标是否合格
    
 { 
    
 	if((x<0)||(y<0)||(z<0)){  
    
 		return 0;  
    
 	}  
    
 	if((x>=l)||(y>=n)||(z>=m)){  
    
 		return 0;  
    
 	}  
    
 	return 1;  
    
 }
    
 void BFS(int x,int y,int z)   //搜索
    
 {
    
     int i;
    
     int sum=1;
    
 	
    
     queue<node>q;
    
     map[x][y][z]=0;       //把这个坐标设为0,下次不再遍历
    
     struct node u;
    
 	
    
     u.x=x;u.y=y;u.z=z;
    
     q.push(u);        //入队列
    
 	struct node v;
    
     while(!q.empty())       //同最大黑区域,队列不为空时一直弹出判断
    
     {
    
     u=q.front();
    
     q.pop();
    
     
    
     for(i=0;i<6;i++)   //六个方向搜索
    
     {
    
         v.x=u.x+d[i][0];
    
 			v.y=u.y+d[i][1];
    
         v.z=u.z+d[i][2];
    
         if(judge(v.x,v.y,v.z)&&map[v.x][v.y][v.z])  //当要搜索的方向坐标符合题意,并且为1时入队列
    
         {
    
 				map[v.x][v.y][v.z]=0;
    
             sum++;
    
             q.push(v);
    
         }
    
     }
    
     }
    
 	if(sum>=t){    //当这次的肿瘤体积大于等于要求的t,答案++
    
 		ans+=sum;
    
 	}
    
 }
    
 int main()
    
 {
    
     scanf("%d%d%d%d",&n,&m,&l,&t);
    
 	int i,j,k;
    
 	for(i=0;i<l;i++){
    
 		for(j=0;j<n;j++){
    
 			for(k=0;k<m;k++)
    
 			{
    
 				scanf("%d",&map[i][j][k]);     //输入数据
    
 			}
    
 		}
    
 	}
    
 	for(i=0;i<l;i++){
    
 		for(j=0;j<n;j++){
    
 			for(k=0;k<m;k++){
    
 				if(map[i][j][k]==1){   //枚举遍历,当是1的时候搜索
    
 					BFS(i,j,k);
    
 				}
    
 			}
    
 		}
    
 	}
    
 	cout<<ans<<endl; 
    
     return 0;
    
 }
    
 双击代码区域可以查看未格式化的原始代码
    
  
    
  
    
 查看提交
    
    
    
    

全部评论 (0)

还没有任何评论哟~