epifilename = 'r100815_epi' + '.stl'
epipdata = vtk_py.readSTL(epifilename)
#stop
rotatedepipdata = vtk_py.rotatePData_w_axis(epipdata, angle, raxis)
rotatedendopdata = vtk_py.rotatePData_w_axis(endopdata, angle, raxis)
height = min(rotatedendopdata.GetBounds()[5],
rotatedepipdata.GetBounds()[5]) - zoffset
clipped_endo, clipped_epi = vtk_py.clipSurfacesForCutLVMesh(rotatedendopdata,
rotatedepipdata,
height,
verbose=True)
vtk_py.writeSTL(clipped_epi, "clipped_epi_tmp.stl")
vtk_py.writeSTL(clipped_endo, "clipped_endo_tmp.stl")
filename = 'New_mesh'
vtk_py.createLVmesh(
filename, 0.3, "clipped_epi_tmp.stl", "clipped_endo_tmp.stl"
) #(filename, 0.45, "clipped_epi_tmp.stl", "clipped_endo_tmp.stl") #
#os.remove("clipped_epi_tmp.stl")
#os.remove("clipped_endo_tmp.stl")
#os.remove("LVtemp.geo.bak")
ugridfilename = filename + '.vtk'
endopdata = vtk_py.readSTL(endofilename)
#epifilename = 'EPI-LV' + '.stl'
epipdata = vtk_py.readSTL(epifilename)
rotatedepipdata = vtk_py.rotatePData_w_axis(epipdata, angle, raxis)
rotatedendopdata = vtk_py.rotatePData_w_axis(endopdata, angle, raxis)
height = min(rotatedendopdata.GetBounds()[5], rotatedepipdata.GetBounds()[5]) - zoffset
clipped_endo, clipped_epi = vtk_py.clipSurfacesForCutLVMesh(rotatedendopdata, rotatedepipdata, height, verbose=True)
#vtk_py.writeSTL(clipped_epi, "clipped_epi_tmp.stl")
vtk_py.writeSTL(clipped_epi, output_dir + "clipped_epi_tmp.stl")
#vtk_py.writeSTL(clipped_endo, "clipped_endo_tmp.stl")
vtk_py.writeSTL(clipped_endo, output_dir + "clipped_endo_tmp.stl")
#filename = 'New_mesh'
filename = output_info["mesh_name"]
vtk_py.createLVmesh(filename, 0.5, "clipped_epi_tmp.stl", "clipped_endo_tmp.stl")#(filename, 0.45, "clipped_epi_tmp.stl", "clipped_endo_tmp.stl") #
#os.remove("clipped_epi_tmp.stl")
#os.remove("clipped_endo_tmp.stl")
# Info for clipping top of volumes #Laxis = np.array([-0.0760398272351652, -0.0945530105122101, 0.992611541781136]) Laxis = ([0, 0, 1]) #Laxis = clip_info["Laxis"] #zoffset = clip_info["z_offset"] zoffset = 0.0 ################################################################### endofilename = 'endo_stl_scaled.stl' epifilename = 'epi_stl_scaled.stl' angle = np.arccos(Laxis[2]) raxis = np.cross(Laxis, np.array([0, 0, 1])) raxis = raxis / np.linalg.norm(raxis) #endofilename = 'ENDO-LV' + '.stl' """endopdata = vtk_py.readSTL(endofilename) #epifilename = 'EPI-LV' + '.stl' epipdata = vtk_py.readSTL(epifilename) rotatedepipdata = vtk_py.rotatePData_w_axis(epipdata, angle, raxis) rotatedendopdata = vtk_py.rotatePData_w_axis(endopdata, angle, raxis) height = min(rotatedendopdata.GetBounds()[5], rotatedepipdata.GetBounds()[5]) - zoffset clipped_endo, clipped_epi = vtk_py.clipSurfacesForCutLVMesh(rotatedendopdata, rotatedepipdata, height, verbose=True) #vtk_py.writeSTL(clipped_epi, "clipped_epi_tmp.stl") vtk_py.writeSTL(clipped_epi,"clipped_epi_tmp.stl")