Daniel Stevenson1, Grant Searchfield2, Xun Xu1. 1. Department of Mechanical Engineering, University of Auckland, New Zealand. 2. Audiology Section, University of Auckland, New Zealand g.searchfield@auckland.ac.nz.
Abstract
The main aim of this study was to investigate the shape variation in the human ear canal and the effects of venting on the spatial design of hearing aids. The second aim was to determine the design feasibility of a multiple venting hearing aid and assess the flexibility of design provided. A statistical shape model based on principal component analysis was created from a dataset of 60 left and 49 right ears. The modal variations of these models were then examined to determine the narrowest portion of the ear canal likely to limit effective venting. Finally, 3D models of two hearing aid shells, one with multiple 0.4-mm vents and the other with a single large 3-mm vent were created. Results showed that more than 50% of the shape variation in the human ear canal can be described by the first three modes of the statistical shape model developed. The narrowest predicted variation of this model had a minimum area of 36.4 mm(2), and the mean ear shape was found to have a minimum area of approximately 48 mm(2). It is estimated that even with a conservative vent packing of 0.4, multiple venting equivalent to at least a single 2-mm vent is achievable. The predicted variation in the human ear canal provides adequate physical space for a feasible multiple vented hearing aid shell. Furthermore, as multiple small vents are able to fit in around other components in the hearing aid design, certain design flexibility is provided by this venting approach.
The main aim of this study was to investigate the shape variation in the human ear canal and the effects of venting on the spatial design of hearing aids. The second aim was to determine the design feasibility of a multiple venting hearing aid and assess the flexibility of design provided. A statistical shape model based on principal component analysis was created from a dataset of 60 left and 49 right ears. The modal variations of these models were then examined to determine the narrowest portion of the ear canal likely to limit effective venting. Finally, 3D models of two hearing aid shells, one with multiple 0.4-mm vents and the other with a single large 3-mm vent were created. Results showed that more than 50% of the shape variation in the human ear canal can be described by the first three modes of the statistical shape model developed. The narrowest predicted variation of this model had a minimum area of 36.4 mm(2), and the mean ear shape was found to have a minimum area of approximately 48 mm(2). It is estimated that even with a conservative vent packing of 0.4, multiple venting equivalent to at least a single 2-mm vent is achievable. The predicted variation in the human ear canal provides adequate physical space for a feasible multiple vented hearing aid shell. Furthermore, as multiple small vents are able to fit in around other components in the hearing aid design, certain design flexibility is provided by this venting approach.
Venting is incorporated into hearing aid design to perform pressure equalization and reduce occlusion effects. Traditional venting of a single parallel vent channel, running from the base plate to the end plate of the hearing aid, is still widely used in hearing aid fitting. Vent diameters vary but normally range between 1 and 3 mm. A number of studies have been conducted looking at the acoustics of various hearing aid venting styles including traditional parallel vents, reverse horn vents, Y vents, and open fitting hearing aids (Kiessling, Brenner, Jespersen, Groth, & Jensen, 2005; Kuk, Keenan, & Lau, 2005, 2009; Kuk, Keenan, Lau, Dinulescu, et al., 2005; Stuart, Allen, Downs, & Carpenter, 1999; Studebaker & Cox, 1977). Most studies focus on the acoustic properties of these vents with little comment about the spatial interactions within the hearing aid or earmold. However, some vents such as diagonal or Y vents and external trench vents have been recommended as particularly useful for those with narrow ear canals (Dillon, 2001; Studebaker & Cox, 1977; Valente, Hosford-Dunn, & Roeser, 2008). Diagonal vents branch diagonally off the sound bore in hearing aid earmolds rather than running parallel to it. This has been shown to limit the high frequency gain available to the hearing aid (Studebaker & Cox, 1977). As such, diagonal vents have been recommended to be used only in cases where there is not enough space for other venting styles.The authors have proposed an alternative venting configuration that uses multiple small vents as opposed to a single large vent. An earlier publication by the authors looking at the acoustics of multiple vented earmolds concluded that multiple vented earmolds are capable of achieving similar feedback and occlusion performance but are likely to require a greater area to achieve this venting (Stevenson, Searchfield, Dodd, & Xu, 2013).This article explores the spatial constraints of a universal in-the-ear hearing aid and the feasibility of a multiple vented configuration within the predicted ear canal shape variations. Understanding of shape and spatial requirements is critical in the design of hearing aids that are comfortable as well as providing good retention and acoustic seal with the ear canal. To take advantage of a specific venting configuration, nonreproducible slit leakages around the ear must be reduced. Unlike customized hearing aids, when designing for a universal style device or a one-size-fits-many device, exact representation of each ear canal cannot be used. As such, one must obtain an understanding of the variances within the target population and design for these variances.
Shape Data Modeling and Analysis
Statistical Shape Model
A statistical shape model has been used to define the shape variations of the ear canal within the population. The statistical shape model used is based on Cootes, Edwards, and Taylor (1998, 2001), Hutton, Buxton, and Hammond (2001), and Paulsen, Larsen, Nielsen, Laugesen, and Ersbøll (2002). Principal component analysis (PCA) of a set of thin plate spline-wrapped ear impression scans was used to determine the main modes of variation within the dataset.A set of 60 left ears and 49 right ears has been used in the model to assess expected shape variation within the ear canal. On each of the scans in the dataset, 15 landmark points were identified. Figure 1 shows a sample landmarked shape. Numbers mark the superior, inferior, anterior, and posterior points on the first (1–4) and second bends (5–8) identified as the points of maximum curvature of the ear canal. Landmark point 9 is located at the crus of the helix where the antihelix meets the helix. Point 10 relates to the point where the helix meets the plane of the concha, and point 11 is the point directly opposite this on the edge of the concha. Finally, points 12–15 are evenly placed along the edge of the concha defined as the area of maximum curvature along the antihelical path.
Figure 1.
Mesh showing 15 landmark locations on sample earmold scan, Xs representing landmark points hidden from view.
Mesh showing 15 landmark locations on sample earmold scan, Xs representing landmark points hidden from view.A template mesh is then wrapped to each shape instance using the landmarked points through a thin plate spline technique. Wrapped data are then aligned in 3D space using Procrustes alignment, and PCA is performed. PCA gives the eigenvectors describing the modes of variation and the eigenvalues describing the relative importance of each mode. For a review of PCA, Procrustes alignment, and thin plate spline techniques, see Bookstein (1989), Cootes, Taylor, Cooper, and Graham (1995), Hutton et al. (2001), Jolliffe (1986), Lorenz and Krahnstöver (2000).
Minimum Cross-sectional Area
Given that a hearing aid pressure equalizing vent is required to provide a path for sound from one side of a hearing aid to the other, the spatially limiting factor is likely to occur at the narrowest point in the ear canal. We will call this point the minimum cross-sectional area. To determine the minimum cross-sectional area, a series of cross sections were taken normal to a central axis of the ear canal. The ear canal central axis has been defined by determining the shortest path along the medial sheet between a central point at the top of the ear canal and at the entrance of the ear canal.The medial sheet forms a part of the medial axis transform. The medial axis transform is a commonly used 3D shape descriptor for shape recognition, modeling, and manipulation (Crimi et al., 2011; Gerig, Styner, Shenton, & Lieberman, 2001; Styner & Gerig, 2001a, 2001b). Shapes are represented by a set of maximal medial balls within the interior of the solid (M-rep) as opposed to the points on the boundary (B-rep; Crimi et al., 2011; Pizer, Thall, & Chen, 1999; Styner & Gerig, 2001b). Maximal medial balls are the largest sphere that does not contain points on the surface of the shape and are not contained completely within another medial ball. The center of the medial ball is a point that has two or more closest points on the surface. The medial sheet is the union of the center of these medial balls within the shape forming a skeletal sheet structure within the shape.An approximation of the medial sheet has been obtained through the use of a 3D Voronoi diagram or Voronoi skeleton (Fabbri, Estrozi, & Costa, 2002; Ogniewicz & Ilg, 1992; Styner & Gerig, 2001b). The Voronoi diagram divides the shape into a set of cells with the edges of each cell containing the points equidistant between two shape vertices. The vertices of these Voronoi cells contained within the shape represent the center of the medial balls of the shape. However, Voronoi vertices are also often formed near the surface of the shape far from the medial axis or medial sheet even for very smooth scans. As such, the vertices must be pruned to include only those near the medial axis or center of the shape.Pruning removes sheet branches of medial skeletons not indicative of the skeleton structure of the shape whilst maintaining the topology of the original shape. A number of different skeleton pruning or thinning algorithms exist for reducing the complexity of medial axis shape representations based on medial or object angle (Dey & Zhao, 2003; Miklos, Giesen, & Pauly, 2010), segmentation by shape significance (Beristain, Graña, & Gonzalez, 2012; Tam & Heidrich, 2003), and the so-called λ-medial axis method (Chaussard, Couprie, & Talbot, 2009, 2011; Chazal & Lieutier, 2005). The method selected is similar to that of λ-medial axis that prunes medial balls based on the size of the circumradius of the medial boundary values or closest surface points. In this study, the medial balls were pruned based purely on the radius of the medial balls in question, which is the distance of the medial ball from the surface. Medial balls with radius (r) less than or equal to approximately 2 or 3 mm were pruned removing small medial balls describing detailed surface information.Pruning small medial balls related to surface detail provides a good approximation of the medial axis. For all shapes analyzed, the medial sheet was viewed to assess whether an appropriate r was used and that the pruned medial sheet provided a good approximation of the true medial axis. Figure 2 shows the medial sheets generated using r of (a) 0 and (b) 3 mm for mean left ear shape. As r increased, more of the vertices closer to the surface are removed, and the triangulation on the medial sheet improves, particularly for the ear canal, which is the area of interest. These small radius vertices describe the fine surface detail of the shape not needed in the medial axis. The method reduced the medial sheet within the ear scan to a thin sheet within the center of the ear canal in the area of interest in the scan (Figure 2b).
Figure 2.
Plots showing medial sheet for the mean left ear generated from medial balls with radius r greater than or equal to (a) 0 mm and (b) 3 mm.
Plots showing medial sheet for the mean left ear generated from medial balls with radius r greater than or equal to (a) 0 mm and (b) 3 mm.A fast marching algorithm that uses a finite difference approximation was used to determine the shortest path along the medial sheet between the two central points. Planar cross sections were taken approximately normal to the ear canal at both the entrance to the ear canal and just beyond the second bend. The centers of gravity of these sections were then taken, and the closest points to these on the medial sheet were used to determine the two end points of the centerline ear canal axis. The fast marching algorithm was then used to determine the shortest path between these two end points along the medial sheet.The fast marching algorithm has been used for a number of other medical applications involving tubular structures such as vessels, veins, and airways (Antiga, Ene-Iordache, & Remuzzi, 2003; Cárdenes, Bogunovic, & Frangi, 2010; Rouchdy & Cohen, 2011). As the fast marching algorithm allows for calculation of the shortest paths across faces and not just along edges, it was found to provide a smoother centerline and to lead to less variability in curve normals compared with other shortest path methods such as the Dijkstra algorithm. For example, Figure 3 shows the fast marching shortest path on the left mean ear medial sheet.
Figure 3.
Plot of shortest path between two points on the mean left ear medial sheet using the fast marching algorithm.
Plot of shortest path between two points on the mean left ear medial sheet using the fast marching algorithm.The Creo 3D solid modeling package (PTC Creo, Needham, MA, USA) was used to create a spline defining the path of fast marching shortest path, and a series of cross sections of the ear canal were taken normal to the spline curve. Eleven evenly spaced cross sections were taken from the bottom of the centerline to the top with distance between the cross sections ranging from 0.88 to 1.78 mm along the centerline curve. From these, a cross-sectional area profile of each studied shape was determined. The cross section approximating the point of the second and first bends was identified and the relative centerline distance between these two features determined.A closed cross-sectional area could not always be achieved. The plane normal to the centerline curve would sometimes section the scan mesh outside of the ear canal into the concha and outside the ear itself. This occurred more in wider canal shapes that had a sharp first bend of the ear canal.
Results
Shape Model
The modes of variation of the left and right data models against the percentage of variation in the training dataset explained by each mode are shown in Figure 4. The first mode explains 23% of the variation in the left ear shapes and 32% of variation in the right ears. More than 50% of the variation in both datasets is explained by the first three modes of variation. This shows similar results to that observed by Paulsen et al. (2002), who used scans of 29 ear impressions to train a deformable shape model to aid in automated hearing aid componentry placement (approximately 26% on the first mode and 51% in the first three modes). That is, doubling the dataset size had little effect on the relative importance of the first few modes. This suggests that these modes are a good indication of variation in the population as a whole. Beyond the first nine modes, the variation explained by each mode is less than what you would expect from random data, and the usefulness of these modes to describe real statistical shape variation is questionable.
Figure 4.
Plot of the percentage of total variation explained by each mode of variation in the left ear and right ear training dataset, as well as results expected from randomizing the dataset.
Plot of the percentage of total variation explained by each mode of variation in the left ear and right ear training dataset, as well as results expected from randomizing the dataset.By varying the weights of the eigenvectors, new shape instances created from the model are well defined and similar to those used in the training set. Furthermore, data within the training set can be reproduced by the model accurately using a limited number of modes. Pairwise plots of the first five eigenvector weights b used to recreate the shape instances in the training dataset indicate that we can treat these modes as independent of one another. No statistically significant correlation existed between any of the first 48 most important modes of the ear shape models. That is, the probability p of getting a correlation between any two modes as large as the observed value by random chance, when the true correlation is zero is greater than 0.05.
Cross-sectional Area
The cross-sectional area against distance from first ear canal bend for the ±3 standard deviations on the first three modes is shown in Figure 5. Highlighted is the cross-sectional area profile of the mean left and right ear (solid gray and black lines, respectively) and the narrowest left and right ear shape instances (dashed gray and black lines, respectively). Both of these instances occurred on the ±3 standard deviations of mode 1, and the widest ear canals were also observed on this mode. The mean distance along the centerline to the second bend is shown. The modal effects on the distance between first and second bend and the cross-sectional area at each mode are shown in Table 1.
Figure 5.
Scatter plot of the cross-sectional area against distance from first ear canal bend along ear canal centerline. Solid gray and black lines are the cross-sectional profile of the mean left and right ear shapes, respectively. Dashed gray and black lines are the cross-sectional profile of the narrowest left and right ear studied.
Table 1.
Cross-sectional Areas and Distance Between the First and Second Bends of the ±3 Standard Deviation Outputs on the First Three Modes of the Left and Right Ear Shape Models.
Measure
Ears
M
Mode 1
Mode 2
Mode 3
+3 SD
−3 SD
+3 SD
−3 SD
+3 SD
−3 SD
Distance between first and second bend (mm)
Left
7.5
5.4
8.9
8.4
6.9
6.3
6.7
Right
8.8
9.7
6.1
7.1
8.3
7.1
6.9
Cross-sectional area at first bend (mm2)
Left
69.8
–
46.3
83.7
56.3
60.6
71.5
Right
78.6
40.3
–
75.4
68.6
91.8
56.7
Cross-sectional area at second bend (mm2)
Left
54.9
59.5
40.7
65.4
44.8
53.7
58.3
Right
53.6
38.7
63.7
44.9
63.2
56.7
52.3
Note. Dashes represent where a valid cross section could not be obtained.
Scatter plot of the cross-sectional area against distance from first ear canal bend along ear canal centerline. Solid gray and black lines are the cross-sectional profile of the mean left and right ear shapes, respectively. Dashed gray and black lines are the cross-sectional profile of the narrowest left and right ear studied.Cross-sectional Areas and Distance Between the First and Second Bends of the ±3 Standard Deviation Outputs on the First Three Modes of the Left and Right Ear Shape Models.Note. Dashes represent where a valid cross section could not be obtained.There is a narrowing down the ear canal toward the second bend on all the shapes studied (Figure 5). The area of the ear canal toward the end ranged from 36.4 to 65.4 mm2. Both mean left and right ears were narrowest at the end, with an area of approximately 48 mm2.
Estimated Available Vents
The cross section of the receiver element is also likely to reduce the usable cross-sectional area at the ear canal’s narrowest point. For this reason, the available cross-sectional area was reduced by 7.29 mm2, corresponding to a receiver with dimension matching that of an Oticon Digifocus II CIC hearing aid (2.7 × 2.7 × 5.9 mm; Smørum, Denmark). Knowing the minimum cross-sectional area available for vents, we can approximate how many small circular vents would fit by estimating the relative packing density of the circles.It has been shown that the most efficient packing of circles in a plane is in a hexagonal lattice, achieving a packing density of approximately 0.906 (Chang & Wang, 2010; Thue, 1910; Tóth, 1943). Optimal packing within bounded areas is seen to be less area efficient depending on the bounding area and the size of the circles. Exact mathematical solutions exist to the circle packing problem for small numbers of circles bounded by circles and squares (Fodor, 2000, 2003; Graham, 1968; Melissen, 1994; Nurmela & Östergård, 1999; Peikert, Würtz, Monagan, & Groot, 1992; Pirl, 1969). A referenced summary of the best-known solutions (although not always proven optimal) for packing of N equal-sized circles in circles and squares, up to N = 1,500 and N = 10,000, respectively, can be found on Specht’s Web site www.packomania.com. Solutions for a variety of other geometrical domains are also available in the literature (Birgin, Bustamante, Callisaya, & Martínez, 2013; Birgin & Gentil, 2010; Galiev & Lisafina, 2013; Specht, 2013). For a review of literature on circle packing problems, see Hifi and M’Hallah (2009). Of particular interest to this research, given the approximate shape of a human ear canal, is Birgin et al.’s (2013) work on packing circles into ellipses. Birgin et al. found lower bound packing density of ellipses ranged from 0.44 to 0.83 with varying ellipse dimensions. Although the ear canal is not perfectly elliptical and the hearing aid receiver provides a packing obstacle, a conservative packing density of approximately 0.4 seems reasonable.To test this assumption, a hexagonal lattice pattern fill of 0.4-mm vents (outside diameter 0.8 mm) was applied to a narrowest ear canal cross section of 36.4 mm2 (Figure 6). A total of 30 vents were found to fit using the pattern given a packing density of approximately 0.5. Packing densities of larger vents were seen to be lower with 2- and 3-mm vents only achieving packing densities of approximately 0.46 and 0.31, respectively, on the same cross section.
Figure 6.
Sketched hexagonal packed 0.4-mm vents on narrowest ear canal cross section, square representing the hearing aid receiver element.
Sketched hexagonal packed 0.4-mm vents on narrowest ear canal cross section, square representing the hearing aid receiver element.The estimated number of vents N of varying diameters that could fit into the narrowest ear assuming a given packing density from 0.3 to 0.7 is shown in Table 2. Calculations assume that vents are spaced by at least 0.4 mm. Also shown for comparison is an equivalent vent diameter d for a traditional single vent with equivalent area. We see that even if we assume a conservative packing density of 0.4, the estimated total vent area available will match that of a traditional 2-mm vent.
Table 2.
Estimated Number of Possible Vents N of Varying Vent Diameters Assuming a Given Circle Packing Density.
Vent diameter
Packing density
0.3
0.4
0.5
0.6
0.7
N
deq
N
deq
N
deq
N
deq
N
deq
3 mm
0
0
1
3.0
1
3.0
1
3.0
2
4.2
2 mm
1
2.0
2
2.8
3
3.5
3
3.5
4
4.0
1 mm
5
2.2
7
2.6
9
3.0
11
3.3
13
3.6
0.8 mm
7
2.1
10
2.5
12
2.8
15
3.1
18
3.4
0.6 mm
11
2.0
14
2.2
18
2.5
22
2.8
25
3.0
0.5 mm
13
1.8
18
2.1
22
2.3
27
2.6
32
2.8
0.4 mm
17
1.6
23
1.9
28
2.1
34
2.3
40
2.5
Note. Area equivalent single vent diameter (d, in mm) is also shown.
Estimated Number of Possible Vents N of Varying Vent Diameters Assuming a Given Circle Packing Density.Note. Area equivalent single vent diameter (d, in mm) is also shown.
Spatial Design Flexibility
To determine the spatial design flexibility provided by a multivented hearing aid, 3D computer-aided design models of two in-the-canal or minicanal hearing aid models have been created. The hearing aid shell shape used for both models was based on the mean right ear mesh determined by the statistical shape model. Models also include electronic components in gray representing an A10 hearing aid battery as well as the chipset and receiver of an Oticon Digifocus II CIC hearing aid. Model A is a multivented hearing aid incorporating 16, 0.5-mm hearing aid vent paths ranging from 11.8 to 19.3 mm long (Figure 7). Model B is a single-vented hearing aid with a large 3-mm vent, 16.2 mm in length (Figure 8).
Figure 7.
3D model A of right ear ITC hearing aid incorporating 16, 0.5-mm vents.
Figure 8.
3D model B of right ear ITC hearing aid incorporating a single, large 3-mm vent.
3D model A of right ear ITC hearing aid incorporating 16, 0.5-mm vents.3D model B of right ear ITC hearing aid incorporating a single, large 3-mm vent.Figures 9 and 10 show views of the models A and B with the end and base plates removed. There is potentially more space for venting within the multivented model. However, 16 vents (equivalent to a 2-mm single vent) were chosen, and there is a clear view of the vent layout. The large vent used in model B causes the other componentry of the hearing aid to be located more toward the right of the model (Figures 7 to 10). The size of the vent was largely set by how compact the other componentry elements could be positioned, as a single large area is required. In contrast, we see that the multiple vented earmold has the smaller vents located in the gaps around the outside of the model. The multiple vents fit in around the important electronic componentry more easily and are able to fill the free space in the shell the larger vents are not. This provides considerable flexibility of location, shape, and size of the other hearing aid commentary within the hearing aid shell. This is particularly important as new technologies such as inductive power charging are beginning to be implemented in hearing aid products.
Figure 9.
Views of the (a) top and (b) bottom of the hearing aid model A with the base and end plates removed.
Figure 10.
Views of the (a) top and (b) bottom of the hearing aid model B with the base and end plates removed.
Views of the (a) top and (b) bottom of the hearing aid model A with the base and end plates removed.Views of the (a) top and (b) bottom of the hearing aid model B with the base and end plates removed.
Manufacturability and Maintenance
Figures 9 and 10 show the multiple vents as tubes separate from the hearing aid outer shell. This could be achieved through a series of tube vent paths (i.e., made of tubing) placed within the hearing aid shell. Other options include forming the vent paths within the hearing aid shell using 3D printing technique. This is not likely to be a major barrier because currently many of the hearing aid manufacturers are moving toward 3D printing. The maintenance of a multivented device presents a challenge in terms of keeping vent paths clear. Small vent paths are likely to become blocked more easily by wax buildup in the ear. Product design may have to include some form of wax filter or wax guard and simple vent path cleaning tools provided.
Conclusions
Investigations into the cross-sectional profile of the narrowest sections of the human ear canal shapes identified by the statistical shape model have shown that multiple venting is spatially feasible within the expected variances of the human ear. It is expected that the maximum total vent area available to a multiple vented hearing aid is at least as much as a traditional 2-mm vent. This has been validated in 3D space using 3D solid modeling of a multiple vented right ear hearing aid.Comparative 3D model representations of two identical hearing aid shells fitted with a single large 3-mm and 16, 0.5-mm multiple vents have demonstrated the potential spatial design flexibility provided by a multiple vented hearing aid. Venting within the multiple vented hearing aid was able to easily fit around other electronic components within the hearing aid shell. This increased the design flexibility around the shape and placement of the other elements within the hearing aid shell. It may also be used alongside other technologies to improve flexibility of design in hearing aids such as inductive power charging without providing physical limitations on their implementation. Traditional hearing aid shell manufacture combined with multiple tube vents or 3D printing have been identified as possible manufacturing methods for creating multiple vented hearing aid shells. The major challenge faced by a multiple vented hearing aid has been identified as avoiding vent blockages. As such, wax blockage prevention measures such as wax guards as well as appropriate vent cleaning tools have been recommended.
Authors: Jürgen Kiessling; Barbara Brenner; Charlotte Thunberg Jespersen; Jennifer Groth; Ole Dyrlund Jensen Journal: J Am Acad Audiol Date: 2005-04 Impact factor: 1.664
Authors: Alessandro Crimi; Martin Lillholm; Mads Nielsen; Anarta Ghosh; Marleen de Bruijne; Erik B Dam; Jon Sporring Journal: IEEE Trans Med Imaging Date: 2011-03-22 Impact factor: 10.048