Skip to content
GitLab
Select Git revision
  • zonalstats_histogram default
Find file BlameHistoryPermalink
DimensionalityReductionModel.txx 4.14 KiB
1 
/*=========================================================================
2
3 
  Program:   ORFEO Toolbox
4 
  Language:  C++
5 
  Date:      $Date$
6 
  Version:   $Revision$
7
8
9 
  Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
10 
  See OTBCopyright.txt for details.
11
12
13 
     This software is distributed WITHOUT ANY WARRANTY; without even
14 
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
15 
     PURPOSE.  See the above copyright notices for more information.
16
17 
=========================================================================*/
18 
#ifndef DimensionalityReductionModel_txx
19 
#define DimensionalityReductionModel_txx
20
21 
#ifdef _OPENMP
22 
 # include <omp.h>
23 
#endif
24
25 
#include "DimensionalityReductionModel.h"
26
27 
#include "itkMultiThreader.h"
28
29 
namespace otb
30 
{
31
32 
template <class TInputValue, class TOutputValue>
33 
DimensionalityReductionModel<TInputValue,TOutputValue>
34 
::DimensionalityReductionModel()
35 
{ this->m_IsDoPredictBatchMultiThreaded=false;}
36
37
38 
template <class TInputValue, class TOutputValue>
39 
DimensionalityReductionModel<TInputValue,TOutputValue>
40 
::~DimensionalityReductionModel()
41 
{}
42
43
44
45 
template <class TInputValue, class TOutputValue>
46 
typename DimensionalityReductionModel<TInputValue,TOutputValue>
47 
::TargetSampleType
48 
DimensionalityReductionModel<TInputValue,TOutputValue>
49 
::Predict(const InputSampleType& input) const
50 
{
51 
  // Call protected specialization entry point
52 
  return this->DoPredict(input);
53 
}
54
55
56 
template <class TInputValue, class TOutputValue>
57 
typename DimensionalityReductionModel<TInputValue,TOutputValue>
58 
::TargetListSampleType::Pointer
59 
DimensionalityReductionModel<TInputValue,TOutputValue>
60 
::PredictBatch(const InputListSampleType * input) const
61 
{
62 
  typename TargetListSampleType::Pointer targets = TargetListSampleType::New();
63 
  targets->Resize(input->Size());
64
65
66 
  if(m_IsDoPredictBatchMultiThreaded)
67 
    {
68 
    // Simply calls DoPredictBatch
69 
    this->DoPredictBatch(input,0,input->Size(),targets);
70 
    return targets;
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
} else { #ifdef _OPENMP // OpenMP threading here unsigned int nb_threads(0), threadId(0), nb_batches(0); #pragma omp parallel shared(nb_threads,nb_batches) private(threadId) { // Get number of threads configured with ITK omp_set_num_threads(itk::MultiThreader::GetGlobalDefaultNumberOfThreads()); nb_threads = omp_get_num_threads(); threadId = omp_get_thread_num(); nb_batches = std::min(nb_threads,(unsigned int)input->Size()); // Ensure that we do not spawn unncessary threads if(threadId<nb_batches) { unsigned int batch_size = ((unsigned int)input->Size()/nb_batches); unsigned int batch_start = threadId*batch_size; if(threadId == nb_threads-1) { batch_size+=input->Size()%nb_batches; } this->DoPredictBatch(input,batch_start,batch_size,targets); } } #else this->DoPredictBatch(input,0,input->Size(),targets); #endif return targets; } } template <class TInputValue, class TOutputValue> void DimensionalityReductionModel<TInputValue,TOutputValue> ::DoPredictBatch(const InputListSampleType * input, const unsigned int & startIndex, const unsigned int & size, TargetListSampleType * targets) const { assert(input != ITK_NULLPTR); assert(targets != ITK_NULLPTR); assert(input->Size()==targets->Size()&&"Input sample list and target label list do not have the same size."); assert(((quality==ITK_NULLPTR)||(quality->Size()==input->Size()))&&"Quality samples list is not null and does not have the same size as input samples list"); if(startIndex+size>input->Size()) { itkExceptionMacro(<<"requested range ["<<startIndex<<", "<<startIndex+size<<"[ partially outside input sample list range.[0,"<<input->Size()<<"["); } for(unsigned int id = startIndex;id<startIndex+size;++id) { const TargetSampleType target = this->DoPredict(input->GetMeasurementVector(id)); targets->SetMeasurementVector(id,target); } } template <class TInputValue, class TOutputValue> void DimensionalityReductionModel<TInputValue,TOutputValue> ::PrintSelf(std::ostream& os, itk::Indent indent) const { // Call superclass implementation Superclass::PrintSelf(os,indent); }
141142143144
} #endif