Mercurial > hg > orthanc-stone

comparison Framework/Volumes/VolumeReslicer.cpp @ 742:fa5febe0f0c2

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
moved OrientedBoundingBox in the Volumes folder
author Sebastien Jodogne <s.jodogne@gmail.com>
date Wed, 22 May 2019 09:41:03 +0200
parents c3bbb130abc4
children cd13a062c9bd
comparison
equal deleted inserted replaced
741:c1d6a566dfd3 742:fa5febe0f0c2
228 228
229 public: 229 public:
230 FastRowIterator(const Orthanc::ImageAccessor& slice, 230 FastRowIterator(const Orthanc::ImageAccessor& slice,
231 const Extent2D& extent, 231 const Extent2D& extent,
232 const CoordinateSystem3D& plane, 232 const CoordinateSystem3D& plane,
233 const OrientedBoundingBox& box, 233 const OrientedVolumeBoundingBox& box,
234 unsigned int y) 234 unsigned int y)
235 { 235 {
236 const double width = static_cast<double>(slice.GetWidth()); 236 const double width = static_cast<double>(slice.GetWidth());
237 const double height = static_cast<double>(slice.GetHeight()); 237 const double height = static_cast<double>(slice.GetHeight());
238 assert(y < height); 238 assert(y < height);
283 { 283 {
284 private: 284 private:
285 const Orthanc::ImageAccessor& slice_; 285 const Orthanc::ImageAccessor& slice_;
286 const Extent2D& extent_; 286 const Extent2D& extent_;
287 const CoordinateSystem3D& plane_; 287 const CoordinateSystem3D& plane_;
288 const OrientedBoundingBox& box_; 288 const OrientedVolumeBoundingBox& box_;
289 unsigned int x_; 289 unsigned int x_;
290 unsigned int y_; 290 unsigned int y_;
291 291
292 public: 292 public:
293 SlowRowIterator(const Orthanc::ImageAccessor& slice, 293 SlowRowIterator(const Orthanc::ImageAccessor& slice,
294 const Extent2D& extent, 294 const Extent2D& extent,
295 const CoordinateSystem3D& plane, 295 const CoordinateSystem3D& plane,
296 const OrientedBoundingBox& box, 296 const OrientedVolumeBoundingBox& box,
297 unsigned int y) : 297 unsigned int y) :
298 slice_(slice), 298 slice_(slice),
299 extent_(extent), 299 extent_(extent),
300 plane_(plane), 300 plane_(plane),
301 box_(box), 301 box_(box),
340 TransferFunction Function> 340 TransferFunction Function>
341 static void ProcessImage(Orthanc::ImageAccessor& slice, 341 static void ProcessImage(Orthanc::ImageAccessor& slice,
342 const Extent2D& extent, 342 const Extent2D& extent,
343 const ImageBuffer3D& source, 343 const ImageBuffer3D& source,
344 const CoordinateSystem3D& plane, 344 const CoordinateSystem3D& plane,
345 const OrientedBoundingBox& box, 345 const OrientedVolumeBoundingBox& box,
346 float scaling, 346 float scaling,
347 float offset) 347 float offset)
348 { 348 {
349 typedef PixelShader<InputFormat, OutputFormat, Interpolation, Function> Shader; 349 typedef PixelShader<InputFormat, OutputFormat, Interpolation, Function> Shader;
350 350
384 Orthanc::PixelFormat OutputFormat> 384 Orthanc::PixelFormat OutputFormat>
385 static void ProcessImage(Orthanc::ImageAccessor& slice, 385 static void ProcessImage(Orthanc::ImageAccessor& slice,
386 const Extent2D& extent, 386 const Extent2D& extent,
387 const ImageBuffer3D& source, 387 const ImageBuffer3D& source,
388 const CoordinateSystem3D& plane, 388 const CoordinateSystem3D& plane,
389 const OrientedBoundingBox& box, 389 const OrientedVolumeBoundingBox& box,
390 ImageInterpolation interpolation, 390 ImageInterpolation interpolation,
391 bool hasLinearFunction, 391 bool hasLinearFunction,
392 float scaling, 392 float scaling,
393 float offset) 393 float offset)
394 { 394 {
450 template <typename RowIterator> 450 template <typename RowIterator>
451 static void ProcessImage(Orthanc::ImageAccessor& slice, 451 static void ProcessImage(Orthanc::ImageAccessor& slice,
452 const Extent2D& extent, 452 const Extent2D& extent,
453 const ImageBuffer3D& source, 453 const ImageBuffer3D& source,
454 const CoordinateSystem3D& plane, 454 const CoordinateSystem3D& plane,
455 const OrientedBoundingBox& box, 455 const OrientedVolumeBoundingBox& box,
456 ImageInterpolation interpolation, 456 ImageInterpolation interpolation,
457 bool hasLinearFunction, 457 bool hasLinearFunction,
458 float scaling, 458 float scaling,
459 float offset) 459 float offset)
460 { 460 {
499 499
500 500
501 501
502 void VolumeReslicer::CheckIterators(const ImageBuffer3D& source, 502 void VolumeReslicer::CheckIterators(const ImageBuffer3D& source,
503 const CoordinateSystem3D& plane, 503 const CoordinateSystem3D& plane,
504 const OrientedBoundingBox& box) const 504 const OrientedVolumeBoundingBox& box) const
505 { 505 {
506 for (unsigned int y = 0; y < slice_->GetHeight(); y++) 506 for (unsigned int y = 0; y < slice_->GetHeight(); y++)
507 { 507 {
508 FastRowIterator fast(*slice_, extent_, plane, box, y); 508 FastRowIterator fast(*slice_, extent_, plane, box, y);
509 SlowRowIterator slow(*slice_, extent_, plane, box, y); 509 SlowRowIterator slow(*slice_, extent_, plane, box, y);
783 // Firstly, compute the intersection of the source volumetric 783 // Firstly, compute the intersection of the source volumetric
784 // image with the reslicing plane. This leads to a polygon with 3 784 // image with the reslicing plane. This leads to a polygon with 3
785 // to 6 vertices. We compute the extent of the intersection 785 // to 6 vertices. We compute the extent of the intersection
786 // polygon, with respect to the coordinate system of the reslicing 786 // polygon, with respect to the coordinate system of the reslicing
787 // plane. 787 // plane.
788 OrientedBoundingBox box(geometry); 788 OrientedVolumeBoundingBox box(geometry);
789 789
790 if (!box.ComputeExtent(extent_, plane)) 790 if (!box.ComputeExtent(extent_, plane))
791 { 791 {
792 // The plane does not intersect with the bounding box of the volume 792 // The plane does not intersect with the bounding box of the volume
793 slice_.reset(new Orthanc::Image(outputFormat_, 0, 0, false)); 793 slice_.reset(new Orthanc::Image(outputFormat_, 0, 0, false));