我有一个 fortran 代码来生成二进制 VTK 格式的网格。此代码生成一个二进制 VTK 文件,如下所示:
# vtk DataFile Version 3.0
vtk output
BINARY
DATASET RECTILINEAR_GRID
DIMENSIONS 2 2 1
X_COORDINATES 2 float
^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@
Y_COORDINATES 2 float
^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@
Z_COORDINATES 1 float
^@^@^@^@^@^@^@^@
当我尝试使用 ParaView 打开它时,它崩溃并显示以下错误消息:
错误:在
/home/user/OpenFOAM/ThirdParty-2.3.0/ParaView-4.1.0/VTK/IO/Legacy/vtkRectilinearGridReader.cxx,
第 311 行 vtkRecti LinearGridReader (0x379f4b0):无法识别的关键字:
...
如果我用 ASCII 编写上述文件,那么它可以在 ParaView 中正确打开,这是打开的 ASCII 类似文件:
# vtk DataFile Version 3.0
vtk output
ASCII
DATASET RECTILINEAR_GRID
DIMENSIONS 2 2 1
X_COORDINATES 2 float
0 1
Y_COORDINATES 2 float
0 1
Z_COORDINATES 1 float
0
我使用以下代码生成二进制 VTK 格式的网格(我给出了最小工作代码):
program VTKBinary
implicit none
real*8 :: x(2) = (0., 1.)
real*8 :: y(2) = (0., 1.)
real*8 :: z(1) = (0.)
character :: buffer*80, lf*1, str1*8, str2*8, str3*8
integer :: ivtk = 9, int
lf = char(10) ! line feed character
open(unit=ivtk,file='test_bin.vtk',form='binary',convert='BIG_ENDIAN')
buffer = '# vtk DataFile Version 3.0'//lf ; write(ivtk) trim(buffer)
buffer = 'vtk output'//lf ; write(ivtk) trim(buffer)
buffer = 'BINARY'//lf ; write(ivtk) trim(buffer)
buffer = 'DATASET RECTILINEAR_GRID'//lf ; write(ivtk) trim(buffer)
! WRITE GRID
write(str1(1:8),'(i8)') size(x)
write(str2(1:8),'(i8)') size(y)
write(str3(1:8),'(i8)') size(z)
buffer = 'DIMENSIONS '//str1//str2//str3//lf ; write(ivtk) trim(buffer)
buffer = 'X_COORDINATES '//str1//' float'//lf ; write(ivtk) trim(buffer)
write(ivtk) x
buffer = lf//'Y_COORDINATES '//str2//' float'//lf ; write(ivtk) trim(buffer)
write(ivtk) y
buffer = lf//'Z_COORDINATES '//str3//' float'//lf ; write(ivtk) trim(buffer)
write(ivtk) z
close(ivtk)
end program VTKBinary
二进制 VTK 文件有什么问题?为什么它退出?