Math and Annotations

Some examples require external input files. Before you start, please follow the link in the Example Scripts section to download the zip file with model and result files.

Example 01 - Curves from Files, Datums, Notes, Legends and PNG, XYDATA Export

The code loads two curves in an XY plot using the CurveXY class, applies SAE filtering and vertical trimming.

Two notes are added by the Note class, dynamically positioned at the maximum of the curves using a Templex expression to get the index of the curve maximum.

In addition, two vertical datum lines are added by the DatumVertical class, positioned similar to the notes at the maximum of the curves using a Templex expression.

The area between the two curves is shaded using the curve attribute shadeStyle = “betweenCurves” on the second curve.

Finally, a JPEG image in a given width and height is exported using the CaptureImageTool class, in addition with a xy data ascii file containing the curve values via the ExportCurves class.

import hw
import hw.hg as hg
import os

ALTAIR_HOME = os.path.abspath(os.environ["ALTAIR_HOME"])
plotFile = os.path.join(
    ALTAIR_HOME, "demos", "mv_hv_hg", "plotting", "madymo", "ANGACC"
)
pngFile = os.path.join(
    os.path.abspath(os.path.dirname(__file__)), "HG_Image_Export.png"
)
xyDataFile = os.path.join(
    os.path.abspath(os.path.dirname(__file__)), "HG_Curve_Export.xydata"
)

# New session
ses = hw.Session()
ses.new()
win = ses.get(hw.Window)
win.type = "plot"

# Define image and curve export tools
captureImage = hw.CaptureImageTool(
    file=pngFile,
    type="png",
    width=2000,
    height=1000,
)
exportCurves = hg.ExportCurves(file=xyDataFile, format="xyData", range="all")

# Layout and channel settings
colorList = [(255, 0, 0), (0, 0, 255)]
requestList = ["50th% Hybrid3   - LOWER TORSO", "50th% Hybrid3   - UPPER TORSO"]
autoPosList = ["right", "left"]

# Loop over curves, notes and datums
for color, request, autoPos in zip(colorList, requestList, autoPosList):

    # Add curves
    cu = hg.CurveXY(
        xFile=plotFile,
        xSource="file",
        xDataType="Time",
        xRequest="Time",
        xComponent="Time",
        yFile=plotFile,
        ySource="file",
        yDataType="Angular Acceleration",
        yRequest=request,
        yComponent="Res. ang. acc.",
        lineColor=color,
        label="{y.HWRequest} - {y.HWComponent}",
    )

    # Add curve math
    cu.mathAdd("SAE", padding=1, direction=1, filterClass=15)
    cu.mathAdd("trim", start=20, end=130)

    cuId = str(cu.id)

    # Add note attached to curve math
    no = hg.Note(
        attachTo="curve",
        attachmentCurve=cu.id,
        attachmentPoint=f'indexofmax(c{cuId}.y)',
        text="Max = {y} at {x}",
        autoPosition=True,
        autoPositionStrategy=autoPos,
        curveColor=True,
    )

    # Add vertical datum
    dv = hg.DatumVertical(
        position=f'c{cuId}.x[indexofmax(c{cuId}.y)]',
        label="Max " + cu.yRequest,
        labelLocation="upperCenter",
        lineColor=color,
    )

# Modify legend
legend = ses.get(hg.Legend)
legFont = hw.Font(size=14)
legend.placement = "insideRightBottomCorner"
legend.font = legFont

# Shade area between 2 curves
cu.setAttributes(
    shadeArea=True,
    shadeStyle="betweenCurves",
    shadeSecondCurve=1,
    shadeColor="Blue",
    shadeAlpha=0.2,
)

win.update()

# Export image and modified curve values in XYDATA format
captureImage.capture()
exportCurves.export()

Figure 1. Exported Image of script ‘Curves from Files, Datums, Notes, Legends and PNG, XYDATA Export’

Example 02 - evalHWC() for Curve Import and Math Operations

The main focus of this example is to execute slightly modified recorded HWC commands using the the evalHWC() wrapper class in Python.

The script executes the following steps:

1. Starting a new session and switching to HyperGraph

2. Plotting two curves using the “Create Curves by File” GUI in a two window side-by-side layout

3. Applying Stackmath operations for “trim”, “SAE” and “resample” to all curves in all windows

import hw
import hw.hg as hg
import os

ALTAIR_HOME = os.path.abspath(os.environ["ALTAIR_HOME"])
plotFile = os.path.join(
    ALTAIR_HOME, "demos", "mv_hv_hg", "plotting", "madymo", "ANGACC"
).replace('\\','/')

# New session
hw.evalHWC("delete session")
# Set Client Type to HyperGraph
hw.evalHWC('hwd window type="HyperGraph 2D"')
# Load curves with build plots module
hw.evalHWC("xy option layout destination=oneplotperrequest")
hw.evalHWC("xy option layout layout=1")
hw.evalHWC(
    f'xy load file="{plotFile}" ydatatype= "Angular Acceleration" \
    yrequest= "50th% Hybrid3   - LOWER TORSO" , "50th% Hybrid3   - UPPER TORSO" ycomponent= all'
)
hw.evalHWC(
    'xy curve stackadd range="w:all i:all" operation=trim start=20 end=130 interpolate=false'
)
# Add curve math for all curves in current plot
hw.evalHWC(
    'xy curve stackadd range="w:all i:all" operation=SAE padding=1 samples=10 direction=3 filterClass=60'
)
hw.evalHWC(
    'xy curve stackadd range="w:all i:all" operation=resample type="Resample Fixed Points" fixnum=100.0'
)

Figure 2. Output of ‘evalHWC() for Curve Import and Math Operations’