## UlRe-NeRF: 3D Ultrasound Imaging through Neural Rendering with Ultrasound Reflection Direction Parameterization

Ziwen Guo, Zi Fang, and Zhuang Fu*

Shanghai Jiao Tong University miixnon1@sjtu.edu.cn

Abstract. Three-dimensional ultrasound imaging is a critical technology widely used in medical diagnostics. However, traditional 3D ultrasound imaging methods have limitations such as fixed resolution, low storage efficiency, and insufficient contextual connectivity, leading to poor performance in handling complex artifacts and reflection characteristics. Recently, techniques based on NeRF (Neural Radiance Fields) have made significant progress in view synthesis and 3D reconstruction, but there remains a research gap in high-quality ultrasound imaging. To address these issues, we propose a new model, UlRe-NeRF, which combines implicit neural networks and explicit ultrasound volume rendering into an ultrasound neural rendering architecture. This model incorporates reflection direction parameterization and harmonic encoding, using a directional MLP module to generate view-dependent high-frequency reflection intensity estimates, and a spatial MLP module to produce the medium's physical property parameters. These parameters are used in the volume rendering process to accurately reproduce the propagation and reflection behavior of ultrasound waves in the medium. Experimental results demonstrate that the UlRe-NeRF model significantly enhances the realism and accuracy of high-fidelity ultrasound image reconstruction, especially in handling complex medium structures.

Keywords: Ultrasound imaging · Implicit Neural Networks · Ultrasound Volume Rendering

## 1 Introduction

Ultrasound (US) is a highly valued imaging modality in medicine due to its costeffectiveness, safety, non-invasiveness, and ability to provide real-time imaging. Unlike MRI and CT, US is portable, radiation-free, and user-friendly, making it ideal for various clinical applications, including intraoperative imaging and surgical guidance. Its affordability and efficiency have led to its growing popularity among clinicians and researchers. Compared to conventional two-dimensional (2D) ultrasound, three-dimensional (3D) ultrasound offers significant advantages. 3D ultrasound allows direct visualization of 3D anatomy, enhancing diagnostic accuracy and reducing dependence on operator skill, while providing

more precise measurements of organ volumes and other quantitative properties, enabling repeatable and consistent imaging for follow-up studies [6]. It also offers the flexibility to view anatomical structures from any orientation, overcoming the limitations imposed by patient positioning in 2D ultrasound.

Three-dimensional (3D) ultrasound imaging involves several key steps to realize view synthesis of novel 2D ultrasound images [7] [12]. The first step is the data acquisition phase, where a series of 2D ultrasound images I i ∈ J along with corresponding image plane pose T I i ∈ T are captured using a 1D ultrasound probe with tracking sensors attached. The next step is the volumetric representation phase, where a 3D representation ˆ = V F θ ( I, T ) is obtained based on the 2D images and their corresponding poses, aiming to preserve as much of the information as possible. The final step is the rendering phase, for any novel image plane pose T l low ∈ T , a synthesized 2D ultrasound image ˆ I new = G µ ( ˆ V , T new ) can be rendered. A high-quality 3D scene reconstruction should achieve synthesized views that are as close to the actual views as possible.

Much past decades, some researches about explicit volumetric representation have recently been developed [12]. Each slice is treated as a slice point cloud P i ∈ P , where each point has corresponding spatial coordinates x i and color c i . In the volumetric representation phase, an explicit 3D voxel grid ˆ V is constructed. For each voxel x voxel ∈ ˆ V , the corresponding color ˆ( c x voxel ) is obtained by using a parameterized interpolation algorithm F θ ( {·} ) to gather the spatial and color information of the point cloud points contained in its spatial neighborhood N δ ( x voxel ) ∩ P .

However, those explicit volumetric representation methods have some flaws, leading to that synthetic images do not conform to physical laws:

- 1. The unreasonable Lambertian assumption, that the same spatial point appears the same from different viewpoints, causes difficulty in reproducing physical phenomena about viewpoints like refraction and mirror image artifacts.
- 2. Spatial neighborhood interpolation methods have a very small receptive field, which weaken the contextual connections between spatially distant points, making it hard to reproduce acoustic shadowing and enhancement caused by different tissue impedances.
- 3. Explicit discrete voxel grids are inefficient in terms of storage and have a fixed resolution, which depends on the grid and is often very low. Querying the values between grid points requires linear interpolation. Furthermore, the discretization makes differentiation difficult.

In the literature, as summarized in Section 2.1, researchers attempt to address this issue through implicit representation methods. In explicit representation, the color value of a query point is obtained by searching the spatial neighborhood in the dataset and interpolating F θ : { ( N δ ( x ) ∩ P } → ) ˆ c . In implicit representation, uses a continuous coordinate-based MLP (Multi-Layer Perceptron) F θ : ( x, d ) → ˆ f to query multidimensional representations ˆ f from threedimensional coordinates x and two-dimensional directions d .

Fig. 1. UlRe-NeRF for High-Fidelity 3D Ultrasound Imaging. This figure illustrates our proposed ultrasound neural rendering architecture, which integrates implicit neural representation (b) and volume rendering (c) to enhance 3D ultrasound imaging. The system comprises two main modules: the spatial MLP and the directional MLP. The spatial MLP predicts physical property parameters, such as attenuation coefficients, reflectivity, scattering density, and amplitude, through high-dimensional position encoding. This simulates the propagation and reflection behavior of ultrasound waves within the medium. The directional MLP processes the input ray directions, generating view-dependent intensity estimates that capture complex reflection phenomena.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00860v3-figure-f53d8a3d-f1f5-4c74-a4c9-ca332aba3d4b.jpg)

Neural Radiance Fields (NeRF) have significantly advanced the general application of implicit representation [24]. NeRF uses implicit representation, specifically a coordinate-based MLP supervised to represent shape or appearance, combined with explicit rendering through a fixed differentiable rendering engine, to achieve novel view synthesis. The strong physically motivated inductive biases act as regularizers, ensuring that the derived implicit representation closely mimics the real world, which leads to better generalization.

Based on the foundational concepts of NeRF and incorporating additional inductive biases, numerous NeRF variants have been proposed to expand their application scenarios, such as ultrasound, CT, and MRI, or to enhance performance, such as improving imaging quality and enabling sparse data training.

