Skip to content

Contour map improvements#321

Open
magnesj wants to merge 2 commits intodevfrom
contour-map-01
Open

Contour map improvements#321
magnesj wants to merge 2 commits intodevfrom
contour-map-01

Conversation

@magnesj
Copy link
Owner

@magnesj magnesj commented Aug 29, 2025 •
edited by qodo-code-review bot
Loading

PR Type

Enhancement, Other

Description

  • Add texture-based rendering for contour maps

  • Migrate obsolete color fields to main color

  • Update version to 2025.09.0-RC_1

  • Add surfio license information

Diagram Walkthrough

flowchart LR
  A["Contour Map"] --> B["Texture Image"]
  A --> C["Triangles with Color"]
  B --> D["Textured Quad Rendering"]
  C --> E["Traditional Triangle Rendering"]
  F["Version Update"] --> G["2025.09.0-RC_1"]
  H["Color Migration"] --> I["Main Color Field"]
Loading

File Walkthrough
Relevant files
Enhancement
6 files
RivContourMapProjectionPartMgr.cpp
Add texture-based quad rendering for contour maps               
+219/-0 
RimContourMapProjection.cpp
Add UI controls for texture and triangle rendering             
+30/-7   
RimEclipseContourMapView.cpp
Implement conditional rendering based on display options 
+24/-6   
RiaEnsembleNameTools.cpp
Remove underscores from ensemble short names                         
+1/-0     
RivContourMapProjectionPartMgr.h
Add texture rendering method declarations                               
+16/-0   
RimContourMapProjection.h
Add texture and triangle display option fields                     
+4/-0     
Bug fix
3 files
RimEnsembleCurveSet.cpp
Migrate obsolete color field to main color                             
+11/-0   
RimEnsembleCurveFilter.cpp
Set default filter mode for older projects                             
+5/-0     
RimEnsembleCurveSet.h
Add obsolete color field for migration                                     
+2/-0     
Configuration changes
2 files
ResInsightVersion.cmake
Update version to 2025.09.0-RC_1                                                 
+5/-24   
ri-objectKeywords.txt
Add project file keywords for version 2025.09                       
+524/-0 
Documentation
1 files
LicenseInformation.txt
Add MIT license information for surfio                                     
+26/-1   
Additional files
4 files
ri-fieldKeywords.txt +6433/-0
RiaGuiApplication.cpp +0/-1     
RiaPreferencesSystem.cpp +0/-12   
RiaPreferencesSystem.h +0/-3     

magnesj added 2 commits August 28, 2025 20:14
@magnesj
Prepare 2025.09.0-RC_1
35d1279 
* Migrates obsolete color field to main color
* Sets filter mode for older project files
* Adds surfio license information
* Removes log to file preference
* Refines short name generation.
@magnesj
Prototype for use of texture on quad
eec01bd
@qodo-code-review
Copy link

qodo-code-review bot commented Aug 29, 2025

PR Reviewer Guide 🔍

Here are some key observations to aid the review process:

⏱️ Estimated effort to review: 3 🔵🔵🔵⚪⚪
🧪 No relevant tests
🔒 No security concerns identified
⚡ Recommended focus areas for review

Z-depth Logic

The z-value for the textured quad is computed as z = originZ + 0.3 * originZ, which scales with originZ and might be behind/above unintended geometry or vary per dataset. Consider a fixed offset in display space or deriving from grid thickness; verify ordering with other parts and transparency sorting.

 
double zValueForContourMap = contourMapGrid.origin3d().z();
zValueForContourMap += 0.3 * zValueForContourMap;

