開放式 VDB 是奧斯卡獎獲獎的稀疏動態(tài)卷的行業(yè)標準庫。在整個視覺效果行業(yè)中，它被用于模擬和渲染水、火、煙、云和大量其他依賴于稀疏體積數(shù)據(jù)的效果。該庫包括一個分層的、動態(tài)的數(shù)據(jù)結(jié)構(gòu)和一套工具，用于高效地存儲和操作三維網(wǎng)格上離散的稀疏體數(shù)據(jù)。庫由 學(xué)院軟件基金會（ ASWF ） 維護。有關(guān)詳細信息，請參見 VDB ：具有動態(tài)拓撲的高分辨率稀疏卷 。

盡管 OpenVDB 提供了性能優(yōu)勢，但它的設(shè)計并沒有考慮到 GPUs 。它對幾個外部庫的依賴使得利用 GPUs 上的 VDB 數(shù)據(jù)變得很麻煩，這正是本文主題的動機。我們將向您介紹 NanoVDB 庫，并提供一些如何在光線跟蹤和碰撞檢測上下文中使用它的示例。

NanoVDB 簡介

最初在 NVIDIA 開發(fā)的 NanoVDB 庫是一個 ASWF OpenVDB 項目的新增功能 。它提供了一個與 OpenVDB 的核心數(shù)據(jù)結(jié)構(gòu)完全兼容的簡化表示，具有在 NanoVDB 和 OpenVDB 數(shù)據(jù)結(jié)構(gòu)之間來回轉(zhuǎn)換、創(chuàng)建和可視化數(shù)據(jù)的功能。

圖 1 OpenVDB 和 NanoVDB 數(shù)據(jù)結(jié)構(gòu)的圖示。

NanoVDB 采用了 VDB 樹結(jié)構(gòu)的壓縮、線性化、只讀表示（圖 1 ），這使得它適合于樹層次結(jié)構(gòu)的快速傳輸和快速、無指針遍歷。為了提高效率，數(shù)據(jù)流經(jīng)過調(diào)整，可以在 GPUs 和 CPU 上使用。

創(chuàng)建 NanoVDB 網(wǎng)格

盡管 NanoVDB 網(wǎng)格是一種只讀數(shù)據(jù)結(jié)構(gòu)，但該庫包含生成或加載數(shù)據(jù)的功能。

所有的 OpenVDB 網(wǎng)格類 – LevelSets 、 FogVolumes 、 PointIndexGrids 和 PointDataGrids ——都支持 NanoVDB 表示，并且可以直接從 OpenVDB 文件（即 。 vdb 系統(tǒng) 文件）加載。還可以將數(shù)據(jù)加載或保存到 NanoVDB 自己的文件格式中或從中保存，該格式本質(zhì)上是其內(nèi)存流的一個轉(zhuǎn)儲，其中包含用于高效檢查的附加元數(shù)據(jù)。

以下代碼示例從 OpenVDB 文件轉(zhuǎn)換：

以下代碼示例從 OpenVDB 文件轉(zhuǎn)換：

openvdb::io::File file(fileName);
auto vdbGrid = file.readGrid(gridName);
auto handle = nanovdb::openToNanoVDB(vdbGrid);

雖然從現(xiàn)有的 OpenVDB 數(shù)據(jù)加載是典型的用例，但是附帶的網(wǎng)格生成器工具允許您直接在內(nèi)存中構(gòu)建 NanoVDB 網(wǎng)格。提供了一些簡單原語的函數(shù)來幫助您入門：

// generate a sparse narrow-band level set (i.e. truncated signed distance field) representation of a sphere.
auto handle = nanovdb::createLevelSetSphere(50, nanovdb::Vec3f(0));

下面的示例顯示了如何使用 lambda 函數(shù)生成小而密集的體積（圖 2 ）：

nanovdb::GridBuilder builder(0);
auto op = [](const nanovdb::Coord& ijk) -> float { return menger(nanovdb::Vec3f(ijk) * 0.01f);
};
builder(op, nanovdb::CoordBBox(nanovdb::Coord(-100), nanovdb::Coord(100)));
// create a FogVolume grid called "menger" with voxel-size 1
auto handle = builder.getHandle<>(1.0, nanovdb::Vec3d(0), "menger", nanovdb::GridClass::FogVolume);

網(wǎng)格控制柄

網(wǎng)格句柄 是一個簡單的類，它擁有它分配的緩沖區(qū)的所有權(quán)，允許網(wǎng)格的范圍劃分（ RAII ）。

它還用于封裝不透明的網(wǎng)格數(shù)據(jù)。盡管網(wǎng)格數(shù)據(jù)本身是以數(shù)據(jù)類型（如 浮動 為模板的），但句柄提供了一種方便的方法來訪問網(wǎng)格的元數(shù)據(jù)，而不必知道網(wǎng)格的數(shù)據(jù)類型 MIG 是什么。這很有用，因為您可以純粹從句柄確定 GridType 。

下面的代碼示例驗證是否有包含級別集函數(shù)的 32 位浮點網(wǎng)格：

