使用Houdini输出四面体网格并输出tetgen格式

news2024/11/28 4:53:22

我们的目标是从houdini输出生成的四面体，希望是tetgen格式的。

众所周知，houdini是不能直接输出四面体的。

有三方案去解决：

  1. 输出点云ply文件，然后利用tetgen生成网格。
  2. 输出Hounidi内置的.geo格式文件，然后写个脚本去解析json，因为这个文件就是个json。
  3. 直接从Houdini中利用Python节点输出tetgen格式。

我探索并试验了以上所有三种方案。优缺点如下：

第一种方案的缺点是四面体是不可控的，因为是tetgen现生成的。

第二种方案的缺点是你要保证geo里面没有多余的数据。例如还存储的颜色或uv等信息，或者houdini的prim上还存了其他信息，就会导致解析失败。当然我们可以使用clean节点清除掉这些多余信息。我会把脚本放到文末。

第三种方案是最好的。可控性最好。因此后面我们会说这种方法。

tetgen的数据格式

我们首先要介绍一下tetgen的数据格式：

我们需要其中三种：
.node结尾代表顶点位置
.ele结尾代表四面体编号
.face结尾代表三角形编号
在这里插入图片描述
如图，其中第一个数代表点/面/单元的数量

ele和face中的顶点编号对应的都是node中点前面那个编号。

例如ele中每一行为：
当前单元的编号 四面体第1个顶点的编号 四面体第2个顶点的编号 四面体第3个顶点的编号 四面体第4个顶点的编号

Houdini中使用python节点写出

在这里插入图片描述

各个节点的作用如下所示：

  • file: 读入obj
  • clean: 清理多余信息（非必要）
  • polyreduce: 简化模型，缩小顶点数（非必要）
  • tetconform: 生成四面体
  • split: 分出表面三角形和四面体
  • null：无作用的节点，只是为了占位
  • python node(write_faces): 写出表面三角形编号
  • python node(write_nodes):写出所有顶点位置
  • python node(write_eles)：写出四面体编号

脚本的内容在附录。

在这里插入图片描述

在这里，我们要讲解一下houdini中的数据。

分为四种：

  • points 几何点的信息，包括点的位置
  • vertices: 拓扑顶点信息。例如可以存顶点的编号
  • primitives: 图元信息。可以存储例如该四面体的体积大小等
  • detail: 整个几何体的信息。

一个常见的误区是混淆points和vertices。points完全是空间中真实存在的一个点。具有位置速度等信息。但是vertices可以认为是对顶点编号的reference。例如一个正方体的角点，可以被三个面同时共享。他都是同一个几何点，但是却有三个不同的vertices归属于不同的面。这样的好处是保证了唯一性：一个point只对应一个vertex，一个vertex只被一个primitive所包含。

这里要注意的是，vertices中存的是什么完全取决于图元是什么。假如是个四面体，就可以是四个顶点编号。假如是三角形，就是三角形三个点编号。

在这里插入图片描述

正是由于三角形与四面体都被存在primitive中，所以我们才用split将其分开。方面后面输出。

注意在tetconform中勾选add surface triangles才会输出三角面。

在这里插入图片描述

最后稍微讲解下python节点中的脚本。

请注意，houdini中的python是完全面向对象的，因此万物皆为对象。

我们仅以write_nodes为例。API请查阅Houdini的官方文档。

import hou # houdini包
geo = hou.pwd().geometry() #当前python节点的第一个输入端口所对应的对象，是个SOP对象
print(geo)

import os
path = hou.hipFile.path()  #hip文件所在的位置
path = os.path.dirname(path) + "/models/bunny2000_try.node"
print(f"path is {path}")

pts = geo.points() # 获取SOP对象上的point对象列表（是个list）
numpts = len(pts)
print("numpts:", numpts)

f = open(path, 'w')
f.write(str(numpts)+"  3  0  0\n")
for i in range(numpts):
    pt = pts[i] 
    pos = pt.position() #获取point对象的position属性，就是位置
    f.write("   "+str(i)+"    "+str(pos[0])+"  "+str(pos[1])+"  "+str(pos[2])+"\n")
f.close()

完毕。

（非必要内容）在taichi中ggui显示

请见
https://github.com/chunleili/learn-meshtaichi
中的tut03

其中主要就是多写了directly_import_surf这个函数而已。

网格文件也请见这里。

结果如图

在这里插入图片描述

附录1：Houdini中python节点的脚本内容

