Python---VTK---vtkplotter之load()、loadImageData()

Python---VTK---vtkplotter之load()、loadImageData()

1、load()

（1）原型：

def load(inputobj, c=None, alpha=1, threshold=False, spacing=(), unpack=True):

（2）Example

:Examples:
    .. code-block:: python

        from vtkplotter import datadir, load, show

        # Return a Mesh object
        g = load(datadir+'250.vtk')
        show(g)

        # Return a list of 2 meshes
        g = load([datadir+'250.vtk', datadir+'270.vtk'])
        show(g)

        # Return a list of meshes by reading all files in a directory
        # (if directory contains DICOM files then a Volume is returned)
        g = load(datadir+'timecourse1d/')
        show(g)

        # Return a Volume. Color/Opacity transfer functions can be specified too.
        g = load(datadir+'embryo.slc')
        g.c(['y','lb','w']).alpha((0.0, 0.4, 0.9, 1))
        show(g)

        # Return a Mesh from a SLC volume with automatic thresholding
        g = load(datadir+'embryo.slc', threshold=True)
        show(g)
"""

from vtkplotter import vtkio,load,show,datadir

    g=vtkio.load(r"C:\Users\ZhangRR\Desktop\CorArtSeg\Rotterdam\raw\CASD13\transform\dataset08\17AHA.vtk")

show(g)

 

 2、loadImageData()

(1) 原型

def loadImageData(filename, spacing=()):

 

（2） Example 

    from vtkplotter import vtkio

    pio = PathIO()
    pth = pio.ReadFile(r'..\SimVascularTools\seg12.pth')

    ds_dir = r'D:\__Dataset__\centerline\Rotterdam\transform\dataset08'
    vol_dir = pjoin(ds_dir, 'dataset08.nii.gz')
    pth_dir = r'..\SimVascularTools\seg12.pth'

    cmpr = CurvedMPR(vol_dir)
    cmpr.add_path_from_file(pth_dir)

    StrenchedVol = cmpr.GetCurvedMPR(size=10., spacing=0.1)

    save_demo(StrenchedVol)


    # 使用vtkio
    Image = vtkio.loadImageData(vol_dir)
    #Image.SetOrigin() 此處VTK有小bug，需要手動設置origin
    reslice = GetSlicevtkImage(pth.GetPathPoints()[0], Image, 10., 0.1)
    # 變成numpy數據
    reslice = vtkImageData2Numpy(reslice)
    show_image(reslice)