An inverse scattering approach for geometric body generation: a machine learning perspective

In this paper, we are concerned with the 2D and 3D geometric shape generation by prescribing a set of characteristic values of a specific geometric body. One of the major motivations of our study is the 3D human body generation in various applications. We develop a novel method that can generate the desired body with customized characteristic values. The proposed method follows a machine-learning flavour that generates the inferred geometric body with the input characteristic parameters from a training dataset. One of the critical ingredients and novelties of our method is the borrowing of inverse scattering techniques in the theory of wave propagation to the body generation. This is done by establishing a delicate one-to-one correspondence between a geometric body and the far-field pattern of a source scattering problem governed by the Helmholtz system. It in turn enables us to establish a one-to-one correspondence between the geometric body space and the function space defined by the far-field patterns. Hence, the far-field patterns can act as the shape generators. The shape generation with prescribed characteristic parameters is achieved by first manipulating the shape generators and then reconstructing the corresponding geometric body from the obtained shape generator by a stable multiple-frequency Fourier method. Our method is easy to implement and produces more efficient and stable body generations. We provide both theoretical analysis and extensive numerical experiments for the proposed method. The study is the first attempt to introduce inverse scattering approaches in combination with machine learning to the geometric body generation and it opens up many opportunities for further developments.


Introduction
With the rapid technological advancement today, the access to realistic 3D human shapes is of great importance in both computer vision and graphics, and has various applications in different industries including virtual game design, film making, bioinformatics [36], healthcare [35], and especially, those related to garment design.Some applications involve fitting predictions, virtual try-on simulations [13,16,23,31] or size recommendations [5][6][7], that help to recommend relevant clothing which would fit specific occasions or fashion trends for online customers.Such applications require a critical ingredient on digital transformation from humans bodies to digital 3D shapes, such that the shapes maintain some of the main features from human bodies.
The traditional approaches to access reliable digital information of a human body are through laser range scanners [1], stereo reconstruction [18,22,30] or structured light methods for 3D sensing [11,20,24].However, considering the cost of data storage, network transmission and expensive scanning equipment, it is rather unpractical to scan individuals for each application.Hence many studies have been done to generate 3D human shapes based on partial input information.These prior systems can be mainly classified into three types: marker-based systems, silhouette-based systems and measurement-based systems.Markerbased system estimates dynamic 3D human body shapes by capturing a sparse set of marker positions.These techneqiues proceed by using a single static scan or multiple scans and a marker motion capture sequence of the person [4].For the static case, silhouette-based system estimates human body shapes based on a set of input images by fitting the silhouette in each view [3,8,9,15].Some apporaches also combine with machine learning that build a correlation between a training dataset of 3D body shapes and a set of 2D images, and then predict a shape based on the correlation [14].
Although marker-based systems and silhouette-based systems could yield satisfactory reconstructions on 3D human body shapes under tight dresses or naked human shapes, most of the schemes are so computationally expensive and the results are easily affected if heavy or loose clothes are worn.To overcome these difficulties, a great deal of efforts have been devoted to the investigation of simple and fast measurement-based systems [17,29,32].Typically, one considers the landmarks or circumferences from the human structures at specific locations as characteristic values.Since such characteristic values are linear or curvilinear, they are relatively invariant to articulation changes than those silhouettes measurements.If the set of characteristic values is well selected, one can achieve meaningful estimation for both global and local body shapes.Kart et al.(2011) built a system which only requires to input some personal information, such as weight, height and age as well as a 2D photograph.The decision algorithm then determines the human shape according to the measurements and the body mass index (BMI) [12].Seo et al. [32] presented a human body generation apporach by taking the anthropometric measurements, e.g.stature, crotch length, arm length, neck girth, chest/bust girth, underbust girth, waist girth and hip girth as input.They derived the relationship between the input characteristic values and the preprocessing database of 3D scanned data of human body models by using radial basis interpolation.At run-time, the system generates new human body shapes from the user input characteristic values by fitting the template model onto each scanned data.
In this paper, we develop a completely novel methodology for the geometric body generation, which fulfils the following two basic requirements: (i) the geometric body generation is automatically determined by the input characteristic sets; (ii) the predicted geometric shape fits for all input characteristic values and moreover it can well approximate the exact geometric body possessing the aforesaid characteristic values.The proposed method follows a machine-learning flavour that generates the inferred geometric body with the customized characteristic parameters from a training dataset.The training dataset consists of some preprocessed body shapes associated with appropriately sampled characteristic parameters.One of the critical ingredients and novelties of our method is the borrowing of inverse scattering techniques in the theory of wave propagation to the body generation.This is done by establishing a delicate one-to-one correspondence between a geometric body and the far-field pattern of a source scattering problem governed by the Helmholtz system.It in turn enables us to establish a one-to-one correspondence between the geometric body space and the function space defined by the far-field patterns.Hence, the far-field patterns can act as the shape generators.The shape generation with prescribed characteristic parameters is achieved by first manipulating the shape generators in the function space and then reconstructing the corresponding geometric body from the obtained shape generator by a stable multiple-frequency Fourier method.The proposed method is in sharp difference from the existing methodologies in the literature, which usually treat the human body as a suitable Riemannian manifold and the generation is based on non-Euclidean approximation and interpolation.In fact, in all of the literature mentioned earlier on manifold learning of body generation, one typically uses Principal Component Analysis (PCA) or Principal Geodesic Analysis (PGA).PCA and PGA are used for optimal reduction of the data and thus efficient deformation by computing statistics on Euclidean manifolds or non-Euclidean manifolds can be achieved; see [10,32] and the references therein for more relevant discussion.In our new approach, the shape generator enables us to train the learning dataset via the algebraic operations in the shape space directly without dealing with the deformation of the manifold meshes between geometric shapes.
The rest of the paper is organized as follows.In section 2, we provide rigorous mathematical formulations of characteristic values and shape space.Section 3 introduces the notion of shape generator via the inverse source scattering associated with the Helmholtz system.In Section 4, we present the mathematical setup of the geometric body generation from a machine-learning perspective.Section 5 is devoted to the development of the new method for the shape generation.In Section 6, we present several two-and three-dimensional numerical examples to show the effectiveness and efficiency of our method.The paper is concluded in Section 7 with some relevant discussion.