write_faces

import hou
geo = hou.pwd().geometry()
print(geo)

import os
path = hou.hipFile.path() 
path = os.path.dirname(path) + "/models/bunny2000_try.face"
print(f"path is {path}")


# write surface triangles
tris = geo.prims()
num_tris = len(tris)
f = open(path[:-4]+"face", 'w')
f.write(str(num_tris)+"  0\n")
for i in range(num_tris):
    tri = tris[i].points()
    f.write("   "+str(i)+"    "+str(tri[0].number())+" "+str(tri[1].number())+" "+str(tri[2].number())+ "  -1" +"\n")
f.close()

write_nodes

import hou
geo = hou.pwd().geometry()
print(geo)

import os
path = hou.hipFile.path() 
path = os.path.dirname(path) + "/models/bunny2000_try.node"
print(f"path is {path}")

pts = geo.points()
numpts = len(pts)
print("numpts:", numpts)

f = open(path, 'w')
f.write(str(numpts)+"  3  0  0\n")
for i in range(numpts):
    pt = pts[i]
    pos = pt.position()
    f.write("   "+str(i)+"    "+str(pos[0])+"  "+str(pos[1])+"  "+str(pos[2])+"\n")
f.close()

write_eles

import hou
geo = hou.pwd().geometry()
print(geo)

import os
path = hou.hipFile.path() 
path = os.path.dirname(path) + "/models/bunny2000_try.ele"
print(f"path is {path}")


eles = geo.prims()
num_eles = len(eles)
print("num_eles:", num_eles)
f1 = open(path[:-4]+".ele", 'w')
f1.write(str(num_eles)+"  4  0\n")
for i in range(num_eles):
    ele = eles[i].points()
    f1.write("   "+str(i)+"    "+str(ele[0].number())+"  "+str(ele[1].number())+"  "+str(ele[2].number())+"  "+str(ele[3].number())+"\n")
f1.close()

附录2： houdni的geo文件解析转换为tetgen格式四面体的脚本

import os
import json

def read_geo(from_path):
    with open(from_path,'r') as f:
        data=json.load(f)

    # 读取顶点个数等信息
    pointcount=data[5]  # 点个数
    vertexcount=data[7] 
    primitivecount=data[9] # 四面体个数

    # 读取四面体的索引
    topology = data[13]
    pointref = topology[1]
    tet_indices = pointref[1]# 四面体的索引，是一个一维数组

    # 读取顶点的位置
    attributes = data[15]
    pointattributes = attributes[1]
    positions = pointattributes[0][1][7][5]

    return tet_indices,positions, pointcount,vertexcount,primitivecount


def write_tetgen(tet_indices,positions, pointcount, primitivecount,to_path, gen_face=False):
    # 写入tetgen的node文件(也就是顶点的位置)
    node_file = to_path+".node"
    if(os.path.exists(node_file)):
        print("remove file: "+node_file)
        os.remove(node_file)
    with open(node_file,'w') as f:
        f.write(str(pointcount)+"  3  0  0\n")
        for i in range(pointcount):
            f.write("   "+str(i)+"    "+str(positions[i][0])+"  "+str(positions[i][1])+"  "+str(positions[i][2])+"\n")

    # 写入tetgen的ele文件(也就是四面体的索引)
    ele_file = to_path+".ele"
    if(os.path.exists(ele_file)):
        print("remove file: "+ele_file)
        os.remove(ele_file)
    with open(ele_file,'w') as f:
        f.write(str(primitivecount)+"  4  0\n")
        for i in range(primitivecount):
            f.write("   "+str(i)+"    "+str(tet_indices[i*4])+"  "+str(tet_indices[i*4+1])+"  "+str(tet_indices[i*4+2])+"  "+str(tet_indices[i*4+3])+"\n")

    # 写入tetgen的face文件(也就是三角面的索引)
    face_file = to_path+".face"
    if(os.path.exists(face_file)):
        print("remove file: "+face_file)
        os.remove(face_file)

    if(gen_face):
        # 由于本身没有三角面，所以如果想生成face，就自己遍历一遍
        facecount = 0
        for i in range(primitivecount):
            facecount += 4
        with open(face_file,'w') as f:
            f.write(str(facecount)+" 0\n")
            face_i = 0
            for i in range(primitivecount):
                f.write("    "+str(face_i)+"    " + str(tet_indices[i*4])+"    "+str(tet_indices[i*4+2])+"    "+str(tet_indices[i*4+1])+"  -1\n")
                face_i += 1
                f.write("    "+str(face_i)+"    " + str(tet_indices[i*4])+"    "+str(tet_indices[i*4+3])+"    "+str(tet_indices[i*4+2])+"  -1\n")
                face_i += 1
                f.write("    "+str(face_i)+"    " + str(tet_indices[i*4])+"    "+str(tet_indices[i*4+1])+"    "+str(tet_indices[i*4+3])+"  -1\n")
                face_i += 1
                f.write("    "+str(face_i)+"    " + str(tet_indices[i*4+1])+"    "+str(tet_indices[i*4+2])+"    "+str(tet_indices[i*4+3])+"  -1\n")
                face_i += 1
    
    print("\n\nwrite tetgen file success! \nnode file: "+node_file+"\nele file: "+ele_file)