cvf::Vec3dArray displayCoords;
displayCoords.reserve( 4 );
for ( int i = 0; i < 4; i++ )
{
    auto displayCoord = displayCoordTransform->transformToDisplayCoord( domainCoords[i] );
    displayCoord.z()  = zValueForContourMap;

    displayCoords.add( displayCoord );
Transparency Handling

Textured rendering uses blending and marks the part as TransparentSeismic while also setting local 'transparent' flag always true. Verify alpha semantics in texture generation (infinite -> alpha 0) and ensure correct depth write/state to avoid halo/sorting artifacts; consider enabling/disabling depth writes or using pre-multiplied alpha consistently.

 
    bool transparent = true;
    auto part        = createSingleTexturedQuadPart( displayCoords, textureImage, transparent );
    part->setSourceInfo( new RivObjectSourceInfo( m_pdmObject.p() ) );
    part->setPriority( RivPartPriority::BaseLevel );

    model->addPart( part.p() );
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
cvf::ref<cvf::Vec2fArray> RivContourMapProjectionPartMgr::createTextureCoords( const std::vector<double>& values,
                                                                               cvf::ScalarMapper*         scalarMapper ) const
{
    cvf::ref<cvf::Vec2fArray> textureCoords = new cvf::Vec2fArray( values.size() );

#pragma omp parallel for
    for ( int i = 0; i < (int)values.size(); ++i )
    {
        if ( values[i] != std::numeric_limits<double>::infinity() )
        {
            cvf::Vec2f textureCoord = scalarMapper->mapToTextureCoord( values[i] );
            textureCoord.y()        = 0.0;
            ( *textureCoords )[i]   = textureCoord;
        }
        else
        {
            ( *textureCoords )[i] = cvf::Vec2f( 0.0, 1.0 );
        }
    }
    return textureCoords;
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RivContourMapProjectionPartMgr::appendContourLinesToModel( const cvf::Camera*                camera,
                                                                cvf::ModelBasicList*              model,
                                                                const caf::DisplayCoordTransform* displayCoordTransform,
                                                                const std::vector<RigContourPolygonsTools::ContourPolygons>& contourLinePolygons,
                                                                const RigContourMapGrid&          contourMapGrid,
                                                                cvf::ScalarMapper*                mapper,
                                                                bool                              showContourLines,
                                                                bool                              showContourLabels,
                                                                RiaNumberFormat::NumberFormatType numberFormat,
                                                                int                               precision )
{
    if ( showContourLines )
    {
        std::vector<std::vector<cvf::BoundingBox>> labelBBoxes;
        std::vector<cvf::ref<cvf::Drawable>>       labelDrawables;

        if ( showContourLabels )
        {
            labelDrawables =
                createContourLabels( camera, displayCoordTransform, &labelBBoxes, contourLinePolygons, contourMapGrid, mapper, numberFormat, precision );
        }

        std::vector<std::vector<cvf::ref<cvf::Drawable>>> contourDrawablesForAllLevels =
            createContourPolygons( displayCoordTransform, labelBBoxes, contourLinePolygons, mapper, contourMapGrid );

        std::vector<double> tickValues;
        mapper->majorTickValues( &tickValues );

        for ( size_t i = 0; i < contourDrawablesForAllLevels.size(); ++i )
        {
            std::vector<cvf::ref<cvf::Drawable>> contourDrawables = contourDrawablesForAllLevels[i];

            cvf::Color3f backgroundColor( mapper->mapToColor( tickValues[i] ) );
            cvf::Color3f lineColor = RiaColorTools::contrastColor( backgroundColor );

            for ( cvf::ref<cvf::Drawable> contourDrawable : contourDrawables )
            {
                if ( contourDrawable.notNull() && contourDrawable->boundingBox().isValid() )
                {
                    caf::MeshEffectGenerator meshEffectGen( lineColor );
                    meshEffectGen.setLineWidth( 1.0f );
                    caf::MeshEffectGenerator::createAndConfigurePolygonOffsetRenderState( caf::PO_1 );

                    cvf::ref<cvf::Effect> effect = meshEffectGen.generateCachedEffect();

                    cvf::ref<cvf::Part> part = new cvf::Part;
                    part->setName( "RivContourMapProjectionPartMgr::contourDrawable_mesh" );
                    part->setDrawable( contourDrawable.p() );
                    part->setEffect( effect.p() );
                    part->setPriority( RivPartPriority::MeshLines );
                    part->setSourceInfo( new RivMeshLinesSourceInfo( m_pdmObject.p() ) );

                    model->addPart( part.p() );
                }
            }
        }

        if ( showContourLabels )
        {
            for ( auto labelDrawableRef : labelDrawables )
            {
                cvf::ref<cvf::Part> part = new cvf::Part;
                part->setName( "RivContourMapProjectionPartMgr::labelDrawableRef" );
                part->setDrawable( labelDrawableRef.p() );
                part->setEffect( m_labelEffect.p() );
                part->setPriority( RivPartPriority::Text );
                part->setSourceInfo( new RivMeshLinesSourceInfo( m_pdmObject.p() ) );
                model->addPart( part.p() );
            }
        }
    }
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
cvf::ref<cvf::DrawableText> RivContourMapProjectionPartMgr::createTextLabel( const cvf::Color3f& textColor, const cvf::Color3f& backgroundColor )
{
    auto font = RiaFontCache::getFont( RiaFontCache::FontSize::FONT_SIZE_10 );

    cvf::ref<cvf::DrawableText> labelDrawable = new cvf::DrawableText();
    labelDrawable->setFont( font.p() );
    labelDrawable->setCheckPosVisible( false );
    labelDrawable->setUseDepthBuffer( true );
    labelDrawable->setDrawBorder( false );
    labelDrawable->setDrawBackground( false );
    labelDrawable->setBackgroundColor( backgroundColor );
    labelDrawable->setVerticalAlignment( cvf::TextDrawer::CENTER );
    labelDrawable->setTextColor( textColor );
    labelDrawable->setBorderColor( textColor );

    return labelDrawable;
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
cvf::ref<cvf::Part> RivContourMapProjectionPartMgr::createProjectionMapPart( const caf::DisplayCoordTransform* displayCoordTransform,
                                                                             const std::vector<cvf::Vec4d>&    vertices,
                                                                             const RigContourMapGrid&          contourMapGrid,
                                                                             const cvf::Color3f&               backgroundColor,
                                                                             cvf::ScalarMapper*                scalarMapper ) const
{
    if ( vertices.size() < 3u )
    {
        return cvf::ref<cvf::Part>();
    }
    cvf::ref<cvf::Vec3fArray> vertexArray = new cvf::Vec3fArray( vertices.size() );
    cvf::ref<cvf::UIntArray>  faceList    = new cvf::UIntArray( vertices.size() );
    std::vector<double>       values( vertices.size() );
    for ( uint i = 0; i < vertices.size(); ++i )
    {
        cvf::Vec3d globalVertex = cvf::Vec3d( vertices[i].x(), vertices[i].y(), vertices[i].z() ) + contourMapGrid.origin3d();
        cvf::Vec3f displayVertexPos( displayCoordTransform->transformToDisplayCoord( globalVertex ) );
        ( *vertexArray )[i] = displayVertexPos;
        ( *faceList )[i]    = i;
        values[i]           = vertices[i].w();
    }

    cvf::ref<cvf::PrimitiveSetIndexedUInt> indexUInt = new cvf::PrimitiveSetIndexedUInt( cvf::PrimitiveType::PT_TRIANGLES, faceList.p() );

    cvf::ref<cvf::DrawableGeo> geo = new cvf::DrawableGeo;
    geo->addPrimitiveSet( indexUInt.p() );
    geo->setVertexArray( vertexArray.p() );

    cvf::ref<cvf::Part> part = new cvf::Part;
    part->setName( "RivContourMapProjectionPartMgr::createProjectionMapPart" );
    part->setDrawable( geo.p() );

    cvf::ref<cvf::Vec2fArray> textureCoords = createTextureCoords( values, scalarMapper );
    RivScalarMapperUtils::applyTextureResultsToPart( part.p(), textureCoords.p(), scalarMapper, 1.0f, caf::FC_NONE, true, backgroundColor );

    part->setSourceInfo( new RivObjectSourceInfo( m_pdmObject.p() ) );
    part->setPriority( RivPartPriority::BaseLevel );
    return part;
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::vector<std::vector<cvf::ref<cvf::Drawable>>>
    RivContourMapProjectionPartMgr::createContourPolygons( const caf::DisplayCoordTransform*                 displayCoordTransform,
                                                           const std::vector<std::vector<cvf::BoundingBox>>& labelBBoxes,
                                                           const std::vector<RigContourPolygonsTools::ContourPolygons>& contourLinePolygons,
                                                           cvf::ScalarMapper*                                           scalarMapper,
                                                           const RigContourMapGrid& contourMapGrid ) const
{
    std::vector<double> tickValues;
    scalarMapper->majorTickValues( &tickValues );

    std::vector<std::vector<cvf::ref<cvf::Drawable>>> contourDrawablesForAllLevels;
    contourDrawablesForAllLevels.resize( tickValues.size() );

    for ( size_t i = 1; i < contourLinePolygons.size(); ++i )
    {
        std::vector<cvf::ref<cvf::Drawable>> contourDrawables;

        for ( size_t j = 0; j < contourLinePolygons[i].size(); ++j )
        {
            if ( contourLinePolygons[i][j].vertices.empty() ) continue;

            // cvf::String::number does not allow precision on 'g' formats, so use Qt.
            QString     qLabelText = QString::number( contourLinePolygons[i][j].value, 'g', 2 );
            cvf::String labelText  = cvfqt::Utils::toString( qLabelText );

            size_t nVertices = contourLinePolygons[i][j].vertices.size();

            std::vector<cvf::Vec3f> displayLines;
            displayLines.reserve( nVertices * 2 );
            for ( size_t v = 0; v < nVertices; ++v )
            {
                cvf::Vec3d globalVertex1  = contourLinePolygons[i][j].vertices[v] + contourMapGrid.origin3d();
                cvf::Vec3d displayVertex1 = displayCoordTransform->transformToDisplayCoord( globalVertex1 );

                cvf::Vec3d globalVertex2;
                if ( v < nVertices - 1 )
                    globalVertex2 = contourLinePolygons[i][j].vertices[v + 1] + contourMapGrid.origin3d();
                else
                    globalVertex2 = contourLinePolygons[i][j].vertices[0] + contourMapGrid.origin3d();

                cvf::Vec3d displayVertex2 = displayCoordTransform->transformToDisplayCoord( globalVertex2 );

                cvf::BoundingBox lineBBox;
                lineBBox.add( displayVertex1 );
                lineBBox.add( displayVertex2 );

                bool addOriginalSegment = true;
                if ( !labelBBoxes.empty() )
                {
                    for ( const cvf::BoundingBox& existingBBox : labelBBoxes[i] )
                    {
                        if ( lineBBox.intersects( existingBBox ) )
                        {
                            if ( existingBBox.contains( displayVertex1 ) && existingBBox.contains( displayVertex2 ) )
                            {
                                addOriginalSegment = false;
                            }
                            else
                            {
                                cvf::Vec3d dir = displayVertex2 - displayVertex1;

                                cvf::Ray ray;
                                ray.setOrigin( displayVertex1 );
                                ray.setDirection( dir.getNormalized() );
                                ray.setMaximumDistance( dir.length() );

                                if ( !existingBBox.contains( displayVertex1 ) )
                                {
                                    cvf::Vec3d intersection;
                                    bool       hit = ray.boxIntersect( existingBBox, &intersection );
                                    if ( hit )
                                    {
                                        displayLines.push_back( cvf::Vec3f( displayVertex1 ) );
                                        displayLines.push_back( cvf::Vec3f( intersection ) );
                                        addOriginalSegment = false;
                                    }
                                }

                                if ( !existingBBox.contains( displayVertex2 ) )
                                {
                                    ray.setOrigin( displayVertex2 );
                                    ray.setDirection( -ray.direction() );
                                    cvf::Vec3d intersection;
                                    bool       hit = ray.boxIntersect( existingBBox, &intersection );
                                    if ( hit )
                                    {
                                        displayLines.push_back( cvf::Vec3f( intersection ) );
                                        displayLines.push_back( cvf::Vec3f( displayVertex2 ) );
                                        addOriginalSegment = false;
                                    }
                                }
                            }
                        }
                    }
                }
                if ( addOriginalSegment )
                {
                    displayLines.push_back( cvf::Vec3f( displayVertex1 ) );
                    displayLines.push_back( cvf::Vec3f( displayVertex2 ) );
                }
            }

            cvf::ref<cvf::Vec3fArray> vertexArray = new cvf::Vec3fArray( displayLines );

            std::vector<cvf::uint> indices;
            indices.reserve( vertexArray->size() );
            for ( cvf::uint k = 0; k < vertexArray->size(); ++k )
            {
                indices.push_back( k );
            }

            cvf::ref<cvf::PrimitiveSetIndexedUInt> indexedUInt = new cvf::PrimitiveSetIndexedUInt( cvf::PrimitiveType::PT_LINES );
            cvf::ref<cvf::UIntArray>               indexArray  = new cvf::UIntArray( indices );
            indexedUInt->setIndices( indexArray.p() );

            cvf::ref<cvf::DrawableGeo> geo = new cvf::DrawableGeo;

            geo->addPrimitiveSet( indexedUInt.p() );
            geo->setVertexArray( vertexArray.p() );
            contourDrawables.push_back( geo );
        }
        if ( !contourDrawables.empty() )
        {
            contourDrawablesForAllLevels[i] = contourDrawables;
        }
    }
    return contourDrawablesForAllLevels;
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::vector<cvf::ref<cvf::Drawable>>
    RivContourMapProjectionPartMgr::createContourLabels( const cvf::Camera*                                           camera,
                                                         const caf::DisplayCoordTransform*                            displayCoordTransform,
                                                         std::vector<std::vector<cvf::BoundingBox>>*                  labelBBoxes,
                                                         const std::vector<RigContourPolygonsTools::ContourPolygons>& contourLinePolygons,
                                                         const RigContourMapGrid&                                     contourMapGrid,
                                                         const cvf::ScalarMapper*                                     scalarMapper,
                                                         RiaNumberFormat::NumberFormatType                            numberFormat,
                                                         int                                                          precision ) const
{
    CVF_ASSERT( camera && displayCoordTransform && labelBBoxes );

    std::vector<cvf::ref<cvf::Drawable>> labelDrawables;
    labelBBoxes->clear();
    labelBBoxes->resize( contourLinePolygons.size() );

    if ( !camera->viewport() || camera->viewport()->width() == 0 ) return labelDrawables;

    if ( scalarMapper == nullptr ) return labelDrawables;

    if ( dynamic_cast<const caf::CategoryMapper*>( scalarMapper ) != nullptr ) return labelDrawables;

    std::vector<double> tickValues;
    scalarMapper->majorTickValues( &tickValues );

    const RigContourPolygonsTools::ContourPolygons* previousLevel = nullptr;
    for ( int64_t i = (int64_t)contourLinePolygons.size() - 1; i > 0; --i )
    {
        cvf::Color3f                backgroundColor( scalarMapper->mapToColor( tickValues[i] ) );
        cvf::Color3f                textColor = RiaColorTools::contrastColor( backgroundColor );
        cvf::ref<cvf::DrawableText> label     = createTextLabel( textColor, backgroundColor );

        for ( size_t j = 0; j < contourLinePolygons[i].size(); ++j )
        {
            if ( contourLinePolygons[i][j].vertices.empty() ) continue;

            QString     qLabelText = RiaNumberFormat::valueToText( contourLinePolygons[i][j].value, numberFormat, precision );
            cvf::String labelText  = cvfqt::Utils::toString( qLabelText );

            size_t nVertices              = contourLinePolygons[i][j].vertices.size();
            size_t nLabels                = nVertices;
            double distanceSinceLastLabel = std::numeric_limits<double>::infinity();
            for ( size_t l = 0; l < nLabels; ++l )
            {
                size_t nVertex    = ( nVertices * l ) / nLabels;
                size_t nextVertex = ( nVertex + 1 ) % nVertices;

                const cvf::Vec3d& localVertex1 = contourLinePolygons[i][j].vertices[nVertex];
                const cvf::Vec3d& localVertex2 = contourLinePolygons[i][j].vertices[nextVertex];

                cvf::Vec3d lineCenter = ( localVertex1 + localVertex2 ) * 0.5;
                double     tolerance  = 1.0e-2 * contourMapGrid.sampleSpacing();

                if ( previousLevel && lineOverlapsWithPreviousContourLevel( lineCenter, *previousLevel, tolerance ) )
                {
                    continue;
                }

                cvf::Vec3d globalVertex1 = localVertex1 + contourMapGrid.origin3d();
                cvf::Vec3d globalVertex2 = localVertex2 + contourMapGrid.origin3d();

                cvf::Vec3d globalVertex = 0.5 * ( globalVertex1 + globalVertex2 );

                cvf::Vec3d segment       = globalVertex2 - globalVertex1;
                cvf::Vec3d displayVertex = displayCoordTransform->transformToDisplayCoord( globalVertex );
                cvf::Vec3d windowVertex;
                camera->project( displayVertex, &windowVertex );
                CVF_ASSERT( !windowVertex.isUndefined() );
                displayVertex.z() += 10.0f;
                cvf::BoundingBox windowBBox = label->textBoundingBox( labelText, cvf::Vec3f::ZERO, cvf::Vec3f( segment.getNormalized() ) );
                cvf::Vec3d       displayBBoxMin, displayBBoxMax;
                camera->unproject( windowBBox.min() + windowVertex, &displayBBoxMin );
                camera->unproject( windowBBox.max() + windowVertex, &displayBBoxMax );

                CVF_ASSERT( !displayBBoxMin.isUndefined() );
                CVF_ASSERT( !displayBBoxMax.isUndefined() );

                cvf::BoundingBox displayBBox( displayBBoxMin - cvf::Vec3d::Z_AXIS * 20.0, displayBBoxMax + cvf::Vec3d::Z_AXIS * 20.0 );

                cvf::Vec3d currentExtent = displayBBoxMax - displayBBoxMin;

                bool overlaps = false;
                if ( distanceSinceLastLabel < currentExtent.length() * 10.0 )
                {
                    overlaps = true;
                }

                if ( !overlaps )
                {
                    for ( auto boxVector : *labelBBoxes )
                    {
                        for ( const cvf::BoundingBox& existingBBox : boxVector )
                        {
                            // Assert on invalid bounding box seen on Linux
                            if ( !displayBBox.isValid() || !existingBBox.isValid() ) continue;

                            double dist = ( displayBBox.center() - existingBBox.center() ).length();
                            if ( dist < segment.length() || existingBBox.intersects( displayBBox ) )
                            {
                                overlaps = true;
                                break;
                            }
                        }
                    }
                }

                if ( !overlaps )
                {
                    cvf::Vec3f displayVertexV( displayVertex );
                    CVF_ASSERT( !displayVertex.isUndefined() );
                    label->addText( labelText, displayVertexV, cvf::Vec3f( segment.getNormalized() ) );
                    labelBBoxes->at( i ).push_back( displayBBox );
                    distanceSinceLastLabel = 0.0;
                }
                else
                {
                    distanceSinceLastLabel += segment.length();
                }
            }
        }
        if ( label->numberOfTexts() != 0u )
        {
            labelDrawables.push_back( label );
        }

        previousLevel = &contourLinePolygons[i];
    }
    return labelDrawables;
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
cvf::ref<cvf::DrawableGeo> RivContourMapProjectionPartMgr::createPickPointVisDrawable( const caf::DisplayCoordTransform* displayCoordTransform,
                                                                                       const cvf::Vec2d&        pickPoint,
                                                                                       const RigContourMapGrid& contourMapGrid ) const
{
    std::vector<cvf::Vec3d> pickPointPolygon = RigContourPolygonsTools::generatePickPointPolygon( pickPoint, contourMapGrid );

    if ( pickPointPolygon.empty() )
    {
        return nullptr;
    }
    cvf::ref<cvf::Vec3fArray> vertexArray = new cvf::Vec3fArray( pickPointPolygon.size() );

    for ( size_t i = 0; i < pickPointPolygon.size(); ++i )
    {
        cvf::Vec3d globalPoint = pickPointPolygon[i] + contourMapGrid.origin3d();
        cvf::Vec3f displayPoint( displayCoordTransform->transformToDisplayCoord( globalPoint ) );
        ( *vertexArray )[i] = displayPoint;
    }

    cvf::ref<cvf::DrawableGeo> geo = nullptr;
    if ( vertexArray->size() > 0u )
    {
        std::vector<cvf::uint> indices;
        indices.reserve( vertexArray->size() );
        for ( cvf::uint j = 0; j < vertexArray->size(); ++j )
        {
            indices.push_back( j );
        }

        cvf::ref<cvf::PrimitiveSetIndexedUInt> indexedUInt = new cvf::PrimitiveSetIndexedUInt( cvf::PrimitiveType::PT_LINES );
        cvf::ref<cvf::UIntArray>               indexArray  = new cvf::UIntArray( indices );
        indexedUInt->setIndices( indexArray.p() );

        geo = new cvf::DrawableGeo;

        geo->addPrimitiveSet( indexedUInt.p() );
        geo->setVertexArray( vertexArray.p() );
    }
    return geo;
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
bool RivContourMapProjectionPartMgr::lineOverlapsWithPreviousContourLevel( const cvf::Vec3d&                               lineCenter,
                                                                           const RigContourPolygonsTools::ContourPolygons& previousLevel,
                                                                           double                                          tolerance )
{
    return RigContourPolygonsTools::lineOverlapsWithContourPolygons( lineCenter, previousLevel, tolerance );
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
cvf::ref<cvf::DrawableGeo> RivContourMapProjectionPartMgr::createXYPlaneQuadGeoWithTexCoords( const cvf::Vec3dArray& cornerPoints )
{
    cvf::ref<cvf::Vec3fArray> vertices = new cvf::Vec3fArray;
    vertices->reserve( 4 );

    for ( const auto& v : cornerPoints )
    {
        vertices->add( cvf::Vec3f( v ) );
    }

    cvf::ref<cvf::Vec2fArray> texCoords = new cvf::Vec2fArray;
    texCoords->reserve( 4 );
    texCoords->add( cvf::Vec2f( 0, 0 ) );
    texCoords->add( cvf::Vec2f( 1, 0 ) );
    texCoords->add( cvf::Vec2f( 1, 1 ) );
    texCoords->add( cvf::Vec2f( 0, 1 ) );

    cvf::ref<cvf::DrawableGeo> geo = new cvf::DrawableGeo;
    geo->setVertexArray( vertices.p() );
    geo->setTextureCoordArray( texCoords.p() );

    cvf::ref<cvf::UIntArray> indices = new cvf::UIntArray;
    indices->reserve( 6 );

    for ( uint i : { 0, 1, 2, 0, 2, 3 } )
    {
        indices->add( i );
    }

    cvf::ref<cvf::PrimitiveSetIndexedUInt> primSet = new cvf::PrimitiveSetIndexedUInt( cvf::PT_TRIANGLES );
    primSet->setIndices( indices.p() );
    geo->addPrimitiveSet( primSet.p() );

    geo->computeNormals();

    return geo;
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
cvf::ref<cvf::Part> RivContourMapProjectionPartMgr::createSingleTexturedQuadPart( const cvf::Vec3dArray&      cornerPoints,
                                                                                  cvf::ref<cvf::TextureImage> image,
                                                                                  bool                        transparent )
{
    cvf::ref<cvf::Part> part = new cvf::Part;

    cvf::ref<cvf::DrawableGeo> geo = createXYPlaneQuadGeoWithTexCoords( cornerPoints );
    geo->computeNormals();

    cvf::ref<cvf::Effect> eff           = new cvf::Effect;
    bool                  surfaceEffect = false;
    if ( surfaceEffect )
    {
        caf::SurfaceEffectGenerator geometryEffgen( cvf::Color3f::YELLOW, caf::PO_1 );
        eff = geometryEffgen.generateCachedEffect();
    }
    else
    {
        cvf::ref<cvf::Texture> texture = new cvf::Texture( image.p() );
        cvf::ref<cvf::Sampler> sampler = new cvf::Sampler;
        sampler->setMinFilter( cvf::Sampler::LINEAR );
        sampler->setMagFilter( cvf::Sampler::NEAREST );
        sampler->setWrapModeS( cvf::Sampler::CLAMP_TO_EDGE );
        sampler->setWrapModeT( cvf::Sampler::CLAMP_TO_EDGE );

        cvf::ref<cvf::RenderStateTextureBindings> textureBindings = new cvf::RenderStateTextureBindings;
        textureBindings->addBinding( texture.p(), sampler.p(), "u_texture2D" );

        eff->setRenderState( textureBindings.p() );
        eff->setShaderProgram( m_textureShaderProg.p() );

        if ( transparent )
        {
            part->setPriority( RivPartPriority::PartType::TransparentSeismic );
            cvf::ref<cvf::RenderStateBlending> blending = new cvf::RenderStateBlending;
            blending->configureTransparencyBlending();
            eff->setRenderState( blending.p() );
        }
    }
Dual Rendering Toggle

Both textured image and triangle rendering can be enabled simultaneously; confirm intended combined output and performance impact. If mutually exclusive is desired, enforce UI or logic to avoid duplicating the same data render.

 
if ( m_contourMapProjection->showTrianglesWithColor() )
{
    m_contourMapProjectionPartMgr->appendProjectionToModel( contourMapProjectionModelBasicList.p(),
                                                            transForm.p(),
                                                            m_contourMapProjection->trianglesWithVertexValues(),
                                                            *m_contourMapProjection->mapGrid(),
                                                            backgroundColor(),
                                                            m_contourMapProjection->legendConfig()->scalarMapper() );
}

if ( m_contourMapProjection->showImage() )
{
    /*
                    auto    image    = RicCreateContourMapPolygonTools::createImage( m_contourMapProjection->mapProjection(),
                                                                               m_contourMapProjection->legendConfig()->scalarMapper()
       ); QString fileName = "f:/scratch/contour-map.png"; image.save( fileName );
    */

    m_contourMapProjectionPartMgr->appendProjectionAsTexturedQuad( contourMapProjectionModelBasicList.p(),
                                                                   transForm.p(),
                                                                   m_contourMapProjection->legendConfig()->scalarMapper(),
                                                                   *m_contourMapProjection->mapProjection(),
                                                                   *m_contourMapProjection->mapGrid() );
}

@qodo-code-review
Copy link

qodo-code-review bot commented Aug 29, 2025

PR Code Suggestions ✨

Explore these optional code suggestions:

CategorySuggestion                                                                                                                                    Impact
Possible issue
Validate indices and values

Guard against out-of-range indices and non-finite values to prevent memory
access and NaN color mapping issues. Skip or clamp invalid indices and treat
non-finite values as fully transparent.

ApplicationLibCode/ModelVisualization/RivContourMapProjectionPartMgr.cpp [784-803]

 auto dataValues = contourMapProjection->aggregatedVertexResultsFiltered();
+const size_t valueCount = dataValues.size();
 
 for ( int y = 0; y < height; ++y )
 {
     for ( int x = 0; x < width; ++x )
     {
-        size_t index = contourMapProjection->vertexIndex( x, y );
-
-        double valueAtVertex = dataValues[index];
-        auto   color         = scalarMapper->mapToColor( valueAtVertex );
-
-        int transparency = 0;
-        if ( valueAtVertex != std::numeric_limits<double>::infinity() )
+        const size_t index = contourMapProjection->vertexIndex( x, y );
+        if ( index >= valueCount )
         {
-            transparency = 255;
+            textureImage->setPixel( x, y, cvf::Color4ub( 0, 0, 0, 0 ) );
+            continue;
         }
 
-        textureImage->setPixel( x, y, cvf::Color4ub( color, transparency ) );
+        const double valueAtVertex = dataValues[index];
+        if ( !std::isfinite( valueAtVertex ) )
+        {
+            textureImage->setPixel( x, y, cvf::Color4ub( 0, 0, 0, 0 ) );
+            continue;
+        }
+
+        const auto color = scalarMapper->mapToColor( valueAtVertex );
+        textureImage->setPixel( x, y, cvf::Color4ub( color, 255 ) );
     }
 }
  • Apply / Chat
Suggestion importance[1-10]: 8

__

Why: The suggestion correctly points out that the existing code does not check for out-of-bounds access on dataValues or handle non-finite values, which could lead to crashes or undefined behavior.

Medium
Use stable Z offset

Avoid scaling Z by itself; if Z is near zero this collapses the offset, and if
negative it may move the quad in the wrong direction. Use a fixed positive
offset relative to grid thickness or a small epsilon to ensure consistent
layering above the surface.

ApplicationLibCode/ModelVisualization/RivContourMapProjectionPartMgr.cpp [139-140]

 double zValueForContourMap = contourMapGrid.origin3d().z();
-zValueForContourMap += 0.3 * zValueForContourMap;
+const double zOffset       = std::max(1e-3, 0.001 * std::abs(contourMapGrid.expandedBoundingBox().size().z()));
+zValueForContourMap += zOffset;
  • Apply / Chat
Suggestion importance[1-10]: 7

__

Why: The suggestion correctly identifies that scaling the z-offset by zValueForContourMap can lead to issues if the value is zero or negative, and proposes a more robust method to ensure a consistent positive offset.

Medium
General
Improve texture filtering

Enable proper minification filtering with mipmaps to avoid aliasing and
shimmering when the quad is viewed at a distance. Generate mipmaps on the
texture and use trilinear filtering for better visual quality.

ApplicationLibCode/ModelVisualization/RivContourMapProjectionPartMgr.cpp [709-714]

 cvf::ref<cvf::Texture> texture = new cvf::Texture( image.p() );
+texture->setGenerateMipMap( true );
 cvf::ref<cvf::Sampler> sampler = new cvf::Sampler;
-sampler->setMinFilter( cvf::Sampler::LINEAR );
-sampler->setMagFilter( cvf::Sampler::NEAREST );
+sampler->setMinFilter( cvf::Sampler::LINEAR_MIPMAP_LINEAR );
+sampler->setMagFilter( cvf::Sampler::LINEAR );
 sampler->setWrapModeS( cvf::Sampler::CLAMP_TO_EDGE );
 sampler->setWrapModeT( cvf::Sampler::CLAMP_TO_EDGE );
  • Apply / Chat
Suggestion importance[1-10]: 5

__

Why: The suggestion correctly proposes using mipmapping to improve texture rendering quality and avoid aliasing, which is a good practice for visual fidelity, though not a critical bug fix.

Low
Add terminating newline

Ensure there is a terminating newline at the end of the file to comply with
common tooling and POSIX expectations. This avoids diff noise and issues with
license scanners.

ApplicationLibCode/Adm/LicenseInformation.txt [780-784]

+IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
 
+===============================================================================
+ Notice for surfio
+===============================================================================
  • Apply / Chat
Suggestion importance[1-10]: 2

__

Why: The suggestion to add a terminating newline is a valid convention, but its impact is minor, and the provided code snippets are incorrect and do not reflect the proposed change at the end of the file.

Low
  • More

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Reviewers

No reviews

Assignees

No one assigned

Labels

Review effort 3/5

Projects

None yet

Milestone

No milestone

Development

Successfully merging this pull request may close these issues.

1 participant

@magnesj