76 lines
2.9 KiB
C
76 lines
2.9 KiB
C
|
//
|
||
|
|
// FastDTW.h
|
||
|
|
// FastDTW-x
|
||
|
|
//
|
||
|
|
// Created by Melo Yao on 12/9/13.
|
||
|
|
// Copyright (c) 2013 melo.yao. All rights reserved.
|
||
|
|
//
|
||
|
|
|
||
|
|
#ifndef __FastDTW_x__FastDTW__
|
||
|
|
#define __FastDTW_x__FastDTW__
|
||
|
|
|
||
|
|
#include "Foundation.h"
|
||
|
|
#include "DTW.h"
|
||
|
|
#include "PAA.h"
|
||
|
|
#include "ExpandedResWindow.h"
|
||
|
|
|
||
|
|
//#include "TimeWarpInfo.h"
|
||
|
|
//#include "TimeSeries.h"
|
||
|
|
|
||
|
|
FD_NS_START
|
||
|
|
|
||
|
|
namespace FAST {
|
||
|
|
|
||
|
|
extern const JInt DEFAULT_SEARCH_RADIUS;
|
||
|
|
|
||
|
|
template <typename ValueType,JInt nDimension, typename DistanceFunction>
|
||
|
|
TimeWarpInfo<ValueType> getWarpInfoBetween(TimeSeries<ValueType,nDimension> const& tsI, TimeSeries<ValueType,nDimension> const& tsJ, JInt searchRadius, DistanceFunction const& distFn)
|
||
|
|
{
|
||
|
|
if (searchRadius < 0) {
|
||
|
|
searchRadius = 0;
|
||
|
|
}
|
||
|
|
JInt minTSsize = searchRadius + 2;
|
||
|
|
if (tsI.size() <= minTSsize || tsJ.size()<=minTSsize) {
|
||
|
|
return STRI::getWarpInfoBetween(tsI, tsJ, distFn);
|
||
|
|
}
|
||
|
|
else
|
||
|
|
{
|
||
|
|
JDouble resolutionFactor = 2.0;
|
||
|
|
PAA<ValueType,nDimension> shrunkI(tsI,(JInt)(tsI.size()/resolutionFactor));
|
||
|
|
PAA<ValueType,nDimension> shrunkJ(tsJ,(JInt)(tsJ.size()/resolutionFactor));
|
||
|
|
// Determine the search window that constrains the area of the cost matrix that will be evaluated based on
|
||
|
|
// the warp path found at the previous resolution (smaller time series).
|
||
|
|
TimeWarpInfo<ValueType> warpInfo = getWarpInfoBetween(shrunkI, shrunkJ, searchRadius, distFn);
|
||
|
|
ExpandedResWindow window(tsI, tsJ, shrunkI, shrunkJ,
|
||
|
|
*(warpInfo.getPath()),
|
||
|
|
searchRadius);
|
||
|
|
return STRI::getWarpInfoBetween(tsI, tsJ, window, distFn);
|
||
|
|
}
|
||
|
|
|
||
|
|
}
|
||
|
|
|
||
|
|
template <typename ValueType,JInt nDimension, typename DistanceFunction>
|
||
|
|
inline ValueType getWarpDistBetween(TimeSeries<ValueType,nDimension> const& tsI,TimeSeries<ValueType,nDimension> const& tsJ, DistanceFunction const& distFn)
|
||
|
|
{
|
||
|
|
return getWarpInfoBetween(tsI, tsJ, DEFAULT_SEARCH_RADIUS, distFn).getDistance();
|
||
|
|
}
|
||
|
|
|
||
|
|
template <typename ValueType,JInt nDimension, typename DistanceFunction>
|
||
|
|
inline TimeWarpInfo<ValueType> getWarpInfoBetween(TimeSeries<ValueType,nDimension> const& tsI,TimeSeries<ValueType,nDimension> const& tsJ,DistanceFunction const& distFn)
|
||
|
|
{
|
||
|
|
return getWarpInfoBetween(tsI, tsJ, DEFAULT_SEARCH_RADIUS, distFn);
|
||
|
|
}
|
||
|
|
|
||
|
|
template <typename ValueType,JInt nDimension, typename DistanceFunction>
|
||
|
|
inline ValueType getWarpDistBetween(TimeSeries<ValueType,nDimension> const& tsI,TimeSeries<ValueType,nDimension> const& tsJ,
|
||
|
|
JInt searchRadius,DistanceFunction const& distFn)
|
||
|
|
{
|
||
|
|
return getWarpInfoBetween(tsI, tsJ, searchRadius, distFn).getDistance();
|
||
|
|
}
|
||
|
|
|
||
|
|
|
||
|
|
|
||
|
|
}
|
||
|
|
FD_NS_END
|
||
|
|
#endif /* defined(__FastDTW_x__FastDTW__) */
|