if __name__ == '__main__':
    from_path="models/bunny1000_dilate/bunny1000_dilate.geo"
    to_path=from_path[:-4]

    tet_indices,positions, pointcount,vertexcount,primitivecount = read_geo(from_path)

    write_tetgen(tet_indices,positions, pointcount,primitivecount,to_path, gen_face=True)

附录3：在太极ggui中显示

learn-meshtaichi tut03

import taichi as ti
import meshtaichi_patcher as Patcher

ti.init()

# CAUTION: 我们只加了这一个函数, 其他的基本不变。这个就是用来读取face文件的
def directly_import_surf():
    import numpy as np
    import os
    pwd = os.getcwd().replace("\\", "/")
    face_file_name = pwd + "/models/bunny_tet/bunny_tet.face"
    # print("face_file_name: ", face_file_name)
    with open(face_file_name, 'r') as f:
        lines = f.readlines()
        NF = int(lines[0].split()[0])
        face_indices = np.zeros((NF, 3), dtype=np.int32)
        for i in range(NF):
            face_indices[i] = np.array(lines[i + 1].split()[1:-1], dtype=np.int32)
    return face_indices.flatten()
armadillo_surf_indices = directly_import_surf()

# 读入四面体网格
def init_tet_mesh(model_name):
    #基本与上面一样，只是多了一个CV关系，表示通过一个cell可以找到它的四个顶点
    theMesh = Patcher.load_mesh(model_name, relations=["CV"])
    theMesh.verts.place({'x' : ti.math.vec3})
    theMesh.verts.x.from_numpy(theMesh.get_position_as_numpy())
    display_indices = ti.field(ti.u32, shape = len(armadillo_surf_indices))
    display_indices.from_numpy(armadillo_surf_indices) #这里直接读入了face文件
    return theMesh, display_indices


model_name = "models/bunny_tet/bunny_tet.node"
armadillo, armadillo_indices = init_tet_mesh(model_name)
armadillo_indices.to_numpy()

window = ti.ui.Window("taichimesh", (1024, 1024))
canvas = window.get_canvas()
scene = ti.ui.Scene()
camera = ti.ui.Camera()
camera.up(0, 1, 0)
camera.fov(75)
camera.position(4.5,4.5,0.6)
camera.lookat(3.8, 3.8, 0.5)
camera.fov(75)

frame = 0
paused = ti.field(int, shape=())
paused[None] = 1
while window.running:
    # 用下面这段代码，通过提前设置一个paused变量，我们就可以在运行的时候按空格暂停和继续了！
    for e in window.get_events(ti.ui.PRESS):
        if e.key == ti.ui.SPACE:
            paused[None] = not paused[None]
            print("paused:", paused[None])
    if not paused[None]:
        # substep()
        print(f"frame: {frame}")
        frame += 1
    # 我们可以通过下面的代码来查看相机的位置和lookat，这样我们就能知道怎么调整相机的位置了
    # print("camera.curr_position",camera.curr_position)
    # print("camera.curr_lookat",camera.curr_lookat)

    # movement_speed=0.05表示移动速度，hold_key=ti.ui.RMB表示按住右键可以移动视角
    # wasdqe可以移动相机
    camera.track_user_inputs(window, movement_speed=0.05, hold_key=ti.ui.RMB)
    scene.set_camera(camera)
    
    scene.mesh(armadillo.verts.x, armadillo_indices, color = (0.5,0.5,0.5))

    scene.point_light(pos=(0.5, 1.5, 0.5), color=(1, 1, 1))
    scene.ambient_light((0.5,0.5,0.5))

    canvas.scene(scene)

    window.show()

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.coloradmin.cn/o/402399.html

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈，一经查实，立即删除！

