[Paraview] Slices generation with pvpython

[Utkarsh Ayachit]

It's not clear what's doing on with out a sample dataset to try this
out with. Can you share a dataset that will help reproduce the
problem?

Utkarsh

On Wed, Apr 8, 2015 at 4:47 AM, Borchi Leonardo
<Leonardo.Borchi at ducati.com> wrote:
> Hi everybody,
>
>
>
> I have a problem when I tried to implement an automated procedure for
> postprocessing analysis.
>
> I have wrote and run (from command prompt) the following script with
> pvpython.exe.
>
>
>
> If the variable numberofslice is <= 20 the script run well, but if I try to
> generate 40 slides for example, after the twenty one, for a couple of slide,
> the image is blank. After that, the following slides are generated in a
> correct manner.
>
>
>
> I have tried to change the py with different type of algorithm ( for, while,
> if ) but the problem remain.
>
>
>
> Have you any ideas for resolve this strange issue??
>
>
>
> Thanks for your support
>
>
>
> Leonardo
>
>
>
> #Import the simple module from the paraview
>
> from paraview.simple import *
>
>
>
> #Disable automatic camera reset on 'Show'
>
> paraview.simple._DisableFirstRenderCameraReset()
>
>
>
> #Create a new 'OpenFOAMReader'
>
> case_foam = OpenFOAMReader(FileName='\\\\SERVERCORSE\\RUN_051\\case.foam')
>
> (xmin,xmax,ymin,ymax,zmin,zmax) = case_foam.GetDataInformation().GetBounds()
>
> xmax = float(2.000)
>
> xmin = float(-0.200)
>
> ymax = float(0.250)
>
> ymin = float(-0.250)
>
> zmax = float(1.350)
>
> zmin = float(0.000)
>
> numberofslice = 20
>
> deltax = (xmax-xmin)/numberofslice
>
> deltay = (ymax-ymin)/numberofslice
>
> deltaz = (zmax-zmin)/numberofslice
>
> case_foam.CellArrays = ['U', 'p']
>
> case_foam.MeshRegions = ['internalMesh']
>
>
>
> #Use node discretization instead of cell one.
>
> pointData = CellDatatoPointData(Input=case_foam)
>
>
>
> #Calculate Cp and Cg coefficients
>
> Cp = Calculator(Input=pointData)
>
> Cp.Function = 'p/(0.5*50*50)'
>
> Cp.ResultArrayName = 'Cp'
>
> Cg = Calculator(Input=Cp)
>
> Cg.Function = '(p+(0.5*(U.U)))/(0.5*50*50)'
>
> Cg.ResultArrayName = 'Cg'
>
>
>
> renderView_ = CreateRenderView()
>
> renderView_.ViewSize = [1920,1080]
>
> renderView_.InteractionMode ='2D'
>
>
>
> # Z SLICE
>
> #renderView_.CameraPosition = [1.104649161719515, 0.1483969716170493,
> 70.80596909030284]
>
> #renderView_.CameraFocalPoint = [1.104649161719515, 0.1483969716170493,
> 0.30000001192092896]
>
> #renderView_.CameraViewUp = [0.0, 1.0, 0.0]
>
> #renderView_.CameraParallelScale = 0.8642834397172128
>
>
>
> # Y SLICE
>
> #renderView_.CameraPosition = [0.76, -70.50, 0.85]
>
> #renderView_.CameraFocalPoint = [0.76, 0.0, 0.85]
>
> #renderView_.CameraViewUp = [0.0, 0.0, 1.0]
>
> #renderView_.CameraParallelScale = 1.25
>
>
>
> # X SLICE
>
> renderView_.CameraPosition = [16.89230484541337, 0.022500131777517626,
> 0.8954779729169607]
>
> renderView_.CameraFocalPoint = [0.5, 0.022500131777517626,
> 0.8954779729169607]
>
> renderView_.CameraViewUp = [0.0, 0.0, 1.0]
>
> renderView_.CameraParallelScale = 0.9233222262058975
>
>
>
> cgCM_ = GetColorTransferFunction('Cg')
>
> cgCM_.RGBPoints = [-1.00, 0.0, 0.0, 1.0, 1.20, 1.0, 0.0, 0.0]
>
> cgCM_.ColorSpace = 'HSV'
>
> cgCM_.NanColor = [0.498039, 0.498039, 0.498039]
>
> cgCM_.ScalarRangeInitialized = 1.0
>
>
>
> cgOM_ = GetOpacityTransferFunction('Cg')
>
> cgOM_.Points = [-1.00, -0.5, 0.0, 0.5, 1.0]
>
> cgOM_.ScalarRangeInitialized = 1
>
>
>
> cgCMColorBar_ = GetScalarBar(cgCM_, renderView_)
>
>
>
> # Z SLICE
>
> #cgCMColorBar_.Position = [0.2924509803921572, 0.8569873380156294]
>
> #cgCMColorBar_.Position2 = [0.4299999999999994, 0.1200000000000001]
>
>
>
> # Y SLICE
>
> #cgCMColorBar_.Position = [0.30, 0.025]
>
> #cgCMColorBar_.Position2 = [0.45, 0.025]
>
>
>
> # X SLICE
>
> cgCMColorBar_.Position = [0.2675490196078431, 0.8524428726877044]
>
> cgCMColorBar_.Position2 = [0.42999999999999955, 0.11999999999999988]
>
>
>
> cgCMColorBar_.Orientation = 'Horizontal'
>
> cgCMColorBar_.Title = 'Cg'
>
> cgCMColorBar_.ComponentTitle = ''
>
> cgCMColorBar_.TextPosition = 'Ticks left/bottom, annotations right/top'
>
> cgCMColorBar_.RangeLabelFormat = '%4.3f'
>
> cgCMColorBar_.NumberOfLabels = 2
>
> cgCMColorBar_.TitleBold = 1
>
> cgCMColorBar_.TitleItalic = 1
>
> cgCMColorBar_.LabelFontFamily = 'Courier'
>
> cgCMColorBar_.LabelBold = 1
>
> cgCM_.LockScalarRange = 1
>
> cgCM_.RescaleTransferFunction(-0.3, 1.2)
>
> cgOM_.RescaleTransferFunction(-0.3, 1.2)
>
>
>
> coord = xmin
>
> iteraz = 0
>
>
>
> while coord < xmax:
>
>
>
>                          slice_ = Slice(Input=Cg)
>
>                          slice_.SliceType = 'Plane'
>
>                          slice_.SliceType.Origin = [coord, 0.000, 0.000]
>
> #                       slice_.SliceType.Origin = [0.000, coord, 0.000]
>
> #                       slice_.SliceType.Origin = [0.000, 0.000, coord]
>
>                          slice_.SliceType.Normal = [1.0, 0.0, 0.0]
>
> #                       slice_.SliceType.Normal = [0.0, 1.0, 0.0]
>
> #                       slice_.SliceType.Normal = [0.0, 0.0, 1.0]
>
>                          slice_.SliceOffsetValues = [0.0]
>
>
>
>                          slice_display_ = Show(slice_, renderView_)
>
>                          slice_display_.ColorArrayName = ['POINTS', 'Cg']
>
>                          cgCM_.RescaleTransferFunction(-0.3, 1.2)
>
>                          cgOM_.RescaleTransferFunction(-0.3, 1.2)
>
>                          slice_display_.LookupTable = cgCM_
>
>                          slice_display_.SetScalarBarVisibility(renderView_,
> True)
>
>
>
>                          name_file =
> '\\\\SERVER\\RUN_051\\postProcessing\\00%4.0d_img_%4.3f.png' % (iteraz,
> coord)
>
>
>
>                          Render(view=renderView_)
>
>
>
>                          SaveScreenshot(name_file, magnification=1,
> quality=100, view=renderView_)
>
>
>
>                          coord = coord + deltax
>
>                          iteraz = iteraz + 1
>
>
>
>                          Delete(slice_)
>
>                          Delete(slice_display_)
>
>
>
>                          del slice_
>
>                          del slice_display_
>
>
>
> Delete(renderView_)
>
> del renderView_
>
>
> _______________________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Please keep messages on-topic and check the ParaView Wiki at:
> http://paraview.org/Wiki/ParaView
>
> Search the list archives at: http://markmail.org/search/?q=ParaView
>
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/mailman/listinfo/paraview
>