Ultra-NeRF is the first work to bring NeRF into the field of ultrasound [29]. It builds on the ultrasound volume rendering method by Salehi et al. and combines a coordinate-based MLP that maps from 3D coordinates to five physical parameters, achieving high-quality synthetic ultrasound images [18]. However, it is still view-independent, so cannot achieve ray interactions and the Fresnel Effect, and it does not allow rendering complex ultrasound artifacts, such as reverberations. To address these issues, we drew inspiration from Ref-NeRF and proposed a new model that considers both reflection intensity and volume rendering (Fig. 1(b) and 1(c)). Specifically, we present a ultrasound neural rendering architecture based on implicit neural networks and volume rendering, optimizing the rendering process to better simulate the physical properties, artifacts, and deformations of the medium. We call this model UlRe-NeRF, which can accurately reproduce the propagation and reflection behavior of ultrasound waves

within the medium, as depicted in Fig. 1. In summary, our main contributions are as follows:

High-fidelity 3D Ultrasound Imaging We introduce a novel model that synthesizes precise ultrasound images by learning view-dependent appearance and geometry from multiple ultrasound scans. This method significantly enhances image quality, particularly in capturing fine details and complex tissue structures. It not only improves diagnostic accuracy but also reduces reliance on operator expertise.

Enhanced Implicit Representation Network We designed an implicit representation network that more accurately models the propagation characteristics of ultrasound waves. By incorporating advanced techniques such as reflection direction parameterization and Reflective Harmonic Encoding (RHE), the network better captures high-frequency information and complex tissue reflections. These improvements enable our method to more precisely reconstruct ultrasound images when dealing with intricate tissue structures.

Ultrasound Volume Rendering Building on the enhanced implicit representation network, we developed a complex rendering model that integrates reflection intensity and volume rendering. This model accurately simulates the medium's physical properties, artifacts, and deformations, and correctly handles light interactions and Fresnel effects within the medium. This advancement is crucial for achieving realistic ultrasound image synthesis, effectively overcoming the limitations of existing methods in simulating complex physical phenomena.

In conclusion, UlRe-NeRF has made significant strides in the field of 3D ultrasound imaging, excelling in both image quality and physical realism. It also shows broad potential for applications in medical imaging. Our research represents an important advancement in ultrasound image reconstruction technology and lays a solid foundation for future research and applications.

## 2 Related Works

## 2.1 Variants of NeRF for quality

Neural Radiance Fields (NeRF) uses implicit representation and explicit rendering to achieve novel view synthesis [14]. Given any input 3D position x and view direction d , NeRF uses a spatial MLP to output the density τ ( x ) and bottleneck vector b x ( ) from 3D position x , and then use a second directional MLP to outputs the color of light c x, d ( ) emitted by a particle at 3D position x and direction d (see Fig. 3 for a visualization). The densities { τ i ( x ) } and colors { c i ( x, d ) } queried along the camera ray x i = o + t d, d i i = d originating at o with direction

d , are alpha composited [11] to obtain the color of the pixel corresponding to the ray.

$$\hat { C } ( o, d ) & = \sum _ { i } w _ { i } c _ { i }, w _ { i } \\ & = e ^ { - \sum _ { j < i } \tau _ { i } ( t _ { j + 1 } - t _ { j } ) } \left ( 1 - e ^ { - \tau _ { i } ( t _ { j + 1 } - t _ { j } ) } \right )$$

Numerous subsequent NeRF variants focus on improving the core component of NeRF to improve imaging quality. MLP tends to learn low-frequency functions, limiting its ability to represent complex geometry and texture [17]. To address this issue, NeRF maps the Cartesian coordinates to a higher dimensional space using sinusoidal positional encoding [14] [23]. Mip-NeRF [2] improved volume rendering by using cone tracing instead of the ray tracing in standard NeRF. They introduced Integrated Positional Encoding (IPE), where a cone cast from the camera's center through the pixel was approximated by a multivariate Gaussian. By IPE over a conical frustum for each quadrature segment along the ray, Mip-NeRF encoded the scene at multiple scales, preventing aliasing when rendering from different positions or resolutions. However, due to differences in imaging principles, these neural volumetric scene representation methods designed for natural images are not well-suited for reconstructing the complex biological tissues in ultrasound images. Consequently, their view synthesis performance falls short of meeting the requirements for 3D ultrasound reconstruction.

Ref-NeRF [25] enhances NeRF by structuring view-dependent appearance into pre-filtered incident radiance, diffuse color, material roughness, and specular tint. This decomposition method allows for more accurate rendering of specular highlights and reflections. To achieve this, view direction reflections are reparameterized as a function of reflection direction, significantly improving the interpolation of reflective appearances. The specific formula is as follows:

$$r = 2 ( \mathbf v \cdot \mathbf n ) \mathbf n - \mathbf v$$

Additionally, Ref-NeRF introduces Integrated Directional Encoding (IDE) to effectively represent the material radiance function as it varies with roughness. IDE encodes the reflection direction using a convolution of spherical harmonics and von Mises-Fisher distribution:

$$I D E ( { \mathbf r } ) = \sum _ { l = 0 } ^ { L } \sum _ { m = - l } ^ { l } A _ { l m } \cdot Y _ { l m } ( { \mathbf r } )$$

This method captures finer high-frequency information in reflection directions and shares lighting information across different roughness regions, significantly enhancing rendering realism and accuracy. By modeling reflection directions based on computer graphics principles, it provides a viable solution for accurately capturing ultrasound wave reflections and scattering at tissue interfaces, thereby improving image detail and clarity.

In summary, although NeRF may not be directly suitable for 3D ultrasound reconstruction due to differences in data characteristics and imaging principles,

we can leverage the capability of implicit neural representations to reconstruct high-quality images from sparse viewpoints. By combining volume rendering, we can develop new models tailored for ultrasound imaging, addressing the reflections and refractions of ultrasound waves at various tissue interfaces, thus generating more realistic ultrasound images.

## 2.2 Implicit neural representations in medical imaging

Implicit neural models have been widely used in medical imaging, particularly in solving inverse imaging problems [8]. These problems involve learning the structure of an object (such as an organ of interest in medical cases) from observations or measurements. By using INRs, it becomes possible to synthesize novel high-resolution scans from sparsely-sampled low-resolution data with minimal memory requirements.