Preliminary knowledge on shape manifold theory
In this section, we present some preliminary knowledge on the shape manifold theory that shall be needed in our subsequent study of body generation.Generally speaking, a geometric shape or a geometric body is a topological n-manifold, n ∈ N, equipped with certain shape descriptors, which give the full information to describe the geometric shape.We call such shape descriptors as characteristic values.We have the following formal definition.Definition 2.1.Let D be a topological n-manifold with n ∈ N. Let Λ D := {λ (j) } j∈C be a set of parameters associated with D that are invariant with respect to isometric deformations and are independent to the parametrizations of D. Here, the cardinality C might be finite or infinite.Λ D is said to be a characteristic set of D if it uniquely determines D. D and its characteristic set Λ D , written as the (D, Λ D ) is referred to as a geometric shape or a geometric body.
Clearly, Definition 2.1 includes much general geometric objects.However, for the present study, we are mainly concerned with the case that D can be embedded into R d , d = 2, 3, as a bounded domain.That means, we exclude some interesting cases such as D is a Riemannian surface with boundary in R 3 .Nevertheless, our study is general enough to include the human body as a specific case.
In Definition 2.1, the set of characteristic values is typically a set of measurements which gives a systematic characterization of the size, shape and composition of a geometric object for us to determine the shape of the object.For example, when considering a rectangular object, once can introduce a set of characteristic values containing its height, width and length, which provide all details to determine a unique rectangular shape.Expanding the same idea to human body shapes, one could also use characteristic sets to represent them.There are many different ways to represent a human shape.We would try to group those characteristic values into four main catagories, including Eucidean distance, geodesic distance, circumference and ratio.The Eucidean or geodisic distance is linear or curvilinear distance between two points on the human model, such as stature, crotch length, arm length, shoulder breadth etc..The circumference can be computed by the horizontal girth of the body, such as neck girth, chest/bust girth, under-bust girth, waist girth, hip girth, etc..The ratio can be information of weight, Body Mass Index, muscle and fat rate.In spite of the above characteristic values, one can also consider some pure measurements such as age or gender as characteristic values.
The full set of characteristic values gives the complete information of a geometric shape without lossing any information.It is easy to imagine that the cardinality of a set of characteristic values depends on the complexity of a shape.Hence, the number of characteristic values required can be considered as the dimensionality of the geometric shape.For those complicated objects, like human shapes, it may require infinite set of characteristic values for accurate formulations.Due to practical reasons, one can consider a reasonable truncation of an infinite characteristic set into a finite one for a complicated geometric shape.In doing so, we can consider our study in the following product space where D is composed of all the bounded domains in R d and V is an M -dimensional vector space containing the characteristic values.In fact, in the present study, the characteristic values are usually real numbers and one can take V = R M with M ∈ N. S is referred to as the geometric shape space.According to the (approximate) one-to-one correspondence between a geometric shape and its characteristic values in Definition 2.1, we readily see that all the shape information can be obtained by a single point of this M -dimensional vector space V .By adjusting the characteristic values, we can obtain new geometric shapes and this is a key ingredient in our human body interpolation.

