Mercurial > hg > orthanc-stone
view OrthancStone/Resources/Graveyard/RTStructTentativeReimplementation-BGO/DicomStructureSetSlicer2.cpp @ 2036:444527d34647 deep-learning
sanity check in textures
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Sat, 28 Jan 2023 12:05:23 +0100 |
parents | affde38b84de |
children | 07964689cb0b |
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-2022 Osimis S.A., Belgium * Copyright (C) 2021-2022 Sebastien Jodogne, ICTEAM UCLouvain, Belgium * * This program is free software: you can redistribute it and/or * modify it under the terms of the GNU Lesser 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 * Lesser General Public License for more details. * * You should have received a copy of the GNU Lesser General Public * License along with this program. If not, see * <http://www.gnu.org/licenses/>. **/ #ifdef BGO_ENABLE_DICOMSTRUCTURESETLOADER2 #include "DicomStructureSetSlicer2.h" #include "../Toolbox/GeometryToolbox.h" #include "../Volumes/IVolumeSlicer.h" #include "../Scene2D/PolylineSceneLayer.h" namespace OrthancStone { DicomStructureSetSlicer2::DicomStructureSetSlicer2(boost::shared_ptr<DicomStructureSet2> structureSet) : structureSet_(structureSet) {} IVolumeSlicer::IExtractedSlice* DicomStructureSetSlicer2::ExtractSlice(const CoordinateSystem3D& cuttingPlane) { // revision is always the same, hence 0 return new DicomStructureSetSlice2(structureSet_, 0, cuttingPlane); } DicomStructureSetSlice2::DicomStructureSetSlice2( boost::weak_ptr<DicomStructureSet2> structureSet, uint64_t revision, const CoordinateSystem3D& cuttingPlane) : structureSet_(structureSet.lock()) , isValid_(false) { bool opposite = false; if (structureSet_->GetStructuresCount() == 0) { isValid_ = false; } else { // some structures seen in real life have no polygons. We must be // careful bool found = false; size_t curStructure = 0; while (!found && curStructure < structureSet_->GetStructuresCount()) { if (structureSet_->GetStructure(curStructure).IsValid()) { found = true; const Vector normal = structureSet_->GetStructure(0).GetNormal(); isValid_ = ( GeometryToolbox::IsParallelOrOpposite(opposite, normal, cuttingPlane.GetNormal()) || GeometryToolbox::IsParallelOrOpposite(opposite, normal, cuttingPlane.GetAxisX()) || GeometryToolbox::IsParallelOrOpposite(opposite, normal, cuttingPlane.GetAxisY())); } } } } ISceneLayer* DicomStructureSetSlice2::CreateSceneLayer( const ILayerStyleConfigurator* configurator, const CoordinateSystem3D& cuttingPlane) { assert(isValid_); std::unique_ptr<PolylineSceneLayer> layer(new PolylineSceneLayer); layer->SetThickness(2); // thickness of the on-screen line for (size_t i = 0; i < structureSet_->GetStructuresCount(); i++) { const DicomStructure2& structure = structureSet_->GetStructure(i); if (structure.IsValid()) { const Color& color = structure.GetColor(); std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments; if (structure.Project(segments, cuttingPlane)) { for (size_t j = 0; j < segments.size(); j++) { PolylineSceneLayer::Chain chain; chain.resize(2); chain[0] = ScenePoint2D(segments[j].first.GetX(), segments[j].first.GetY()); chain[1] = ScenePoint2D(segments[j].second.GetX(), segments[j].second.GetY()); layer->AddChain(chain, false /* NOT closed */, color); } } } } return layer.release(); } } #endif // BGO_ENABLE_DICOMSTRUCTURESETLOADER2