cxDataReaderWriter.cpp

Go to the documentation of this file.

/*=========================================================================
This file is part of CustusX, an Image Guided Therapy Application.

Copyright (c) 2008-2014, SINTEF Department of Medical Technology
All rights reserved.

Redistribution and use in source and binary forms, with or without 
modification, are permitted provided that the following conditions are met:

1. Redistributions of source code must retain the above copyright notice, 
   this list of conditions and the following disclaimer.

2. Redistributions in binary form must reproduce the above copyright notice, 
   this list of conditions and the following disclaimer in the documentation 
   and/or other materials provided with the distribution.

3. Neither the name of the copyright holder nor the names of its contributors 
   may be used to endorse or promote products derived from this software 
   without specific prior written permission.

THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" 
AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 
IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE 
DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE 
FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 
DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR 
SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER 
CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, 
OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE 
OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
=========================================================================*/


#include "cxDataReaderWriter.h"

#include <sstream>
#include <vtkImageData.h>
#include <vtkMetaImageReader.h>
#include <vtkSmartPointer.h>
#include <vtkMetaImageWriter.h>

#include <vtkPolyData.h>
#include <vtkPolyDataReader.h>
#include <vtkPolyDataWriter.h>
#include <vtkSTLReader.h>
#include <vtkImageChangeInformation.h>
#include "vtkMINCImageReader.h"
#include "vtkTransform.h"
#include "vtkCommand.h"
#include <vtkPNGReader.h>

#include <QtCore>
#include <QDomDocument>
#include <QFileInfo>
#include <QFile>
#include <QTextStream>
#include <QDir>

#include "cxTransform3D.h"
#include "cxRegistrationTransform.h"
#include "cxLogger.h"
#include "cxTypeConversions.h"
#include "cxUtilHelpers.h"
#include "cxVideoSource.h"
#include "cxCustomMetaImage.h"
#include "cxConfig.h"

//#include "cxPointMetric.h"
//#include "cxDistanceMetric.h"
//#include "cxPlaneMetric.h"
//#include "cxAngleMetric.h"
//#include "cxShapedMetric.h"
//#include "cxSphereMetric.h"
//#include "cxFrameMetric.h"
//#include "cxToolMetric.h"

#include "cxImageLUT2D.h"
#include "cxImageTF3D.h"


typedef vtkSmartPointer<class vtkPNGReader> vtkPNGReaderPtr;

namespace cx
{

//---------------------------------------------------------
StaticMutexVtkLocker::StaticMutexVtkLocker()
{
/*      if (!mMutex)
                mMutex.reset(new QMutex(QMutex::Recursive));

        mMutex->lock();*/
}
StaticMutexVtkLocker::~StaticMutexVtkLocker()
{
//      mMutex->unlock();
}
boost::shared_ptr<QMutex> StaticMutexVtkLocker::mMutex;
//---------------------------------------------------------

class ErrorObserver: public vtkCommand
{
public:
        ErrorObserver()
        {
        }
        static ErrorObserver* New()
        {
                return new ErrorObserver;
        }
        virtual void Execute(vtkObject* caller, unsigned long, void* text)
        {
                mMessage = QString(reinterpret_cast<char*> (text));
        }
        QString mMessage;