相关文章

[Java Web]Request对象 | 超1w字带你熟悉Servlet中的request请求

[Java Web]Request对象 | 超1w字带你熟悉Servlet中的request请求

⭐作者介绍：大二本科网络工程专业在读，持续学习Java，输出优质文章 ⭐所属专栏：Java Web ⭐如果觉得文章写的不错，欢迎点个关注😉有写的不好的地方也欢迎指正，一同进步😁 目录 Reque…
阅读更多...
Codeforces Round 857 (Div. 2)【A-C】

Codeforces Round 857 (Div. 2)【A-C】

文章目录A. Likes【贪心、模拟】B. Settlement of Guinea Pigs【贪心】C. The Very Beautiful Blanket【构造、观察】链接传送门A. Likes【贪心、模拟】 分析 为了使得当前时间点赞的尽可能大，那么前面的赞的数目也要尽可能大，所以前面把能赞的都要先赞…
阅读更多...
2-8 SpringCloud快速开发入门: Eureka 服务注册中心自我保护机制

2-8 SpringCloud快速开发入门: Eureka 服务注册中心自我保护机制

接上一章节Eureka 注册中心高可用集群搭建，这里讲讲Eureka 服务注册中心自我保护机制 Eureka 服务注册中心自我保护机制 自我保护机制是 Eureka 注册中心的重要特性，当 Eureka 注册中心进入自我保护模式时，在 Eureka Server 首页会输出如下警…
阅读更多...
Python JS逆向篇（一）

Python JS逆向篇（一）

Python JS逆向篇（一）效果实现思路最后一步逆向 p.a.HmacSHA256(t, s["a"].state.commonStore.cupid_sign_key)JS实现py实现（先苦后甜）逆向主题：51job请求头headers中携带的sign参数。 （注&#x…
阅读更多...
Windows基于Nginx搭建RTMP流媒体服务器（附带所有组件下载地址及验证方法）

Windows基于Nginx搭建RTMP流媒体服务器（附带所有组件下载地址及验证方法）

RTMP服务时常用于直播时提供拉流推流传输数据的一种服务。前段时间由于朋友想搭建一套直播时提供稳定数据传输的服务器，所以就研究了一下如何搭建及使用。 1、下载nginx 首先我们要知道一般nginx不能直接配置rtmp服务，在Windows系统上需要特殊nginx版本…
阅读更多...
centos8 安装 pcs pacemaker

centos8 安装 pcs pacemaker

一、背景 在centos-8中安装pcs、pacemaker会显示找不到源 （yum install pcs pacemaker 也是一样的） 通过搜索引擎，有说：dnf config-manager --set-enable HighAvailability 也有的说：执行dnf update 也有的说执行 dn…
阅读更多...
AB测试——流程介绍（设计实验）

AB测试——流程介绍（设计实验）

前言： 作为AB测试的学习记录，接上文内容， 本文继续介绍假设建立和实验设计部分，包括实验对象、样本量计算（显著性水平、统计功效及最小可检测效应）、实验周期。 相关文章： AB测试——原理介绍 A…
阅读更多...
【PyTorch】教程：torch.nn.PReLU

【PyTorch】教程：torch.nn.PReLU

torch.nn.PReLU 原型 CLASS torch.nn.PReLU(num_parameters1, init0.25, deviceNone, dtypeNone) 参数 num_parameters ([int]) – 需要学习的 aaa 的数量，尽管作为输入，只有两个值是合法的，1 或者 输入的通道数，默认为 1 。ini…
阅读更多...
各种光照模型的shader实现

各种光照模型的shader实现

大家好，我是阿赵。 这里打算给大家介绍一些常用的光照模型的shader实现方法。虽然这些光照模型很多都会在各大引擎内置，一般不需要自己写。但我觉得学习一下，首先对了解渲染原理有帮助，然后对写一些复合效果的shader时&#xff0c…
阅读更多...
java安全编码规范考试

java安全编码规范考试

