rename("orig");
setBatchMode(true);
max_xy = 5;
max_z = 5;
newImage("res","32-bit",max_xy,max_z,1);
for(i=1;i<=max_xy;i++){
for(j=1;j<=max_z;j++){
	selectImage("orig");
	run("3D Morpho Gray", "image=orig kernel_image=orig kernel_type=Spherical operation=open radius_x="+i+" radius_y="+i+" radius_z="+j+" nb_iteration=1 white_object");
	mm = 0;
	for(s=1;s<=nSlices;s++){
		setSlice(s);
		getStatistics(area,mean);
		mm = mm + mean;
	}
	close();
	selectImage("res");
	setPixel(i,j,mm);
}}
setBatchMode(false);
