关于torch-points3d可视化的问题--安装环境

#!/usr/bin/env bash

PYTHON=${PYTHON:-"python"}
#cuda 10.1

$PYTHON -mpip install ./torch-1.6.0+cu101-cp37-cp37m-linux_x86_64.whl
$PYTHON -mpip install ./torchvision-0.7.0+cu101-cp37-cp37m-linux_x86_64.whl
$PYTHON -mpip install ./torch_cluster-1.5.8+cu101-cp37-cp37m-linux_x86_64.whl
$PYTHON -mpip install ./torch_scatter-2.0.5+cu101-cp37-cp37m-linux_x86_64.whl
$PYTHON -mpip install ./torch_sparse-0.6.8+cu101-cp37-cp37m-linux_x86_64.whl
$PYTHON -mpip install ./torch_spline_conv-1.2.0+cu101-cp37-cp37m-linux_x86_64.whl

$PYTHON -mpip install torch_geometric==1.7.0 -i https://pypi.tuna.tsinghua.edu.cn/simple

$PYTHON -mpip install torch_points3d==1.2.0 -i https://pypi.tuna.tsinghua.edu.cn/simple

参考:

http://ziyedy.top/page/plyfile.html

https://www.pythonf.cn/read/69331

https://blog.csdn.net/weixin_44389243/article/details/106746240

https://blog.csdn.net/guaiderzhu1314/article/details/105749413/

输出的ply无法在cloudcompare上打开

我切换成pcd进行打开

# python3.8
# -*- coding: utf-8 -*-
# ---
# @Software: PyCharm
# @File: vis_ply.py
# @Author: ---
# @Institution: BeiJing, China
# @E-mail: lgdyangninghua@163.com
# @Site: 
# @Time: 4月 26, 2021
# ---

import os
import numpy as np
import open3d as o3d
from plyfile import PlyData, PlyElement

HEADER = '''\
# .PCD v0.7 - Point Cloud Data file format
VERSION 0.7
FIELDS x y z label object
SIZE 4 4 4 4 4
TYPE F F F I I
COUNT 1 1 1 1 1
WIDTH {}
HEIGHT 1
VIEWPOINT 0 0 0 1 0 0 0
POINTS {}
DATA ascii
'''

def write_pcd(points, save_pcd_path):
    n = len(points)
    lines = []
    for i in range(n):
        x, y, z, label, object = points[i]
        lines.append('{:.6f} {:.6f} {:.6f} {} {}'.format(x, y, z, label, object))
    with open(save_pcd_path, 'w') as f:
        f.write(HEADER.format(n, n))
        f.write('\n'.join(lines))

def write_pcd(points, save_pcd_path):
    with open(save_pcd_path, 'w') as f:
        f.write(HEADER.format(len(points), len(points)) + '\n')
        np.savetxt(f, points, delimiter=' ', fmt='%f %f %f %d %d')


def read_ply(filename):
    """ 读取坐标,返回值为 n * 3 """
    plydata = PlyData.read(filename)
    x = plydata.elements[0].data['x']
    y = plydata.elements[0].data['y']
    z = plydata.elements[0].data['z']
    p = plydata.elements[0].data['p'] #pred
    l = plydata.elements[0].data['l'] #label
    point = np.array([[_x, _y, _z, _p, _l] for _x, _y, _z, _p, _l in zip(x, y, z, p, l)])
    return point

