k-Wave Toolbox

makeSphericalSection

Create a binary map of a sphere segment within a 3D grid

Syntax

ss = makeSphericalSection(radius, height)
ss = makeSphericalSection(radius, height, width)
ss = makeSphericalSection(radius, height, width, plot_section)
ss = makeSphericalSection(radius, height, [], plot_section)
ss = makeSphericalSection(radius, height, width, plot_section, binary)
ss = makeSphericalSection(radius, height, [], [] binary)

[ss, dist_map] = makeSphericalSection(radius, height)
[ss, dist_map] = makeSphericalSection(radius, height, width)
[ss, dist_map] = makeSphericalSection(radius, height, [], plot_section)
[ss, dist_map] = makeSphericalSection(radius, height, width, plot_section)
[ss, dist_map] = makeSphericalSection(radius, height, width, plot_section, binary)
[ss, dist_map] = makeSphericalSection(radius, height, [], [] binary)

Description

makeSphericalSection creates a binary map of a section of a spherical surface within a three-dimensional matrix. The sphere is created using an extension of the midpoint circle algorithm. A single grid point is taken as the sphere center so the total diameter will always be an odd number of grid points. The sphere is then truncated based on the values for height and width (a diagram of the input sizes is given below). The face of the spherical section faces in the positive x-direction and the optional width parameter truncates the size in the y-direction.

If the optional input parameter binary is set to false (the default), the section map is returned as a double precision matrix. If it is set to true, the map is returned as a logical matrix. The average distance between each grid point in the spherical section and its contiguous neighbours can also be returned. This is given as a ratio compared to the average neighbour distance for a flat surface.

Examples

makeSphericalSection(25, 6, 15, true);

makeSphericalSection(50, 20, [], true);

Inputs

`radius`	radius of curvature [grid points]
`height`	transducer height [grid points]

Optional Inputs

`width`	section width (must be specified as an odd number) [grid points]
`plot_section`	Boolean controlling whether the spherical section is plotted using `voxelPlot` (default = `false`)
`binary`	Boolean controlling whether the spherical section is returned as a double precision matrix (`false`) or a logical matrix (`true`) (default = `false`)

Outputs

`ss`	binary matrix containing spherical section
`dist_map`	ratio of average neighbour distance for each grid point within the spherical section compared to a flat surface