std::shared_ptr<GLTextureRGBA> Application::loadTextureRGBA( const std::string& partialFile) { ImageRGBA image; ReadImageFromFile(context->getFullPath(partialFile), image); return std::shared_ptr<GLTextureRGBA>( new GLTextureRGBA(image, true, context)); }
void buildImageIntensityModelOnROI(const char* referenceFilename, const char* maskFilename, const char* dir, const char* outputImageFilename) { typedef itk::PCAModelBuilder<RepresenterType> ModelBuilderType; typedef itk::StatisticalModel<RepresenterType> StatisticalModelType; typedef std::vector<std::string> StringVectorType; typedef itk::DataManager<RepresenterType> DataManagerType; RepresenterType::Pointer representer = RepresenterType::New(); typedef itk::ImageFileReader< ImageType > MaskReaderType; MaskReaderType::Pointer maskReader = MaskReaderType::New(); maskReader->SetFileName( maskFilename ); maskReader->Update(); representer->SetReference( ReadImageFromFile(referenceFilename), maskReader->GetOutput() ); StringVectorType filenames; getdir(dir, filenames, ".vtk"); DataManagerType::Pointer dataManager = DataManagerType::New(); dataManager->SetRepresenter(representer); for (StringVectorType::const_iterator it = filenames.begin(); it != filenames.end(); it++) { std::string fullpath = (std::string(dir) + "/") + *it; dataManager->AddDataset( ReadImageFromFile(fullpath), fullpath.c_str()); } ModelBuilderType::Pointer pcaModelBuilder = ModelBuilderType::New(); StatisticalModelType::Pointer model = pcaModelBuilder->BuildNewModel(dataManager->GetSampleDataStructure(), 0); std::cout<<"dimensionality of the data: "<<model->GetDomain().GetNumberOfPoints()<<", dimension of the images: "<<(*dataManager->GetSampleDataStructure().begin())->GetSample()->GetLargestPossibleRegion().GetNumberOfPixels()<<std::endl; std::cout<<"writing the mean sample to a png file..."<<std::endl; typedef itk::ImageFileWriter< ImageType > ImageWriterType; ImageWriterType::Pointer writer = ImageWriterType::New(); writer->SetFileName( outputImageFilename ); writer->SetInput(model->DrawSample()); writer->Update(); }