Mercurial > hg > orthanc-stone
view OrthancStone/Sources/Volumes/DicomVolumeImageMPRSlicer.cpp @ 1768:226718777702
fix DicomVolumeImageMPRSlicer::Slice::CreateSceneLayer() for opposite normals
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Tue, 11 May 2021 17:18:39 +0200 |
parents | 9ac2a65d4172 |
children | a217140dd41a |
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-2021 Osimis S.A., 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/>. **/ #include "DicomVolumeImageMPRSlicer.h" #include "../StoneException.h" #include "../Toolbox/ImageToolbox.h" #include <Logging.h> #include <OrthancException.h> //#include <Images/PngWriter.h> #include <Images/JpegWriter.h> namespace OrthancStone { void DicomVolumeImageMPRSlicer::Slice::CheckValid() const { if (!valid_) { LOG(ERROR) << "DicomVolumeImageMPRSlicer::Slice::CheckValid(): (!valid_)"; throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); } } DicomVolumeImageMPRSlicer::Slice::Slice(const DicomVolumeImage& volume, const CoordinateSystem3D& cuttingPlane) : volume_(volume), revision_(volume_.GetRevision()) { valid_ = (volume_.HasDicomParameters() && volume_.GetGeometry().DetectSlice(projection_, sliceIndex_, cuttingPlane)); } VolumeProjection DicomVolumeImageMPRSlicer::Slice::GetProjection() const { CheckValid(); return projection_; } unsigned int DicomVolumeImageMPRSlicer::Slice::GetSliceIndex() const { CheckValid(); return sliceIndex_; } ISceneLayer* DicomVolumeImageMPRSlicer::Slice::CreateSceneLayer(const ILayerStyleConfigurator* configurator, const CoordinateSystem3D& cuttingPlane) { CheckValid(); if (configurator == NULL) { throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer, "A style configurator is mandatory for textures"); } std::unique_ptr<TextureBaseSceneLayer> texture; { const DicomInstanceParameters& parameters = volume_.GetDicomParameters(); ImageBuffer3D::SliceReader reader(volume_.GetPixelData(), projection_, sliceIndex_); texture.reset(dynamic_cast<TextureBaseSceneLayer*> (configurator->CreateTextureFromDicom(reader.GetAccessor(), parameters))); } const CoordinateSystem3D& system = volume_.GetGeometry().GetProjectionGeometry(projection_); /** * TODO => There was a shift of (0.5, 0.5) introduced by * TextureBaseSceneLayer::GetTransform(). Is it an error? **/ Vector pixelSpacing = volume_.GetGeometry().GetVoxelDimensions(projection_); double x0, y0, x1, y1, x2, y2; cuttingPlane.ProjectPoint(x0, y0, system.GetOrigin()); cuttingPlane.ProjectPoint(x1, y1, system.GetOrigin() + system.GetAxisX() * pixelSpacing[0]); cuttingPlane.ProjectPoint(x2, y2, system.GetOrigin() + system.GetAxisY() * pixelSpacing[1]); /** A = [ a11 a12 ; a21 a22 ] (1) A * (0 ; 0) + (b1 ; b2) = (x0 ; y0) (2) A * (1 ; 0) + (b1 ; b2) = (x1 ; y1) (3) A * (0 ; 1) + (b1 ; b2) = (x2 ; y2) (2-1) A * (1 ; 0) = (x1 - x0 ; y1 - y0) <=> (a11 ; a21) = (x1 - x0 ; y1 - y0) (3-1) A * (0 ; 1) = (x2 - x0 ; y2 - y0) <=> (a12 ; a22) = (x2 - x0 ; y2 - y0) **/ Matrix m(3, 3); m(0, 0) = x1 - x0; // a11 m(0, 1) = x2 - x0; // a12 m(0, 2) = x0; // b1 m(1, 0) = y1 - y0; // a21 m(1, 1) = y2 - y0; // a22 m(1, 2) = y0; // b2 m(2, 0) = 0; m(2, 1) = 0; m(2, 2) = 1; texture->SetTransform(AffineTransform2D(m)); return texture.release(); } DicomVolumeImageMPRSlicer::~DicomVolumeImageMPRSlicer() { LOG(TRACE) << "DicomVolumeImageMPRSlicer::~DicomVolumeImageMPRSlicer()"; } IVolumeSlicer::IExtractedSlice* DicomVolumeImageMPRSlicer::ExtractSlice(const CoordinateSystem3D& cuttingPlane) { if (volume_->HasGeometry()) { return new Slice(*volume_, cuttingPlane); } else { return new IVolumeSlicer::InvalidSlice; } } }