again

layerfMRI · Apr 24, 2018 · 62bbeae · 62bbeae
1 parent 504a685
commit 62bbeae
Show file tree

Hide file tree

Showing 8 changed files with 501 additions and 1 deletion.
diff --git a/.DS_Store b/.DS_Store
diff --git a/LN_3DCOLUMNS.cpp b/LN_3DCOLUMNS.cpp
diff --git a/LN_3DGROW_LAYERS.cpp b/LN_3DGROW_LAYERS.cpp
@@ -1,6 +1,9 @@
 
 // Ausf¸hren mit . ./layers border_example_resized.nii brain_maskexample_resized.nii 0
 
+
+// @ Renzo to do check the vincinity loops. the "+1" are missing and the z_dimanetions is limieted by REAN dimention????
+
 #include <stdio.h>
 #include "nifti2_io.h"
 //#include "nifti2.h"
@@ -538,7 +541,7 @@ cout << " correct for pytagoras error .... " << endl;
 
 	    	for(int iy_i=max(0,(*(WMkoordy2_data   +  nxy*iz + nx*ix  + iy  ))-grow_vinc);     iy_i<min((*(WMkoordy2_data   +  nxy*iz + nx*ix  + iy  ))+grow_vinc,sizePhase); ++iy_i){
 	    	  for(int ix_i=max(0,(*(WMkoordx2_data   +  nxy*iz + nx*ix  + iy  ))-grow_vinc);   ix_i<min((*(WMkoordx2_data   +  nxy*iz + nx*ix  + iy  ))+grow_vinc,sizeRead);  ++ix_i){
-	    	   for(int iz_i=max(0,(*(WMkoordz2_data   +  nxy*iz + nx*ix  + iy  ))-grow_vinc); iz_i<min((*(WMkoordz2_data  +  nxy*iz + nx*ix  + iy  ))+grow_vinc,sizeRead);  ++iz_i){
+	    	   for(int iz_i=max(0,(*(WMkoordz2_data   +  nxy*iz + nx*ix  + iy  ))-grow_vinc); iz_i<min((*(WMkoordz2_data  +  nxy*iz + nx*ix  + iy  ))+grow_vinc,sizeSlice);  ++iz_i){
 			    if (*(nim_input_data  + nxy*iz_i + nx*ix_i  + iy_i  )  == 2){
 
 		  		 dist_i = dist((float)ix,(float)iy,(float)iz,(float)ix_i,(float)iy_i,(float)iz_i,dX,dY,dZ); 

diff --git a/LN_FIX_RIM.cpp b/LN_FIX_RIM.cpp
diff --git a/LN_FLOAT_ME.cpp b/LN_FLOAT_ME.cpp