Some research has been conducted on reconstructing various medical images using innovative approaches.For MRI modality, IREM [28] and ArSSR [27] leverage implicit continuous representations to recover high-resolution MRI images from low-resolution ones, achieving significant improvements in anatomical detail and image contrast.For CT modality, NeRP [21] reconstructs CT images from sparsely sampled measurements and incorporates longitudinal prior knowledge into deep learning frameworks as network parameters.For X-ray modality, Med-NeRF [4] explores the inference of anatomical 3D structures from a few or a single-view X-ray by combining GRAF [19] with discriminator architecture adapted for CT-projections. In the field of ultrasound, ImplicitVol [30] and work by Li et al. [10] are the early explorations utilizing implicit representations and neutral rendering for high-quality 3D ultrasound reconstruction. But they pay little attention to the unique features of ultrasound, lacking some strong physically motivated inductive biases as regularizers. This results in poor image quality and defects that violate physical laws.

To address these issues, Ultra-NeRF [29] propose an implicit neural representation with physics-based rendering for US imaging. It builds on the ultrasound volume rendering method by Salehi et al. [18] and combines a coordinatebased MLP that maps from 3D coordinates to five physical parameters, achieving high-quality synthetic ultrasound images. However, this approach is still viewindependent, preventing it from achieving ray interactions and the Fresnel Effect. It also does not allow for rendering complex ultrasound artifacts, such as reverberations. Our work draws on this combination of traditional ray tracing and physics-based volume rendering techniques, achieving better reconstruction results for 3D ultrasound scenes through implicit neural methods.

## 3 Methods

NeRF and its early variants perform poorly in parameterizing reflection directions and handling high-frequency information, especially in complex media structures and high-fidelity image rendering [5]. These models typically fail to

Fig. 2. Ultrasound Wave Interaction and Reflection Parameterization at the Medium Interface. (a) Shows the reflection and refraction phenomena when the ultrasound wave I i impinges on the medium interface, describing the basic physical processes and corresponding formulas. (b) Describes the process of calculating the specular reflection direction ˆ I rs by adjusting the macroscale normal vector ̂ n to the microfacet normal vector ̂ n ′ , highlighting the impact of the microfacet structure on reflection characteristics.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00860v3-figure-4ba4253f-09c8-4a51-bf01-0aab58ed2a7c.jpg)

accurately capture complex reflection phenomena, resulting in suboptimal rendering quality. Ultra-NeRF, the first to introduce NeRF concepts into ultrasound imaging, also faces limitations in handling light interactions and Fresnel effects, further affecting its performance in rendering complex artifacts and reflection characteristics [29]. To address these issues, we proposed a new model, UlRe-NeRF, which considers both reflection intensity and ultrasound volume rendering. Specifically, we designed an ultrasound neural rendering architecture based on implicit neural networks and volume rendering, optimizing the rendering process (Figure 3) to better simulate the medium's physical characteristics, artifacts, and deformations, accurately reproducing the propagation and reflection behavior of ultrasound waves within the medium. Our approach provides a more effective and precise solution for ultrasound imaging.

## 3.1 Ultrasound Reflection Direction Parameterization

To better simulate the reflection characteristics of ultrasound waves at the medium interface, we propose a new method for parameterizing reflection directions, as shown in Figure 2(a). The reflection and refraction properties of ultrasound waves, treated as rays, can be calculated using the laws of wave physics [20]. When an ultrasound wave I i impinges on the medium interface at a certain angle, reflection and refraction occur. The intensities of reflection and refraction can be calculated using the following formulas:

$$I _ { r } = I _ { i } \left ( \frac { Z _ { 2 } - Z _ { 1 } } { Z _ { 2 } + Z _ { 1 } } \right ) ^ { 2 }, \ \ I _ { t } = I _ { i } \left ( \frac { 4 Z _ { 2 } Z _ { 1 } } { Z _ { 2 } + Z _ { 1 } } \right ) ^ { 2 }.$$

To overcome the limitations of traditional methods in handling complex interfaces, we introduced the Phong model from the field of computer graphics

to accurately simulate specular reflection phenomena [16]. This model analyzes reflected light intensity by considering ambient light, diffuse reflection, and specular reflection:

$$\tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t {  \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t {  \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t { \tt t \tt                                                                                                                                                                                                         }                                                                                                                                                                                                       
                                                                                                                                                                                                        |                                                                                                                                                                                                        -                                                                                                                                                                                                        +                                                                                                                                                                                                       .                                                                                                                                                                                                        )                                                                                                                                                                                                        (                                                                                                                                                                                                        1                                                                                                                                                                                                        0                                                                                                                                                                                                        2                                                                                                                                                                                                        3                                                                                                                                                                                                        4                                                                                                                                                                                                        6                                                                                                                                                                                                        7$$

The Phong model's Bidirectional Reflectance Distribution Function (BRDF) exhibits rotational symmetry with respect to the reflection angle, thus satisfying f ( I i , I o ) = p ( I r · I o ) . When ignoring complex optical phenomena such as interreflections and self-occlusions, the radiance in the view direction can be determined solely as a function of the specular reflection direction ˆ I rs :

$$L _ { \text{out} } \, \left ( \hat { I } _ { o } \right ) \, \in \, \int L _ { \text{in} } \, \left ( \hat { I } _ { i } \right ) \, p \left ( \hat { I } _ { r s } \cdot \hat { I } _ { i } \right ) d \hat { I } _ { i } = F \left ( \hat { I } _ { r s } \right ) \, \text{ \quad \ \ } \, ( 6 )$$

Therefore, we parameterize the relationship between the specular reflection direction ˆ I rs and the viewing direction ˆ I o :

$$\hat { I } _ { r s } = 2 \left ( \hat { I } _ { o } \cdot \hat { \mathbf n } \right ) \hat { \mathbf n } - \hat { I } _ { o }$$

Where ˆ I o = -ˆ d represents the unit vector from the reflection point on the boundary to the transducer, and ̂ n is the normal vector at that point. In practice, the medium interface is usually not a perfect mirror but consists of many microfacets. This results in the macroscale normal vector ̂ n being a perturbed normal vector ̂ n ′ , also known as the microfacet normal. As shown in Figure 2(b), the corrected normal vector ̂ n ′ is calculated as follows:

$$\widehat { n } ^ { \prime } = \widehat { n } + \delta \left ( \hat { I } _ { o } - \widehat { n } \left ( \widehat { n } \cdot \hat { I } _ { o } \right ) \right )$$

Where δ represents the roughness at that location. By accounting for the microfacet structure of the medium, a more accurate specular reflection direction can be calculated:

$$\hat { I } _ { r s } = 2 \left ( \hat { \mathbf n } ^ { \prime } \cdot \hat { I } _ { o } \right ) \hat { \mathbf n } ^ { \prime } - \hat { I } _ { o }.$$

By parameterizing the specular reflection direction ˆ I rs , we input it into a multilayer perceptron (MLP), training the model to output the integrated BRDF as a function of ˆ I rs . This approach not only captures the reflection characteristics of ultrasound waves on various surfaces and materials but also allows for accurate simulation of complex scenes through model training. We also supplement the output of diffuse reflection using a spatial MLP, providing a more comprehensive description of the interaction between sound waves and interfaces. To better adapt to the actual conditions of ultrasound propagation, we introduce the angle between the normal vector and the outgoing direction ̂ n ′ · ˆ I o as an additional input parameter, enabling the model to dynamically adjust reflection characteristics.

## 3.2 Reflective Harmonic Encoding

In the traditional NeRF framework, positional encoding methods have limitations when handling high-frequency information, leading to an inability to accurately reproduce high-frequency details in the propagation of ultrasound waves [9] [14]. This poses challenges in simulating the behavior of ultrasound waves in various media, especially in biological tissues with complex structures and heterogeneous properties, where traditional methods fail to accurately capture reflection direction information. To address this issue, we drew inspiration from the Integrated Directional Encoding (IDE) method in Mip-NeRF [2]and applied it to ultrasound imaging, calling it Reflective Harmonic Encoding (RHE). RHE uses spherical harmonics to encode high-frequency information of reflection directions, making it particularly suitable for biological tissues with complex characteristics [1]. It finely encodes reflection directions, improving reconstruction accuracy. Specifically, we use the von Mises-Fisher (vMF) distribution to define the distribution of reflection directions and encode these directions using a set of spherical harmonics. First, we define the distribution of reflection directions on the unit sphere using the von Mises-Fisher distribution, centered at the reflection vector ˆ I r , with the concentration parameter κ defined as the reciprocal of roughness ρ , κ = 1 /ρ . Next, we represent the encoding of reflection directions using a weighted sum of spherical harmonic functions { Y m ℓ } . The encoding form under the vMF distribution is defined as:

$$\text{RHE} \left ( \hat { \boldsymbol I } _ { r }, \kappa \right ) = \sum _ { ( \ell, m ) \in \mathcal { L } } A _ { \ell } ( \kappa ) Y _ { \ell } ^ { m } \left ( \hat { \boldsymbol I } _ { r } \right ) \quad \text{(10)}$$

Where L represents the set of spherical harmonic functions. The weighting function A ℓ ( κ ) depends on the concentration parameter κ and can be approximated by the following formula:

$$A _ { \ell } ( \kappa ) \approx \exp \left ( - \frac { \ell ( \ell + 1 ) } { 2 \kappa } \right ).$$

This encoding method allows us to more accurately capture and express information about ultrasound wave reflections. Through RHE, we can finely encode these reflection directions, significantly improving the accuracy of boundary reconstructions. Especially when dealing with reflections in complex tissue structures and rich details, RHE provides higher resolution and accuracy.

## 3.3 Sine Activation Function

Sinusoidal representation networks (SIRENs) [22] utilize sine as the activation function in MLPs to model high-frequency information accurately. They introduced an initialization scheme to prevent vanishing gradients in periodic activation functions, showing that SIRENs fit complex signals and their derivatives robustly. This provides a solid foundation for using the Sine activation function in ultrasound image reconstruction.The periodicity and continuity of the Sine

activation function enable it to capture high-frequency components and subtle variations effectively. Its expression is:

$$\text{Sine } ( x ) = \sin ( x )$$

Compared to traditional activation functions like ReLU and Sigmoid, which are sensitive to initial parameters and unstable with complex signals, the Sine activation function is more stable and robust over a wide range of parameters [15]. In the Ultra-NeRF model, results were initially unstable and highly dependent on network configuration. By introducing the Sine activation function, our model demonstrated more consistent performance across different initial configurations, reducing the impact of initial settings on the final imaging outcome. This improvement enhances the model's robustness and stability, making it better suited for various configurations and complex ultrasound imaging scenarios.

## 3.4 Ultrasound Neural Rendering Architecture

To address the complexity of ultrasound scenes and accurately reconstruct ultrasound characteristics within the medium, we designed an ultrasound neural rendering architecture based on implicit neural networks and volume rendering (as shown in Figure 3). This architecture integrates two main modules, the directional MLP and the spatial MLP, and performs volume rendering based on ray tracing and physical principles, significantly improving the reconstruction accuracy of ultrasound images.

The primary function of the directional MLP is to process ray direction information. Specifically, it receives the ray direction ˆ d as input and generates view-dependent reflection intensity information. By generating the reflection intensity estimate ̂ C ref , it ensures that the reflection characteristics under each view are accurately captured, enhancing the ability to simulate complex reflection properties. The formula is as follows:

$$\dot { C } _ { r e f } = \mu \left ( \mathbf c _ { d } + \mathbf c _ { s } \right )$$

Where ̂ C ref is the reflection intensity estimate, c s is the specular reflection information encoded by direction, c d is the diffuse reflection information produced during the interaction of ultrasound waves with the interface, which is only position-dependent, and µ is a scaling function that limits the output value to the range [0,1].

The spatial MLP, on the other hand, is responsible for predicting a set of physical property parameters, including the attenuation coefficient α , reflectivity β , scattering density ρ s , and amplitude ϕ . We apply Fourier feature positional encoding to map the spatial position x into a high-dimensional space [23], capturing high-frequency details in ultrasound propagation. This information is then fed into the spatial MLP to generate the aforementioned physical property parameters, accurately reconstructing the ultrasound characteristics of the medium. During the rendering process, these physical properties are used to simulate

Fig. 3. Ultrasound Neural Rendering Architecture of UlRe-NeRF. This figure illustrates our proposed ultrasound neural rendering architecture. The view-dependent reflection intensity information and physical property parameters generated by the directional MLP and spatial MLP modules are combined in the rendering section, achieving high-fidelity simulation of ultrasound wave propagation and reflection behavior.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00860v3-figure-600ba113-f9a2-4f22-81b4-b1b0eb1ffd9c.jpg)

the propagation and reflection behavior of ultrasound waves in different media, resulting in the rendering outcome ̂ C render . Finally, the reflection information generated by the directional MLP is combined with the rendering outcome to produce the final ultrasound imaging output. By improving the network architecture and rendering process, UlRe-NeRF not only accurately captures and simulates the propagation and reflection behavior of ultrasound waves in different media but also dynamically adapts to view changes, enhancing the detail and accuracy of the reconstructed images. This approach excels in handling complex ultrasound imaging effects and enhances the model's robustness and adaptability, achieving high-fidelity simulation of ultrasound wave propagation and reflection.

## 3.5 US volume rendering

The traditional Neural Radiance Fields (NeRF) method generates images by accumulating point data in 3D space along the line of sight. However, this approach has significant limitations when applied to ultrasound imaging. The physical phenomena involved in ultrasound imaging are more complex, requiring consideration of wave effects and acoustic impedance matching, which are not factors in natural image rendering. Therefore, a more specialized approach is necessary. We improved the Ultra-NeRF method by proposing a volume rendering technique that utilizes the physical properties of the medium. This new method combines parameters generated by directional MLP and spatial MLP and incorporates classical ultrasound imaging theory and ray tracing technology [3],

thereby achieving more accurate simulations of ultrasound wave propagation and reflection to produce high-fidelity rendered images.

Specifically, we synthesize views by querying the ultrasound characteristics of the medium along the path of the sound wave: attenuation coefficient α , reflectivity β , scattering density ρ s , and amplitude ϕ . As shown in Figure 3, we treat ultrasound waves as rays emitted from the transducer, simulating their propagation paths through different media. For each scanline r , we define the ultrasound intensity C render ( r, t ) at a distance t as the sum of the reflected energy R r, t ( ) and backscattered energy B r, t ( ) , describing the energy changes of the ultrasound wave during propagation:

$$C _ { r e n d e r } \left ( r, t \right ) = R ( r, t ) + B ( r, t )$$

Based on the Beer-Lambert law, we introduce the concept of residual energy I ( r, t ) , whose variation is calculated using the following formula:

$$I ( r, t ) = I _ { 0 } \cdot \exp \left ( - \int _ { 0 } ^ { t - 1 } ( \alpha \cdot f \cdot \Delta t ) d \tau \right ) \cdot \prod _ { n = 0 } ^ { t - 1 } ( 1 - \beta ( r, n ) ) \quad \ \ ( 1 5 )$$

Where α is the attenuation coefficient of the medium, f is the frequency of the ultrasound wave, and β r, n ( ) is the reflection coefficient of the medium in the tissue. This formula describes the exponential decay of the signal as it propagates through the medium.

The reflected energy R r, t ( ) refers to the energy generated at the medium interface due to reflection. The calculation of R r, t ( ) is based on the residual energy β r, t ( ) , the reflection coefficient , and a predefined two-dimensional point spread function (PSF):

$$R ( r, t ) = I ( r, t ) \cdot \beta ( r, t ) \otimes \text{PSF} ( r ),$$

In our model, the point spread function (PSF) is simplified to include only the axial and lateral directions, excluding the elevation direction. The PSF function is approximated as a two-dimensional Gaussian function modulated by a cosine function, specifically:

$$P S F ( x, y ) = \exp \left ( - \frac { 1 } { 2 } \left ( \frac { x ^ { 2 } } { \sigma _ { x } ^ { 2 } } + \frac { y ^ { 2 } } { \sigma _ { y } ^ { 2 } } \right ) \right ) \cdot \cos ( 2 \pi f x )$$

Where σ x and σ y represent the standard deviations in the axial and lateral directions, respectively, and f is the frequency of the ultrasound wave.

Backscattered energy B r, t ( ) refers to the energy generated by scattering, which consists of the residual energy I ( r, t ) and the two-dimensional distribution of scatter points T r, t ( ) :

$$B ( r, t ) = I ( r, t ) \cdot S ( r, t ) \otimes \text{PSF} ( r ),$$

Here, the two-dimensional distribution of scatter points S r, t ( ) is defined as:

$$S ( r, t ) = \eta ( r, t ) \cdot \phi ( r, t )$$

where the function η r, t ( ) indicates the intensity of the scatterers. To add a degree of smoothness while maintaining boundary clarity, we use the Beta distribution to parameterize the scattering density ρ s , which provides better continuity compared to the Bernoulli distribution. The intensity of the scatter points is determined by their amplitude ϕ .

By integrating these energy estimates, we obtain the ultrasound intensity at different positions along each scanline. By normalizing and combining the intensity values ̂ C render generated by the rendering process with the reflection intensity ̂ C ref produced by the directional MLP, we get the final rendered result:

$$\widehat { C } _ { \text{final} } \, = \, \mathcal { G } \left ( \widehat { C } _ { \text{render} } \, + \widehat { C } _ { \text{ref} } \, \right ).$$

This improved method not only addresses the shortcomings of Ultra-NeRF in handling complex artifacts and reflection characteristics but also enhances rendering consistency and accuracy across different views and complex scenes. Through this new approach, UlRe-NeRF can more accurately simulate the propagation and reflection behavior of ultrasound waves in various media, significantly improving the physical realism and visual quality of the images.

## 4 Experiments

## 4.1 Experimental Setup

Our experiments were conducted on a computer equipped with an NVIDIA RTX 3070 GPU and 32GB of RAM. We utilized the Ultra-NeRF open-source dataset to simulate liver ultrasound images generated from CT scans, following the image generation methods outlined in Ultra-NeRF's work and using the ImFusion tool. The dataset includes seven scans: six at different oblique angles and one vertical scan. Each scan contains 200 2D ultrasound images along with their tracking information. The oblique scans provide images from different perspectives and directions, with frames containing occlusions caused by scanning between ribs. We used four oblique scans (800 frames) for training and one vertical scan plus two oblique scans (600 frames in total) for testing. We employed the same quantitative metrics as previous view synthesis studies [2] [14] [31]. Performance evaluation metrics include Peak Signal-to-Noise Ratio (PSNR), Structural Similarity Index (SSIM) [26], Learned Perceptual Image Patch Similarity (LPIPS) [32], and Mean Squared Error (MSE). Among these, LPIPS focuses on the global structure and noise levels of the images and performs exceptionally well in measuring the similarity of images with significant differences and noticeable disturbances, being highly correlated with human perceptual preferences [13].

## 4.2 Ablation Study Results of the Neural Network Architecture

Although NeRF and Mip-NeRF (the previous best-performing techniques) have shown outstanding performance in view synthesis and image reconstruction, their

(a)

(b)

Fig. 4. Our model demonstrates accurate and high-fidelity performance, which highlights two different frames of ultrasound images. (a) Overall image comparison; (b) Detailed section comparison.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00860v3-figure-6d300015-d438-4314-8ef1-9fb810160589.jpg)

