+
+
+//--------------------------------------------------------------------
+template<class ImageType>
+bool AnatomicalFeatureDatabase::
+CheckImage(std::string tag)
+{
+ if (m_MapOfTag.find(tag) == m_MapOfTag.end()) {
+ return false;
+ }
+ else {
+ typename ImageType::Pointer image;
+ if (m_MapOfImage[tag]) {
+ image = static_cast<ImageType *>(m_MapOfImage[tag]);
+ }
+ else {
+ std::string s = m_MapOfTag[tag];
+ // Read the file
+ std::string n = GetPath()+"/"+s;
+ // image = readImage<ImageType>();
+ typedef itk::ImageFileReader<ImageType> ReaderType;
+ typename ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName(n.c_str());
+ try {
+ reader->Update();
+ } catch(itk::ExceptionObject & err) {
+ return false;
+ }
+ }
+ }
+ return true;
+}
+//--------------------------------------------------------------------
+
+
+//--------------------------------------------------------------------
+template<class ImageType>
+void AnatomicalFeatureDatabase::
+SetImage(TagType tag, std::string f, typename ImageType::Pointer image, bool write)
+{
+ SetImageFilename(tag, f);
+ m_MapOfImage[tag] = &(*image);
+ // I add a reference count because the cache is not a smartpointer
+ image->SetReferenceCount(image->GetReferenceCount()+1);
+ if (write) {
+ writeImage<ImageType>(image, f);
+ }
+}
+//--------------------------------------------------------------------
+
+
+//--------------------------------------------------------------------
+template<class ImageType>
+void AnatomicalFeatureDatabase::
+ReleaseImage(std::string tag)
+{
+ if (m_MapOfTag.find(tag) == m_MapOfTag.end()) {
+ clitkExceptionMacro("Could not find the tag <" << tag << "> of type Image Filename in the DB");
+ }
+ else {
+ typename ImageType::Pointer image = GetImage<ImageType>(tag);
+ image->SetReferenceCount(image->GetReferenceCount()-1);
+ m_MapOfImage.erase(tag);
+ }
+}
+//--------------------------------------------------------------------