const nanovdb::GridMetaData* metadata = handle.gridMetaData();
if (!metadata->isLevelSet() || !metadata->gridType() == GridType::Float) throw std::runtime_error("Not the right stuff!");

網(wǎng)格緩沖區(qū)

NanoVDB 被設(shè)計成支持許多不同的平臺， CPU ， CUDA 甚至圖形 api 。為了實現(xiàn)這一點，數(shù)據(jù)結(jié)構(gòu)被存儲在一個平坦的連續(xù)內(nèi)存緩沖區(qū)中。

使這個緩沖區(qū)駐留在 CUDA 設(shè)備上很簡單。為了完全控制，您可以使用 CUDA api 分配設(shè)備內(nèi)存，然后將句柄的數(shù)據(jù)上載到其中。

void* d_gridData;
cudaMalloc(&d_gridData, handle.size());
cudaMemcpy(d_gridData, handle.data(), handle.size(), cudaMemcpyHostToDevice);
const nanovdb::FloatGrid* d_grid = reinterpret_cast(d_gridData);

NanoVDB 的 GridHandle 模板化在緩沖區(qū)類型上，緩沖區(qū)類型是其內(nèi)存分配的包裝器。它默認為使用主機系統(tǒng)內(nèi)存的主機緩沖區(qū)；然而， NanoVDB 還提供了CUDA 緩沖器，以便輕松創(chuàng)建 CUDA 設(shè)備緩沖區(qū)。

auto handle = nanovdb::openToNanoVDB<nanovdb::CudaDeviceBuffer>(vdbGrid);
handle->deviceUpload();
const nanovdb::FloatGrid* grid = handle->deviceGrid();

將數(shù)據(jù)流解釋為納米網(wǎng)格類型后，可以使用這些方法訪問網(wǎng)格中的數(shù)據(jù)。有關(guān)更多詳細信息，請參閱相關(guān) API 的文檔。本質(zhì)上，它反映了 OpenVDB 中只讀方法的基本子集。

auto hostOrDeviceOp = [grid] __host__ __device__ (nanovdb::Coord ijk) -> float {
  // Note that ReadAccessor (see below) should be used for performance.
  return grid->tree().getValue(ijk);
};

可以構(gòu)造自定義緩沖區(qū)來處理不同的內(nèi)存空間。有關(guān)創(chuàng)建可與圖形 API 交互操作的緩沖區(qū)的示例的更多信息，請參閱存儲庫中的示例。

致使

由于 NanoVDB 網(wǎng)格提供了一個緊湊的只讀 VDB 樹，因此它們很適合渲染任務(wù)。光線將 VDB 網(wǎng)格跟蹤到圖像中。使用每線程一條光線，并使用一個自定義的 雷吉諾 functor 生成光線，該函數(shù)接受像素偏移并創(chuàng)建世界空間光線。完整的代碼在存儲庫示例中可用。

考慮到沿射線采樣具有空間相干性這一事實，可以使用 讀寫器 來加速采樣。當光線向前移動時，這會緩存樹遍歷堆棧，從而允許自底向上的樹遍歷，這比傳統(tǒng)的自上而下遍歷要快得多，后者涉及相對較慢的無界根節(jié)點。

auto renderTransmittanceOp = [image, grid, w, h, rayGenOp, imageOp, dt] __host__ __device__ (int i) {
    nanovdb::Ray wRay = rayGenOp(i, w, h);
    // transform the ray to the grid's index-space...
    nanovdb::Ray iRay = wRay.worldToIndexF(*grid);
    // clip to bounds.
    if (iRay.clip(grid->tree().bbox()) == false) {
        imageOp(image, i, w, h, 1.0f);
        return;
    }
    // get an accessor.
    auto acc = grid->tree().getAccessor();
    // integrate along ray interval...
    float transmittance = 1.0f;
    for (float t = iRay.t0(); t < iRay.t1(); t+=dt) {
        float sigma = acc.getValue(nanovdb::Coord::Floor(iRay(t)));
        transmittance *= 1.0f - sigma * dt;
    }
    imageOp(image, i, w, h, transmittance );
};

由于光線與水平集網(wǎng)格相交是一項常見任務(wù)， NanoVDB 實現(xiàn)了一個零交叉功能，并使用分層 DDA （ HDDA ）作為沿光線的根搜索的一部分來清空空間跳躍（圖 5 ）。有關(guān) HDDA 的更多信息，請參閱 OpenVDB 中高效光線行進的分層數(shù)字微分分析儀 。下面是代碼示例：

...
    auto acc = grid->tree().getAccessor();
    // intersect with zero level-set...
    float iT0;
    nanovdb::Coord ijk;
    float v;
    if (nanovdb::ZeroCrossing(iRay, acc, ijk, v, iT0)) { 
        // convert intersection distance (iT0) to world-space
        float wT0 = iT0 * grid->voxelSize();
        imageOp(image, i, w, h, wT0);
    } else {
        imageOp(image, i, w, h, 0.0f);
    }
...

碰撞檢測

