OpenCASCADE Color Scale

最新推荐文章于 2024-09-08 13:09:05 发布

weixin_34082854

最新推荐文章于 2024-09-08 13:09:05 发布

阅读量250

点赞数

文章标签： c/c++

原文链接：http://www.cnblogs.com/opencascade/p/3983651.html

版权

OpenCASCADE Color Scale

eryar@163.com

Abstract. The color scale is a specialized label object that displays a color map and an accompanying numerical scale for color mapped or contour data plots. As the geometry modeling kernel of SALOME, OpenCASCADE provides the color scale function. The paper focus on the usage of color scale in OpenCASCADE.

Key Words. OpenCASCADE, Color Scale, 应力云图, 颜色映射表

1. Introduction

应力云图是一种应用广泛的标量场可视化方法，也经常用于表达矢量场或张量场的一个分量。云图在颜色与标量数据之间建立某种确定的映射关系，在计算机屏幕上绘制一个颜色离散变化的图像来显示科学计算结果。云图将标量场定义域内的数据值映射成不同的颜色，从而通过颜色变化来反映数据场中数据的变化规律。

很多数值分析软件后处理的结果都有绘制云图的功能。如下图所示为Abaqus软件产生的云图：

Figure 1.1 A Colored Stress Patterns by Abaqus

OpenCASCADE作为SALOME开源数值分析软件的几何造型内核，也提供了云图可视化的功能。本文主要介绍如何使用OpenCASCADE来显示模型的应力云图中的颜色映射表(Color Scale)。

2. Color Scale

云图绘制依赖于颜色集合与标量场数值集合之间的一一映射，即颜色映射表(Color Scale)。颜色映射表为区域填充时建立数值与颜色的映射关系。如下图所示：

Figure 2.1 Color Scale

作为分析对比参照标准，用户可以根据颜色线性表对应的数值范围，来判断某区域内数值分布规律。颜色线性表不仅应用于云图绘制，也应用于其他计算可视化算法中。

颜色线性表可以定义成不同的形式，多以对比鲜明的色彩作为线性表段颜色，在两个对比色之间采用过渡颜色。

在OpenCASCADE中颜色映射表的绘制是由Viewer来实现的。其Tcl命令为vcolorscale，如下所示：

Figure 2.2 Color Scale Tcl Command: vcolorscale

默认的颜色映射表的显示效果如下图所示：

Figure 2.3 Color Scale in Draw Test Harness

3. Tcl Test

OpenCASCADE基于Tcl/Tk的Draw Test Harness环境很方便测试一些想法。现在在Draw Test Harness中来显示一个完整的云图，如下图所示：

Figure 3.1 Color Scale in Draw Test Harness

实现上图的Tcl脚本代码如下所示：

#
#    Copyright (c) 2014 eryar All Rights Reserved.
#
#        File    : colorscale.tcl
#        Author  : eryar@163.com
#        Date    : 2014-09-20 18:10
#        Version : 1.0v
#
#    Description : Demonstrate the usage of OpenCASCADE color scale.
#
pload ALL

meshfromstl m data/stl/head.stl

meshcolors m nodaltex 0

# show the color sacle.
vcolorscale

vfit

首先加载所需要的所有模块，再从head.stl中加载网格模型；设置网格顶点颜色后就用vcolorscale命令打开发颜色映射表。

4.Code Analysis

根据Tcl命令找到对应的C++实现代码如下所示：

//=============================================================================
//function : VColorScale
//purpose  : representation color scale
//=============================================================================
#include <V3d_ColorScale.hxx>

static int VColorScale (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
{
  if ( argc != 1 && argc != 4 && argc != 5 && argc != 6 && argc != 8 )
  {
    di << "Usage : " << argv[0] << " [RangeMin = 0 RangeMax = 100 Intervals = 10 HeightFont = 16 Position = Right X = 0 Y = 0]  " << "\n";
    return 1;
  }

  Handle(AIS_InteractiveContext) aContext = ViewerTest::GetAISContext();
  if(aContext.IsNull()) {
    di << argv[0] << " ERROR : use 'vinit' command before " << "\n";
    return -1;
  }

  Standard_Real minRange = 0. , maxRange = 100. ;

  Standard_Integer numIntervals = 10 ;
  Standard_Integer textHeight = 16;
  Aspect_TypeOfColorScalePosition position = Aspect_TOCSP_RIGHT;
  Standard_Real X = 0., Y = 0. ;

  if ( argc < 9 )
  {
     if( argc > 3 )
     {
       minRange = Draw::Atof( argv[1] );
       maxRange = Draw::Atof( argv[2] );
       numIntervals = Draw::Atoi( argv[3] );
     }
     if ( argc > 4 )
       textHeight = Draw::Atoi( argv[4] );
     if ( argc > 5 )
       position = (Aspect_TypeOfColorScalePosition)Draw::Atoi( argv[5] );
     if ( argc > 7 )
     {
       X = Draw::Atof( argv[6] );
       Y = Draw::Atof( argv[7] );
     }
  }
  Handle(V3d_View) curView = ViewerTest::CurrentView( );
  if ( curView.IsNull( ) )
    return 1;
  Handle(Aspect_ColorScale) aCSV = curView->ColorScale( );
  Handle(V3d_ColorScale) aCS = ( Handle( V3d_ColorScale )::DownCast( aCSV ) );
  if( ! aCS.IsNull( ) )
  {
    aCS->SetPosition( X , Y );
    aCS->SetHeight( 0.95) ;
    aCS->SetTextHeight( textHeight );
    aCS->SetRange( minRange , maxRange );
    aCS->SetNumberOfIntervals( numIntervals );
    aCS->SetLabelPosition( position );
    if( !curView->ColorScaleIsDisplayed() )
      curView->ColorScaleDisplay( );
  }
  return 0;
}