ray integration-based rendering method has limitations in ultrasound imaging. The traditional NeRF method determines pixel color values by integrating along a ray, which is suitable for visible light images. However, ultrasound imaging involves interactions of ultrasound waves with different tissues along the propagation path. To address this, we made specific modifications to NeRF and MipNeRF to better handle the characteristics of ultrasound images. Specifically, we used a Directional MLP in the model to directly output simulated ultrasound images, capturing reflection information of ultrasound waves in different directions. Furthermore, we conducted ablation studies on the performance of each model using the ultrasound dataset.

In Figure 4, we compare the performance of different methods in ultrasound imaging. As shown in Figure 4(a), NeRF fails to accurately capture the complex reflection and scattering phenomena in ultrasound images, resulting in blurred image edges. Although Mip-NeRF shows some improvement, it still lacks refinement, with noticeable high-frequency detail and fine feature loss. Ultra-NeRF improves in detail restoration and visual effect but still suffers from noise and blurriness when handling complex tissues, and its robustness is poor, with performance heavily dependent on model initialization. Figure 4(b) shows a detailed comparison, where only our model accurately captures reflection and scattering information, exhibiting higher image clarity and less noise. By introducing the reflection direction parameterization and RHE, our method excels in restoring high-frequency details and edge information, with less noise and overall visual effects closer to real scenes, achieving high-fidelity image reconstruction. This advantage is particularly evident in the LPIPS metric. The quantitative results in Table 1 indicate that our model significantly outperforms others in the LPIPS metric, demonstrating a substantial advantage in perceptual image quality. The various improvement modules described in section 3 each contribute to this per-

