prints out average lengths and torutoisity in the networks

pranathivemuri
1 parent 1b78fd76
Showing 1 changed file with 9 additions and 1 deletions Show diff stats
main.cpp
@@ -249,7 +249,15 @@ int main(int argc, char* argv[])
 		num_nets = 2;										//set the number of networks to two
 		float sigma = args["sigma"].as_float();				//get the sigma value from the user
 		T.load_obj(args.arg(1));							//load the second (test) network
-
+		double avgL, avgT;
+		avgL = GT.Lengths();
+		avgT = GT.Tortuosities();
+		std::cout<<avgL<<"---average segment length in the ground truth network"<<std::endl;
+		std::cout<<avgT<<"---average segment tortuosity in the ground truth network"<<std::endl;
+		avgL = T.Lengths();
+		avgT = T.Tortuosities();
+		std::cout<<avgL<<"---average segment length in the truth network"<<std::endl;
+		std::cout<<avgT<<"---average segment tortuosity in the truth network"<<std::endl;
 		GT = GT.resample(resample_rate * sigma);			//resample both networks based on the sigma value
 		T = T.resample(resample_rate * sigma);