Merge lp:~forsslundsystems/forssim/statistics into lp:forssim

[jocke]

[Martin Flodin (mflodin)]

Looks good. For future development the text should be changed, but as this was only a proof of concept it does what it should. I'll merge it into the main branch.

/Martin

review: Approve

=== modified file 'FsWisdom/ApplicationNode.cpp'
--- FsWisdom/ApplicationNode.cpp	2009-05-06 09:26:07 +0000
+++ FsWisdom/ApplicationNode.cpp	2009-05-11 08:37:54 +0000
@@ -291,6 +291,7 @@
 		this->setOverrideCursor(QCursor(Qt::WaitCursor));
 		volumeModel = static_cast <VolumeModel*> (glWindow->getSceneNode("VM"));
 		volumeModel->saveTimeStamps();
+		volumeModel->saveStatistics();
 		this->restoreOverrideCursor();
 		// falls down into no
 	case QMessageBox::No:

=== modified file 'src/VolumeModel.cpp'
--- src/VolumeModel.cpp	2009-05-04 08:19:33 +0000
+++ src/VolumeModel.cpp	2009-05-11 08:37:54 +0000
@@ -836,6 +836,78 @@
 	timeStamps->saveDataToNrrdFile(filePath,noOfDimensions,dimensionLengths,spacings);
 }
 
+void VolumeModel::saveStatistics()
+{
+
+	LoggingNode* loggingNode = static_cast< LoggingNode * >(loggingNodeField->getValue() );
+	string userLoggingPath = loggingNode->userLoggingPathField->getValue();
+
+	MaterialSegmentationModel* materialSegmentationNode =
+		static_cast< MaterialSegmentationModel * >(materialSegmentationModelField->getValue());
+
+    const float mm3PerM3 = 1000000.0;
+	// construct string with current date and time
+	time_t rawTime = time(0);
+	struct tm* timeInfo;
+	char strBuffer [80];
+	timeInfo = localtime ( &rawTime );
+	strftime (strBuffer,80,"%y%m%d%H%M%S",timeInfo);
+	string timeString(strBuffer);
+
+	// construct file name
+	string filePath = userLoggingPath+"/";
+	filePath.append(timeString).append("_stat").append(".txt");
+
+	//creating file
+	ofstream file;
+	file.open(filePath.c_str(),ios::out);
+
+	/*Data*/
+	//Total numbers of voxels and segments
+	int totVoxels = height*width*depth;
+	file<<"Total numbers of voxels: " << totVoxels << " (including Air)" << endl;
+	file << "Number of Segments: "<<materialSegmentationNode->getNoOfSegments()<<endl;
+
+	//Total time on case
+	double timeSpendBoring = 0;
+	if(timestamp_map.size()>1)
+		timeSpendBoring = (*(timestamp_map.rbegin())).first - (*(timestamp_map.begin())).first;
+	file << "Time spent boring: "<<timeSpendBoring<<" secs"<<endl<<endl;
+
+	//Number of Forbidden voxels
+	ForbiddenSegmentationModel* forbiddenSegmentationModel = static_cast< ForbiddenSegmentationModel * >(forbiddenSegmentationModelField->getValue());
+
+	//Number of voxels drilled in each segment
+	int noOfSegments = materialSegmentationNode->getNoOfSegments();
+	int voxelsEachSegArray[noOfSegments];
+
+	for (int i = 0; i < noOfSegments; i++)
+	{
+		voxelsEachSegArray[i] = 0;
+	}
+
+	multimap<double, Vec3i>::iterator it;
+	for ( it=timestamp_map.begin() ; it !=timestamp_map.end(); it++ )
+	{
+		int temp = getMaterialSegmentationId((*it).second);
+		voxelsEachSegArray[temp] = voxelsEachSegArray[temp] + 1;
+	}
+	file << " Voxels drilled in " << endl;
+	for (int i = 1; i < noOfSegments; i++) //Start at 1 since 0=Air
+	{
+		file << "Segment " << i << " - " << materialSegmentationNode->getName(i) << ": " << voxelsEachSegArray[i] << endl;
+	}
+
+	//Forbidden segments
+	file << endl;
+	file << "(TODO These values are not the same as displayed in GUI. Check code...) Number of voxels drilled in forbidden areas" <<endl;
+	file << "Wisdom Backside: " << forbiddenSegmentationModel->noOfVoxelsBoredByUser->getValue()[0] << endl;
+	file << "Other teeth: " << forbiddenSegmentationModel->noOfVoxelsBoredByUser->getValue()[1] <<endl<<endl;
+
+	file.close();
+	cout << "Saving statistics to file " << filePath << endl;
+}
+
 /****************************************************************/
 
 H3D::Vec3f VolumeModel::getGlobalCoordinate()

=== modified file 'src/VolumeModel.h'
--- src/VolumeModel.h	2009-03-31 09:45:57 +0000
+++ src/VolumeModel.h	2009-05-11 08:37:54 +0000
@@ -324,6 +324,11 @@
 	void saveTimeStamps();
 
 	/**
+	 * Save drilling statistics. Is done on Exit.
+	 */
+	void saveStatistics();
+
+	/**
 	 * Remove a partition giving its name
 	 */
 	void removePartitionByName(string partitionName);