Table 1. Baseline Comparisons and Ablation Study

|                        |   LPIPS↓ |   SSIM↑ |   PSNR↑ |    MSE↓ |
|------------------------|----------|---------|---------|---------|
| NeRF                   |    0.565 |   0.57  |   29.66 | 0.00143 |
| Mip-NeRF               |    0.529 |   0.561 |   29.41 | 0.00149 |
| Ultra-NeRF             |    0.321 |   0.449 |   20.51 | 0.0585  |
| Ours, no reflection    |    0.553 |   0.486 |   22.84 | 0.00857 |
| Ours, concat. viewdir  |    0.566 |   0.48  |   22.4  | 0.0101  |
| Ours, dotprod          |    0.551 |   0.493 |   23.1  | 0.0083  |
| Ours, no roughness     |    0.261 |   0.544 |   28.71 | 0.00178 |
| Ours, no diffuse color |    0.552 |   0.494 |   23.13 | 0.00826 |
| Ours, ReLu             |    0.547 |   0.585 |   29.67 | 0.00141 |
| Ours, PE               |    0.553 |   0.494 |   23.12 | 0.00829 |
| Ours, ReLu. PE         |    0.552 |   0.586 |   29.71 | 0.0014  |
| Ours                   |    0.25  |   0.546 |   28.75 | 0.00178 |

formance. When using view direction instead of reflection direction as input to the Directional MLP ("no reflection"), the model's reconstruction metrics notably decrease, highlighting the benefit of parameterizing reflective radiance. Adding view direction as an additional input to the Directional MLP ("concat. viewdir") severely impacts performance, underscoring the importance of parameterizing specular appearance as a function of view direction. Moreover, not feeding ̂ n ′ · ˆ I o into the Directional MLP leads to significant degradation in rendering. Using fixed RHE weight parameters slightly reduces performance, demonstrating the robustness of spherical harmonic encoding. Removing diffuse color ("no diffuse color") greatly affects performance, further proving the crucial role of diffuse reflection in a complete reflection model. Finally, replacing the Sine activation function with the ReLU activation function used in NeRF also impacts performance. Deleting RHE structural components (roughness, diffuse color, etc.) and replacing them with NeRF's standard positional encoding ("PE") shows a performance drop, further validating the effectiveness of our proposed improvements. The combination of ReLU activation function and positional encoding ("ReLu. PE") also shows a decline in performance. Our improvements significantly enhance detail restoration and overall image quality by better capturing high-frequency details and accurately representing ultrasound wave reflection information. Qualitative and quantitative analyses clearly show the superiority of our model in ultrasound image reconstruction tasks.

## 4.3 Ablation Study Results of the Rendering Components

To demonstrate the effectiveness of UlRe-NeRF's rendering method in ultrasound image reconstruction, we designed various ablation experiments comparing it to Ultra-NeRF. These experiments included replacing Ultra-NeRF's rendering method with ours and using Ultra-NeRF's rendering method in our model. Table 2 shows the quantitative comparison results of different models across several metrics. It is evident that UlRe-NeRF exhibits significant advantages in all

Table 2. Rendering Method Comparisons and Ablation Study

|                              |   LPIPS↓ |   SSIM↑ |   PSNR↑ |    MSE↓ |
|------------------------------|----------|---------|---------|---------|
| Ultra-NeRF                   |    0.321 |   0.449 |   20.51 | 0.0585  |
| Ultra-NeRF, our rendering    |    0.314 |   0.498 |   22.31 | 0.0475  |
| Ours, Ultra-NeRF's rendering |    0.361 |   0.409 |   24.97 | 0.0042  |
| Ours                         |    0.25  |   0.546 |   28.75 | 0.00178 |

Ground truth

Ultra-NeRF New

Ultra-NeRF

Ours

Ours Orig.

LPIPS /PSNR

0.314/22.31

0.321/20.51

0.250/28.75

0.361/24.97

Ļ

Ĺ

Fig. 5. Visual Comparison of Different Rendering Methods. Ours Orig. indicates the use of Ultra-NeRF's rendering method, while Ultra-NeRF New indicates the use of our rendering method.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00860v3-figure-c44209ac-9853-4770-9958-646b5582225d.jpg)

metrics, especially in significantly reducing the LPIPS value. This improvement is reflected in the visual quality of the images, with a noticeable suppression of noise and artifacts, resulting in clearer and more realistic images. In contrast, using Ultra-NeRF's original rendering method results in significant noise and detail loss in the images. Additionally, when Ultra-NeRF adopts our method, there is a significant improvement in image quality and metrics, indicating that our method better reconstructs ultrasound signals and effectively reduces noise interference. Quantitative analysis of the results shows that our method can more accurately capture and reconstruct details when handling complex tissue structures, significantly enhancing image clarity and realism.

Figure 5 presents the visual effects of four different rendering methods in ultrasound images. Images generated by Ultra-NeRF exhibit high levels of noise and shadow issues, particularly with blurred edges and details, leading to overall lower image quality. When Ultra-NeRF's rendering method is replaced with ours, there is a reduction in shadows and an improvement in detail clarity, although some noise remains. In contrast, our volume rendering method performs best in detail restoration and noise reduction, significantly enhancing the overall image quality. This improvement in visual effects is not only reflected in the clarity of details and edges but also in the overall contrast and visual consistency of the images. Using our rendering method, ultrasound images can better preserve

