package require vtk package require vtkinteraction # create pipeline # vtkPLOT3DReader pl3d pl3d SetXYZFileName "$VTK_DATA_ROOT/Data/combxyz.bin" pl3d SetQFileName "$VTK_DATA_ROOT/Data/combq.bin" pl3d SetScalarFunctionNumber 100 pl3d SetVectorFunctionNumber 202 pl3d Update # create three line probes vtkLineSource line line SetResolution 30 vtkTransform transL1 transL1 Translate 3.7 0.0 28.37 transL1 Scale 5 5 5 transL1 RotateY 90 vtkTransformPolyDataFilter tf tf SetInput [line GetOutput] tf SetTransform transL1 vtkProbeFilter probe probe SetInput [tf GetOutput] probe SetSource [pl3d GetOutput] vtkTransform transL2 transL2 Translate 9.2 0.0 31.20 transL2 Scale 5 5 5 transL2 RotateY 90 vtkTransformPolyDataFilter tf2 tf2 SetInput [line GetOutput] tf2 SetTransform transL2 vtkProbeFilter probe2 probe2 SetInput [tf2 GetOutput] probe2 SetSource [pl3d GetOutput] vtkTransform transL3 transL3 Translate 13.27 0.0 33.40 transL3 Scale 4.5 4.5 4.5 transL3 RotateY 90 vtkTransformPolyDataFilter tf3 tf3 SetInput [line GetOutput] tf3 SetTransform transL3 vtkProbeFilter probe3 probe3 SetInput [tf3 GetOutput] probe3 SetSource [pl3d GetOutput] vtkAppendPolyData appendF appendF AddInput [probe GetPolyDataOutput] appendF AddInput [probe2 GetPolyDataOutput] appendF AddInput [probe3 GetPolyDataOutput] vtkTubeFilter tuber tuber SetInput [appendF GetOutput] tuber SetRadius 0.1 vtkPolyDataMapper lineMapper lineMapper SetInput [tuber GetOutput] vtkActor lineActor lineActor SetMapper lineMapper # probe the line and plot it vtkGlyphSource2D triangle triangle SetGlyphTypeToTriangle vtkGlyphSource2D cross cross SetGlyphTypeToCross vtkXYPlotActor xyplot xyplot AddInput [probe GetOutput] xyplot AddInput [probe2 GetOutput] xyplot AddInput [probe3 GetOutput] [xyplot GetPositionCoordinate] SetValue 0.0 0.67 0 [xyplot GetPosition2Coordinate] SetValue 1.0 0.33 0;#relative to Position xyplot SetXValuesToArcLength xyplot SetNumberOfXLabels 6 xyplot SetTitle "Pressure vs. Arc Length (Zoomed View)" xyplot SetXTitle "" xyplot SetYTitle "P" xyplot SetXRange .1 .35 xyplot SetYRange .2 .4 [xyplot GetProperty] SetColor 0 0 0 xyplot PlotLinesOn [xyplot GetProperty] SetLineWidth 1 xyplot PlotPointsOn [xyplot GetProperty] SetPointSize 3 xyplot LegendOn xyplot SetPlotSymbol 2 [triangle GetOutput] xyplot SetPlotColor 2 0 0 1 xyplot SetGlyphSize 0.025 vtkSphereSource vertexGlyph vtkXYPlotActor xyplot2 xyplot2 AddInput [probe GetOutput] xyplot2 AddInput [probe2 GetOutput] xyplot2 AddInput [probe3 GetOutput] [xyplot2 GetPositionCoordinate] SetValue 0.00 0.33 0 [xyplot2 GetPosition2Coordinate] SetValue 1.0 0.33 0;#relative to Position xyplot2 SetXValuesToNormalizedArcLength xyplot2 SetNumberOfXLabels 6 xyplot2 SetXRange 0.2 1.0 xyplot2 SetTitle "VTK Quality vs. Alcohol Consumption" xyplot2 SetXTitle "" xyplot2 SetYTitle "A" xyplot2 PlotPointsOn xyplot2 PlotLinesOff xyplot2 LegendOn xyplot2 SetLegendPosition 0.4 0.6 xyplot2 SetLegendPosition2 0.40 0.25 [xyplot2 GetProperty] SetColor 1 0 0 [xyplot2 GetProperty] SetPointSize 2 xyplot2 SetPlotSymbol 0 [vertexGlyph GetOutput] xyplot2 SetPlotLabel 0 "Ken's Mudslide Consumption" xyplot2 SetPlotSymbol 1 [cross GetOutput] xyplot2 SetPlotColor 1 1 0 1 xyplot2 SetPlotLabel 1 "Bill's Beer Consumption" xyplot2 SetPlotSymbol 2 [triangle GetOutput] xyplot2 SetPlotColor 2 0 0 1 xyplot2 SetPlotLabel 2 "VTK Quality" #Okay exercise data object stuff vtkDataSetToDataObjectFilter ds2do ds2do SetInput [probe GetOutput] vtkDataSetToDataObjectFilter ds2do2 ds2do2 SetInput [probe GetOutput] vtkDataSetToDataObjectFilter ds2do3 ds2do3 SetInput [probe GetOutput] vtkXYPlotActor xyplot3 xyplot3 AddDataObjectInput [ds2do GetOutput] xyplot3 SetDataObjectXComponent 0 2 xyplot3 SetDataObjectYComponent 0 5 xyplot3 SetPlotColor 0 1 0 0 xyplot3 SetPlotLabel 0 "Mx" xyplot3 AddDataObjectInput [ds2do2 GetOutput] xyplot3 SetDataObjectXComponent 1 2 xyplot3 SetDataObjectYComponent 1 6 xyplot3 SetPlotColor 1 0 1 0 xyplot3 SetPlotLabel 1 "My" xyplot3 AddDataObjectInput [ds2do3 GetOutput] xyplot3 SetDataObjectXComponent 2 2 xyplot3 SetDataObjectYComponent 2 7 xyplot3 SetPlotColor 2 0 0 1 xyplot3 SetPlotLabel 2 "Mz" [xyplot3 GetPositionCoordinate] SetValue 0.0 0.0 0 [xyplot3 GetPosition2Coordinate] SetValue 1.0 0.33 0;#relative to Position xyplot3 SetXValuesToIndex xyplot3 SetNumberOfXLabels 6 xyplot3 SetTitle "Momentum Component vs. Point Id" xyplot3 SetXTitle "Point Id" xyplot3 SetYTitle "M" xyplot3 PlotPointsOn [xyplot3 GetProperty] SetColor 0 0 1 [xyplot3 GetProperty] SetPointSize 3 xyplot3 LegendOn xyplot3 SetLegendPosition 0.8 0.35 xyplot3 SetLegendPosition2 0.20 0.20 # draw an outline vtkStructuredGridOutlineFilter outline outline SetInput [pl3d GetOutput] vtkPolyDataMapper outlineMapper outlineMapper SetInput [outline GetOutput] vtkActor outlineActor outlineActor SetMapper outlineMapper [outlineActor GetProperty] SetColor 0 0 0 # Create graphics stuff # vtkRenderer ren1 vtkRenderer ren2 vtkRenderWindow renWin renWin AddRenderer ren1 renWin AddRenderer ren2 vtkRenderWindowInteractor iren iren SetRenderWindow renWin ren1 SetBackground 0.6784 0.8471 0.9020 ren1 SetViewport 0 0 .5 1 ren1 AddActor outlineActor ren1 AddActor lineActor ren2 SetBackground 1 1 1 ren2 SetViewport 0.5 0.0 1.0 1.0 ren2 AddActor2D xyplot ren2 AddActor2D xyplot2 ren2 AddActor2D xyplot3 renWin SetSize 500 250 set cam1 [ren1 GetActiveCamera] $cam1 SetClippingRange 3.95297 100 $cam1 SetFocalPoint 8.88908 0.595038 29.3342 $cam1 SetPosition -12.3332 31.7479 41.2387 $cam1 SetViewUp 0.060772 -0.319905 0.945498 iren Initialize renWin Render # render the image # iren AddObserver UserEvent {wm deiconify .vtkInteract} # prevent the tk window from showing up then start the event loop wm withdraw .