mitkTimeHelper.h
Go to the documentation of this file.
00001 /*========================================================================= 00002 00003 Program: Medical Imaging & Interaction Toolkit 00004 Language: C++ 00005 Date: $Date$ 00006 Version: $Revision$ 00007 00008 Copyright (c) German Cancer Research Center, Division of Medical and 00009 Biological Informatics. All rights reserved. 00010 See MITKCopyright.txt or https://www.mitk.org/copyright.html for details. 00011 00012 This software is distributed WITHOUT ANY WARRANTY; without even 00013 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 00014 PURPOSE. See the above copyright notices for more information. 00015 00016 =========================================================================*/ 00017 00018 00019 #ifndef MITKTIMEHELPER_H_HEADER_INCLUDED_C1C2FCD2 00020 #define MITKTIMEHELPER_H_HEADER_INCLUDED_C1C2FCD2 00021 00022 namespace mitk 00023 { 00024 00025 //## @brief convert the start- and end-index-time of output-region in 00026 //## start- and end-index-time of input-region via millisecond-time 00027 template <class TOutputRegion, class TInputRegion> 00028 void ITK_EXPORT GenerateTimeInInputRegion(const mitk::TimeSlicedGeometry *outputTimeGeometry, const TOutputRegion& outputRegion, const mitk::TimeSlicedGeometry *inputTimeGeometry, TInputRegion& inputRegion) 00029 { 00030 assert(outputTimeGeometry!=NULL); 00031 assert(inputTimeGeometry!=NULL); 00032 00033 // convert the start-index-time of output in start-index-time of input via millisecond-time 00034 ScalarType timeInMS = outputTimeGeometry->TimeStepToMS(outputRegion.GetIndex(3)); 00035 int timestep = inputTimeGeometry->MSToTimeStep( timeInMS ); 00036 if( ( timeInMS > ScalarTypeNumericTraits::NonpositiveMin() ) && ( inputTimeGeometry->IsValidTime( timestep ) ) ) 00037 inputRegion.SetIndex( 3, timestep ); 00038 else 00039 inputRegion.SetIndex( 3, 0 ); 00040 // convert the end-index-time of output in end-index-time of input via millisecond-time 00041 timeInMS = outputTimeGeometry->TimeStepToMS(outputRegion.GetIndex(3)+outputRegion.GetSize(3)-1); 00042 timestep = inputTimeGeometry->MSToTimeStep( timeInMS ); 00043 if( ( timeInMS > ScalarTypeNumericTraits::NonpositiveMin() ) && ( outputTimeGeometry->IsValidTime( timestep ) ) ) 00044 inputRegion.SetSize( 3, timestep - inputRegion.GetIndex(3) + 1 ); 00045 else 00046 inputRegion.SetSize( 3, 1 ); 00047 } 00048 00049 //##Documentation 00050 //## @brief convert the start- and end-index-time of output in 00051 //## start- and end-index-time of input1 and input2 via millisecond-time 00052 template <class TOutputData, class TInputData> 00053 void ITK_EXPORT GenerateTimeInInputRegion(const TOutputData* output, TInputData* input) 00054 { 00055 assert(output!=NULL); 00056 assert(input!=NULL); 00057 00058 const typename TOutputData::RegionType& outputRegion = output->GetRequestedRegion(); 00059 typename TInputData::RegionType inputRegion; 00060 00061 if(outputRegion.GetSize(3)<1) 00062 { 00063 typename TInputData::RegionType::SizeType inputsize; 00064 inputsize.Fill(0); 00065 inputRegion.SetSize(inputsize); 00066 input->SetRequestedRegion( &inputRegion ); 00067 } 00068 00069 // convert the start-index-time of output in start-index-time of input via millisecond-time 00070 inputRegion = input->GetRequestedRegion(); 00071 GenerateTimeInInputRegion(output->GetTimeSlicedGeometry(), outputRegion, input->GetTimeSlicedGeometry(), inputRegion); 00072 input->SetRequestedRegion( &inputRegion ); 00073 } 00074 00075 } // end namespace mitk 00076 00077 //#ifndef ITK_MANUAL_INSTANTIATION 00078 //#include "mitkTimeHelper.txx" 00079 #include "MitkExtExports.h" 00080 //#endif 00081 00082 #endif // MITKTIMEHELPER_H_HEADER_INCLUDED_C1C2FCD2