碰撞檢測和解決是 NanoVDB 的另一項任務(wù)，因為它們通常需要有效地查找實體碰撞對象的有符號距離值。窄帶電平集表示非常理想，因為它們用符號對內(nèi)部/外部拓撲信息（碰撞檢測所需）進行了緊湊編碼。此外，最近點變換（沖突解決所需的）很容易從水平集函數(shù)的梯度計算。

下面的代碼示例是一個用于處理沖突的函數(shù)。使用 讀寫器 是很有用的，因為用于沖突解決的梯度計算涉及到同一空間附近的多個提取。

auto collisionOp = [grid, positions, velocities, dt] __host__ __device__ (int i) {
    nanovdb::Vec3f wPos = positions[i];
    nanovdb::Vec3f wVel = velocities[i];
    nanovdb::Vec3f wNextPos = wPos + wVel * dt;
    // transform the position to a custom space...
    nanovdb::Vec3f iNextPos = grid.worldToIndexF(wNextPos);
    // the grid index coordinate.
    nanovdb::Coord ijk = nanovdb::Coord::Floor(iNextPos);
    // get an accessor.
    auto acc = grid->tree().getAccessor();
    if (tree.isActive(ijk)) { // are you inside the narrow band?
        float wDistance = acc.getValue(ijk);
        if (wDistance <= 0) { // are you inside the levelset?
            // get the normal for collision resolution.
            nanovdb::Vec3f normal(wDistance);
            ijk[0] += 1;
            normal[0] += acc.getValue(ijk);
            ijk[0] -= 1;
            ijk[1] += 1;
            normal[1] += acc.getValue(ijk);
            ijk[1] -= 1;
            ijk[2] += 1;
            normal[2] += acc.getValue(ijk);
            normal.normalize();
            
            // handle collision response with the surface.
            collisionResponse(wPos, wNextPos, normal, wDistance, wNextPos, wNextVel);
        }
    }
    positions[i] = wNextPos;
    velocities[i] = wNextVel;
};

同樣，完整的代碼可以在存儲庫中找到。

基準

NanoVDB 被開發(fā)成在主機和設(shè)備上同樣運行良好。使用 modernCUDA 中的擴展 lambda 支持，您可以輕松地在兩個平臺上運行相同的代碼。

本節(jié)包括比較英特爾線程構(gòu)建塊和 CPU CUDA 上光線跟蹤和碰撞檢測性能的基準測試。計時以毫秒為單位，與 NVIDIA NVIDIA 8000 相比，是在 Xeon E5-2696 v4 x2 –（ 88 個 CPU 線程）上生成的。使用的 FogVolume 是兔子云， LevelSet 是 dragon 數(shù)據(jù)集。兩者都可以從 OpenVDB 網(wǎng)站下載。渲染的分辨率為 1024×1024 。這次碰撞試驗?zāi)M了一億顆彈道粒子。

雖然基準測試（圖 6 和下表）顯示了 NanoVDB 高效表示加速 CPU 上 OpenVDB 的好處，但它真正突出了使用 GPU 對 VDB 數(shù)據(jù)進行只讀訪問以進行碰撞檢測和光線跟蹤的好處。

結(jié)論

NanoVDB 是一個小而強大的庫，它通過使用 GPUs 來加速某些 OpenVDB 應(yīng)用程序。開源存儲庫現(xiàn)在可用了！要下載源代碼、構(gòu)建示例

關(guān)于作者

Wil Braithwaite 在倫敦和洛杉磯的工作室工作了 15 年的視覺特效。他的職位包括研究、技術(shù)指導(dǎo)、合成、 CG 監(jiān)督和 MOCAP 監(jiān)督。他開創(chuàng)了圖形硬件在 VFX 管道中的應(yīng)用，在 NVIDIA 擔任高級應(yīng)用工程師，專門從事咨詢、培訓(xùn)和使用 NVIDIA 技術(shù)協(xié)助 VFX 工作室項目的開發(fā)。

Ken Museth 是模擬技術(shù)的高級主管，并于 2020 年初加入 NVIDIA ，當時他發(fā)起了 NanoVDB 的開發(fā)。他以前在開發(fā)虛擬現(xiàn)實技術(shù)的時候，一直致力于虛擬現(xiàn)實的開發(fā)。他是 VDB 的創(chuàng)建者和 OpenVDB 的首席架構(gòu)師，也是其技術(shù)指導(dǎo)委員會的主席。此外，肯在 SpaceX 公司工作了六年，對新的猛禽火箭發(fā)動機進行大規(guī)模流體動力學(xué)模擬。在 2017 年加入 Weta 之前，他在夢工廠動畫和數(shù)字領(lǐng)域工作了 10 年，在此之前，他曾在加州理工學(xué)院和林科平大學(xué)擔任研究員和全職教授。他擁有哥本哈根大學(xué)量子動力學(xué)博士學(xué)位，并獲得電影藝術(shù)與科學(xué)學(xué)院頒發(fā)的技術(shù)成就獎?？鲜?SIGGRAPH 2020 技術(shù)論文委員會成員。

審核編輯：郭婷