Shape generators via inverse source scattering
In the previous section, we introduce the important notion of shape space for our study.We proceed to introduce another critical ingredient, shape generator, for our subsequent study of the geometric body generation.In fact, the generation of a new geometric shape shall be based on algebraic interpolation of exemplar models from the shape space.If the algebraic operations are to be conducted directly in the shape space, dealing with geometric deformations of manifolds, one would certainly encounter very complicated and tedious calculations and manipulations because of the lack of global parametrizations for the non-Eucidean shapes involved.The shape generator can overcome this challenge by bridging the geometric shape space and the function space.To that end, we next introduce the inverse scattering problem in finding an active source from its generated far-field pattern.
Let f : R d → C be a function having a compact support, f = χ D ϕ, where D ⊂ R d is a bounded domain and ϕ ∈ L ∞ (R d ).The set D is the external shape of f while ϕ describes the intensity of the source at various points in D. We assume that ϕ and D do not depend on the wavenumber k ∈ R + .In other words we are considering monochomatic scattering.The source f produces a scattered wave u ∈ H 2 loc (R d ) given by the unique solution to where r = |x| for x ∈ R d .The limit in (3.1) is known as the Sommerfeld radiation condition which characterizes the outgoing nature of the radiating wave.By the limiting absorption principle (cf.[21]), the solution to (3.1) can be computed as follows, where signifies the Fourier transform of f .Inverting the Fourier transform in (3.2), one has the following integral representation, where H (1) is the first-kind Hankel function of order (d−2)/2.Stationary phase applied to (3.4) yields that where The far-field pattern of u is given by then it is known on the whole set by analytic continuation.The inverse source scattering problem is concerned with the recovery of According to our discussion above, without loss of generality, we always assume that Σ = S d−1 × R + in what follows.The inverse source problem arises in a variety of important applications including detection of hazardous chemicals, medical imaging, photoacoustic and thermoacoustic tomography, brain imaging, artificial intelligence in gesture computing and others.We refer to two recent articles [25,26] by two of the authors of this article for some recent developments on the inverse source problem.
Next, let us consider a specific case by assuming a source supported in a domain D with a constant density 1.Then clearly by (3.6), there is a one-to-one correspondence between D and u ∞ (Σ; (3.7) Based on (3.7), we next introduce Definition 3.1.For a geometric shape (D, Λ D ) ∈ S, defined via the Helmholtz system (3.1) is called a shape generator for D.
Remark 3.1.By Definition 3.1, a geometric body D can be completely determined by a shape generator u ∞ (Σ).Since u ∞ (Σ) is from a function space, this paves the way for the new body generation through function interpolations.
Remark 3.2.By (3.6), we know the far-field pattern is actually the Fourier transform of the source density up a dimensional constant.However, introducing the shape generator via the inverse scattering approach shall provide more physical insights in our study, and moreover it enables us to borrow ideas from the inverse scattering literature of recovering the geometric shape D from the associated far-field pattern.This also paves the way of extending the idea by using other inverse scattering models that have such one-to-one correspondence between geometric shapes and far-field patterns; see more relevant discussion in Section 7.

Mathematical setup for the geometric body generation
In this section, we introduce the mathematical formulation of the geometric body generation for our study from a machine learning perspective.For a geometric shape (D, Λ D ) with the associated shape generator u ∞ (Σ; D), the pair of the high dimensional variables, written as {(Λ D , u ∞ (Σ; D))}, is referred to as an input-output pair.Let {(Λ D i , u ∞ (Σ; D i ))} i∈N with N = {0, 1, . . ., N pair } be a set of input-output pairs associated with the characteristic sets {Λ D i } i∈N .Here the input characteristic sets Λ D i are introduced as with C = {1, . . ., M }.In (4.1), the notation λ represents the characteristic value of the i-th geometric shape D i in the j-th direction of its characteristic set.Here, the cardinality N is finite.The training dataset of the geometric body generation is introduced to be The training dataset consists of certain pre-sampled geometric shapes with statistically well selected characteristic values.The corresponding shape generator of a specific body in the training dataset can also be pre-calculated and stored.The main goal of our study is to first infer a learning model from the training dataset, T Z : R M → L 2 (S d−1 × R + ) that fulfils the following requirements: (1) It fits the training data well in the sense that (2) It can be used to infer the shape generator for a given new shape with prescribed characteristic values, namely, and with a statically well selected training dataset, it is justifiable to expect that û∞ (Σ; where u ∞ (Σ; D new ) is the shape generator for D new .If a learning model can be achieved that fulfils the two requirements as described above, then the body generation can be proceeded as follows.For a given new set of characteristic values, one first generates the learned shape generator as in (4.5).By a certain inverse scattering approach, one can then reconstruct the (approximate) shape D new from the corresponding shape generator û∞ (Σ; D new ).In the next section, we shall develop the two critical ingredients in the body generation procedure described above, namely, the learning model and the reconstruction method.To be more definite and specific, we first introduce the following definition from a machine learning perspective.
Let H be a compact subset of L 2 (S d−1 × R + ).T Z ∈ H (with specified coefficients C) is said to be the best fit learning model associated with the training dataset Z if it is the minimizer of the following optimization problem, min According to Definition 4.1, the choice of the learning subspace H plays a critical role.However, we note that the shape generator is actually (real) analytic in all of its arguments.Hence, instead of solving the computationally costly optimization problem (4.8), we can make use of the functional interpolation to produce a well-rounded shape learning model.This is one of the main advantages of introducing the shape generator through the inverse scattering model.In the next section, for a given training dataset as in (4.7), we shall derive a learning model using the cubic B-spline interpolation through the use of the high-dimensional data-points (4.4).For the reconstruction of the approximate body shape from the shape generator obtained through the learning model, we shall make use of a multiple-frequency Fourier method, and it can also produce an efficient and stable recovery.Throughout, we assume that the characteristic values in the training dataset is statistically well selected and it is not the focus of the present article.

A scheme for geometric body generation
In this section, we develop the details of our scheme for the geometric body generation following the general discussion made in the previous section.We first derive the learning model through the functional interpolation of the high-dimensional data in the training dataset.To that end, we present some preliminary knowledge on the cubic B-spline, and we also refer to [2,19,33,34] for more relevant discussion on the cubic B-spline.
define on the intervals [a j , b j ] as M sets of points λ (j) Here, k j is the greatest number of distinct characteristic values in the j-th direction of the characteristic set.We remark that if the characteristic values of the training dataset are all collected in distinct values, then k j is actually the last index of the training dataset, N pair .However, the training dataset might be collected in such a way that some body shapes may possess the same characteristic value in the j-th direction, and hence k j is usually smaller than N pair .
With the above notation, the training dataset stored as the array in (4.2) can be represented as elements on the grid mesh corresponding to the characteristic numbers λ (j) as described above.The interpolation data are the corresponding shape generators and are written as ) where g j = 0, . . ., k j , j = 1, . . ., M. In the subsequent study, we shall stick to the same notation g 1 , g 2 , . . ., g M to represent the linear indexing.The following example demonstrates a real application for the human body generation.
Example 5.1.The training dataset consists of 20 bodies with two characteristic values as consideration, say, height and relative weight.Here, the height and relative weight are the two directions of the grids, λ (height) and λ (weight) .Suppose the heights of the sampled bodies are given by 1.5m, 1.6m, 1.7m, 1.8m, 1.9m and the relative weights of the sampled bodies are given by 60%, 80%, 100%, 120%.Then the first grid ∆ height = {λ Table 1.Training dataset with the given characteristic grids (λ (height) , λ (weight) ).
Definition 5.1.The sets of k j + k, j = 1, . . ., M , B-spline basis functions {B l,k (λ (j) )} k j +k l=1 of degree k of the function space S k (∆ j ) are defined based on concurrent boundary knots vectors with Cox-deBoor recurrence [19], for l = 1, . . ., k j + k, where k = 3 is for the cubic B-spline and λ (j) l are elements of the knot vectors, satisfying the relation λ All methods are in the following using splines with a k + 1 regular knot vector, and the interior knots are the grid points.Based on Definition 5.1, we next introduce a general learning model for the geometric shape generation through the multidimensional cubic Bspline interpolation. (1), . . ., λ (M ) ) for the geometric body generation associated with the sets of the grids {∆ j } j∈C is defined as follows

Body Learning Model I. Given the training dataset
which satisfies the following conditions by (4.4) where c g 1 ,g 2 ,...,g M with g j = 0, . . ., k j , j = 1, . . ., M are the coefficients to be determined from the training dataset Z, B j g j ,3 (λ (j) ), j ∈ C are the B-spline basis functions of degree 3 defined in (5.5), and k j , j = 1, . . ., M is the number of different characteristic values in each direction.
Remark 5.1.In Learning Model I, (5.7) presents a general form of the learning model T Z for the geometric body generation associated with the non-uniform grids {∆ j } j∈C .The learning model eventually generates a B-spline interpolation with the associated spacing for each segment.If the training dataset consists of equidistant grids, we can derive a faster and easier learning model and this shall be provided in the next subsection.

Uniform B-spline.
In this subsection, we derive a learning model for a special case with the training dataset consisting of equidistant grids.By (5.1), for the M set of equidistant grids ∆ j = {λ k j } j∈C with additional conditions λ (j) g j = a j + g j h j , h j = b j − a j k j , , g j = 1, . . ., k j , (5.9) the B-spline basis function β k (t) of degree k is a symmetrical, bell-shaped function constructed from k + 1 times self-convolution of the β 0 (t) basis function of degree zero which is a centered rectangle around origin [34] (5.11) The centered symmetric B-spline of degree k has an explicit expression [33] where the function x + is defined as follows x, for x > 0, 0, otherwise.(5.13)In this paper, we are particular intereted in the cubic B-spline.By (5.12), the closed-form representation of the cubic B-spline basis function can be also expressed as elsewhere, (5.14) which is used for preforming the interpolation.Then we choose the interpolation kernels to be L g j (λ (j) ) = β 3 λ (j) , k j = 1, . . ., N pair + 3, (5.15) as the basis of S 3 (∆ j ) in the λ (j) -direction such that L = {L 1 , L 2 , . . ., L k j +3 } is a basis of the (k j + 3)-dimensional space S 3 (∆) and hence, the basis of the j∈C (k j + 3)-dimensional space S 3 (∆ 1 , . . ., ∆ M ) in the directions of λ (1) , . . ., λ (M ) is given by (5.16) Based on (5.16),(5.1)and (5.9), we next introduce the learning model for the uniform case.
Body Learning Model II.Given the training dataset 1) , . . ., λ (M ) ) for the geometric body generation associated with the sets of equidistent grids {∆ j } j∈C defined in (5.9) is defined as follows which is required to satisfy the following conditions where c g 1 ,g 2 ,...,g M with g j = 0, . . ., k j , j = 1, . . ., M. are the coefficients to be determined from the training dataset Z, L j g j (λ (j) ) are B-spline basis functions of degree 3 defined in (5.15) and k j , j = 1, . . ., M is the number of different characteristic values in each direction.5.3.Natural Spline.In Learning Models I and II, the learning functionals for the nonuniform and uniform case are respectively considered in the j∈C (k j + 3)-dimensional space S 3 (∆ 1 , . . ., ∆ M ).
j∈C (k j + 3) interpolation conditions are required to determine the coefficients c i 1 ,i 2 ,...,i M in the training models.However, there are only j∈C (k j + 1) shape generators to specify j∈C (k j + 1) conditions in (5.8) or (5.18).To obtain a unique correlation between the characteristic values and the shape generator, we need to add j∈C 2 = 2 M conditions, which define the second-order derivatives of the spline function at the boundary a j and b j to be equal to 0 and lead to a natural spline.
where the coefficients c g 1 ,g 2 ,...,g M in (5.19) could be determined by solving the natural spline problem.
5.5.Reconstruction.In this subsection, we briefly outline the Fourier method for the reconstruction of geometry shape D new by using the shape generators u ∞ (Σ; D new ).Define the periodic Sobolev space by where σ ≥ 1, ξ ∈ Z d and g(ξ) denote the Fourier coefficients of g.Suppose that f ∈ H σ (R d ) has a compact support in domain V 0 = (−a/2, a/2) d , (a > 0), then the Fourier transform of f is represented by where the overbar stands for the complex conjugate and the Fourier basis functions are given by Definition 5.2 (Admissible wavenumbers and observation directions).Let µ be a sufficiently small positive constant such that 0 < µ < 1 and (µ, 0, 0), d = 3, then the admissible set of wavenumbers is defined by correspondingly, the admissiable set of observation directions is given by where Due to supp S ⊂⊂ V 0 , for ξ ∈ Z 3 \{0}, the far-field pattern defined in (3.6) can be written as where k ∈ K and x ∈ X depend on ξ.Combining (5.20) and (5.21), one has For ξ = 0, using the Fourier expansion of f , we derive that which implies Therefore, the Fourier method is to approximate f by a truncated Fourier expansion where N ∈ N + denotes the truncation order and the Fourier coefficients are given by (5.22) and (5.23).Hence the domain D new is determined since the set D new is the external shape of f .Next, we investigate the stability of the proposed Fourier method.In practical computation, there exists some noise between the shape generators u ∞ (Σ; D new ) and the predictions û∞ (Σ; D new ), which satisfies where δ > 0 denotes the noise level.Noting that (x, k) ∈ X×K ⊂ Σ, then the approximation of f from predicted shape generators is given by where Theorem 5.1.Let f be a compactly supported function in H σ (R d ), σ ≥ 1, with suppf ⊂⊂ V 0 , then we have the following estimate where d ≤ τ ∈ R + and C is a constant which depends on f, a, d, µ.
Let N = τ δ − 2 2+d , τ ≥ d , here [X] denotes the largest integer that is smaller than X + 1.From definition 5.2, the truncated wavenumbers and observation directions can be written as Thus, the truncated Fourier expansion of f from the predictions {û ∞ (Σ; D newq )} q∈Q takes the form where (5.31) 5.6.Summary.Motivated by above discussion, we are ready to present our novel modeling methodology for geometric shape in R d , d = 2, 3, see Algorithm 1.
Algorithm 1 Inverse-scattering-based geometric body generation scheme 1: Select the parameter N , the set of admissible wavenumbers K N and the set of admissible observation directions X N .
the coefficients c g 1 ,...,g M of the learning model T Z by solving the problem of the natural spline interpolation.3: Given the characteristic sets Λ Dnew , predict the new shape generators {û ∞ (x, k)} for x ∈ X N , k ∈ K N with the use of the learning model T Z .4: Compute the Fourier coefficients f δ 0 and f δ ξ defined in (5.31) for 1 ≤ |ξ| ∞ ≤ N t .5: Select a sampling mesh T h in a region V 0 .For each sampling point z j ∈ T , calculate the imaging function f δ N defined in (5.30).D new is obtained as the external shape of f δ N .

