void mitk::SegmentationInterpolationController::SetChangedSlice( const Image* sliceDiff, unsigned int sliceDimension, unsigned int sliceIndex, unsigned int timeStep )
{
if ( !sliceDiff ) return;
if ( sliceDimension > 2 ) return;
if ( timeStep >= m_SegmentationCountInSlice.size() ) return;
if ( sliceIndex >= m_SegmentationCountInSlice[timeStep][sliceDimension].size() ) return;
unsigned int dim0(0);
unsigned int dim1(1);
// determine the other two dimensions
switch (sliceDimension)
{
default:
case 2:
dim0 = 0; dim1 = 1; break;
case 1:
dim0 = 0; dim1 = 2; break;
case 0:
dim0 = 1; dim1 = 2; break;
}
//mitkIpPicDescriptor* rawSlice = const_cast<Image*>(sliceDiff)->GetSliceData()->GetPicDescriptor(); // we promise not to change anything!
unsigned char* rawSlice = (unsigned char*) const_cast<Image*>(sliceDiff)->GetData();
if (!rawSlice) return;
AccessFixedDimensionByItk_1( sliceDiff, ScanChangedSlice, 2, SetChangedSliceOptions(sliceDimension, sliceIndex, dim0, dim1, timeStep, rawSlice) );
//PrintStatus();
Modified();
}
void mitk::SegmentationInterpolationController::ScanWholeVolume( itk::Image<DATATYPE, 3>*, const Image* volume, unsigned int timeStep )
{
if (!volume) return;
if ( timeStep >= m_SegmentationCountInSlice.size() ) return;
for (unsigned int slice = 0; slice < volume->GetDimension(2); ++slice)
{
DATATYPE* rawVolume = static_cast<DATATYPE*>( const_cast<Image*>(volume)->GetVolumeData(timeStep)->GetData() ); // we again promise not to change anything, we'll just count
//DATATYPE* rawSlice = static_cast<DATATYPE*>( volume->GetSliceData(slice)->GetData() ); // TODO THIS wouldn't work. Did I mess up with some internal mitk::Image data structure?
DATATYPE* rawSlice = rawVolume + ( volume->GetDimension(0) * volume->GetDimension(1) * slice );
ScanChangedSlice<DATATYPE>( NULL, SetChangedSliceOptions(2, slice, 0, 1, timeStep, rawSlice) );
}
}
void mitk::SegmentationInterpolationController::ScanWholeVolume( const itk::Image<DATATYPE, 3>*, const Image* volume, unsigned int timeStep )
{
if (!volume) return;
if ( timeStep >= m_SegmentationCountInSlice.size() ) return;
ImageReadAccessor readAccess(volume, volume->GetVolumeData(timeStep));
for (unsigned int slice = 0; slice < volume->GetDimension(2); ++slice)
{
const DATATYPE* rawVolume = static_cast<const DATATYPE*>( readAccess.GetData() ); // we again promise not to change anything, we'll just count
const DATATYPE* rawSlice = rawVolume + ( volume->GetDimension(0) * volume->GetDimension(1) * slice );
ScanChangedSlice<DATATYPE>( nullptr, SetChangedSliceOptions(2, slice, 0, 1, timeStep, rawSlice) );
}
}
void mitk::SliceBasedInterpolationController::SetChangedSlice(const Image *slice,
unsigned int sliceDimension,
unsigned int sliceIndex,
unsigned int timeStep)
{
if (!slice)
return;
if (slice->GetDimension() != 2)
return;
if (sliceDimension > 2)
return;
if (m_WorkingImage.IsNull())
return;
// check if the number of labels has changed
int numberOfLabels = m_WorkingImage->GetNumberOfLabels();
if (m_LabelCountInSlice[0][0][0].size() != numberOfLabels)
return;
unsigned int dim0(0);
unsigned int dim1(1);
// determine the other two dimensions
switch (sliceDimension)
{
default:
case 2:
dim0 = 0;
dim1 = 1;
break;
case 1:
dim0 = 0;
dim1 = 2;
break;
case 0:
dim0 = 1;
dim1 = 2;
break;
}
AccessFixedDimensionByItk_1(
slice, ScanSliceITKProcessing, 2, SetChangedSliceOptions(sliceDimension, sliceIndex, dim0, dim1, timeStep));
// this->Modified();
}
