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 ...
随机推荐
- Delphi多线程下的ADO编程
前言: 几个月前接到一个任务:将一后台程序访问数据库的方式从BDE改为ADO,原因是由于业务量的增加,通过BDE不论是向数据库写入数据还是从数据库中读出数据的速度都变得无法忍受,大家都知道ADO在数据 ...
- Windows 10开发基础——指针事件和操作事件(一)
主要内容: 1.指针事件 2.操作事件 1.指针事件 指针事件由各种活动输入源引发,包括触摸.触摸板.笔和鼠标(它们替代传统的鼠标事件).指针事件基于单一输入点(手指.笔尖.鼠标光标),但不支持基于速 ...
- Java Date Calendar DateFormat Details
From https://www.ntu.edu.sg/home/ehchua/programming/java/DateTimeCalendar.html Date and Time - Creat ...
- <%@ Application Codebehind="Global.asax.cs" Inherits="XXX.MvcApplication" Language="C#" %>
<%@ Application Codebehind="Global.asax.cs" Inherits="XXX.MvcApplication" Lan ...
- SetForegroundWindow API函数还不够(好多好多解决方案,真是奇思妙想)
好多好多解决方案: var Input: TInput; begin ZeroMemory(@Input, SizeOf(Input)); SendInput(, Input, SizeOf(Inpu ...
- Ruby元编程:执行某个目录下的全部测试用例
目前手里有个测试项目各个feature的测试用例都放在对应的子目录下,虽然有自动化测试框架的帮助执行起来很方便,但是偶尔也有需要在本地执行某个feature的全部测试用例集合.因为本人对shell脚本 ...
- 一、OpenScenGraph环境搭建
1.OpenSceneGraph 3.4.0 网址 http://www.openscenegraph.org/ github源码地址 https://github.com/openscenegr ...
- SYN6109型 NTP网络子钟
SYN6109型 NTP网络子钟 产品概述 SYN6109型NTP网络子钟是由西安同步电子科技有限公司精心设计.自行研发生产的一套通过网口与母钟连接的子钟,接收母钟发送来的时间信息(信息内容:年.月. ...
- vue实现Excel文件的上传与下载
一.前言项目中使用到比较多的关于Excel的前端上传与下载,整理出来,以便后续使用或分析他人. 1.前端vue:模板下载与导入Excel 导入Excel封装了子组件,点击导入按钮可调用子组件,打开文件 ...
- hgoi#20190516
T1-Buying A House 给你一个长度为n的序列a,给你目标房子m,最多花的钱k 如果a[i]为0,这座房子无法购买,否则可以购买,求能买的距离目标房子最近的房子,输出最小距离 两座房子之间 ...