Numerical examples
In this section, several numerical examples are conducted to show that the proposed method is effective and efficient.
The proposed algorithm is implemented by using Matlab 2016.The shape generator {u ∞ (x, k; D i )} i∈N is obtained by solving the direct problem of (3.1).To avoid the inverse crime, we use the quadratic finite elements on a truncated spherical domain enclosed by a PML layer.The mesh of the forward solver is successively refined till the relative error of the successive measured scattered data is below 0.1%.Then artificial shape generators are generated by applying the Kirchhoff integral formula to the scattered data.Thus the training dataset is given by Next, we present the implementation of interpolation.The characteristic value set is given by } has M variables.For a fixed wavenumber k j , we use cubic spline interpolation to obtain the coefficients c i 1 ...,i M from the characteristic value Λ D i and the shape generator u ∞ (x j , k j ; D i ).Therefore, given a new characteristic value Λ Dnew = {λ Dnew }, we obtain the predicted shape generators {û Finally, we specify details of recovering the geometry.As discussed above, reconstructing the geometry shape is equally to reconstructing the source function f .In the discrete formula, the domain V 0 is divided into a uniform mesh with size 100×100 in two dimensions and size 100 × 100 × 100 in three dimensions.Further, the approximated Fourier series f δ N are computed at the mesh nodes T j , j = 1, • • • , 100 3 in (5.30).Thus, the geometry shape D new is approximated by the boundary of the imaging results f δ N .6.1.Kite shaped Experiments.In the first example, we aim to reconstruct a kite shaped domain with scale changing.The kite shaped domain is parameterized by    3, where the black dotted lines denote the shadows of the exact cube boundary.Due to discontinuities of the source, there is Gibbs phenomena on the boundary of the rectangular solids.On the whole, given the characteristic values, our proposed method is valid for determining the geometry shape.This experiments consider two characteristic values, i.e., height and relative weight.Define the exact weight by EW and standard weight by SW , then the relative weight RW is calculated by Here the height of the body is given by 1.5m, 1.6m, 1.7m, 1.8m, 1.9m and the relative weight of the body is given by 60%, 80%, 100%, 120%, 140%.Some human body shapes in the training dataset are presented in Figure 4.In addition, we choose two characteristic values of human body which are not covered by the training data.The first body's height is 1.55m and the relative weight is 130%.The second body's height is 1.85m and the relative weight is 110%.Figure 5(a) and Figure 6(a) present the exact body shape with the given characteristic value.

Concluding remarks
In this paper, we develop a machine-learning method in generating a geometric body shape through prescribing a set of characteristic values of the body.The generation is mainly based on a given training dataset consisting of certain pre-selected body shapes with statistically well-sampled characteristic values.A major novelty and critical ingredient of our study is the borrowing of inverse scattering techniques in the theory of wave propagation to the geometric shape generation.We introduce the notion of shape generator which establishes a one-to-one correspondence between the geometric shape space and the function space consisting of the multiple-frequency far-field patterns associated with the time-harmonic source scattering problem.The shape generator plays an intermediate role in the geometric shape generation.First, the training dataset of geometric shapes is converted into a subset of the function space consisting of the corresponding shape generators.Then a learning model is derived through a functional interpolation of the aforementioned shape generators.For a given set of characteristic values, one then uses the learning model to obtain the shape generator of the underlying geometric body and finally reconstructs it through a multiple-frequency Fourier method.the left column is 60%; the center column is 100%, the right column is 140%; height: the top row is 1.50m, the center row is 1.70m, the bottom row is 1.90m.
To our best knowledge, the present study is the first attempt to introduce inverse scattering approaches in combination with machine learning to the geometric body generation and it opens up many opportunities for further developments.For example, in the current article, the shape generator is introduced through an inverse source scattering model where we make use of the one-to-one correspondence between a geometric shape and the multiple-frequency far-field pattern associated with a compactly-supported acoustic source.One may consider to introduce the shape generator through other inverse scattering models, e.g. the inverse acoustic obstacle scattering model (cf.[28]) or the inverse electromagnetic scattering model (cf.[27]).In doing so, one may achieve other geometric shape generation schemes that are suitable for different applications.

Definition 4 . 1 .
(Body Learning Model) Given a training dataset

5. 1 .
Preliminary knowledge on the cubic B-spline.Consider the training dataset (4.2).Let the M sets of unique grids in the directions of {λ (j) } j∈C

5. 4 .
Prediction on shape generator.With the Learning Models I and II established in the previous subsections, for an input new set of characteristic values Λ Dnew associated with a new geometric body D new , the unknown shape generator can be generated as follows,

6. 4 .
Body Experiments.In last example, we consider a challenging case and verify the proposed method by using a set of synthetic experiments on 3D human body shape.The training dataset consists of 25 bodies which are generated by the MakeHuman 1.1.1soft.

Figure 5 (
b) andFigure 6(b)  show the prediction of the human body shape with the given characteristic value.The results show that our method is efficient to predict the human body shape.

Figure 4 .
Figure 4. Isosurface plots of some training body data.Relative weight:the left column is 60%; the center column is 100%, the right column is 140%; height: the top row is 1.50m, the center row is 1.70m, the bottom row is 1.90m.
Rectangular Solid Experiments.In the third example, we verify the proposed method by using a set of artificial experiments on rectangular solid.The training dataset consists of 125 rectangular solids with different height, width and length.Here the height, width and length are uniformly distributed on [1, 2] with 5 amounts, i.e., 1, 1.25, 1.5, 1.75, 2. Here, We consider four sets of different height, width and length of rectangular solids which are not covered by the training data.The imaging results with different characteristic values are shown in Figure ). multi-domain with different characteristic values via contour plots, where the black dotted lines denote the exact boundary.It demonstrates very good imaging performance of the approach.6.3.