OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

#5539 Support ensembles in Analysis plots

Browse Source
This commit is contained in:
Jacob Støren
2020-02-26 13:51:57 +01:00
parent d540737a53
commit fd4cbc6c87
1 changed files with 97 additions and 78 deletions
Download Patch File Download Diff File Expand all files Collapse all files

175
ApplicationCode/ProjectDataModel/AnalysisPlots/RimAnalysisPlot.cpp
View File

@@ -644,7 +644,7 @@ void RimAnalysisPlot::onLoadDataAndUpdate()
/// ///
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
QString assignGroupingText( RimAnalysisPlot::SortGroupType sortGroup, QString assignGroupingText( RimAnalysisPlot::SortGroupType sortGroup,
const RimAnalysisPlotDataEntry* dataEntry, const RiaSummaryCurveDefinition dataEntry,
const QString& timestepString ) const QString& timestepString )
{ {
QString groupingText; QString groupingText;
@@ -653,29 +653,29 @@ QString assignGroupingText( RimAnalysisPlot::SortGroupType sortGroup,
{ {
case RimAnalysisPlot::SUMMARY_ITEM: case RimAnalysisPlot::SUMMARY_ITEM:
{ {
RifEclipseSummaryAddress addr = dataEntry->summaryAddress(); RifEclipseSummaryAddress addr = dataEntry.summaryAddress();
groupingText = QString::fromStdString( addr.itemUiText() ); groupingText = QString::fromStdString( addr.itemUiText() );
} }
break; break;
case RimAnalysisPlot::CASE: case RimAnalysisPlot::CASE:
{ {
if ( dataEntry->summaryCase() ) if ( dataEntry.summaryCase() )
{ {
groupingText = dataEntry->summaryCase()->displayCaseName(); groupingText = dataEntry.summaryCase()->displayCaseName();
} }
} }
break; break;
case RimAnalysisPlot::ENSEMBLE: case RimAnalysisPlot::ENSEMBLE:
{ {
if ( dataEntry->ensemble() ) if ( dataEntry.ensemble() )
{ {
groupingText = dataEntry->ensemble()->name(); groupingText = dataEntry.ensemble()->name();
} }
} }
break; break;
case RimAnalysisPlot::QUANTITY: case RimAnalysisPlot::QUANTITY:
{ {
RifEclipseSummaryAddress addr = dataEntry->summaryAddress(); RifEclipseSummaryAddress addr = dataEntry.summaryAddress();
groupingText = QString::fromStdString( addr.quantityName() ); groupingText = QString::fromStdString( addr.quantityName() );
} }
@@ -708,94 +708,113 @@ void RimAnalysisPlot::addDataToChartBuilder( RiuGroupedBarChartBuilder& chartBui
for ( const RimAnalysisPlotDataEntry* dataEntry : m_data ) for ( const RimAnalysisPlotDataEntry* dataEntry : m_data )
{ {
if ( !dataEntry->summaryCase() ) continue; // Todo, ensembles RiaSummaryCurveDefinition orgBarDataEntry = dataEntry->curveDefinition();
RifSummaryReaderInterface* reader = dataEntry->summaryCase()->summaryReader(); std::vector<RiaSummaryCurveDefinition> barDataDefinitions;
if ( !reader ) continue;
// Todo: if ( orgBarDataEntry.isEnsembleCurve() )
// If is RimGridSummaryCase and using summary item ans legend and summary items are wells, then:
/// use color from eclCase->defaultWellColor( wellName );
const std::vector<time_t>& timesteps = reader->timeSteps( dataEntry->summaryAddress() );
std::vector<double> values;
reader->values( dataEntry->summaryAddress(), &values );
// Find selected timestep indices
std::vector<int> selectedTimestepIndices;
for ( time_t tt : selectedTimesteps )
{ {
for ( int timestepIdx = 0; static_cast<unsigned>( timestepIdx ) < timesteps.size(); ++timestepIdx ) std::vector<RimSummaryCase*> sumCases = orgBarDataEntry.ensemble()->allSummaryCases();
for ( auto sumCase : sumCases )
{ {
if ( timesteps[timestepIdx] == tt ) barDataDefinitions.push_back(
{ RiaSummaryCurveDefinition( sumCase, orgBarDataEntry.summaryAddress(), orgBarDataEntry.ensemble() ) );
selectedTimestepIndices.push_back( timestepIdx );
break;
}
} }
} }
else
for ( int timestepIdx : selectedTimestepIndices )
{ {
double sortValue = std::numeric_limits<double>::infinity(); barDataDefinitions.push_back( orgBarDataEntry );
}
QDateTime dateTime = RiaQDateTimeTools::fromTime_t( timesteps[timestepIdx] ); for ( const RiaSummaryCurveDefinition& curveDef : barDataDefinitions )
QString formatString = RiaQDateTimeTools::createTimeFormatStringFromDates( {dateTime} ); {
RifSummaryReaderInterface* reader = curveDef.summaryCase()->summaryReader();
if ( !reader ) continue;
QString timestepString = dateTime.toString( formatString ); // Todo:
// If is RimGridSummaryCase and using summary item as legend and summary items are wells, then:
/// use color from eclCase->defaultWellColor( wellName );
QString majorText = assignGroupingText( m_majorGroupType(), dataEntry, timestepString ); const std::vector<time_t>& timesteps = reader->timeSteps( curveDef.summaryAddress() );
QString medText = assignGroupingText( m_mediumGroupType(), dataEntry, timestepString );
QString minText = assignGroupingText( m_minorGroupType(), dataEntry, timestepString );
QString legendText = assignGroupingText( m_sortGroupForLegend(), dataEntry, timestepString );
double value = values[timestepIdx]; std::vector<double> values;
reader->values( curveDef.summaryAddress(), &values );
switch ( m_valueSortOperation() ) // Find selected timestep indices
std::vector<int> selectedTimestepIndices;
for ( time_t tt : selectedTimesteps )
{ {
case VALUE: for ( int timestepIdx = 0; static_cast<unsigned>( timestepIdx ) < timesteps.size(); ++timestepIdx )
sortValue = value;
break;
case ABS_VALUE:
sortValue = fabs( value );
break;
}
QString barText;
QString separator = " ";
if ( m_useBarText() )
{
if ( m_useQuantityInBarText )
{ {
barText += QString::fromStdString( dataEntry->summaryAddress().quantityName() ) + separator; if ( timesteps[timestepIdx] == tt )
} {
selectedTimestepIndices.push_back( timestepIdx );
if ( m_useSummaryItemInBarText ) break;
{ }
barText += QString::fromStdString( dataEntry->summaryAddress().itemUiText() ) + separator;
}
if ( m_useCaseInBarText && dataEntry->summaryCase() )
{
barText += dataEntry->summaryCase()->displayCaseName() + separator;
}
if ( m_useEnsembleInBarText && dataEntry->ensemble() )
{
barText += dataEntry->ensemble()->name() + separator;
}
if ( m_useTimeStepInBarText )
{
barText += timestepString + separator;
} }
} }
chartBuilder.addBarEntry( majorText, medText, minText, sortValue, legendText, barText, values[timestepIdx] ); for ( int timestepIdx : selectedTimestepIndices )
{
double sortValue = std::numeric_limits<double>::infinity();
QDateTime dateTime = RiaQDateTimeTools::fromTime_t( timesteps[timestepIdx] );
QString formatString = RiaQDateTimeTools::createTimeFormatStringFromDates( {dateTime} );
QString timestepString = dateTime.toString( formatString );
QString majorText = assignGroupingText( m_majorGroupType(), curveDef, timestepString );
QString medText = assignGroupingText( m_mediumGroupType(), curveDef, timestepString );
QString minText = assignGroupingText( m_minorGroupType(), curveDef, timestepString );
QString legendText = assignGroupingText( m_sortGroupForLegend(), curveDef, timestepString );
double value = values[timestepIdx];
switch ( m_valueSortOperation() )
{
case VALUE:
sortValue = value;
break;
case ABS_VALUE:
sortValue = fabs( value );
break;
}
QString barText;
QString separator = " ";
if ( m_useBarText() )
{
if ( m_useQuantityInBarText )
{
barText += QString::fromStdString( curveDef.summaryAddress().quantityName() ) + separator;
}
if ( m_useSummaryItemInBarText )
{
barText += QString::fromStdString( curveDef.summaryAddress().itemUiText() ) + separator;
}
if ( m_useCaseInBarText && curveDef.summaryCase() )
{
barText += curveDef.summaryCase()->displayCaseName() + separator;
}
if ( m_useEnsembleInBarText && curveDef.ensemble() )
{
barText += curveDef.ensemble()->name() + separator;
}
if ( m_useTimeStepInBarText )
{
barText += timestepString + separator;
}
}
chartBuilder.addBarEntry( majorText, medText, minText, sortValue, legendText, barText, values[timestepIdx] );
}
} }
} }
} }