java安全编码规范考试 整理不易，收点币！！ 安全编码规范考试.md 下面对zip文件的安全解压缩描述，错误的是 A.zip文件解压时，可以使用entry.getSize(）对解压缩文件进行文件大小判断 B.zip文件解压时&…
阅读更多...
Spring Cloud （ Eureka集群的搭建 ）

Spring Cloud （ Eureka集群的搭建 ）

操作步骤： 添加主机映射创建Eureka服务 导入依赖编写启动类编写yml复制服务实例测试一、添加主机映射 以 Windows系统为例，如果要构建集群，需要修改 hosts 文件，为其添加主机名的映射。 打开C:\Windows\System32\drivers\etc\h…
阅读更多...
DJ1-4 计算机网络和因特网

DJ1-4 计算机网络和因特网

目录 一、协议层及其服务模型 ISO/OSI 七层参考模型 TCP/IP 参考模型 1. 网际协议栈（protocol stack） 2. 分层：逻辑通信 3. 协议分层与数据 二、攻击威胁下的网络 1. 植入恶意软件 2. 攻击服务器和网络基础设施 3. 嗅探分组 4. 伪…
阅读更多...
为什么机器人操作系统ROS1、ROS2这么难学

为什么机器人操作系统ROS1、ROS2这么难学

AI回复：ROS1、ROS2是机器人操作系统，涉及到很多复杂的概念和技术，需要有一定的编程和机器人知识基础才能学习。此外，ROS1和ROS2的架构和设计也有很大的不同，需要花费一定的时间和精力去学习和适应。但是，一…
阅读更多...
DHCP的配置

DHCP的配置

实验目的熟悉DHCP的应用场景掌握DHCP的配置方法实验拓扑DHCP的配置如图15-2所示： 图15-2：DHCP的配置 实验步骤配置IP地址<Huawei>system-view Enter system view, return user view with Ctrl+Z.
阅读更多...
新手入门python实现神经网络，超级简单!

新手入门python实现神经网络，超级简单!

前言&#xff1a; 这篇文章完全是为新手准备的。我们会通过用Python从头实现一个神经网络来理解神经网络的原理。 文章目录神经元1、一个简单的例子2、编码一个神经元把神经元组装成网络1、例子&#xff1a;前馈2、编码神经网络&#xff1a;前馈训练神经网络 第一部分1、损失2、…
阅读更多...
C51---串口发送指令,控制LED灯亮灭

C51---串口发送指令,控制LED灯亮灭

1.Code: #include "reg52.h" #include "intrins.h" sfr AUXR 0x8E; sbit D5 P3^7; void UartInit(void) //9600bps11.0592MHz { //PCON & 0x7F; //波特率不倍速 AUXR 0x01; SCON 0x50; //8位数据,可变波…
阅读更多...
Spark-RDD 转换算子（双 Value 类型、Key - Value 类型）

Spark-RDD 转换算子（双 Value 类型、Key - Value 类型）

双 Value 类型 1、intersection&#xff08;交集&#xff09; 2、union&#xff08;并集&#xff09; 3、subtract&#xff08;差集&#xff09; 4、zip&#xff08;拉链&#xff09; Key - Value 类型 1、partitionBy 2、reduceByKey 3、groupByKey 4、aggregateByK…
阅读更多...
FinOps首次超越安全成为企业头等大事｜云计算趋势报告

FinOps首次超越安全成为企业头等大事｜云计算趋势报告

随着云计算在过去十年中的广泛应用&#xff0c;云计算用户所面临的一个持续不变的趋势是&#xff1a;安全一直是用户面临的首要挑战。然而&#xff0c;这种情况正在发生转变。 知名IT软件企业 Flexera 对云计算决策者进行年度调研已经持续12年&#xff0c;而今年安全问题首次…
阅读更多...
3.初识Vue

3.初识Vue

目录 1 vue 浏览器调试工具 1.1 安装 1.2 配置 2 数据驱动视图与双向数据绑定 3 简单使用 3.1 下载 3.2 将信息渲染到DOM上 4 使用vue浏览器调试工具 5 vue指令 1 vue 浏览器调试工具 chrome可能是我浏览器的原因&#xff0c;装上用不了&#xff0c;我们使…
阅读更多...
javaWeb核心05-FilterListenerAjax

javaWeb核心05-FilterListenerAjax

文章目录Filter&Listener&Ajax1&#xff0c;Filter1.1 Filter概述1.2 Filter快速入门1.2.1 开发步骤1.2.2 代码演示1.3 Filter执行流程1.4 Filter拦截路径配置1.5 过滤器链1.5.1 概述1.5.2 代码演示1.5.3 问题1.6 案例1.6.1 需求1.6.2 分析1.6.3 代码实现1.6.3.1 创建F…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023