Mercurial > hg > orthanc-stone
view Applications/Samples/SingleVolumeApplication.h @ 1327:4f8db2d202c8 broker
OrthancSeriesProgressiveLoader now has two modes that
can be selected at object creation :
- progressive (will first load jpeg50, then jpeg90 then PAM)
- non-progressive (will directly load PAM (uncompressed))
Please note that the slice loading order remains dynamic
and depending upon the slice that the client code wishes
to extract from the volume.
author | Benjamin Golinvaux <bgo@osimis.io> |
---|---|
date | Wed, 25 Mar 2020 14:34:27 +0100 |
parents | 8a0a62189f46 |
children |
line wrap: on
line source
/** * Stone of Orthanc * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics * Department, University Hospital of Liege, Belgium * Copyright (C) 2017-2020 Osimis S.A., Belgium * * This program is free software: you can redistribute it and/or * modify it under the terms of the GNU Affero General Public License * as published by the Free Software Foundation, either version 3 of * the License, or (at your option) any later version. * * This program is distributed in the hope that it will be useful, but * WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU * Affero General Public License for more details. * * You should have received a copy of the GNU Affero General Public License * along with this program. If not, see <http://www.gnu.org/licenses/>. **/ #pragma once #include "SampleApplicationBase.h" #include "../../Framework/dev.h" #include "../../Framework/Layers/LineMeasureTracker.h" #include "../../Framework/Layers/CircleMeasureTracker.h" #include <Core/Toolbox.h> #include <Core/Logging.h> #include <Plugins/Samples/Common/OrthancHttpConnection.h> // TODO REMOVE #include "../../Framework/Layers/DicomStructureSetSlicer.h" // TODO REMOVE #include "../../Framework/Toolbox/MessagingToolbox.h" // TODO REMOVE namespace OrthancStone { namespace Samples { class SingleVolumeApplication : public SampleApplicationBase { private: class Interactor : public VolumeImageInteractor { private: SliceViewerWidget& widget_; size_t layer_; protected: virtual void NotifySliceContentChange(const ISlicedVolume& volume, const size_t& sliceIndex, const Slice& slice) { const OrthancVolumeImage& image = dynamic_cast<const OrthancVolumeImage&>(volume); RenderStyle s = widget_.GetLayerStyle(layer_); if (image.FitWindowingToRange(s, slice.GetConverter())) { //printf("Windowing: %f => %f\n", s.customWindowCenter_, s.customWindowWidth_); widget_.SetLayerStyle(layer_, s); } } virtual void MouseOver(CairoContext& context, WorldSceneWidget& widget, const ViewportGeometry& view, double x, double y, IStatusBar* statusBar) { const SliceViewerWidget& w = dynamic_cast<const SliceViewerWidget&>(widget); Vector p = w.GetSlice().MapSliceToWorldCoordinates(x, y); printf("%f %f %f\n", p[0], p[1], p[2]); } public: Interactor(OrthancVolumeImage& volume, SliceViewerWidget& widget, VolumeProjection projection, size_t layer) : VolumeImageInteractor(volume, widget, projection), widget_(widget), layer_(layer) { } }; public: virtual void DeclareStartupOptions(boost::program_options::options_description& options) { boost::program_options::options_description generic("Sample options"); generic.add_options() ("series", boost::program_options::value<std::string>(), "Orthanc ID of the series") ("instance", boost::program_options::value<std::string>(), "Orthanc ID of a multi-frame instance that describes a 3D volume") ("threads", boost::program_options::value<unsigned int>()->default_value(3), "Number of download threads") ("projection", boost::program_options::value<std::string>()->default_value("axial"), "Projection of interest (can be axial, sagittal or coronal)") ("reverse", boost::program_options::value<bool>()->default_value(false), "Reverse the normal direction of the volume") ; options.add(generic); } virtual void Initialize(IStatusBar& statusBar, const boost::program_options::variables_map& parameters) { using namespace OrthancStone; if (parameters.count("series") > 1 || parameters.count("instance") > 1) { LOG(ERROR) << "Only one series or instance is allowed"; throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); } if (parameters.count("series") == 1 && parameters.count("instance") == 1) { LOG(ERROR) << "Cannot specify both a series and an instance"; throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); } std::string series; if (parameters.count("series") == 1) { series = parameters["series"].as<std::string>(); } std::string instance; if (parameters.count("instance") == 1) { instance = parameters["instance"].as<std::string>(); } if (series.empty() && instance.empty()) { LOG(ERROR) << "The series ID or instance ID is missing"; throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); } //unsigned int threads = parameters["threads"].as<unsigned int>(); //bool reverse = parameters["reverse"].as<bool>(); std::string tmp = parameters["projection"].as<std::string>(); Orthanc::Toolbox::ToLowerCase(tmp); VolumeProjection projection; if (tmp == "axial") { projection = VolumeProjection_Axial; } else if (tmp == "sagittal") { projection = VolumeProjection_Sagittal; } else if (tmp == "coronal") { projection = VolumeProjection_Coronal; } else { LOG(ERROR) << "Unknown projection: " << tmp; throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); } std::unique_ptr<SliceViewerWidget> widget(new SliceViewerWidget); #if 1 std::unique_ptr<OrthancVolumeImage> volume(new OrthancVolumeImage(context.GetWebService(), true)); if (series.empty()) { volume->ScheduleLoadInstance(instance); } else { volume->ScheduleLoadSeries(series); } widget->AddLayer(new VolumeImageMPRSlicer(*volume)); context_->AddInteractor(new Interactor(*volume, *widget, projection, 0)); context_->AddSlicedVolume(volume.release()); if (1) { RenderStyle s; //s.drawGrid_ = true; s.alpha_ = 1; s.windowing_ = ImageWindowing_Bone; widget->SetLayerStyle(0, s); } else { RenderStyle s; s.alpha_ = 1; s.applyLut_ = true; s.lut_ = Orthanc::EmbeddedResources::COLORMAP_JET; s.interpolation_ = ImageInterpolation_Bilinear; widget->SetLayerStyle(0, s); } #else std::unique_ptr<OrthancVolumeImage> ct(new OrthancVolumeImage(context_->GetWebService(), false)); //ct->ScheduleLoadSeries("15a6f44a-ac7b88fe-19c462d9-dddd918e-b01550d8"); // 0178023P //ct->ScheduleLoadSeries("dd069910-4f090474-7d2bba07-e5c10783-f9e4fb1d"); //ct->ScheduleLoadSeries("a04ecf01-79b2fc33-58239f7e-ad9db983-28e81afa"); // IBA //ct->ScheduleLoadSeries("03677739-1d8bca40-db1daf59-d74ff548-7f6fc9c0"); // 0522c0001 TCIA ct->ScheduleLoadSeries("295e8a13-dfed1320-ba6aebb2-9a13e20f-1b3eb953"); // Captain std::unique_ptr<OrthancVolumeImage> pet(new OrthancVolumeImage(context_->GetWebService(), true)); //pet->ScheduleLoadSeries("48d2997f-8e25cd81-dd715b64-bd79cdcc-e8fcee53"); // 0178023P //pet->ScheduleLoadSeries("aabad2e7-80702b5d-e599d26c-4f13398e-38d58a9e"); //pet->ScheduleLoadInstance("830a69ff-8e4b5ee3-b7f966c8-bccc20fb-d322dceb"); // IBA 1 //pet->ScheduleLoadInstance("337876a1-a68a9718-f15abccd-38faafa1-b99b496a"); // IBA 2 //pet->ScheduleLoadInstance("830a69ff-8e4b5ee3-b7f966c8-bccc20fb-d322dceb"); // IBA 3 //pet->ScheduleLoadInstance("269f26f4-0c83eeeb-2e67abbd-5467a40f-f1bec90c"); // 0522c0001 TCIA pet->ScheduleLoadInstance("f080888c-0ab7528a-f7d9c28c-84980eb1-ff3b0ae6"); // Captain 1 //pet->ScheduleLoadInstance("4f78055b-6499a2c5-1e089290-394acc05-3ec781c1"); // Captain 2 std::unique_ptr<StructureSetLoader> rtStruct(new StructureSetLoader(context_->GetWebService())); //rtStruct->ScheduleLoadInstance("c2ebc17b-6b3548db-5e5da170-b8ecab71-ea03add3"); // 0178023P //rtStruct->ScheduleLoadInstance("54460695-ba3885ee-ddf61ac0-f028e31d-a6e474d9"); // IBA //rtStruct->ScheduleLoadInstance("17cd032b-ad92a438-ca05f06a-f9e96668-7e3e9e20"); // 0522c0001 TCIA rtStruct->ScheduleLoadInstance("96c889ab-29fe5c54-dda6e66c-3949e4da-58f90d75"); // Captain widget->AddLayer(new VolumeImageMPRSlicer(*ct)); widget->AddLayer(new VolumeImageMPRSlicer(*pet)); widget->AddLayer(new DicomStructureSetSlicer(*rtStruct)); context_->AddInteractor(new Interactor(*pet, *widget, projection, 1)); //context_->AddInteractor(new VolumeImageInteractor(*ct, *widget, projection)); context_->AddSlicedVolume(ct.release()); context_->AddSlicedVolume(pet.release()); context_->AddVolumeLoader(rtStruct.release()); { RenderStyle s; //s.drawGrid_ = true; s.alpha_ = 1; s.windowing_ = ImageWindowing_Bone; widget->SetLayerStyle(0, s); } { RenderStyle s; //s.drawGrid_ = true; s.SetColor(255, 0, 0); // Draw missing PET layer in red s.alpha_ = 0.5; s.applyLut_ = true; s.lut_ = Orthanc::EmbeddedResources::COLORMAP_JET; s.interpolation_ = ImageInterpolation_Bilinear; s.windowing_ = ImageWindowing_Custom; s.customWindowCenter_ = 0; s.customWindowWidth_ = 128; widget->SetLayerStyle(1, s); } #endif statusBar.SetMessage("Use the keys \"b\", \"l\" and \"d\" to change Hounsfield windowing"); statusBar.SetMessage("Use the keys \"t\" to track the (X,Y,Z) mouse coordinates"); statusBar.SetMessage("Use the keys \"m\" to measure distances"); statusBar.SetMessage("Use the keys \"c\" to draw circles"); widget->SetTransmitMouseOver(true); context_->SetCentralWidget(widget.release()); } }; } }