Fig. 6. Optimization Effects on Various Physical Property Parameters. We used normalized values to unify the scale for display in the same dimension.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00860v3-figure-6b509c6b-2269-429b-879d-efa169eace2e.jpg)

artifacts and real details when dealing with complex tissue structures, resulting in high-fidelity reconstructed images.

Figure 6 further illustrates the optimization effects on various physical property parameters, including attenuation coefficient α , reflectivity β , scattering density ρ s , and amplitude ϕ . For the attenuation coefficient α , our method provides a more accurate description of attenuation behavior, thereby improving image contrast and clarity. The optimized reflectivity β significantly enhances the precision of interface reflections, making reflection boundaries sharper and more distinct, ensuring clearer delineation between tissues. The scattering density ρ s optimized using the beta distribution more effectively describes scattering phenomena compared to the Bernoulli distribution, reducing blurriness caused by scattering. Finally, the optimized amplitude ϕ parameter improves signal consistency, resulting in clearer details within the image, especially in areas with complex tissue structures, further enhancing the realism and visual quality of the images.

## 5 Discussion

In this paper, we present UlRe-NeRF, a more effective and precise Neural Rendering method for 3D Ultrasound Imaging. We incorporate additional inductive biases about ultrasound: enhancing the implicit representation network with reflection direction parameterization and Reflective Harmonic Encoding (RHE) to capture complex tissue reflections, and improving explicit volumetric rendering by integrating reflection intensity and volume rendering to simulate the medium's physical properties and light interactions.

Our experiments show that UlRe-NeRF achieves high-fidelity 3D ultrasound imaging synthesis. It effectively learns view-dependent geometry, captures complex tissue structures, simulates various physical properties of the medium, and accurately reproduces the propagation and reflection behavior of ultrasound waves. This study marks an important advance in ultrasound image reconstruction techniques, demonstrating the potential of ultrasound neural rendering architectures based on implicit neural networks and explicit ultrasound volume rendering to improve image quality.

However, the inherent characteristics of ultrasonic imaging limit the usage of hierarchical sampling strategies. Thus, methods specifically designed for ultrasound rendering to enhance sample efficiency are required. Although we have introduced new explicit ultrasound volume rendering techniques to improve contrast and clarity in physical property simulations, further advancements are needed to develop more efficient and precise methods for simulating the propagation and reflection behavior of ultrasound waves. Additionally, in more complex freehand scanning situations, where ultrasound data is sparse, additional regularization for view synthesis from sparse inputs is necessary.

## 6 Conclusion

The UlRe-NeRF model represents a significant advancement in 3D ultrasound imaging by combining implicit neural networks with explicit volumetric rendering. Our approach successfully addresses the limitations of traditional ultrasound imaging methods, providing high-fidelity reconstructions of complex tissue structures and realistic simulations of wave propagation. While our work has made considerable progress, further research is needed to improve sample efficiency and address challenges associated with sparse data in freehand scanning scenarios. We anticipate that this study will inspire continued exploration and refinement of neural rendering techniques for medical applications, contributing to the future development of ultrasound imaging technology.

## References

- 1. Ahrens, J., Thomas, M.R., Tashev, I.: Hrtf magnitude modeling using a nonregularized least-squares fit of spherical harmonics coefficients on incomplete data. In: Proceedings of The 2012 Asia Pacific Signal and Information Processing Association Annual Summit and Conference. pp. 1-5. IEEE (2012)

- 2. Barron, J.T., Mildenhall, B., Tancik, M., Hedman, P., Martin-Brualla, R., Srinivasan, P.P.: Mip-nerf: A multiscale representation for anti-aliasing neural radiance fields. In: Proceedings of the IEEE/CVF international conference on computer vision. pp. 5855-5864 (2021)
- 4. Corona-Figueroa, A., Frawley, J., Bond-Taylor, S., Bethapudi, S., Shum, H.P., Willcocks, C.G.: Mednerf: Medical neural radiance fields for reconstructing 3d-aware ct-projections from a single x-ray. In: 2022 44th annual international conference of the IEEE engineering in medicine &amp; Biology society (EMBC). pp. 3843-3848. IEEE (2022)
- 3. Burger, B., Bettinghausen, S., Radle, M., Hesser, J.: Real-time gpu-based ultrasound simulation using deformable mesh models. IEEE transactions on medical imaging 32 (3), 609-618 (2012)
- 5. Debbagh, M.: Neural radiance fields (nerfs): A review and some recent developments. arXiv preprint arXiv:2305.00375 (2023)
- 7. Hsu, P.W., Prager, R.W., Gee, A.H., Treece, G.M.: Freehand 3d ultrasound calibration: a review. Advanced imaging in biology and medicine: technology, software environments, applications pp. 47-84 (2009)
- 6. Fenster, A., Downey, D.B., Cardinal, H.N.: Three-dimensional ultrasound imaging. Physics in medicine &amp; biology 46 (5), R67 (2001)
- 8. Huang, Q., Zeng, Z.: A review on real-time 3d ultrasound imaging technology. BioMed research international 2017 (1), 6027029 (2017)
- 10. Li, H., Chen, H., Jing, W., Li, Y., Zheng, R.: 3d ultrasound spine imaging with application of neural radiance field method. In: 2021 IEEE International Ultrasonics Symposium (IUS). pp. 1-4. IEEE (2021)
- 9. Huang, X., Li, W., Hu, J., Chen, H., Wang, Y.: Refsr-nerf: Towards high fidelity and super resolution view synthesis. In: Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition. pp. 8244-8253 (2023)
- 11. Max, N.: Optical models for direct volume rendering. IEEE Transactions on Visualization and Computer Graphics 1 (2), 99-108 (1995)
- 13. Mier, J.C., Huang, E., Talebi, H., Yang, F., Milanfar, P.: Deep perceptual image quality assessment for compression. In: 2021 IEEE International Conference on Image Processing (ICIP). pp. 1484-1488. IEEE (2021)
- 12. Mercier, L., Langø, T., Lindseth, F., Collins, D.L.: A review of calibration techniques for freehand 3-d ultrasound systems. Ultrasound in medicine &amp; biology 31 (4), 449-471 (2005)
- 14. Mildenhall, B., Srinivasan, P.P., Tancik, M., Barron, J.T., Ramamoorthi, R., Ng, R.: Nerf: Representing scenes as neural radiance fields for view synthesis. Communications of the ACM 65 (1), 99-106 (2021)
- 16. Phong, B.T.: Illumination for computer generated pictures. In: Seminal graphics: pioneering efforts that shaped the field, pp. 95-101 (1998)
- 15. Parascandolo, G., Huttunen, H., Virtanen, T.: Taming the waves: sine as activation function in deep neural networks (2016)
- 17. Rahaman, N., Baratin, A., Arpit, D., Draxler, F., Lin, M., Hamprecht, F., Bengio, Y., Courville, A.: On the spectral bias of neural networks. In: International conference on machine learning. pp. 5301-5310. PMLR (2019)
- 18. Salehi, M., Ahmadi, S.A., Prevost, R., Navab, N., Wein, W.: Patient-specific 3d ultrasound simulation based on convolutional ray-tracing and appearance optimization. In: Medical Image Computing and Computer-Assisted Intervention-MICCAI 2015: 18th International Conference, Munich, Germany, October 5-9, 2015, Proceedings, Part II 18. pp. 510-518. Springer (2015)