        static bool checkedRead(vtkSmartPointer<vtkAlgorithm> reader, QString filename)
        {
                vtkSmartPointer<ErrorObserver> errorObserver = vtkSmartPointer<ErrorObserver>::New();
                reader->AddObserver("ErrorEvent", errorObserver);

                {
                        StaticMutexVtkLocker lock;
                        reader->Update();
                }
//              ErrorObserver::threadSafeUpdate(reader);

                if (!errorObserver->mMessage.isEmpty())
                {
                        reportError("Load of data [" + filename + "] failed with message:\n"
                                + errorObserver->mMessage);
                        return false;
                }
                return true;
        }
};

//-----
vtkImageDataPtr MetaImageReader::loadVtkImageData(QString filename)
{
        //load the image from file
        vtkMetaImageReaderPtr reader = vtkMetaImageReaderPtr::New();
        reader->SetFileName(cstring_cast(filename));
        reader->ReleaseDataFlagOn();

        if (!ErrorObserver::checkedRead(reader, filename))
                return vtkImageDataPtr();

        vtkImageChangeInformationPtr zeroer = vtkImageChangeInformationPtr::New();
        zeroer->SetInputConnection(reader->GetOutputPort());
        zeroer->SetOutputOrigin(0, 0, 0);
        zeroer->Update();
        return zeroer->GetOutput();
}

bool MetaImageReader::readInto(DataPtr data, QString filename)
{
        return this->readInto(boost::dynamic_pointer_cast<Image>(data), filename);
}
bool MetaImageReader::readInto(ImagePtr image, QString filename)
{
        if (!image)
                return false;

        CustomMetaImagePtr customReader = CustomMetaImage::create(filename);
        Transform3D rMd = customReader->readTransform();

        vtkImageDataPtr raw = this->loadVtkImageData(filename);
        if(!raw)
                return false;

        image->setVtkImageData(raw);
//      ImagePtr image(new Image(uid, raw));

        //  RegistrationTransform regTrans(rMd, QFileInfo(filename).lastModified(), "From MHD file");
        //  image->get_rMd_History()->addRegistration(regTrans);
        image->get_rMd_History()->setRegistration(rMd);
        image->setModality(customReader->readModality());
        image->setImageType(customReader->readImageType());

        bool ok1 = true;
        bool ok2 = true;
        double level = customReader->readKey("WindowLevel").toDouble(&ok1);
        double window = customReader->readKey("WindowWidth").toDouble(&ok2);

        if (ok1 && ok2)
        {
                image->setInitialWindowLevel(window, level);
                image->resetTransferFunctions();
        }

        return true;
}

//-----
DataPtr MetaImageReader::load(const QString& uid, const QString& filename)
{
        ImagePtr image(new Image(uid, vtkImageDataPtr()));
        this->readInto(image, filename);
        return image;
}

void MetaImageReader::saveImage(ImagePtr image, const QString& filename)
{
        vtkMetaImageWriterPtr writer = vtkMetaImageWriterPtr::New();
        writer->SetInputData(image->getBaseVtkImageData());
        writer->SetFileDimensionality(3);
        writer->SetFileName(cstring_cast(filename));
        QDir().mkpath(QFileInfo(filename).path());

//      std::cout << "SAVING MHD COMPRESSED " << filename << std::endl;
//      writer->SetCompression(true);
        writer->SetCompression(false);
//      writer->Update(); // caused writing of (null).0 files - not necessary
        writer->Write();

        writer = 0;

        CustomMetaImagePtr customReader = CustomMetaImage::create(filename);
        customReader->setTransform(image->get_rMd());
        customReader->setModality(image->getModality());
        customReader->setImageType(image->getImageType());
        customReader->setKey("WindowLevel", qstring_cast(image->getInitialWindowLevel()));
        customReader->setKey("WindowWidth", qstring_cast(image->getInitialWindowWidth()));
        customReader->setKey("Creator", QString("CustusX_%1").arg(CustusX_VERSION_STRING));
}



bool PNGImageReader::readInto(DataPtr data, QString filename)
{
        return this->readInto(boost::dynamic_pointer_cast<Image>(data), filename);
}

bool PNGImageReader::readInto(ImagePtr image, QString filename)
{
        if (!image)
                return false;
        vtkImageDataPtr raw = this->loadVtkImageData(filename);
        if(!raw)
                return false;
        image->setVtkImageData(raw);
        return true;
}

DataPtr PNGImageReader::load(const QString& uid, const QString& filename)
{
        ImagePtr image(new Image(uid, vtkImageDataPtr()));
        this->readInto(image, filename);
        return image;
}

vtkImageDataPtr PNGImageReader::loadVtkImageData(QString filename)
{
        vtkPNGReaderPtr pngReader = vtkPNGReaderPtr::New();
        pngReader->SetFileName(filename.toStdString().c_str());
        pngReader->Update();
        return pngReader->GetOutput();
}


bool PolyDataMeshReader::readInto(DataPtr data, QString filename)
{
        return this->readInto(boost::dynamic_pointer_cast<Mesh>(data), filename);
}

bool PolyDataMeshReader::readInto(MeshPtr mesh, QString filename)
{
        if (!mesh)
                return false;
        vtkPolyDataPtr raw = this->loadVtkPolyData(filename);
        if(!raw)
                return false;
        mesh->setVtkPolyData(raw);
        return true;
}
//-----
vtkPolyDataPtr PolyDataMeshReader::loadVtkPolyData(QString fileName)
{
        vtkPolyDataReaderPtr reader = vtkPolyDataReaderPtr::New();
        reader->SetFileName(cstring_cast(fileName));

        if (!ErrorObserver::checkedRead(reader, fileName))
                return vtkPolyDataPtr();

        vtkPolyDataPtr polyData = reader->GetOutput();
        return polyData;
}

DataPtr PolyDataMeshReader::load(const QString& uid, const QString& filename)
{
        MeshPtr mesh(new Mesh(uid));
        this->readInto(mesh, filename);
        return mesh;
}

bool StlMeshReader::readInto(DataPtr data, QString filename)
{
        return this->readInto(boost::dynamic_pointer_cast<Mesh>(data), filename);
}

bool StlMeshReader::readInto(MeshPtr mesh, QString filename)
{
        if (!mesh)
                return false;
        vtkPolyDataPtr raw = this->loadVtkPolyData(filename);
        if(!raw)
                return false;
        mesh->setVtkPolyData(raw);
        return true;
}


vtkPolyDataPtr StlMeshReader::loadVtkPolyData(QString fileName)
{
        vtkSTLReaderPtr reader = vtkSTLReaderPtr::New();
        reader->SetFileName(cstring_cast(fileName));

        if (!ErrorObserver::checkedRead(reader, fileName))
                return vtkPolyDataPtr();

        vtkPolyDataPtr polyData = reader->GetOutput();
        return polyData;
}

DataPtr StlMeshReader::load(const QString& uid, const QString& filename)
{
        MeshPtr mesh(new Mesh(uid));
        this->readInto(mesh, filename);
        return mesh;
}

DataReaderWriter::DataReaderWriter()
{
        mDataReaders.insert(DataReaderPtr(new MetaImageReader()));
        mDataReaders.insert(DataReaderPtr(new PolyDataMeshReader()));
        mDataReaders.insert(DataReaderPtr(new StlMeshReader()));
        mDataReaders.insert(DataReaderPtr(new PNGImageReader()));
}

DataReaderPtr DataReaderWriter::findReader(const QString& path, const QString& type)
{
        for (DataReadersType::iterator iter = mDataReaders.begin(); iter != mDataReaders.end(); ++iter)
        {
                if ((*iter)->canLoad(type, path))
                        return *iter;
        }
        return DataReaderPtr();
}

vtkImageDataPtr DataReaderWriter::loadVtkImageData(QString filename)
{
        DataReaderPtr reader = this->findReader(filename);
        if (reader)
                return reader->loadVtkImageData(filename);
        return vtkImageDataPtr();
}

vtkPolyDataPtr DataReaderWriter::loadVtkPolyData(QString filename)
{
        DataReaderPtr reader = this->findReader(filename);
        if (reader)
                return reader->loadVtkPolyData(filename);
        return vtkPolyDataPtr();
}

QString DataReaderWriter::findDataTypeFromFile(QString filename)
{
        DataReaderPtr reader = this->findReader(filename);
        if (reader)
                return reader->canLoadDataType();
        return "";
}

void DataReaderWriter::readInto(DataPtr data, QString path)
{
        DataReaderPtr reader = this->findReader(path, data->getType());
        if (reader)
                reader->readInto(data, path);

        if(data)
        {
                QFileInfo fileInfo(qstring_cast(path));
                data->setName(changeExtension(fileInfo.fileName(), ""));
                data->setFilename(path); // need path even when not set explicitly: nice for testing
        }

}

} // namespace cx