plot3d网格读取写入与可视化
说明
plot3d格式是NASA制定并大量使用的CFD网格文件格式,在CFD编程过程中经常涉及到。本文利用Python语言编写一个读取plot3d文件,写入plot3d文件的通用类,并利用vtk完成网格的显示。
对于程序的说明
- 使用说明
在主函数中,写入如下代码:
# 新建对象
mesh = Mesh()
# 读文件
fname = 'demo.fmt'
mesh.read(fname)
fname1 = "1111.dat"
mesh.write(fname1)
- 程序思路
定义Mesh类,包含有:
属性:
1. 网格维度: DIM = 3
2. 网格块数目: int BlockNum
3. 网格块坐标数目:np.array([]) Block_list
4. 网格块列表 list(Block)
方法:
1. setBlockNum: 设置网格块数目
2. setBlockList: 设置网格块数目
3. MeshInit: 创建mesh对象,给BlockNum和Block_list赋值
4. setBlock: 设置每一个块的坐标和IBLANK
每个Block具有:
1. id,jd,kd: 三个方向上坐标数目
2. X,Y,Z: 三个方向的坐标分量(一维列阵)
3. IBlank:当前块的IBlank值
源码
# !/usr/bin/env python
# -*- coding:utf-8 -*-
import math
import numpy as np
class Block():
"""
描述Plot3d文件中的一个Block
"""
def __init__(self, ID, JD, KD):
self.id = ID
self.jd = JD
self.kd = KD
self.X = np.zeros(self.id)
self.Y = np.zeros(self.jd)
self.Z = np.zeros(self.kd)
self.IBlank = np.zeros(self.id * self.jd * self.kd)
def setX(self, x_list):
self.X = x_list
def setY(self, y_list):
self.Y = y_list
def setZ(self, z_list):
self.Z = z_list
def setBlank(self, blank_list):
self.IBlank = blank_list
def setXYZ(self, x_list, y_list, z_list):
self.X = x_list
self.Y = y_list
self.Z = z_list
class Mesh():
"""
描述一个plot3d网格文件
"""
def __init__(self):
self.DIM = 3
self.blockNum = 1
self.block_list = []
self.BLOCK = []
def setBlockNum(self, block_num):
self.blockNum = block_num
def setBlockList(self, block_list):
self.block_list = block_list
def MeshInit(self, block_num, block_list):
self.setBlockNum(block_num)
self.setBlockList(block_list)
for i in range(self.blockNum):
bl = Block(*self.block_list[3 * i:3 * (i + 1)])
self.BLOCK.append(bl)
def setBlock(self, i, coordinates_list, iblank_list):
self.BLOCK[i - 1].setXYZ(*coordinates_list)
self.BLOCK[i - 1].setBlank(iblank_list)
def check(self):
if len(self.block_list) != self.blockNum * self.DIM:
print("The Block Num and the Block list is not consistent.")
for i in range(self.blockNum):
block = self.BLOCK[i]
ID = block.id
JD = block.jd
KD = block.kd
pointNum = ID * JD * KD
flag = len(block.X) == pointNum and len(block.Y) == pointNum
flag = flag and len(block.Z) == pointNum
if not flag:
print("The point number is not consistent.")
def write(self, fname):
self.check()
data = ''
try:
with open(fname, 'w') as f:
# 写入Block数目
line = '\t{:>8d}\n'
data += line.format(self.blockNum)
# f.write(line.format(self.blockNum))
line = '\t'
for i in range(len(self.block_list)):
line = line + '{' + str(i) + ':>8d}\t'
if (i + 1) % 6 == 0:
line += '\n\t'
line += '\n'
data += line.format(*self.block_list)
# f.write(line.format(*self.block_list))
for block in self.BLOCK:
# 循环写入X,Y,Z坐标
for xi in [block.X, block.Y, block.Z]:
line = ' '
for i in range(len(xi)):
line = line + '{' + str(i) + ': 7E} '
if (i + 1) % 5 == 0:
line += '\n '
line += '\n'
data += line.format(*xi)
# f.write(line.format(*xi))
# 写入iBlank
line = '\t'
for i in range(len(block.IBlank)):
line = line + '{' + str(i) + ': 9d}' + '\t'
if (i + 1) % 6 == 0:
line += '\n\t'
line += '\n'
data += line.format(*block.IBlank)
# f.write(line.format(*block.IBlank))
f.write(data)
except (BaseException):
print("写出文件时发生错误!请检查您的输出数据结构。")
exit(1)
def read(self, fname):
try:
with open(fname, 'r') as f:
# 读取块号
line = f.readline().strip()
self.blockNum = int(line)
# 读取block_list
num = math.ceil((self.DIM * self.blockNum) / 6)
line = ''
for i in range(num):
line += f.readline()
line = line.split()
block_list = np.array([int(i) for i in line])
self.block_list = block_list
# 初始化Mesh
self.MeshInit(self.blockNum, self.block_list)
# 读取每个块信息
for ibl in range(self.blockNum):
bl = self.BLOCK[ibl]
# 获取block_list信息
bl.id, bl.jd, bl.kd = self.block_list[self.DIM *
ibl:self.DIM *
(ibl + 1)]
point_num = bl.id * bl.jd * bl.kd
num = math.ceil(point_num / 5)
# 获得x
line = ''
for i in range(num):
line += f.readline()
line = line.split()
x_list = np.array([float(im) for im in line])
# 获得y
line = ''
for i in range(num):
line += f.readline()
line = line.split()
y_list = np.array([float(im) for im in line])
# 获得z
line = ''
for i in range(num):
line += f.readline()
line = line.split()
z_list = np.array([float(im) for im in line])
# 获得iblank
line = ''
num = math.ceil(bl.id * bl.jd * bl.kd / 6)
for i in range(num):
line += f.readline()
line = line.split()
iblank_list = np.array([int(im) for im in line])
self.setBlock(ibl, [x_list, y_list, z_list], iblank_list)
except (BaseException):
pass
self.check()
def visual(self):
"""
显示网格信息
"""
pass
if __name__ == '__main__':
mesh = Mesh()
# 读文件
fname = 'demo.fmt'
mesh.read(fname)
fname1 = "1111.dat"
mesh.write(fname1)
plot3d网格读取写入与可视化的更多相关文章
- OpenMesh 读写网格控制(读取写入纹理坐标,法向等)
OpenMesh读取网格默认是不自动读取obj网格中的法向,纹理坐标等信息的,写入网格同样也是.所以要读取(或写入)这些信息需要修改默认的选项. 先看一下其读写网格的函数 template<cl ...
- 用多态来实现U盘,Mp3,移动硬盘和电脑的对接,读取写入数据。
using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.T ...
- Java里的IO流里的 ObjectInputStream 的读取\写入!
各位好!!我又来了!!今天遇见了一个小问题!! IO流里的对象读取总是出错!各种报错!!神烦啊!!百思不得其解啊!然后就上网百度!找了好久终于让我找到了!下面就让我来说一说! 当你用IO流里的对象流写 ...
- 关于Python中读取写入文件并进行文件与用户交互的操作
一.提前知识点 在python中是同样和其他语言一样可以进行文件的读取写入操作,值得注意的是,Python中打开文件读取的方式有几种,分别是以下几种: f = open('username.txt') ...
- JXL读取,写入Excel
JXL读取,写入Excel2003 相关阅读:poi 读写excel2003:http://www.cnblogs.com/gavinYang/p/3576739.htmlpoi 读写excel200 ...
- J2SE 8的输入输出--读取/写入文本文件和读取/写入二进制数据
读取/写入文本文件 // 1. 文本输入 // (1) 短小文本直接转入字符串 String string = new String(Files.readAllBytes(Paths.get(&quo ...
- Golang 读取写入Etcd
http://www.yunweipai.com/archives/8131.html 新一代etcd:etcd3 etcd是一个高可用的 Key/Value 存储系统,主要用于分享配置和服务发现.简 ...
- Java poi读取,写入Excel2003
Java poi读取,写入Excel2003 相关阅读:poi读写Excel2007:http://www.cnblogs.com/gavinYang/p/3576741.htmljxl读写excel ...
- Java poi读取,写入Excel2007
Java poi读取,写入Excel2007 相关阅读:poi读写Excel2003:http://www.cnblogs.com/gavinYang/p/3576739.htmljxl读写excel ...
随机推荐
- GetDateFormat和GetTimeFormat两个API
https://msdn.microsoft.com/en-us/library/windows/desktop/dd318086(v=vs.85).aspx
- LIBCMTD.lib(exe_winmain.obj) : error LNK2019: 无法解析的外部符号 _WinMain@16,该符号在函数 "int __cdecl invoke_main(void)" (?invoke_main@@YAHXZ) 中被引用
这个问题是没找到程序入口 在网上查这个问题,一般都是说两条: (若是win32程序) 一是在项目属性\CC++\预处理器\预处理器定义\里添加 _WINDOWS 一是在项目属性\链接\系统 里选择 窗 ...
- LINQ查询表达式---------group子句
LINQ查询表达式---------group子句 LINQ表达式必须以from子句开头,以select或group子句结束.使用guoup子句来返回元素分组后的结果.group 子句返回一个 IGr ...
- Android零基础入门第25节:最简单最常用的LinearLayout线性布局
原文:Android零基础入门第25节:最简单最常用的LinearLayout线性布局 良好的布局设计对于UI界面至关重要,在前面也简单介绍过,目前Android中的布局主要有6种,创建的布局文件默认 ...
- Android Studio 添加 Genymotion插件
原文:Android Studio 添加 Genymotion插件 1.下载Genymotion:官网地址,必须先注册才能下载,下载带有VirtualBox的版本 2.安装:安装时会连VirtualB ...
- UBUNTU 16.04 + CUDA8.0 + CUDNN6.0 + OPENCV3.2 + MKL +CAFFE + tensorflow
首先说一下自己机子的配置 CPU:Intel(R) Core(TM) i5-5600 CUP @3.20GHz *4 GPU : GTX 1060 OS : 64bit Ubuntu16.04LTS ...
- MySQL InnoDB缓冲池(Buffer Pool)
InnoDB缓冲池并不仅仅缓存索引,它还会缓存行的数据.自适应哈希索引.插入缓冲(Insert Buffer).锁,以及其他内部数据结构. InnoDB还使用缓冲池来帮助延迟写入,这样就能合并多个写入 ...
- 今天想安装 windowsl ive 提示安装失败 错误码
Windows Live installation error: OnCatalogResult: 0x80072ee6 看了了这个老兄的回答,试了试,果然OK,谢谢@普洛提亚从这里下载安装包,然后安 ...
- C++ 王者归来:对编程语言的需求总结为四个:效率,灵活,抽象,生产率(C++玩的是前三个,Java和C#玩的是后两个)
Why C++ ? 王者归来(转载) 因为又有人邀请我去Quora的C2C网站去回答问题去了,这回是 关于 @laiyonghao 的这篇有点争议的博文<2012 不宜进入的三个技术点>A ...
- Windows服务器下的IIS和Apache性能比较
目前最流行的建立网站的服务工具就要属Apache与IIS了.那么他们之间到底哪个性能更好呢?到底哪个工具才是最适合我们的呢?最近我也对这方面的问题进行了一番研究. 如果是基于Linux平台的话,那不必 ...