- 19. Schwarz, K., Liao, Y., Niemeyer, M., Geiger, A.: Graf: Generative radiance fields for 3d-aware image synthesis. Advances in Neural Information Processing Systems 33 , 20154-20166 (2020)
- 21. Shen, L., Pauly, J., Xing, L.: Nerp: implicit neural representation learning with prior embedding for sparsely sampled image reconstruction. IEEE Transactions on Neural Networks and Learning Systems 35 (1), 770-782 (2022)
- 20. Shams, R., Hartley, R., Navab, N.: Real-time simulation of medical ultrasound from ct images. In: Medical Image Computing and Computer-Assisted InterventionMICCAI 2008: 11th International Conference, New York, NY, USA, September 6-10, 2008, Proceedings, Part II 11. pp. 734-741. Springer (2008)
- 22. Sitzmann, V., Martel, J., Bergman, A., Lindell, D., Wetzstein, G.: Implicit neural representations with periodic activation functions. Advances in neural information processing systems 33 , 7462-7473 (2020)
- 24. Tewari, A., Thies, J., Mildenhall, B., Srinivasan, P., Tretschk, E., Yifan, W., Lassner, C., Sitzmann, V., Martin-Brualla, R., Lombardi, S., et al.: Advances in neural rendering. In: Computer Graphics Forum. vol. 41, pp. 703-735. Wiley Online Library (2022)
- 23. Tancik, M., Srinivasan, P., Mildenhall, B., Fridovich-Keil, S., Raghavan, N., Singhal, U., Ramamoorthi, R., Barron, J., Ng, R.: Fourier features let networks learn high frequency functions in low dimensional domains. Advances in neural information processing systems 33 , 7537-7547 (2020)
- 25. Verbin, D., Hedman, P., Mildenhall, B., Zickler, T., Barron, J.T., Srinivasan, P.P.: Ref-nerf: Structured view-dependent appearance for neural radiance fields. In: 2022 IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR). pp. 5481-5490. IEEE (2022)
- 27. Wu, Q., Li, Y., Sun, Y., Zhou, Y., Wei, H., Yu, J., Zhang, Y.: An arbitrary scale super-resolution approach for 3d mr images via implicit neural representation. IEEE Journal of Biomedical and Health Informatics 27 (2), 1004-1015 (2022)
- 26. Wang, Z., Bovik, A.C., Sheikh, H.R., Simoncelli, E.P.: Image quality assessment: from error visibility to structural similarity. IEEE transactions on image processing 13 (4), 600-612 (2004)
- 28. Wu, Q., Li, Y., Xu, L., Feng, R., Wei, H., Yang, Q., Yu, B., Liu, X., Yu, J., Zhang, Y.: Irem: High-resolution magnetic resonance image reconstruction via implicit neural representation. In: Medical Image Computing and Computer Assisted Intervention-MICCAI 2021: 24th International Conference, Strasbourg, France, September 27-October 1, 2021, Proceedings, Part VI 24. pp. 65-74. Springer (2021)
- 30. Yeung, P.H., Hesse, L., Aliasi, M., Haak, M., Xie, W., Namburete, A.I., et al.: Implicitvol: Sensorless 3d ultrasound reconstruction with deep implicit representation. arXiv preprint arXiv:2109.12108 (2021)
- 29. Wysocki, M., Azampour, M.F., Eilers, C., Busam, B., Salehi, M., Navab, N.: Ultranerf: Neural radiance fields for ultrasound imaging. In: Medical Imaging with Deep Learning. pp. 382-401. PMLR (2024)
- 31. Zhang, K., Luan, F., Wang, Q., Bala, K., Snavely, N.: Physg: Inverse rendering with spherical gaussians for physics-based material editing and relighting. In: Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition. pp. 5453-5462 (2021)
- 32. Zhang, R., Isola, P., Efros, A.A., Shechtman, E., Wang, O.: The unreasonable effectiveness of deep features as a perceptual metric. In: Proceedings of the IEEE conference on computer vision and pattern recognition. pp. 586-595 (2018)
Ziwen Guo Zi Fang Zhuang Fu
UlRe-NeRF: 3D Ultrasound Imaging through Neural Rendering with Ultrasound Reflection Direction Parameterization
Three-dimensional ultrasound imaging is a critical technology widely used in
medical diagnostics. However, traditional 3D ultrasound imaging methods have
limitations such as fixed resolution, low storage efficiency, and insufficient
contextual connectivity, leading to poor performance in handling complex
artifacts and reflection characteristics. Recently, techniques based on NeRF
(Neural Radiance Fields) have made significant progress in view synthesis and
3D reconstruction, but there remains a research gap in high-quality ultrasound
imaging. To address these issues, we propose a new model, UlRe-NeRF, which
combines implicit neural networks and explicit ultrasound volume rendering into
an ultrasound neural rendering architecture. This model incorporates reflection
direction parameterization and harmonic encoding, using a directional MLP
module to generate view-dependent high-frequency reflection intensity
estimates, and a spatial MLP module to produce the medium's physical property
parameters. These parameters are used in the volume rendering process to
accurately reproduce the propagation and reflection behavior of ultrasound
waves in the medium. Experimental results demonstrate that the UlRe-NeRF model
significantly enhances the realism and accuracy of high-fidelity ultrasound
image reconstruction, especially in handling complex medium structures.