# Error opening .vtm files in pyqt6 gui
**URL:** https://discourse.vtk.org/t/error-opening-vtm-files-in-pyqt6-gui/15087
**Category:** Support
**Tags:** python
**Created:** [December 29, 2024, 8:39am UTC](https://discourse.vtk.org/t/error-opening-vtm-files-in-pyqt6-gui/15087 "2024-12-29T08:39:29Z")
**Posts on this page:** 4
**Page:** 1
### Author:  [@Jimmyenricohu](https://discourse.vtk.org/u/Jimmyenricohu)
#### Post date: [December 29, 2024, 8:39am UTC](https://discourse.vtk.org/t/error-opening-vtm-files-in-pyqt6-gui/15087/1 "2024-12-29T08:39:29Z")
I have builded a program for CFD results interpolation by vtk interpolation functions.
It run successfully in jupyter, but while I write a Pyqt6 GUI and packaged by pyinstaller it show the error like below.
I don‘t now if similar problems have occurred before.
---
### Author:  [@mwestphal](https://discourse.vtk.org/u/mwestphal)
#### Post date: [December 30, 2024, 9:04am UTC](https://discourse.vtk.org/t/error-opening-vtm-files-in-pyqt6-gui/15087/2 "2024-12-30T09:04:53Z")
please share your data
---
### Author:  [@Jimmyenricohu](https://discourse.vtk.org/u/Jimmyenricohu)
#### Post date: [January 7, 2025, 1:46pm UTC](https://discourse.vtk.org/t/error-opening-vtm-files-in-pyqt6-gui/15087/3 "2025-01-07T13:46:24Z")
here is my vtk interpolation class
import vtkmodules.all as vtk
import numpy as np
import pandas as pd
import csv
from pandas.core.frame import DataFrame
class sindainterp:
def **init** (self,sourcefilename,sindafilename):
self.sourcefie = sourcefilename[0]
self.Time = sourcefilename[1]
self.Type = sourcefilename[2]
self.sindafie = sindafilename
self.CaseNum = len(self.Time)
if (len(self.Time) == len(self.Type))and (len(self.Type) == len(self.sourcefie)):
print(‘Import Finished!’)
else:
print(‘Please Check Imported Data!’)
```
def read_SourceGrid(self,filename):
reader = vtk.vtkXMLMultiBlockDataReader()
reader.SetFileName(filename)
reader.Update()
# 获取第一个块（只输出单一壁面时只有一个块）
block = reader.GetOutput()
grid = block.GetBlock(0)
c2p = vtk.vtkCellDataToPointData()
c2p.SetInputData(grid)
c2p.Update()
sourcegrid = c2p.GetOutput()
return sourcegrid
def readGrid_fluent(self,fluentfile):
pointsdf = pd.read_csv(fluentfile)
points = pointsdf.iloc[:,1:4].values
temp = pointsdf.iloc[:,4].values
vtk_points = vtk.vtkPoints()
for p in points:
# print(p)
vtk_points.InsertNextPoint(p)
grid = vtk.vtkUnstructuredGrid()
grid.SetPoints(vtk_points)
scalars = vtk.vtkFloatArray()
scalars.SetName("TEMPERATURE")
for i in range(temp.shape[0]):
scalars.InsertNextValue(temp[i])
grid.GetPointData().SetScalars(scalars)
multiblock = vtk.vtkMultiBlockDataSet()
multiblock.SetBlock(0, grid)
sourcegrid = multiblock.GetBlock(0)
return sourcegrid
def read_SindaGrid(self,filename):
Points = []
Elements = []
with open(filename, 'r') as f:
reader = csv.reader(f)
for row in reader:
if row[0] == 'GRID':
Points.append(row)
elif row[0] == 'CTRIA3':
Elements.append(row)
Points = DataFrame(Points,columns=['Type','Node id','Tab', 'x', 'y', 'z']) #.iloc[:,1:].values()
Elements = DataFrame(Elements,columns=['Type','Element id','Tab', 'N1', 'N2', 'N3']) ##.iloc[:,1:].values()
Points['Node id'] = Points['Node id'].astype(int)
Points = Points.sort_values(by = ['Node id'], axis=0, ascending=True, inplace=False)
self.SindaPoints = Points.loc[:,['x','y','z']].astype(float).values
self.SindaElements = Elements.loc[:,['N1','N2','N3']].astype(int).values
vtk_points = vtk.vtkPoints()
for p in self.SindaPoints:
# print(p)
vtk_points.InsertNextPoint(np.array(p))
Sindagrid = vtk.vtkPolyData()
Sindagrid.SetPoints(vtk_points)
# print(Sindagrid)
return Sindagrid
def interpolation(self,kernel):
'''
插值核函数：vtk.vtkShepardKernel()
'''
interpolator = vtk.vtkPointInterpolator()
interpolator.SetInputData(self.sindagrid)
interpolator.SetSourceData(self.sourcegrid)
interpolator.SetKernel(kernel)
interpolator.Update()
self.sindagrid = interpolator.GetOutput()
self.sindagridpoints = self.sindagrid.GetPoints()
self.sindagridpoints_data = self.sindagrid.GetPointData()
# print(self.sindagridpoints)
def Gensindainpfile(self,outputname):
self.outputname = outputname
self.sindagrid = self.read_SindaGrid(self.sindafie)
with open(outputname, 'w', newline='') as f:
f.write('DATA: TEMPERATURE\n')
f.write('UNITS LENGTH meters\n')
f.write('UNITS TIME SECONDS\n')
f.write('UNITS TEMPERATURE K\n')
writer = csv.writer(f,delimiter='\t')
## write Nodes
for i in range(self.SindaPoints.shape[0]):
point = self.SindaPoints[i,:]
f.write('Node\t'+str(i+1)+'\t')
f.write(str(point[0]) + '\t'+str(point[1]) + '\t'+str(point[2]) + '\n')
for i in range(self.SindaElements.shape[0]):
simp = self.SindaElements[i]
simp.sort()
f.write('TRI\t'+str(i+1)+'\t')
f.write(str(simp[0]) + '\t'+str(simp[1]) + '\t'+str(simp[2]) + '\n')
for k in range(len(self.Time)):
if self.Type[k] == 'VTM':
print('case:',k,'VTM')
self.sindagrid = self.read_SindaGrid(self.sindafie)
self.sourcegrid = self.read_SourceGrid(self.sourcefie[k])
# print(self.sourcegrid)
elif self.Type[k] == 'FLUENT':
self.sindagrid = self.read_SindaGrid(self.sindafie)
self.sourcegrid = self.readGrid_fluent(self.sourcefie[k])
print('case:',k,'FLUENT')
# print(self.sourcegrid)
self.interpolation(vtk.vtkShepardKernel())
points = self.sindagridpoints
points_data = self.sindagridpoints_data
print(points_data.GetNumberOfArrays())
f.write('TIME'+' '+str(self.Time[k])+'\n')
for i in range(points.GetNumberOfPoints()):
# if self.Type[k] == 'FLUENT':
# print('i:',i)
for j in range(points_data.GetNumberOfArrays()):
# if self.Type[k] == 'FLUENT':
# print('j:',j)
array = points_data.GetArray(j)
value = array.GetComponent(i, 0) # 假设每个点数据只有一个分量
message = tuple([value]) +tuple(['! Temperature Value Node '])+ tuple([i+1])
writer.writerow(message)
print('Export file ',outputname, ' successfully')
```
---
### Author:  [@Jimmyenricohu](https://discourse.vtk.org/u/Jimmyenricohu)
#### Post date: [January 7, 2025, 1:46pm UTC](https://discourse.vtk.org/t/error-opening-vtm-files-in-pyqt6-gui/15087/4 "2025-01-07T13:46:56Z")
here is the QT GUI code
from PyQt6.QtWidgets import QApplication, QMainWindow, QTableWidgetItem, QMessageBox
from sindainterpolation import Ui\_MainWindow # 导入由pyuic6生成的Python文件
import VTK2SINDA
import tkinter as tk
from tkinter import filedialog
from pathlib import Path
import pandas as pd
class MainWindow(QMainWindow, Ui\_MainWindow):
def **init** (self):
super(). **init** ()
self.setupUi(self)
```
self.getPathButtom.clicked.connect(self.getEngineerPath)
self.TableAdd.clicked.connect(self.TableAddRow)
self.TableDelete.clicked.connect(self.TableDelRow)
self.ImportEnvelopeButtom.clicked.connect(self.readEnvelopeTable)
self.InterpButtom.clicked.connect(self.interpolation)
self.EnvelopeTable.setSortingEnabled(True)
self.msg_box = QMessageBox()
def getEngineerPath(self): # 浏览按钮被点击
# 创建一个Tk根窗口，但不显示
root = tk.Tk()
root.withdraw()
# 打开文件选择对话框，并获取用户选择的文件路径
file_path = filedialog.askopenfilename(
title="选择网格文件",
initialdir="/", # 起始目录
filetypes=[("nastran mesh files", "*.bdf")] # 文件类型过滤
)
# 返回选择的文件路径
directory,filename = self.get_file_path_and_split(file_path)
self.MeshFileName.setText(filename)
self.DirPath.setText(directory)
def get_file_path_and_split(self,file_path):
# 创建Path对象
path = Path(file_path)
# 拆分文件名和路径
directory = path.parent # 获取目录路径
filename = path.name # 获取文件名
return str(directory), filename
def TableAddRow(self):
self.EnvelopeTable.insertRow(self.EnvelopeTable.rowCount())
def TableDelRow(self):
# 获取选中的行号
selected_index = self.EnvelopeTable.selectionModel().selectedRows()
model = self.EnvelopeTable.model()
for index in selected_index:
model.removeRow(index.row())
def readEnvelopeTable(self):
root = tk.Tk()
root.withdraw()
if self.DirPath.toPlainText() is None:
dir = '/'
else:
dir = self.DirPath.toPlainText()
# 打开文件选择对话框，并获取用户选择的文件路径
file_path = filedialog.askopenfilename(
title="选择包线csv",
initialdir= dir, # 起始目录
filetypes=[("包线配置文件", "*.csv")] # 文件类型过滤
)
if file_path == '':
self.msg_box.setWindowTitle("Message")
self.msg_box.setText("未选择合适文件")
return self.msg_box.exec()
data = pd.read_csv(file_path)
self.sourcefile = [data.iloc[:, 0].tolist(), data.iloc[:, 1].tolist(), data.iloc[:, 2].tolist()]
self.EnvelopeTable.clear()
self.EnvelopeTable.setRowCount(data.shape[0])
self.EnvelopeTable.setColumnCount(data.shape[1])
for j in range(data.shape[1]):
for i,cell in enumerate(self.sourcefile[j]):
item = QTableWidgetItem(str(cell))
self.EnvelopeTable.setItem(i, j, item)
def interpolation(self):
self.sindafile = self.DirPath.toPlainText() + '\\' + self.MeshFileName.toPlainText()
self.outputname = self.DirPath.toPlainText() + '\\' + self.OutputFileName.toPlainText()
interpolator = VTK2SINDA.sindainterp(self.sourcefile, self.sindafile)
interpolator.Gensindainpfile(self.outputname)
self.msg_box.setWindowTitle("Message")
self.msg_box.setText("导出成功")
self.msg_box.exec()
```
if **name** == “ **main** ”:
app = QApplication()
window = MainWindow()
window.show()
app.exec()