///////////////////////////////////////////////////////////////////////////////// // // Copyright (C) 2015- Statoil ASA // Copyright (C) 2015- Ceetron Solutions AS // // ResInsight 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. // // ResInsight 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 at // for more details. // ///////////////////////////////////////////////////////////////////////////////// #include "RimWellLogPlot.h" #include "RiaApplication.h" #include "RigWellLogCurveData.h" #include "RimWellAllocationPlot.h" #include "RimWellLogCurve.h" #include "RimWellLogTrack.h" #include "RimWellRftPlot.h" #include "RimWellPltPlot.h" #include "RiuMainPlotWindow.h" #include "RiuWellLogPlot.h" #include "RiuWellLogTrack.h" #include "cvfAssert.h" #include #define RI_LOGPLOT_MINDEPTH_DEFAULT 0.0 #define RI_LOGPLOT_MAXDEPTH_DEFAULT 1000.0 namespace caf { template<> void caf::AppEnum< RimWellLogPlot::DepthTypeEnum >::setUp() { addItem(RimWellLogPlot::MEASURED_DEPTH, "MEASURED_DEPTH", "Measured Depth"); addItem(RimWellLogPlot::TRUE_VERTICAL_DEPTH, "TRUE_VERTICAL_DEPTH", "True Vertical Depth"); addItem(RimWellLogPlot::PSEUDO_LENGTH, "PSEUDO_LENGTH", "Pseudo Length"); addItem(RimWellLogPlot::CONNECTION_NUMBER, "CONNECTION_NUMBER", "Connection Number"); setDefault(RimWellLogPlot::MEASURED_DEPTH); } } // End namespace caf CAF_PDM_SOURCE_INIT(RimWellLogPlot, "WellLogPlot"); //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- RimWellLogPlot::RimWellLogPlot() { CAF_PDM_InitObject("Well Log Plot", ":/WellLogPlot16x16.png", "", ""); m_viewer = nullptr; CAF_PDM_InitField(&m_userName, "PlotDescription", QString("Well Log Plot"),"Name", "", "", ""); caf::AppEnum< RimWellLogPlot::DepthTypeEnum > depthType = MEASURED_DEPTH; CAF_PDM_InitField(&m_depthType, "DepthType", depthType, "Depth type", "", "", ""); caf::AppEnum< RiaDefines::DepthUnitType > depthUnit = RiaDefines::UNIT_METER; CAF_PDM_InitField(&m_depthUnit, "DepthUnit", depthUnit, "Depth unit", "", "", ""); CAF_PDM_InitField(&m_minVisibleDepth, "MinimumDepth", 0.0, "Min", "", "", ""); CAF_PDM_InitField(&m_maxVisibleDepth, "MaximumDepth", 1000.0, "Max", "", "", ""); CAF_PDM_InitField(&m_isAutoScaleDepthEnabled, "AutoScaleDepthEnabled", true, "Auto Scale", "", "", ""); CAF_PDM_InitField(&m_showTrackLegends, "ShowTrackLegends", true, "Show Legends", "", "", ""); CAF_PDM_InitFieldNoDefault(&m_tracks, "Tracks", "", "", "", ""); m_tracks.uiCapability()->setUiHidden(true); m_minAvailableDepth = HUGE_VAL; m_maxAvailableDepth = -HUGE_VAL; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- RimWellLogPlot::~RimWellLogPlot() { removeMdiWindowFromMdiArea(); m_tracks.deleteAllChildObjects(); deleteViewWidget(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::fieldChangedByUi(const caf::PdmFieldHandle* changedField, const QVariant& oldValue, const QVariant& newValue) { RimViewWindow::fieldChangedByUi(changedField, oldValue, newValue); if (changedField == &m_minVisibleDepth || changedField == &m_maxVisibleDepth) { applyDepthZoomFromVisibleDepth(); m_isAutoScaleDepthEnabled = false; } else if (changedField == &m_isAutoScaleDepthEnabled) { updateDepthZoom(); } else if (changedField == &m_userName) { updateMdiWindowTitle(); } if ( changedField == &m_depthType ) { RimWellAllocationPlot* wellAllocPlot; firstAncestorOrThisOfType(wellAllocPlot); if (wellAllocPlot) wellAllocPlot->loadDataAndUpdate(); else if (isRftPlotChild()) rftPlot()->loadDataAndUpdate(); else updateTracks(); } if ( changedField == &m_depthUnit) { updateTracks(); } if ( changedField == &m_showTrackLegends) { updateTracks(); if (m_viewer) m_viewer->updateChildrenLayout(); } RimWellRftPlot* rftPlot(nullptr); this->firstAncestorOrThisOfType(rftPlot); if (rftPlot) { rftPlot->updateConnectedEditors(); } RimWellPltPlot* pltPlot(nullptr); this->firstAncestorOrThisOfType(pltPlot); if (pltPlot) { pltPlot->updateConnectedEditors(); } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- QList RimWellLogPlot::calculateValueOptions(const caf::PdmFieldHandle* fieldNeedingOptions, bool* useOptionsOnly) { QList options; if (fieldNeedingOptions == &m_depthType ) { using DepthAppEnum = caf::AppEnum< DepthTypeEnum >; for (size_t i = 0; i < DepthAppEnum::size(); ++i) { DepthTypeEnum enumVal = DepthAppEnum::fromIndex(i); if (m_disabledDepthTypes.count( enumVal) == 0) { options.push_back(caf::PdmOptionItemInfo(DepthAppEnum::uiText(enumVal), enumVal)); } } } else if ( fieldNeedingOptions == &m_depthUnit) { using UnitAppEnum = caf::AppEnum< RiaDefines::DepthUnitType >; options.push_back(caf::PdmOptionItemInfo(UnitAppEnum::uiText(RiaDefines::UNIT_METER), RiaDefines::UNIT_METER)); options.push_back(caf::PdmOptionItemInfo(UnitAppEnum::uiText(RiaDefines::UNIT_FEET), RiaDefines::UNIT_FEET)); } (*useOptionsOnly) = true; return options; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- QImage RimWellLogPlot::snapshotWindowContent() { QImage image; if (m_viewer) { QPixmap pix = QPixmap::grabWidget(m_viewer); image = pix.toImage(); } return image; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::addTrack(RimWellLogTrack* track) { m_tracks.push_back(track); if (m_viewer) { track->recreateViewer(); m_viewer->addTrackPlot(track->viewer()); } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::insertTrack(RimWellLogTrack* track, size_t index) { m_tracks.insert(index, track); if (m_viewer) { track->recreateViewer(); m_viewer->insertTrackPlot(track->viewer(), index); } updateTrackNames(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::removeTrack(RimWellLogTrack* track) { if (track) { if (m_viewer) m_viewer->removeTrackPlot(track->viewer()); m_tracks.removeChildObject(track); } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::removeTrackByIndex(size_t index) { CVF_ASSERT(index < m_tracks.size()); RimWellLogTrack* track = m_tracks[index]; this->removeTrack(track); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::moveTracks(RimWellLogTrack* insertAfterTrack, const std::vector& tracksToMove) { for (size_t tIdx = 0; tIdx < tracksToMove.size(); tIdx++) { RimWellLogTrack* track = tracksToMove[tIdx]; RimWellLogPlot* wellLogPlot; track->firstAncestorOrThisOfType(wellLogPlot); if (wellLogPlot) { wellLogPlot->removeTrack(track); wellLogPlot->updateTrackNames(); wellLogPlot->updateConnectedEditors(); } } size_t index = m_tracks.index(insertAfterTrack) + 1; for (size_t tIdx = 0; tIdx < tracksToMove.size(); tIdx++) { insertTrack(tracksToMove[tIdx], index + tIdx); } updateTrackNames(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- RimWellLogTrack* RimWellLogPlot::trackByIndex(size_t index) { return m_tracks[index]; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::setDepthZoomByFactorAndCenter(double zoomFactor, double zoomCenter) { double newMinimum = zoomCenter - (zoomCenter - m_minVisibleDepth)*zoomFactor; double newMaximum = zoomCenter + (m_maxVisibleDepth - zoomCenter)*zoomFactor; setDepthZoomMinMax(newMinimum, newMaximum); setDepthAutoZoom(false); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::panDepth(double panFactor) { double delta = panFactor*(m_maxVisibleDepth - m_minVisibleDepth); setDepthZoomMinMax(m_minVisibleDepth + delta, m_maxVisibleDepth + delta); setDepthAutoZoom(false); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::setDepthZoomMinMax(double minimumDepth, double maximumDepth) { m_minVisibleDepth = minimumDepth; m_maxVisibleDepth = maximumDepth; m_minVisibleDepth.uiCapability()->updateConnectedEditors(); m_maxVisibleDepth.uiCapability()->updateConnectedEditors(); applyDepthZoomFromVisibleDepth(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::calculateAvailableDepthRange() { double minDepth = HUGE_VAL; double maxDepth = -HUGE_VAL; for (size_t tIdx = 0; tIdx < m_tracks.size(); tIdx++) { double minTrackDepth = HUGE_VAL; double maxTrackDepth = -HUGE_VAL; if (m_tracks[tIdx]->isVisible()) { m_tracks[tIdx]->availableDepthRange(&minTrackDepth, &maxTrackDepth); if (minTrackDepth < minDepth) { minDepth = minTrackDepth; } if (maxTrackDepth > maxDepth) { maxDepth = maxTrackDepth; } } } m_minAvailableDepth = minDepth; m_maxAvailableDepth = maxDepth; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::availableDepthRange(double* minimumDepth, double* maximumDepth) const { if (hasAvailableDepthRange()) { *minimumDepth = m_minAvailableDepth; *maximumDepth = m_maxAvailableDepth; } else { *minimumDepth = RI_LOGPLOT_MINDEPTH_DEFAULT; *maximumDepth = RI_LOGPLOT_MAXDEPTH_DEFAULT; } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- bool RimWellLogPlot::hasAvailableDepthRange() const { return m_minAvailableDepth < HUGE_VAL && m_maxAvailableDepth > -HUGE_VAL; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::zoomAll() { setDepthAutoZoom(true); updateDepthZoom(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- QWidget* RimWellLogPlot::viewWidget() { return m_viewer; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::setDepthAutoZoom(bool on) { m_isAutoScaleDepthEnabled = on; m_isAutoScaleDepthEnabled.uiCapability()->updateConnectedEditors(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- QString RimWellLogPlot::asciiDataForPlotExport() const { QString out; RimWellAllocationPlot* wellAllocPlot = nullptr; this->firstAncestorOrThisOfType(wellAllocPlot); if (wellAllocPlot) { out += wellAllocPlot->description(); } else { out += description(); } out += "\n"; for (RimWellLogTrack* track : m_tracks) { if (!track->isVisible()) continue; out += "\n" + track->description() + "\n"; std::vector curves = track->curvesVector(); std::vector curveNames; std::vector curveDepths; std::vector > curvesPlotXValues; for (RimWellLogCurve* curve : curves) { if (!curve->isCurveVisible()) continue; const RigWellLogCurveData* curveData = curve->curveData(); if (!curveData) continue; curveNames.push_back(curve->curveName()); if (curveNames.size() == 1) { if (depthType() == TRUE_VERTICAL_DEPTH) { curveDepths = curveData->trueDepthPlotValues(depthUnit()); } else { curveDepths = curveData->measuredDepthPlotValues(depthUnit()); } } std::vector xPlotValues = curveData->xPlotValues(); if (curveDepths.size() != xPlotValues.size() || xPlotValues.empty()) { curveNames.pop_back(); if (curveNames.empty()) { curveDepths.clear(); } continue; } curvesPlotXValues.push_back(xPlotValues); } for (size_t i = 0; i < curveDepths.size(); ++i) { if (i == 0) { if (depthType() == CONNECTION_NUMBER) out += "Connection"; else if (depthType() == MEASURED_DEPTH) out += "MD "; else if (depthType() == PSEUDO_LENGTH) out += "PL "; else if (depthType() == TRUE_VERTICAL_DEPTH) out += "TVD "; for (QString name : curveNames) out += " \t" + name; out += "\n"; } else if (curveDepths[i] == curveDepths[i-1]) { continue; } out += QString::number(curveDepths[i], 'f', 3); for (std::vector plotVector : curvesPlotXValues) { out += " \t" + QString::number(plotVector[i], 'g'); } out += "\n"; } } return out; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- RimWellRftPlot* RimWellLogPlot::rftPlot() const { RimWellRftPlot* rftPlot; firstAncestorOrThisOfType(rftPlot); return rftPlot; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- bool RimWellLogPlot::isRftPlotChild() const { return rftPlot() != nullptr; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- RimWellPltPlot* RimWellLogPlot::pltPlot() const { RimWellPltPlot* pltPlot; firstAncestorOrThisOfType(pltPlot); return pltPlot; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- bool RimWellLogPlot::isPltPlotChild() const { return pltPlot() != nullptr; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::uiOrderingForVisibleDepthRange(caf::PdmUiOrdering& uiOrdering) { caf::PdmUiGroup* gridGroup = uiOrdering.addNewGroup("Visible Depth Range"); gridGroup->add(&m_isAutoScaleDepthEnabled); gridGroup->add(&m_minVisibleDepth); gridGroup->add(&m_maxVisibleDepth); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::uiOrderingForPlot(caf::PdmUiOrdering& uiOrdering) { RimWellRftPlot* rftp; firstAncestorOrThisOfType(rftp); if (!rftp) { uiOrdering.add(&m_depthType); } RimWellAllocationPlot* wap; firstAncestorOrThisOfType(wap); if (!(wap || rftp)) { uiOrdering.add(&m_depthUnit); } uiOrdering.add(&m_showTrackLegends); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::depthZoomMinMax(double* minimumDepth, double* maximumDepth) const { *minimumDepth = m_minVisibleDepth; *maximumDepth = m_maxVisibleDepth; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::defineUiOrdering(QString uiConfigName, caf::PdmUiOrdering& uiOrdering) { uiOrdering.add(&m_userName); uiOrderingForPlot(uiOrdering); uiOrderingForVisibleDepthRange(uiOrdering); uiOrdering.skipRemainingFields(true); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::onLoadDataAndUpdate() { updateMdiWindowVisibility(); updateTracks(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::updateTracks() { if (m_showWindow) { updateDisabledDepthTypes(); for (size_t tIdx = 0; tIdx < m_tracks.size(); ++tIdx) { m_tracks[tIdx]->loadDataAndUpdate(); } calculateAvailableDepthRange(); updateDepthZoom(); } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::updateTrackNames() { for (size_t tIdx = 0; tIdx < m_tracks.size(); tIdx++) { m_tracks[tIdx]->setDescription(QString("Track %1").arg(tIdx + 1)); } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::updateDepthZoom() { if (m_isAutoScaleDepthEnabled) { applyZoomAllDepths(); } else { applyDepthZoomFromVisibleDepth(); } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::applyDepthZoomFromVisibleDepth() { if (m_viewer) { double minDepth = m_minVisibleDepth < HUGE_VAL ? m_minVisibleDepth : RI_LOGPLOT_MINDEPTH_DEFAULT; double maxDepth = m_maxVisibleDepth > -HUGE_VAL ? m_maxVisibleDepth : RI_LOGPLOT_MAXDEPTH_DEFAULT; m_viewer->setDepthZoomAndReplot(minDepth, maxDepth); } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::applyZoomAllDepths() { if (hasAvailableDepthRange()) { setDepthZoomMinMax(m_minAvailableDepth, m_maxAvailableDepth + 0.01*(m_maxAvailableDepth - m_minAvailableDepth)); } else { setDepthZoomMinMax(RI_LOGPLOT_MINDEPTH_DEFAULT, RI_LOGPLOT_MAXDEPTH_DEFAULT); } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::recreateTrackPlots() { CVF_ASSERT(m_viewer); for (size_t tIdx = 0; tIdx < m_tracks.size(); ++tIdx) { m_tracks[tIdx]->recreateViewer(); m_viewer->addTrackPlot(m_tracks[tIdx]->viewer()); } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::detachAllCurves() { for (size_t tIdx = 0; tIdx < m_tracks.size(); ++tIdx) { m_tracks[tIdx]->detachAllCurves(); } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::setDescription(const QString& description) { m_userName = description; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- QString RimWellLogPlot::description() const { return m_userName(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- QWidget* RimWellLogPlot::createViewWidget(QWidget* mainWindowParent) { m_viewer = new RiuWellLogPlot(this, mainWindowParent); recreateTrackPlots(); return m_viewer; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::deleteViewWidget() { detachAllCurves(); if (m_viewer) { m_viewer->deleteLater(); m_viewer = nullptr; } } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- RimWellLogPlot::DepthTypeEnum RimWellLogPlot::depthType() const { return m_depthType.value(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::setDepthType(DepthTypeEnum depthType) { m_depthType = depthType; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- RiaDefines::DepthUnitType RimWellLogPlot::depthUnit() const { return m_depthUnit.value(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- QString RimWellLogPlot::depthPlotTitle() const { QString depthTitle = "Depth"; switch (m_depthType.value()) { case MEASURED_DEPTH: depthTitle = "MD"; break; case TRUE_VERTICAL_DEPTH: depthTitle = "TVD"; break; case PSEUDO_LENGTH: depthTitle = "PL"; break; case CONNECTION_NUMBER: depthTitle = "Connection"; break; } if (m_depthType() == CONNECTION_NUMBER) return depthTitle; if (m_depthUnit == RiaDefines::UNIT_METER) { depthTitle += " [m]"; } else if (m_depthUnit == RiaDefines::UNIT_FEET) { depthTitle += " [ft]"; } else if (m_depthUnit == RiaDefines::UNIT_NONE) { depthTitle += ""; } return depthTitle; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- bool RimWellLogPlot::isTrackLegendsVisible() const { return m_showTrackLegends(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::setTrackLegendsVisible(bool doShow) { m_showTrackLegends = doShow; } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- size_t RimWellLogPlot::trackIndex(RimWellLogTrack* track) { return m_tracks.index(track); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::setDepthUnit(RiaDefines::DepthUnitType depthUnit) { m_depthUnit = depthUnit; updateTracks(); } //-------------------------------------------------------------------------------------------------- /// //-------------------------------------------------------------------------------------------------- void RimWellLogPlot::updateDisabledDepthTypes() { m_disabledDepthTypes.clear(); RimWellAllocationPlot* wap; firstAncestorOrThisOfType(wap); if (wap) { m_disabledDepthTypes.insert(MEASURED_DEPTH); if (m_disabledDepthTypes.count(m_depthType() )) { m_depthType = CONNECTION_NUMBER; } } else if (isRftPlotChild()) { m_disabledDepthTypes.insert(MEASURED_DEPTH); m_disabledDepthTypes.insert(PSEUDO_LENGTH); m_disabledDepthTypes.insert(CONNECTION_NUMBER); } else if (isPltPlotChild()) { m_disabledDepthTypes.insert(TRUE_VERTICAL_DEPTH); m_disabledDepthTypes.insert(CONNECTION_NUMBER); } else { m_disabledDepthTypes.insert(PSEUDO_LENGTH); m_disabledDepthTypes.insert(CONNECTION_NUMBER); } }