#9773 Add method for finding min and max time steps.

OPM · Apr 19, 2023 · aaa45d9 · aaa45d9
1 parent 88042db
commit aaa45d9
Show file tree

Hide file tree

Showing 2 changed files with 62 additions and 11 deletions.
diff --git a/ApplicationLibCode/ProjectDataModel/Summary/RimEnsembleStatisticsCase.cpp b/ApplicationLibCode/ProjectDataModel/Summary/RimEnsembleStatisticsCase.cpp
@@ -137,9 +137,9 @@ void RimEnsembleStatisticsCase::calculate( const std::vector<RimSummaryCase*>& s
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------
-void RimEnsembleStatisticsCase::calculate( const std::vector<RimSummaryCase*> sumCases,
-                                           const RifEclipseSummaryAddress&    inputAddress,
-                                           bool                               includeIncompleteCurves )
+void RimEnsembleStatisticsCase::calculate( const std::vector<RimSummaryCase*>& sumCases,
+                                           const RifEclipseSummaryAddress&     inputAddress,
+                                           bool                                includeIncompleteCurves )
 {
     std::vector<time_t>              allTimeSteps;
     std::vector<std::vector<double>> caseAndTimeStepValues;
@@ -220,9 +220,9 @@ void RimEnsembleStatisticsCase::clearData()
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------
-std::vector<RimSummaryCase*> RimEnsembleStatisticsCase::validSummaryCases( const std::vector<RimSummaryCase*> allSumCases,
-                                                                           const RifEclipseSummaryAddress&    inputAddress,
-                                                                           bool                               includeIncompleteCurves )
+std::vector<RimSummaryCase*> RimEnsembleStatisticsCase::validSummaryCases( const std::vector<RimSummaryCase*>& allSumCases,
+                                                                           const RifEclipseSummaryAddress&     inputAddress,
+                                                                           bool                                includeIncompleteCurves )
 {
     std::vector<RimSummaryCase*>                    validCases;
     std::vector<std::pair<RimSummaryCase*, time_t>> times;
@@ -259,5 +259,54 @@ std::vector<RimSummaryCase*> RimEnsembleStatisticsCase::validSummaryCases( const
 
         validCases.push_back( sumCase );
     }
+
     return validCases;
 }
+
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+std::pair<time_t, time_t> RimEnsembleStatisticsCase::findMinMaxTimeStep( const std::vector<RimSummaryCase*>& sumCases,
+                                                                         const RifEclipseSummaryAddress&     inputAddress )
+{
+    time_t minTimeStep = std::numeric_limits<time_t>::max();
+    time_t maxTimeStep = 0;
+
+    for ( const auto& sumCase : sumCases )
+    {
+        const auto& reader = sumCase->summaryReader();
+        if ( reader )
+        {
+            const std::vector<time_t>& timeSteps = reader->timeSteps( inputAddress );
+            if ( !timeSteps.empty() )
+            {
+                if ( timeSteps.front() < minTimeStep ) minTimeStep = timeSteps.front();
+                if ( timeSteps.back() > maxTimeStep ) maxTimeStep = timeSteps.back();
+            }
+        }
+    }
+
+    return std::make_pair( minTimeStep, maxTimeStep );
+}
+
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+RiaDefines::DateTimePeriod RimEnsembleStatisticsCase::findBestResamplingPeriod( time_t minTimeStep, time_t maxTimeStep )
+{
+    std::vector<RiaDefines::DateTimePeriod> periods = { RiaDefines::DateTimePeriod::DAY,
+                                                        RiaDefines::DateTimePeriod::HOUR,
+                                                        RiaDefines::DateTimePeriod::MINUTE };
+
+    for ( auto p : periods )
+    {
+        size_t numSamples = RiaTimeHistoryCurveResampler::timeStepsFromTimeRange( p, minTimeStep, maxTimeStep ).size();
+        // Resampled data should ideally have at least 100 samples to look good.
+        if ( numSamples > 100 )
+        {
+            return p;
+        }
+    }
+
+    return RiaDefines::DateTimePeriod::DAY;
+}
diff --git a/ApplicationLibCode/ProjectDataModel/Summary/RimEnsembleStatisticsCase.h b/ApplicationLibCode/ProjectDataModel/Summary/RimEnsembleStatisticsCase.h
@@ -55,11 +55,13 @@ class RimEnsembleStatisticsCase : public RimSummaryCase
     RiaDefines::EclipseUnitSystem unitSystem() const;
 
 private:
-    void calculate( const std::vector<RimSummaryCase*> sumCases, const RifEclipseSummaryAddress& inputAddress, bool includeIncompleteCurves );
-    void                         clearData();
-    std::vector<RimSummaryCase*> validSummaryCases( const std::vector<RimSummaryCase*> allSumCases,
-                                                    const RifEclipseSummaryAddress&    inputAddress,
-                                                    bool                               includeIncompleteCurves );
+    void calculate( const std::vector<RimSummaryCase*>& sumCases, const RifEclipseSummaryAddress& inputAddress, bool includeIncompleteCurves );
+    void                                clearData();
+    static std::vector<RimSummaryCase*> validSummaryCases( const std::vector<RimSummaryCase*>& allSumCases,
+                                                           const RifEclipseSummaryAddress&     inputAddress,
+                                                           bool                                includeIncompleteCurves );
+    static std::pair<time_t, time_t>    findMinMaxTimeStep( const std::vector<RimSummaryCase*>& sumCases,
+                                                            const RifEclipseSummaryAddress&     inputAddress );
 
 private:
     std::unique_ptr<RifEnsembleStatisticsReader> m_statisticsReader;