path = "torch-points3d-master/outputs/benchmark/benchmark-pointnet2_charlesssg-20210426_160121/viz/99/train/"
lines = os.listdir(path)
for m_ind, m_val in enumerate(lines):
    if m_val.endswith('ply'):
        point = read_ply(os.path.join(path, m_val))
        xyzpred = point[:, :4]
        h,w = xyzpred.shape
        object = np.ones((h,1), dtype=np.int8)*(-1)
        save_pcd_point = np.hstack((xyzpred, object))
        write_pcd(save_pcd_point, os.path.join(path, m_val).replace('ply', 'pcd'))

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
RandLA-Net是一种基于点云数据的深度学习模型,用于点云分割和场景理解。下面是使用PyTorch实现RandLA-Net的简单步骤: 1. 安装依赖库 在Python环境安装以下库: - PyTorch - NumPy - Open3D - Scikit-learn 其中PyTorch是必须的,其余库是为了可视化和数据预处理。 2. 下载数据集 下载点云数据集,例如S3DIS数据集,该数据集包含了用于建筑物场景的点云数据。可以从官方网站下载数据集。 3. 数据预处理 使用Open3D库读取点云数据并进行预处理。具体来说,可以使用Open3D库将点云数据转换为numpy数组,然后将其分为小的块,以便在GPU上进行训练。 ```python import open3d as o3d import numpy as np import os def load_data(path): pcd = o3d.io.read_point_cloud(path) points = np.asarray(pcd.points) return points def process_data(points, block_size=3.0, stride=1.5): blocks = [] for x in range(0, points.shape[0], stride): for y in range(0, points.shape[1], stride): for z in range(0, points.shape[2], stride): block = points[x:x+block_size, y:y+block_size, z:z+block_size] if block.shape[0] == block_size and block.shape[1] == block_size and block.shape[2] == block_size: blocks.append(block) return np.asarray(blocks) # Example usage points = load_data("data/room1.pcd") blocks = process_data(points) ``` 这将生成大小为3x3x3的块,每个块之间的距离为1.5。 4. 构建模型 RandLA-Net是一个基于点云的分割模型,它使用了局部注意力机制和多层感知器(MLP)。这里给出一个简单的RandLA-Net模型的实现: ```python import torch import torch.nn as nn class RandLANet(nn.Module): def __init__(self, input_channels, num_classes): super(RandLANet, self).__init__() # TODO: Define the model architecture self.conv1 = nn.Conv1d(input_channels, 32, 1) self.conv2 = nn.Conv1d(32, 64, 1) self.conv3 = nn.Conv1d(64, 128, 1) self.conv4 = nn.Conv1d(128, 256, 1) self.conv5 = nn.Conv1d(256, 512, 1) self.mlp1 = nn.Sequential( nn.Linear(512, 256), nn.BatchNorm1d(256), nn.ReLU(), nn.Linear(256, 128), nn.BatchNorm1d(128), nn.ReLU(), nn.Linear(128, num_classes), nn.BatchNorm1d(num_classes) ) def forward(self, x): # TODO: Implement the forward pass x = self.conv1(x) x = self.conv2(x) x = self.conv3(x) x = self.conv4(x) x = self.conv5(x) x = torch.max(x, dim=-1)[0] x = self.mlp1(x) return x ``` 这个模型定义了5个卷积层和一个多层感知器(MLP)。在前向传递过程中,点云数据被送入卷积层,然后通过局部最大池化层进行处理。最后,通过MLP将数据转换为预测的类别。 5. 训练模型 在准备好数据和模型之后,可以使用PyTorch的内置函数训练模型。这里使用交叉熵损失函数和Adam优化器: ```python import torch.optim as optim device = torch.device("cuda" if torch.cuda.is_available() else "cpu") # TODO: Initialize the model model = RandLANet(input_channels=3, num_classes=13).to(device) # TODO: Initialize the optimizer and the loss function optimizer = optim.Adam(model.parameters(), lr=0.001) loss_fn = nn.CrossEntropyLoss() # TODO: Train the model for epoch in range(num_epochs): running_loss = 0.0 for i, batch in enumerate(train_loader): # Move the batch to the GPU batch = batch.to(device) # Zero the gradients optimizer.zero_grad() # Forward pass outputs = model(batch) loss = loss_fn(outputs, batch.labels) # Backward pass and optimization loss.backward() optimizer.step() # Record the loss running_loss += loss.item() # Print the epoch and the loss print('Epoch [%d], Loss: %.4f' % (epoch+1, running_loss / len(train_loader))) ``` 这里使用Adam优化器和交叉熵损失函数进行训练。训练完成后,可以使用预测函数对新数据进行分类: ```python def predict(model, data): with torch.no_grad(): # Move the data to the GPU data = data.to(device) # Make predictions outputs = model(data) _, predicted = torch.max(outputs.data, 1) # Move the predictions back to CPU predicted = predicted.cpu().numpy() return predicted # Example usage data = load_data("data/room2.pcd") data = process_data(data) data = torch.from_numpy(data).float().permute(0, 2, 1) predicted = predict(model, data) ``` 这将返回点云数据的分类预测。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值