(#541) Added feature for export of well log plot curves to a LAS file

Did some refactoring/improvements by introducing new class
RigWellLogCurveData.

ApplicationCode/ReservoirDataModel/CMakeLists_files.cmake

@@ -36,7 +36,7 @@ ${CEE_CURRENT_LIST_DIR}RigResultAccessor2d.h
 ${CEE_CURRENT_LIST_DIR}RigTernaryResultAccessor2d.h
 ${CEE_CURRENT_LIST_DIR}RigNativeStatCalc.h
 ${CEE_CURRENT_LIST_DIR}RigMultipleDatasetStatCalc.h
-${CEE_CURRENT_LIST_DIR}RigWellLogFile.h
+${CEE_CURRENT_LIST_DIR}RigWellLogCurveData.h
 )
 
 set (SOURCE_GROUP_SOURCE_FILES
@@ -67,7 +67,7 @@ ${CEE_CURRENT_LIST_DIR}cvfGeometryTools.cpp
 ${CEE_CURRENT_LIST_DIR}RigTernaryResultAccessor2d.cpp
 ${CEE_CURRENT_LIST_DIR}RigNativeStatCalc.cpp
 ${CEE_CURRENT_LIST_DIR}RigMultipleDatasetStatCalc.cpp
-${CEE_CURRENT_LIST_DIR}RigWellLogFile.cpp
+${CEE_CURRENT_LIST_DIR}RigWellLogCurveData.cpp
 )
 
 list(APPEND CODE_HEADER_FILES

ApplicationCode/ReservoirDataModel/CMakeLists_filesNotToUnitTest.cmake

@@ -7,11 +7,13 @@ endif()
 set (SOURCE_GROUP_HEADER_FILES
     ${CEE_CURRENT_LIST_DIR}RigGeoMechWellLogExtractor.h
     ${CEE_CURRENT_LIST_DIR}RigCaseToCaseCellMapper.h
+    ${CEE_CURRENT_LIST_DIR}RigWellLogFile.h
 )
 
 set (SOURCE_GROUP_SOURCE_FILES
     ${CEE_CURRENT_LIST_DIR}RigGeoMechWellLogExtractor.cpp
     ${CEE_CURRENT_LIST_DIR}RigCaseToCaseCellMapper.cpp
+    ${CEE_CURRENT_LIST_DIR}RigWellLogFile.cpp
 )
 
 list(APPEND CODE_HEADER_FILES

ApplicationCode/ReservoirDataModel/RigWellLogCurveData.cpp

@@ -0,0 +1,187 @@
+/////////////////////////////////////////////////////////////////////////////////
+//
+//  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 <http://www.gnu.org/licenses/gpl.html> 
+//  for more details.
+//
+/////////////////////////////////////////////////////////////////////////////////
+
+#include "RigWellLogCurveData.h"
+
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+RigWellLogCurveData::RigWellLogCurveData()
+{
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+RigWellLogCurveData::~RigWellLogCurveData()
+{
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+void RigWellLogCurveData::setPoints(const std::vector<double>& xValues, const std::vector<double>& yValues)
+{
+    CVF_ASSERT(xValues.size() == yValues.size());
+
+    m_xValues = xValues;
+    m_yValues = yValues;
+    
+    calculateValidPointIntervals();
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+const std::vector<double>& RigWellLogCurveData::xValues() const
+{
+    return m_xValues;
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+const std::vector<double>& RigWellLogCurveData::yValues() const
+{
+    return m_yValues;
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+const std::vector< std::pair<size_t, size_t> >& RigWellLogCurveData::intervals() const
+{
+    return m_validPointIntervals;
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+std::vector<double> RigWellLogCurveData::validXValues() const
+{
+    std::vector<double> filteredValues;
+    addValuesFromIntervals(m_xValues, m_validPointIntervals, &filteredValues);
+
+    return filteredValues;
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+std::vector<double> RigWellLogCurveData::validYValues() const
+{
+    std::vector<double> filteredValues;
+    addValuesFromIntervals(m_yValues, m_validPointIntervals, &filteredValues);
+
+    return filteredValues;
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+std::vector< std::pair<size_t, size_t> > RigWellLogCurveData::validPointsIntervals() const
+{
+    std::vector< std::pair<size_t, size_t> > intervals;
+    filteredIntervals(m_validPointIntervals, &intervals);
+
+    return intervals;
+}
+
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+void RigWellLogCurveData::calculateValidPointIntervals()
+{
+    calculateIntervalsOfValidValues(m_xValues, &m_validPointIntervals);
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+void RigWellLogCurveData::calculateIntervalsOfValidValues(const std::vector<double>& values, std::vector< std::pair<size_t, size_t> >* intervals)
+{
+    CVF_ASSERT(intervals);
+
+    int startIdx = -1;
+    size_t vIdx = 0;
+
+    size_t valueCount = values.size();
+    while (vIdx < valueCount)
+    {
+        double value = values[vIdx];
+        if (value == HUGE_VAL || value == -HUGE_VAL || value != value)
+        {
+            if (startIdx >= 0)
+            {
+                intervals->push_back(std::make_pair(startIdx, vIdx - 1));
+                startIdx = -1;
+            }
+        }
+        else if (startIdx < 0)
+        {
+            startIdx = (int)vIdx;
+        }
+
+        vIdx++;
+    }
+
+    if (startIdx >= 0 && startIdx < ((int)valueCount))
+    {
+        intervals->push_back(std::make_pair(startIdx, valueCount - 1));
+    }
+}
+
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+void RigWellLogCurveData::addValuesFromIntervals(const std::vector<double>& values, const std::vector< std::pair<size_t, size_t> >& intervals, std::vector<double>* filteredValues)
+{
+    CVF_ASSERT(filteredValues);
+
+    for (size_t intIdx = 0; intIdx < intervals.size(); intIdx++)
+    {
+        for (size_t vIdx = intervals[intIdx].first; vIdx <= intervals[intIdx].second; vIdx++)
+        {
+            filteredValues->push_back(values[vIdx]);
+        }
+    }
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+void RigWellLogCurveData::filteredIntervals(const std::vector< std::pair<size_t, size_t> >& intervals, std::vector< std::pair<size_t, size_t> >* fltrIntervals)
+{
+    CVF_ASSERT(fltrIntervals);
+
+    const size_t intervalCount = intervals.size();
+    if (intervalCount < 1) return;
+
+    size_t index = 0;
+    for (size_t intIdx = 0; intIdx < intervalCount; intIdx++)
+    {
+        size_t intervalSize = intervals[intIdx].second - intervals[intIdx].first + 1;
+        fltrIntervals->push_back(std::make_pair(index, index + intervalSize - 1));
+
+        index += intervalSize;
+    }
+}

ApplicationCode/ReservoirDataModel/RigWellLogCurveData.h

@@ -0,0 +1,57 @@
+/////////////////////////////////////////////////////////////////////////////////
+//
+//  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 <http://www.gnu.org/licenses/gpl.html> 
+//  for more details.
+//
+/////////////////////////////////////////////////////////////////////////////////
+
+#pragma once
+
+#include "cvfBase.h"
+#include "cvfObject.h"
+
+#include <vector>
+
+//==================================================================================================
+/// 
+//==================================================================================================
+class RigWellLogCurveData : public cvf::Object
+{
+public:
+    RigWellLogCurveData();
+    virtual ~RigWellLogCurveData();
+
+    void setPoints(const std::vector<double>& xValues, const std::vector<double>& yValues);
+
+    const std::vector<double>&                      xValues() const;
+    const std::vector<double>&                      yValues() const;
+    const std::vector< std::pair<size_t, size_t> >& intervals() const;
+
+    std::vector<double>                             validXValues() const;
+    std::vector<double>                             validYValues() const;
+    std::vector< std::pair<size_t, size_t> >        validPointsIntervals() const;
+
+private:
+    void        calculateValidPointIntervals();
+
+    static void calculateIntervalsOfValidValues(const std::vector<double>& values, std::vector< std::pair<size_t, size_t> >* intervals);
+    static void addValuesFromIntervals(const std::vector<double>& values, const std::vector< std::pair<size_t, size_t> >& intervals, std::vector<double>* filteredValues);
+    static void filteredIntervals(const std::vector< std::pair<size_t, size_t> >& intervals, std::vector< std::pair<size_t, size_t> >* fltrIntervals);
+
+private:
+    std::vector<double>                         m_xValues;
+    std::vector<double>                         m_yValues;
+    std::vector< std::pair<size_t, size_t> >    m_validPointIntervals;
+};

ApplicationCode/ReservoirDataModel/RigWellLogFile.cpp

@@ -19,6 +19,8 @@
 
 #include "RigWellLogFile.h"
 
+#include "RimWellLogPlotCurve.h"
+
 #include "well.hpp"
 #include "laswell.hpp"
 
@@ -234,3 +236,43 @@ QString RigWellLogFile::wellLogChannelUnit(const QString& wellLogChannelName) co
 
     return unit;
 }
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+bool RigWellLogFile::exportToLasFile(const RimWellLogPlotCurve* curve, const QString& fileName)
+{
+    CVF_ASSERT(curve);
+
+    const RigWellLogCurveData* curveData = curve->curveData();
+    if (!curveData)
+    {
+        return false;
+    }
+
+    double minX, maxX;
+    curve->valueRange(&minX, &maxX);
+
+    // Might want to use a different way to find an absent/"null" value, maybe use the default if possible
+    double absentValue = ((size_t) maxX) + 1000;
+
+    std::vector<double> wellLogValues = curveData->xValues();
+    for (size_t vIdx = 0; vIdx < wellLogValues.size(); vIdx++)
+    {
+        double value = wellLogValues[vIdx];
+        if (value == HUGE_VAL || value == -HUGE_VAL || value != value)
+        {
+            wellLogValues[vIdx] = absentValue;
+        }
+    }
+
+    NRLib::LasWell lasFile;
+    lasFile.AddLog("DEPTH", "m", "Depth [m]", curveData->yValues());
+    lasFile.AddLog(curve->name().trimmed().toStdString(), "NO_UNIT", "PARAMETERINFO", wellLogValues);
+    lasFile.SetMissing(absentValue);
+
+    std::vector<std::string> commentHeader;
+    lasFile.WriteToFile(fileName.toStdString(), commentHeader);
+
+    return true;
+}

ApplicationCode/ReservoirDataModel/RigWellLogFile.h

@@ -30,6 +30,7 @@ namespace NRLib
     class Well;
 }
 
+class RimWellLogPlotCurve;
 
 //==================================================================================================
 /// 
@@ -51,6 +52,8 @@ public:
     QString depthUnit() const;
     QString wellLogChannelUnit(const QString& wellLogChannelName) const;
 
+    static bool exportToLasFile(const RimWellLogPlotCurve* curve, const QString& fileName);
+
 private:
     void close();