/
qMRMLSegmentEditorWidget.cxx
3915 lines (3441 loc) · 149 KB
/
qMRMLSegmentEditorWidget.cxx
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
/*==============================================================================
Program: 3D Slicer
Copyright (c) Laboratory for Percutaneous Surgery (PerkLab)
Queen's University, Kingston, ON, Canada. All Rights Reserved.
See COPYRIGHT.txt
or http://www.slicer.org/copyright/copyright.txt for details.
Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.
This file was originally developed by Csaba Pinter, PerkLab, Queen's University
and was supported through the Applied Cancer Research Unit program of Cancer Care
Ontario with funds provided by the Ontario Ministry of Health and Long-Term Care
==============================================================================*/
// Segmentations includes
#include "qMRMLSegmentEditorWidget.h"
#include "qMRMLSortFilterSegmentsProxyModel.h"
#include "ui_qMRMLSegmentEditorWidget.h"
#include "vtkMRMLSegmentationNode.h"
#include "vtkMRMLSegmentationDisplayNode.h"
#include "vtkMRMLSegmentEditorNode.h"
#include "qMRMLSegmentationGeometryDialog.h"
// vtkSegmentationCore Includes
#include "vtkSegmentation.h"
#include "vtkSegmentationHistory.h"
#include "vtkSegment.h"
#include "vtkOrientedImageData.h"
#include "vtkOrientedImageDataResample.h"
#include "vtkSlicerSegmentationsModuleLogic.h"
#include "vtkBinaryLabelmapToClosedSurfaceConversionRule.h"
// Segment editor effects includes
#include "qSlicerSegmentEditorAbstractEffect.h"
#include "qSlicerSegmentEditorAbstractLabelEffect.h"
#include "qSlicerSegmentEditorEffectFactory.h"
// VTK includes
#include <vtkAddonMathUtilities.h>
#include <vtkAlgorithmOutput.h>
#include <vtkCallbackCommand.h>
#include <vtkCollection.h>
#include <vtkDataArray.h>
#include <vtkGeneralTransform.h>
#include <vtkImageThreshold.h>
#include <vtkImageExtractComponents.h>
#include <vtkInteractorObserver.h>
#include <vtkNew.h>
#include <vtkPointData.h>
#include <vtkRenderer.h>
#include <vtkRendererCollection.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkSmartPointer.h>
#include <vtkWeakPointer.h>
// Slicer includes
#include <vtkMRMLSliceLogic.h>
#include <vtkSlicerApplicationLogic.h>
// MRML includes
#include <vtkMRMLCameraDisplayableManager.h>
#include <vtkMRMLCameraWidget.h>
#include <vtkMRMLCrosshairDisplayableManager.h>
#include <vtkMRMLLabelMapVolumeNode.h>
#include <vtkMRMLSelectionNode.h>
#include <vtkMRMLScene.h>
#include <vtkMRMLSliceIntersectionWidget.h>
#include <vtkMRMLSliceNode.h>
#include <vtkMRMLTransformNode.h>
#include <vtkMRMLViewNode.h>
#include <vtkMRMLInteractionNode.h>
#include <vtkMRMLSliceCompositeNode.h>
// Slicer includes
#include <qSlicerApplication.h>
#include <qSlicerLayoutManager.h>
#include "qSlicerModuleManager.h"
#include <qSlicerAbstractModule.h>
#include <qSlicerAbstractModuleWidget.h>
#include <qMRMLSegmentationFileExportWidget.h>
#include <qMRMLSliceWidget.h>
#include <qMRMLSliceView.h>
#include <qMRMLThreeDWidget.h>
#include <qMRMLThreeDView.h>
// Qt includes
#include <QAbstractItemView>
#include <QAction>
#include <QButtonGroup>
#include <QDebug>
#include <QGridLayout>
#include <QInputDialog>
#include <QMainWindow>
#include <QMenu>
#include <QMessageBox>
#include <QPointer>
#include <QScrollArea>
#include <QSettings>
#include <QShortcut>
#include <QTableView>
#include <QToolButton>
#include <QVBoxLayout>
#include <QWidgetAction>
// CTK includes
#include <ctkSliderWidget.h>
#include <ctkCollapsibleButton.h>
static const int BINARY_LABELMAP_SCALAR_TYPE = VTK_UNSIGNED_CHAR;
// static const unsigned char BINARY_LABELMAP_VOXEL_FULL = 1; // unused
static const unsigned char BINARY_LABELMAP_VOXEL_EMPTY = 0;
static const char NULL_EFFECT_NAME[] = "NULL";
//---------------------------------------------------------------------------
class vtkSegmentEditorEventCallbackCommand : public vtkCallbackCommand
{
public:
static vtkSegmentEditorEventCallbackCommand *New()
{
return new vtkSegmentEditorEventCallbackCommand;
}
/// Segment editor widget observing the event
QPointer<qMRMLSegmentEditorWidget> EditorWidget;
/// Slice widget or 3D widget
QPointer<qMRMLWidget> ViewWidget;
};
//-----------------------------------------------------------------------------
struct SegmentEditorEventObservation
{
vtkSmartPointer<vtkSegmentEditorEventCallbackCommand> CallbackCommand;
vtkWeakPointer<vtkObject> ObservedObject;
QVector<int> ObservationTags;
};
//-----------------------------------------------------------------------------
// qMRMLSegmentEditorWidgetPrivate methods
//-----------------------------------------------------------------------------
class qMRMLSegmentEditorWidgetPrivate: public Ui_qMRMLSegmentEditorWidget
{
Q_DECLARE_PUBLIC(qMRMLSegmentEditorWidget);
protected:
qMRMLSegmentEditorWidget* const q_ptr;
public:
qMRMLSegmentEditorWidgetPrivate(qMRMLSegmentEditorWidget& object);
~qMRMLSegmentEditorWidgetPrivate();
void init();
/// Simple mechanism to let the effects know that default modifier labelmap has changed
void notifyEffectsOfReferenceGeometryChange(const std::string& geometry);
/// Simple mechanism to let the effects know that master volume has changed
void notifyEffectsOfMasterVolumeNodeChange();
/// Simple mechanism to let the effects know that layout has changed
void notifyEffectsOfLayoutChange();
/// Select first segment in table view
void selectFirstSegment();
/// Enable or disable effects and their options based on input selection
void updateEffectsEnabledFromMRML();
/// Set cursor for effect. If effect is nullptr then the cursor is reset to default.
void setEffectCursor(qSlicerSegmentEditorAbstractEffect* effect);
/// Updates default modifier labelmap based on reference geometry (to set origin, spacing, and directions)
/// and existing segments (to set extents). If reference geometry conversion parameter is empty
/// then existing segments are used for determining origin, spacing, and directions and the resulting
/// geometry is written to reference geometry conversion parameter.
bool resetModifierLabelmapToDefault();
/// Updates selected segment labelmap in a geometry aligned with default modifierLabelmap.
bool updateSelectedSegmentLabelmap();
/// Updates a resampled master volume in a geometry aligned with default modifierLabelmap.
bool updateAlignedMasterVolume();
/// Updates mask labelmap.
/// Geometry of mask will be the same as current modifierLabelmap.
/// This mask only considers segment-based regions (and ignores masking based on
/// master volume intensity).
bool updateMaskLabelmap();
bool updateReferenceGeometryImage();
static std::string getReferenceImageGeometryFromSegmentation(vtkSegmentation* segmentation);
std::string referenceImageGeometry();
bool segmentationDisplayableInView(vtkMRMLAbstractViewNode* viewNode);
QToolButton* toolButton(qSlicerSegmentEditorAbstractEffect* effect);
/// Return segmentation node's internal labelmap IJK to renderer world coordinate transform.
/// If cannot be retrieved (segmentation is not defined, non-linearly transformed, etc.)
/// then false is returned;
bool segmentationIJKToRAS(vtkMatrix4x4* ijkToRas);
/// Updates surface smoothing factor in segmentation node and updates surface representation
/// if it is enabled.
bool setSurfaceSmoothingFactor(double smoothingFactor);
public:
/// Segment editor parameter set node containing all selections and working images
vtkWeakPointer<vtkMRMLSegmentEditorNode> ParameterSetNode;
vtkWeakPointer<vtkMRMLSegmentationNode> SegmentationNode;
vtkSmartPointer<vtkSegmentationHistory> SegmentationHistory;
vtkWeakPointer<vtkMRMLScalarVolumeNode> MasterVolumeNode;
// Observe InteractionNode to detect when mouse mode is changed
vtkWeakPointer<vtkMRMLInteractionNode> InteractionNode;
/// Lock widget to make segmentation read-only.
// In the future locked state may be read from the Segmentation node.
bool Locked;
/// Ordering of effects
QStringList EffectNameOrder;
bool UnorderedEffectsVisible;
int EffectColumnCount;
/// List of registered effect instances
QList<qSlicerSegmentEditorAbstractEffect*> RegisteredEffects;
/// Active effect
qSlicerSegmentEditorAbstractEffect* ActiveEffect;
/// Last active effect
/// Stored to allow quick toggling between no effect/last active effect.
qSlicerSegmentEditorAbstractEffect* LastActiveEffect;
/// Structure containing necessary objects for each slice and 3D view handling interactions
QVector<SegmentEditorEventObservation> EventObservations;
/// List of view node IDs where custom cursor is set
QSet<QString> CustomCursorInViewNodeIDs;
/// Indicates if views and layouts are observed
/// (essentially, the widget is active).
bool ViewsObserved;
/// List of view node IDs in display nodes, which were specified when views observation was set up.
/// If node IDs change (segmentation node is shown/hidden in a specific view) then view observations has to be refreshed.
QMap<QString, std::vector<std::string> > ObservedViewNodeIDs; // <SegmentationDisplayNodeID, ViewNodeIDs>
bool AutoShowMasterVolumeNode;
/// Button group for the effects
QButtonGroup EffectButtonGroup;
/// These volumes are owned by this widget and a pointer is given to each effect
/// so that they can access and modify it
vtkOrientedImageData* AlignedMasterVolume;
/// Modifier labelmap that is kept in memory to avoid memory reallocations on each editing operation.
/// When update of this labelmap is requested its geometry is reset and its content is cleared.
vtkOrientedImageData* ModifierLabelmap;
vtkOrientedImageData* SelectedSegmentLabelmap;
vtkOrientedImageData* MaskLabelmap;
/// Image that contains reference geometry. Scalars are not allocated.
vtkOrientedImageData* ReferenceGeometryImage;
/// Input data that is used for computing AlignedMasterVolume.
/// It is stored so that it can be determined that the master volume has to be updated
vtkMRMLVolumeNode* AlignedMasterVolumeUpdateMasterVolumeNode;
vtkMRMLTransformNode* AlignedMasterVolumeUpdateMasterVolumeNodeTransform;
vtkMRMLTransformNode* AlignedMasterVolumeUpdateSegmentationNodeTransform;
int MaskModeComboBoxFixedItemsCount;
/// If reference geometry changes compared to this value then we notify effects and
/// set this value to the current value. This allows notifying effects when there is a change.
std::string LastNotifiedReferenceImageGeometry;
QList< QShortcut* > KeyboardShortcuts;
Qt::ToolButtonStyle EffectButtonStyle;
// When segmentation node selector is visible then rotate warning is displayed in that row.
// However, if the node selector is hidden then the rotation warning button would take
// and entire row.
// To prevent this, rotation warning button is displayed in the add/remove/etc. segment
// button row when segmentation node selector is hidden.
// Qt does not have an API to check if a widget is in a layout, therefore we store this
// information in this flag.
bool RotateWarningInNodeSelectorLayout;
QAction* SurfaceSmoothingEnabledAction;
ctkSliderWidget* SurfaceSmoothingSlider;
QString DefaultTerminologyEntrySettingsKey;
QString DefaultTerminologyEntry;
};
//-----------------------------------------------------------------------------
qMRMLSegmentEditorWidgetPrivate::qMRMLSegmentEditorWidgetPrivate(qMRMLSegmentEditorWidget& object)
: q_ptr(&object)
, Locked(false)
, ActiveEffect(nullptr)
, LastActiveEffect(nullptr)
, ViewsObserved(false)
, AutoShowMasterVolumeNode(true)
, AlignedMasterVolume(nullptr)
, ModifierLabelmap(nullptr)
, SelectedSegmentLabelmap(nullptr)
, MaskLabelmap(nullptr)
, ReferenceGeometryImage(nullptr)
, AlignedMasterVolumeUpdateMasterVolumeNode(nullptr)
, AlignedMasterVolumeUpdateMasterVolumeNodeTransform(nullptr)
, AlignedMasterVolumeUpdateSegmentationNodeTransform(nullptr)
, MaskModeComboBoxFixedItemsCount(0)
, EffectButtonStyle(Qt::ToolButtonIconOnly)
, RotateWarningInNodeSelectorLayout(true)
{
this->AlignedMasterVolume = vtkOrientedImageData::New();
this->ModifierLabelmap = vtkOrientedImageData::New();
this->MaskLabelmap = vtkOrientedImageData::New();
this->SelectedSegmentLabelmap = vtkOrientedImageData::New();
this->ReferenceGeometryImage = vtkOrientedImageData::New();
this->SegmentationHistory = vtkSmartPointer<vtkSegmentationHistory>::New();
// Define default effect order
this->EffectNameOrder
// Thresholding is the the starting point for most segmentations
// (it can often create usable segmentation by itself, or used to define intensity range for painting)
<< "Threshold"
// Local painting
<< "Paint" << "Draw" << "Erase" << "Level tracing" << "Grow from seeds" << "Fill between slices"
// Global processing
<< "Margin" << "Hollow" << "Smoothing"
// Global splitting, merging
<< "Scissors" << "Islands" << "Logical operators"
// Operating on volumes
<< "Mask volume";
this->UnorderedEffectsVisible = true;
this->DefaultTerminologyEntrySettingsKey = "Segmentations/DefaultTerminologyEntry";
this->EffectColumnCount = 2;
}
//-----------------------------------------------------------------------------
qMRMLSegmentEditorWidgetPrivate::~qMRMLSegmentEditorWidgetPrivate()
{
Q_Q(qMRMLSegmentEditorWidget);
q->removeViewObservations();
foreach(qSlicerSegmentEditorAbstractEffect* effect, this->RegisteredEffects)
{
delete effect;
}
this->RegisteredEffects.clear();
if (this->AlignedMasterVolume)
{
this->AlignedMasterVolume->Delete();
this->AlignedMasterVolume = nullptr;
}
if (this->ModifierLabelmap)
{
this->ModifierLabelmap->Delete();
this->ModifierLabelmap = nullptr;
}
if (this->MaskLabelmap)
{
this->MaskLabelmap->Delete();
this->MaskLabelmap = nullptr;
}
if (this->SelectedSegmentLabelmap)
{
this->SelectedSegmentLabelmap->Delete();
this->SelectedSegmentLabelmap = nullptr;
}
if (this->ReferenceGeometryImage)
{
this->ReferenceGeometryImage->Delete();
this->ReferenceGeometryImage = nullptr;
}
}
//-----------------------------------------------------------------------------
void qMRMLSegmentEditorWidgetPrivate::init()
{
Q_Q(qMRMLSegmentEditorWidget);
this->setupUi(q);
this->SliceRotateWarningButton->setMaximumHeight(this->SegmentationNodeComboBox->sizeHint().height());
this->SliceRotateWarningButton->setMaximumWidth(this->SegmentationNodeComboBox->sizeHint().height());
this->SpecifyGeometryButton->setMaximumHeight(this->MasterVolumeNodeComboBox->sizeHint().height());
this->SpecifyGeometryButton->setMaximumWidth(this->MasterVolumeNodeComboBox->sizeHint().height());
this->MaskModeComboBox->addItem(qMRMLSegmentEditorWidget::tr("Everywhere"), vtkMRMLSegmentEditorNode::PaintAllowedEverywhere);
this->MaskModeComboBox->addItem(qMRMLSegmentEditorWidget::tr("Inside all segments"), vtkMRMLSegmentEditorNode::PaintAllowedInsideAllSegments);
this->MaskModeComboBox->addItem(qMRMLSegmentEditorWidget::tr("Inside all visible segments"), vtkMRMLSegmentEditorNode::PaintAllowedInsideVisibleSegments);
this->MaskModeComboBox->addItem(qMRMLSegmentEditorWidget::tr("Outside all segments"), vtkMRMLSegmentEditorNode::PaintAllowedOutsideAllSegments);
this->MaskModeComboBox->addItem(qMRMLSegmentEditorWidget::tr("Outside all visible segments"), vtkMRMLSegmentEditorNode::PaintAllowedOutsideVisibleSegments);
this->MaskModeComboBox->insertSeparator(this->MaskModeComboBox->count());
this->MaskModeComboBoxFixedItemsCount = this->MaskModeComboBox->count();
this->OverwriteModeComboBox->addItem(qMRMLSegmentEditorWidget::tr("Overwrite all"), vtkMRMLSegmentEditorNode::OverwriteAllSegments);
this->OverwriteModeComboBox->addItem(qMRMLSegmentEditorWidget::tr("Overwrite visible"), vtkMRMLSegmentEditorNode::OverwriteVisibleSegments);
this->OverwriteModeComboBox->addItem(qMRMLSegmentEditorWidget::tr("Allow overlap"), vtkMRMLSegmentEditorNode::OverwriteNone);
this->SwitchToSegmentationsButton->setIcon(q->style()->standardIcon(QStyle::SP_ArrowRight));
QMenu* segmentationsButtonMenu = new QMenu(qMRMLSegmentEditorWidget::tr("Segmentations"), this->SwitchToSegmentationsButton);
QAction* importExportAction = new QAction("Import/export nodes...", segmentationsButtonMenu);
segmentationsButtonMenu->addAction(importExportAction);
QObject::connect(importExportAction, SIGNAL(triggered()), q, SLOT(onImportExportActionClicked()));
QAction* exportToFileAction = new QAction("Export to files...", segmentationsButtonMenu);
segmentationsButtonMenu->addAction(exportToFileAction);
QObject::connect(exportToFileAction, SIGNAL(triggered()), q, SLOT(onExportToFilesActionClicked()));
this->SwitchToSegmentationsButton->setMenu(segmentationsButtonMenu);
QMenu* show3DButtonMenu = new QMenu(qMRMLSegmentEditorWidget::tr("Show 3D"), this->Show3DButton);
this->SurfaceSmoothingEnabledAction = new QAction(qMRMLSegmentEditorWidget::tr("Surface smoothing"), show3DButtonMenu);
this->SurfaceSmoothingEnabledAction->setToolTip(qMRMLSegmentEditorWidget::tr("Apply smoothing when converting binary lablemap to closed surface representation."));
this->SurfaceSmoothingEnabledAction->setCheckable(true);
show3DButtonMenu->addAction(this->SurfaceSmoothingEnabledAction);
QObject::connect(this->SurfaceSmoothingEnabledAction, SIGNAL(toggled(bool)), q, SLOT(onEnableSurfaceSmoothingToggled(bool)));
QMenu* surfaceSmoothingFactorMenu = new QMenu(qMRMLSegmentEditorWidget::tr("Smoothing factor"), show3DButtonMenu);
surfaceSmoothingFactorMenu->setObjectName("slicerSpacingManualMode");
surfaceSmoothingFactorMenu->setIcon(QIcon(":/Icon/SlicerManualSliceSpacing.png"));
this->SurfaceSmoothingSlider = new ctkSliderWidget(surfaceSmoothingFactorMenu);
this->SurfaceSmoothingSlider->setToolTip(qMRMLSegmentEditorWidget::tr("Higher value means stronger smoothing during closed surface representation conversion."));
this->SurfaceSmoothingSlider->setDecimals(2);
this->SurfaceSmoothingSlider->setRange(0.0, 1.0);
this->SurfaceSmoothingSlider->setSingleStep(0.1);
this->SurfaceSmoothingSlider->setValue(0.5);
this->SurfaceSmoothingSlider->setTracking(false);
QObject::connect(this->SurfaceSmoothingSlider, SIGNAL(valueChanged(double)),
q, SLOT(onSurfaceSmoothingFactorChanged(double)));
QWidgetAction* sliceSpacingAction = new QWidgetAction(surfaceSmoothingFactorMenu);
sliceSpacingAction->setCheckable(true);
sliceSpacingAction->setDefaultWidget(this->SurfaceSmoothingSlider);
surfaceSmoothingFactorMenu->addAction(sliceSpacingAction);
show3DButtonMenu->addMenu(surfaceSmoothingFactorMenu);
this->Show3DButton->setMenu(show3DButtonMenu);
// Make connections
QObject::connect( this->SegmentationNodeComboBox, SIGNAL(currentNodeChanged(vtkMRMLNode*)),
q, SLOT(onSegmentationNodeChanged(vtkMRMLNode*)) );
QObject::connect(this->SliceRotateWarningButton, SIGNAL(clicked()),
q, SLOT(rotateSliceViewsToSegmentation()));
QObject::connect( this->MasterVolumeNodeComboBox, SIGNAL(currentNodeChanged(vtkMRMLNode*)),
q, SLOT(onMasterVolumeNodeChanged(vtkMRMLNode*)) );
QObject::connect( this->SpecifyGeometryButton, SIGNAL(clicked()),
q, SLOT(showSegmentationGeometryDialog()));
QObject::connect( this->SegmentsTableView, SIGNAL(selectionChanged(QItemSelection,QItemSelection)),
q, SLOT(onSegmentSelectionChanged(QItemSelection,QItemSelection)) );
QObject::connect( this->SegmentsTableView, SIGNAL(segmentAboutToBeModified(QString)),
q, SLOT(saveStateForUndo()) );
QObject::connect( this->AddSegmentButton, SIGNAL(clicked()), q, SLOT(onAddSegment()) );
QObject::connect( this->RemoveSegmentButton, SIGNAL(clicked()), q, SLOT(onRemoveSegment()) );
QObject::connect( this->SwitchToSegmentationsButton, SIGNAL(clicked()), q, SLOT(onSwitchToSegmentations()) );
QObject::connect( this->Show3DButton, SIGNAL(toggled(bool)), q, SLOT(onCreateSurfaceToggled(bool)) );
QObject::connect( this->MaskModeComboBox, SIGNAL(currentIndexChanged(int)), q, SLOT(onMaskModeChanged(int)));
QObject::connect( this->MasterVolumeIntensityMaskCheckBox, SIGNAL(toggled(bool)), q, SLOT(onMasterVolumeIntensityMaskChecked(bool)));
QObject::connect( this->MasterVolumeIntensityMaskRangeWidget, SIGNAL(valuesChanged(double,double)), q, SLOT(onMasterVolumeIntensityMaskRangeChanged(double,double)));
QObject::connect( this->OverwriteModeComboBox, SIGNAL(currentIndexChanged(int)), q, SLOT(onOverwriteModeChanged(int)));
QObject::connect( this->UndoButton, SIGNAL(clicked()), q, SLOT(undo()) );
QObject::connect( this->RedoButton, SIGNAL(clicked()), q, SLOT(redo()) );
q->qvtkConnect(this->SegmentationHistory, vtkCommand::ModifiedEvent,
q, SLOT(onSegmentationHistoryChanged()));
// Widget properties
this->SegmentsTableView->setSelectionMode(QAbstractItemView::SingleSelection);
this->SegmentsTableView->setHeaderVisible(true);
this->SegmentsTableView->setVisibilityColumnVisible(true);
this->SegmentsTableView->setColorColumnVisible(true);
this->SegmentsTableView->setOpacityColumnVisible(false);
this->AddSegmentButton->setEnabled(false);
this->RemoveSegmentButton->setEnabled(false);
this->Show3DButton->setEnabled(false);
this->SwitchToSegmentationsButton->setEnabled(false);
this->EffectsGroupBox->setEnabled(false);
this->OptionsGroupBox->setEnabled(false);
this->EffectsGroupBox->setLayout(new QGridLayout(this->EffectsGroupBox));
this->EffectButtonGroup.setExclusive(true);
QObject::connect(&this->EffectButtonGroup, SIGNAL(buttonClicked(QAbstractButton*)), q, SLOT(onEffectButtonClicked(QAbstractButton*) ) );
// Create layout for effect options
QVBoxLayout* layout = new QVBoxLayout(this->EffectsOptionsFrame);
layout->setContentsMargins(0, 0, 0, 0);
layout->setSpacing(0);
// Update effect buttons
q->updateEffectList();
this->OptionsGroupBox->hide();
this->OptionsGroupBox->setTitle("");
this->EffectHelpBrowser->setText("");
this->MaskingGroupBox->hide();
q->updateSliceRotateWarningButtonVisibility();
}
//-----------------------------------------------------------------------------
QToolButton* qMRMLSegmentEditorWidgetPrivate::toolButton(qSlicerSegmentEditorAbstractEffect* effect)
{
QList<QAbstractButton*> effectButtons = this->EffectButtonGroup.buttons();
foreach(QAbstractButton* effectButton, effectButtons)
{
qSlicerSegmentEditorAbstractEffect* foundEffect = qobject_cast<qSlicerSegmentEditorAbstractEffect*>(
effectButton->property("Effect").value<QObject*>());
if (effect == foundEffect)
{
return dynamic_cast<QToolButton*>(effectButton);
}
}
// not found
return nullptr;
}
//-----------------------------------------------------------------------------
void qMRMLSegmentEditorWidgetPrivate::notifyEffectsOfReferenceGeometryChange(const std::string& geometry)
{
if (geometry.compare(this->LastNotifiedReferenceImageGeometry) == 0)
{
// no change
return;
}
this->LastNotifiedReferenceImageGeometry = geometry;
foreach(qSlicerSegmentEditorAbstractEffect* effect, this->RegisteredEffects)
{
effect->referenceGeometryChanged();
}
}
//-----------------------------------------------------------------------------
void qMRMLSegmentEditorWidgetPrivate::notifyEffectsOfMasterVolumeNodeChange()
{
foreach(qSlicerSegmentEditorAbstractEffect* effect, this->RegisteredEffects)
{
effect->masterVolumeNodeChanged();
}
}
//-----------------------------------------------------------------------------
void qMRMLSegmentEditorWidgetPrivate::notifyEffectsOfLayoutChange()
{
foreach(qSlicerSegmentEditorAbstractEffect* effect, this->RegisteredEffects)
{
effect->layoutChanged();
}
}
//-----------------------------------------------------------------------------
bool qMRMLSegmentEditorWidgetPrivate::resetModifierLabelmapToDefault()
{
std::string referenceImageGeometry = this->referenceImageGeometry();
if (referenceImageGeometry.empty())
{
qCritical() << Q_FUNC_INFO << ": Cannot determine default modifierLabelmap geometry";
return false;
}
std::string modifierLabelmapReferenceImageGeometryBaseline = vtkSegmentationConverter::SerializeImageGeometry(this->ModifierLabelmap);
// Set reference geometry to labelmap (origin, spacing, directions, extents) and allocate scalars
vtkNew<vtkMatrix4x4> referenceGeometryMatrix;
int referenceExtent[6] = {0,-1,0,-1,0,-1};
vtkSegmentationConverter::DeserializeImageGeometry(referenceImageGeometry, referenceGeometryMatrix.GetPointer(), referenceExtent);
vtkSegmentationConverter::DeserializeImageGeometry(referenceImageGeometry, this->ModifierLabelmap, true, BINARY_LABELMAP_SCALAR_TYPE, 1);
vtkOrientedImageDataResample::FillImage(this->ModifierLabelmap, BINARY_LABELMAP_VOXEL_EMPTY);
return true;
}
//-----------------------------------------------------------------------------
bool qMRMLSegmentEditorWidgetPrivate::updateSelectedSegmentLabelmap()
{
if (!this->ParameterSetNode)
{
qCritical() << Q_FUNC_INFO << ": Invalid segment editor parameter set node";
return false;
}
vtkMRMLSegmentationNode* segmentationNode = this->ParameterSetNode->GetSegmentationNode();
std::string referenceImageGeometry = this->referenceImageGeometry();
if (!segmentationNode || referenceImageGeometry.empty())
{
return false;
}
const char* selectedSegmentID = this->ParameterSetNode->GetSelectedSegmentID();
if (!selectedSegmentID)
{
qCritical() << Q_FUNC_INFO << ": Invalid segment selection";
return false;
}
vtkSegment* selectedSegment = segmentationNode->GetSegmentation()->GetSegment(selectedSegmentID);
if (selectedSegment == nullptr)
{
qWarning() << Q_FUNC_INFO << " failed: Segment " << selectedSegmentID << " not found in segmentation";
return false;
}
vtkOrientedImageData* segmentLabelmap = vtkOrientedImageData::SafeDownCast(
selectedSegment->GetRepresentation(vtkSegmentationConverter::GetSegmentationBinaryLabelmapRepresentationName()));
if (!segmentLabelmap)
{
qCritical() << Q_FUNC_INFO << ": Failed to get binary labelmap representation in segmentation " << segmentationNode->GetName();
return false;
}
int* extent = segmentLabelmap->GetExtent();
if (extent[0] > extent[1] || extent[2] > extent[3] || extent[4] > extent[5])
{
vtkSegmentationConverter::DeserializeImageGeometry(referenceImageGeometry, this->SelectedSegmentLabelmap, false);
this->SelectedSegmentLabelmap->SetExtent(0, -1, 0, -1, 0, -1);
this->SelectedSegmentLabelmap->AllocateScalars(VTK_UNSIGNED_CHAR, 1);
return true;
}
vtkNew<vtkImageThreshold> threshold;
threshold->SetInputData(segmentLabelmap);
threshold->ThresholdBetween(selectedSegment->GetLabelValue(), selectedSegment->GetLabelValue());
threshold->SetInValue(1);
threshold->SetOutValue(0);
threshold->Update();
vtkSmartPointer<vtkOrientedImageData> thresholdedSegmentLabelmap = vtkSmartPointer<vtkOrientedImageData>::New();
thresholdedSegmentLabelmap->ShallowCopy(threshold->GetOutput());
thresholdedSegmentLabelmap->CopyDirections(segmentLabelmap);
vtkNew<vtkOrientedImageData> referenceImage;
vtkSegmentationConverter::DeserializeImageGeometry(referenceImageGeometry, referenceImage.GetPointer(), false);
vtkOrientedImageDataResample::ResampleOrientedImageToReferenceOrientedImage(
thresholdedSegmentLabelmap, referenceImage.GetPointer(), this->SelectedSegmentLabelmap, /*linearInterpolation=*/false);
return true;
}
//-----------------------------------------------------------------------------
bool qMRMLSegmentEditorWidgetPrivate::updateAlignedMasterVolume()
{
if (!this->ParameterSetNode)
{
qCritical() << Q_FUNC_INFO << ": Invalid segment editor parameter set node";
return false;
}
vtkMRMLSegmentationNode* segmentationNode = this->ParameterSetNode->GetSegmentationNode();
vtkMRMLScalarVolumeNode* masterVolumeNode = this->ParameterSetNode->GetMasterVolumeNode();
std::string referenceImageGeometry = this->referenceImageGeometry();
if (!segmentationNode || !masterVolumeNode || !masterVolumeNode->GetImageData()
|| !masterVolumeNode->GetImageData()->GetPointData() || referenceImageGeometry.empty())
{
return false;
}
vtkNew<vtkOrientedImageData> referenceImage;
vtkSegmentationConverter::DeserializeImageGeometry(referenceImageGeometry, referenceImage.GetPointer(), false);
int* referenceImageExtent = referenceImage->GetExtent();
int* alignedMasterVolumeExtent = this->AlignedMasterVolume->GetExtent();
// If master volume node and transform nodes did not change and the aligned master volume covers the entire reference geometry
// then we don't need to update the aligned master volume.
if (vtkOrientedImageDataResample::DoGeometriesMatch(referenceImage.GetPointer(), this->AlignedMasterVolume)
&& alignedMasterVolumeExtent[0] <= referenceImageExtent[0] && alignedMasterVolumeExtent[1] >= referenceImageExtent[1]
&& alignedMasterVolumeExtent[2] <= referenceImageExtent[2] && alignedMasterVolumeExtent[3] >= referenceImageExtent[3]
&& alignedMasterVolumeExtent[4] <= referenceImageExtent[4] && alignedMasterVolumeExtent[5] >= referenceImageExtent[5]
&& vtkOrientedImageDataResample::DoExtentsMatch(referenceImage.GetPointer(), this->AlignedMasterVolume)
&& this->AlignedMasterVolumeUpdateMasterVolumeNode == masterVolumeNode
&& this->AlignedMasterVolumeUpdateMasterVolumeNodeTransform == masterVolumeNode->GetParentTransformNode()
&& this->AlignedMasterVolumeUpdateSegmentationNodeTransform == segmentationNode->GetParentTransformNode() )
{
// Extents and nodes are matching, check if they have not been modified since the aligned master
// volume generation.
bool updateAlignedMasterVolumeRequired = false;
if (masterVolumeNode->GetMTime() > this->AlignedMasterVolume->GetMTime())
{
updateAlignedMasterVolumeRequired = true;
}
else if (masterVolumeNode->GetParentTransformNode() && masterVolumeNode->GetParentTransformNode()->GetMTime() > this->AlignedMasterVolume->GetMTime())
{
updateAlignedMasterVolumeRequired = true;
}
else if (segmentationNode->GetParentTransformNode() && segmentationNode->GetParentTransformNode()->GetMTime() > this->AlignedMasterVolume->GetMTime())
{
updateAlignedMasterVolumeRequired = true;
}
if (!updateAlignedMasterVolumeRequired)
{
return true;
}
}
// Get a read-only version of masterVolume as a vtkOrientedImageData
vtkNew<vtkOrientedImageData> masterVolume;
if (masterVolumeNode->GetImageData()->GetNumberOfScalarComponents() == 1)
{
masterVolume->vtkImageData::ShallowCopy(masterVolumeNode->GetImageData());
}
else
{
vtkNew<vtkImageExtractComponents> extract;
extract->SetInputData(masterVolumeNode->GetImageData());
extract->Update();
extract->SetComponents(0); // TODO: allow user to configure this
masterVolume->vtkImageData::ShallowCopy(extract->GetOutput());
}
vtkSmartPointer<vtkMatrix4x4> ijkToRasMatrix = vtkSmartPointer<vtkMatrix4x4>::New();
masterVolumeNode->GetIJKToRASMatrix(ijkToRasMatrix);
masterVolume->SetGeometryFromImageToWorldMatrix(ijkToRasMatrix);
vtkNew<vtkGeneralTransform> masterVolumeToSegmentationTransform;
vtkMRMLTransformNode::GetTransformBetweenNodes(masterVolumeNode->GetParentTransformNode(), segmentationNode->GetParentTransformNode(), masterVolumeToSegmentationTransform.GetPointer());
double backgroundValue = masterVolumeNode->GetImageBackgroundScalarComponentAsDouble(0);
vtkOrientedImageDataResample::ResampleOrientedImageToReferenceOrientedImage(masterVolume, referenceImage, this->AlignedMasterVolume,
/*linearInterpolation=*/true, /*padImage=*/false, masterVolumeToSegmentationTransform, backgroundValue);
this->AlignedMasterVolumeUpdateMasterVolumeNode = masterVolumeNode;
this->AlignedMasterVolumeUpdateMasterVolumeNodeTransform = masterVolumeNode->GetParentTransformNode();
this->AlignedMasterVolumeUpdateSegmentationNodeTransform = segmentationNode->GetParentTransformNode();
return true;
}
//-----------------------------------------------------------------------------
bool qMRMLSegmentEditorWidgetPrivate::updateMaskLabelmap()
{
if (!this->ParameterSetNode)
{
qCritical() << Q_FUNC_INFO << ": Invalid segment editor parameter set node";
return false;
}
vtkMRMLSegmentationNode* segmentationNode = this->ParameterSetNode->GetSegmentationNode();
if (!segmentationNode)
{
qCritical() << Q_FUNC_INFO << ": Invalid segmentation node, modifier labelmap, or mask labelmap";
return false;
}
// GenerateEditMask can add intensity range based mask, too. We do not use it here, as currently
// editable intensity range is taken into account in qSlicerSegmentEditorAbstractEffect::modifySelectedSegmentByLabelmap.
// It would simplify implementation if we passed master volume and intensity range to GenerateEditMask here
// and removed intensity range based masking from modifySelectedSegmentByLabelmap.
if (!segmentationNode->GenerateEditMask(this->MaskLabelmap, this->ParameterSetNode->GetMaskMode(), this->ModifierLabelmap,
this->ParameterSetNode->GetSelectedSegmentID() ? this->ParameterSetNode->GetSelectedSegmentID() : "",
this->ParameterSetNode->GetMaskSegmentID() ? this->ParameterSetNode->GetMaskSegmentID() : ""))
{
qCritical() << Q_FUNC_INFO << ": Mask generation failed";
return false;
}
return true;
}
//-----------------------------------------------------------------------------
bool qMRMLSegmentEditorWidgetPrivate::updateReferenceGeometryImage()
{
std::string geometry = this->referenceImageGeometry();
if (geometry.empty())
{
return false;
}
return vtkSegmentationConverter::DeserializeImageGeometry(geometry, this->ReferenceGeometryImage, false /* don't allocate scalars */);
}
//-----------------------------------------------------------------------------
void qMRMLSegmentEditorWidgetPrivate::selectFirstSegment()
{
if (!this->ParameterSetNode)
{
qCritical() << Q_FUNC_INFO << ": Invalid segment editor parameter set node";
return;
}
vtkMRMLSegmentationNode* segmentationNode = this->ParameterSetNode->GetSegmentationNode();
if ( segmentationNode
&& segmentationNode->GetSegmentation()->GetNumberOfSegments() > 0 )
{
std::vector<std::string> segmentIDs;
segmentationNode->GetSegmentation()->GetSegmentIDs(segmentIDs);
QStringList firstSegmentID;
firstSegmentID << QString(segmentIDs[0].c_str());
this->SegmentsTableView->setSelectedSegmentIDs(firstSegmentID);
}
}
//-----------------------------------------------------------------------------
void qMRMLSegmentEditorWidgetPrivate::setEffectCursor(qSlicerSegmentEditorAbstractEffect* effect)
{
qSlicerLayoutManager* layoutManager = qSlicerApplication::application()->layoutManager();
if (!layoutManager)
{
// application is closing
return;
}
// We update the default cursor as well so that if the user hovers the mouse over
// a markup, the cursor shape is then restored to the effect cursor.
foreach(QString sliceViewName, layoutManager->sliceViewNames())
{
qMRMLSliceWidget* sliceWidget = layoutManager->sliceWidget(sliceViewName);
QString viewNodeID = QString::fromStdString(sliceWidget->mrmlSliceNode()->GetID());
if (!this->segmentationDisplayableInView(sliceWidget->mrmlSliceNode()))
{
// segmentation is not displayable in this view anymore
if (!this->CustomCursorInViewNodeIDs.contains(viewNodeID))
{
// we did not use this view previously either, so don't change the cursor there
continue;
}
}
if (effect && effect->showEffectCursorInSliceView())
{
sliceWidget->sliceView()->setViewCursor(effect->createCursor(sliceWidget));
sliceWidget->sliceView()->setDefaultViewCursor(effect->createCursor(sliceWidget));
this->CustomCursorInViewNodeIDs.insert(viewNodeID);
}
else
{
sliceWidget->sliceView()->setViewCursor(QCursor());
sliceWidget->sliceView()->setDefaultViewCursor(QCursor());
this->CustomCursorInViewNodeIDs.remove(viewNodeID);
}
}
for (int threeDViewId = 0; threeDViewId < layoutManager->threeDViewCount(); ++threeDViewId)
{
qMRMLThreeDWidget* threeDWidget = layoutManager->threeDWidget(threeDViewId);
QString viewNodeID = QString::fromStdString(threeDWidget->mrmlViewNode()->GetID());
if (!this->segmentationDisplayableInView(threeDWidget->mrmlViewNode()))
{
// segmentation is not displayable in this view anymore
if (!this->CustomCursorInViewNodeIDs.contains(viewNodeID))
{
// we did not use this view previously either, so don't change the cursor there
continue;
}
}
if (effect && effect->showEffectCursorInThreeDView())
{
threeDWidget->threeDView()->setViewCursor(effect->createCursor(threeDWidget));
threeDWidget->threeDView()->setDefaultViewCursor(effect->createCursor(threeDWidget));
this->CustomCursorInViewNodeIDs.insert(viewNodeID);
}
else
{
threeDWidget->threeDView()->setViewCursor(QCursor());
threeDWidget->threeDView()->setDefaultViewCursor(QCursor());
this->CustomCursorInViewNodeIDs.remove(viewNodeID);
}
}
}
//-----------------------------------------------------------------------------
std::string qMRMLSegmentEditorWidgetPrivate::getReferenceImageGeometryFromSegmentation(vtkSegmentation* segmentation)
{
if (!segmentation)
{
return "";
}
// If "reference image geometry" conversion parameter is set then use that
std::string referenceImageGeometry = segmentation->GetConversionParameter(vtkSegmentationConverter::GetReferenceImageGeometryParameterName());
if (!referenceImageGeometry.empty())
{
// Extend reference image geometry to contain all segments (needed for example for properly handling imported segments
// that do not fit into the reference image geometry)
vtkSmartPointer<vtkOrientedImageData> commonGeometryImage = vtkSmartPointer<vtkOrientedImageData>::New();
vtkSegmentationConverter::DeserializeImageGeometry(referenceImageGeometry, commonGeometryImage, false);
// Determine extent that contains all segments
int commonSegmentExtent[6] = { 0, -1, 0, -1, 0, -1 };
segmentation->DetermineCommonLabelmapExtent(commonSegmentExtent, commonGeometryImage);
if (commonSegmentExtent[0] <= commonSegmentExtent[1]
&& commonSegmentExtent[2] <= commonSegmentExtent[3]
&& commonSegmentExtent[4] <= commonSegmentExtent[5])
{
// Expand commonGeometryExtent as needed to contain commonSegmentExtent
int commonGeometryExtent[6] = { 0, -1, 0, -1, 0, -1 };
commonGeometryImage->GetExtent(commonGeometryExtent);
for (int i = 0; i < 3; i++)
{
commonGeometryExtent[i * 2] = std::min(commonSegmentExtent[i * 2], commonGeometryExtent[i * 2]);
commonGeometryExtent[i * 2 + 1] = std::max(commonSegmentExtent[i * 2 + 1], commonGeometryExtent[i * 2 + 1]);
}
commonGeometryImage->SetExtent(commonGeometryExtent);
referenceImageGeometry = vtkSegmentationConverter::SerializeImageGeometry(commonGeometryImage);
}
// TODO: Use oversampling (if it's 'A' then ignore and changed to 1)
return referenceImageGeometry;
}
if (segmentation->ContainsRepresentation(vtkSegmentationConverter::GetSegmentationBinaryLabelmapRepresentationName()))
{
// If no reference image geometry is specified but there are labels already then determine geometry from that
referenceImageGeometry = segmentation->DetermineCommonLabelmapGeometry();
return referenceImageGeometry;
}
return "";
}
//-----------------------------------------------------------------------------
std::string qMRMLSegmentEditorWidgetPrivate::referenceImageGeometry()
{
if (!this->ParameterSetNode)
{
qCritical() << Q_FUNC_INFO << ": Invalid segment editor parameter set node";
notifyEffectsOfReferenceGeometryChange("");
return "";
}
vtkMRMLSegmentationNode* segmentationNode = this->ParameterSetNode->GetSegmentationNode();
vtkSegmentation* segmentation = segmentationNode ? segmentationNode->GetSegmentation() : nullptr;
if (!segmentationNode || !segmentation)
{
qCritical() << Q_FUNC_INFO << ": Invalid segmentation";
notifyEffectsOfReferenceGeometryChange("");
return "";
}
std::string referenceImageGeometry;
referenceImageGeometry = this->getReferenceImageGeometryFromSegmentation(segmentation);
if (referenceImageGeometry.empty())
{
// If no reference image geometry could be determined then use the master volume's geometry
vtkMRMLScalarVolumeNode* masterVolumeNode = this->ParameterSetNode->GetMasterVolumeNode();
if (!masterVolumeNode)
{
// cannot determine reference geometry
return "";
}
// Update the referenceImageGeometry parameter for next time
segmentationNode->SetReferenceImageGeometryParameterFromVolumeNode(masterVolumeNode);
// Update extents to include all existing segments
referenceImageGeometry = this->getReferenceImageGeometryFromSegmentation(segmentation);
}
notifyEffectsOfReferenceGeometryChange(referenceImageGeometry);
return referenceImageGeometry;
}
//-----------------------------------------------------------------------------
bool qMRMLSegmentEditorWidgetPrivate::segmentationDisplayableInView(vtkMRMLAbstractViewNode* viewNode)
{
if (!viewNode)
{
qWarning() << Q_FUNC_INFO << ": failed. Invalid viewNode.";
return false;
}
if (!this->ParameterSetNode)
{
return false;
}
vtkMRMLSegmentationNode* segmentationNode = this->ParameterSetNode->GetSegmentationNode();
if (!segmentationNode)
{
return false;
}
const char* viewNodeID = viewNode->GetID();
int numberOfDisplayNodes = segmentationNode->GetNumberOfDisplayNodes();
for (int displayNodeIndex = 0; displayNodeIndex < numberOfDisplayNodes; displayNodeIndex++)
{
vtkMRMLDisplayNode* segmentationDisplayNode = segmentationNode->GetNthDisplayNode(displayNodeIndex);
if (segmentationDisplayNode && segmentationDisplayNode->IsDisplayableInView(viewNodeID))
{
return true;
}
}
return false;
}
//-----------------------------------------------------------------------------
bool qMRMLSegmentEditorWidgetPrivate::segmentationIJKToRAS(vtkMatrix4x4* ijkToRas)
{
if (!this->ParameterSetNode)
{
return false;
}
if (!this->updateReferenceGeometryImage())
{
return false;
}
vtkMRMLSegmentationNode* segmentationNode = this->ParameterSetNode->GetSegmentationNode();
if (!segmentationNode || !segmentationNode->GetSegmentation())