source: addons/HeuristicLab.DataImporter/HeuristicLab.DataImporter.Command/TimeSeries/InterpolateMissingValues.cs @ 16567

#region License Information
/* HeuristicLab
 * Copyright (C) 2002-2013 Heuristic and Evolutionary Algorithms Laboratory (HEAL)
 *
 * This file is part of HeuristicLab.
 *
 * HeuristicLab is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * HeuristicLab is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with HeuristicLab. If not, see <http://www.gnu.org/licenses/>.
 */
#endregion

using System;
using System.Collections.Generic;
using System.Linq;
using System.Windows.Forms;
using HeuristicLab.DataImporter.Data;
using HeuristicLab.DataImporter.Data.CommandBase;
using HeuristicLab.DataImporter.Data.Model;
using HeuristicLab.Persistence.Default.CompositeSerializers.Storable;
using HEAL.Attic;

namespace HeuristicLab.DataImporter.Command {
  [StorableType("3975A950-D715-48C8-91FB-C982DA6887CA")]
  [ViewableCommandInfoAttribute("Interpolate Missing Values", 1, ColumnGroupState.Sorted | ColumnGroupState.DoubleColumnSelected, "Time Series",
Position = 2)]
  public class InterpolateMissingValues : FillMissingValueCommandBase {
    private DateTimeColumn columnToSample;

    [StorableConstructor]
    protected InterpolateMissingValues(StorableConstructorFlag _) : base(_) { }
    public InterpolateMissingValues(DataSet dataSet, string columnGroupName, int[] affectedColumns)
      : base(dataSet, columnGroupName, affectedColumns) {
    }

    public override string Description {
      get { return "Interpolate Missing Values"; }
    }

    public override void Execute() {
      base.Execute();
      if (!(ColumnGroup.GetColumn(ColumnGroup.SortedColumnIndexes.ElementAt(0)) is DateTimeColumn)
   || ColumnGroup.SortedColumnIndexes.Count != 1)
        throw new CommandExecutionException("ColumnGroup must be sorted after a datetime column.", this);
      columnToSample = (DateTimeColumn)ColumnGroup.Columns.ElementAt(ColumnGroup.SortedColumnIndexes.ElementAt(0));
      if (columnToSample.SortOrder != SortOrder.Ascending)
        throw new CommandExecutionException("ColumnGroup must be sorted ascending by the datetime column.", this);
      if (columnToSample.ContainsNullValues)
        throw new CommandExecutionException("The datetime column must not contain null values.", this);

      oldSortOrder = ColumnGroup.SortOrdersForColumns.ToList();
      oldSortedColumnIndexes = new List<int>(ColumnGroup.SortedColumnIndexes);
      DoubleColumn column;
      int startIndex = -1;
      int endIndex = -1;
      foreach (int col in AffectedColumns) {
        if (ColumnGroup.GetColumn(col) is DoubleColumn && ColumnGroup.GetColumn(col).ContainsNullValues) {
          column = (DoubleColumn)ColumnGroup.Columns.ElementAt(col);
          nullValues.Add(col, new List<int>());
          startIndex = FindNextNotNullValueIndex(column, startIndex);
          endIndex = FindNextNotNullValueIndex(column, startIndex);
          while (endIndex != -1) {
            if (endIndex - startIndex != 1) {
              nullValues[col].AddRange(Enumerable.Range(startIndex + 1, endIndex - startIndex - 1));
              InterpolateValues(column, startIndex, endIndex);
            }
            startIndex = endIndex;
            endIndex = FindNextNotNullValueIndex(column, endIndex);
          }
          startIndex = -1;
        }
      }
      ColumnGroup.FireChanged();
      columnToSample = null;
      this.ColumnGroup = null;
    }

    private int FindNextNotNullValueIndex(DoubleColumn column, int startIndex) {
      int endIndex = -1;
      for (int i = startIndex + 1; i < column.TotalValuesCount && endIndex == -1; i++)
        if (column.GetValue(i) != null)
          endIndex = i;
      return endIndex;
    }

    private void InterpolateValues(DoubleColumn column, int startIndex, int endIndex) {
      double startValue = (double)column.GetValue(startIndex);
      double endValue = (double)column.GetValue(endIndex);
      DateTime start = (DateTime)columnToSample.GetValue(startIndex);
      DateTime end = (DateTime)columnToSample.GetValue(endIndex);

      for (int i = startIndex + 1; i < endIndex; i++) {
        if (column.GetValue(i) != null)
          throw new CommandExecutionException("Tried to manipulate non null value.", this);
        column.ChangeValue(i, ((double)(((DateTime)columnToSample.GetValue(i) - start).TotalSeconds / (end - start).TotalSeconds)) *
          (endValue - startValue) + startValue);
      }
    }
  }
}