1 Introduction
Reconstruction problems of convex bodies (compact convex sets with nonempty interiors in $\mathbb {R}^n$ ) for some known geometric functions and measures (e.g., area measures, curvatures, brightness function, and lightness functions in this paper) have caught great attention in convex geometric analysis. For instance, a convex body K can be uniquely determined, up to translations, by its surface area measure $S(K,\cdot )$ , which is a spherical Borel measure defined by
where $\mathcal {H}^{n-1}$ denotes the $(n-1)$ -dimensional Hausdorff measure, and $\nu _K: \partial K \to \mathbb {S}^{n-1}$ is the Gauss map defined $\mathcal {H}^{n-1}$ -almost everywhere on the boundary $\partial K.$ Investigation about whether a spherical Borel measure $\mu $ is the surface area measure of a convex body K is the classical Minkowski problem. Its complete solution goes back to the work of Aleksandrov and Fenchel-Jessen (see [Reference Alexandrov1, Reference Huang, Lutwak, Yang and Zhang19, Reference Schneider28] for references). Around 1937, Aleksandrov also showed that a centrally symmetric convex body K is uniquely determined, up to translations, by its brightness function
Independently, Campi [Reference Campi8] and Bourgain and Lindenstrauss [Reference Bourgain and Lindenstrauss6] proved a stability version of this result.
1.1 Reconstruction problems for a geometric function (or measure)
In view of much of previous work, we roughly summarize the reconstruction problems into the following three types:
-
(I) Study whether a geometric function (or measure) can uniquely and stably determine a convex body.
-
(II) Given a function (or measure) $\Xi $ , study both the necessary and sufficient conditions of $\Xi $ so that there is a convex body K whose geometric function (or measure) is the pre-given $\Xi $ .
-
(III) Suppose the geometric function (or measure) of a convex body K is given. Try to build a computational procedure to reconstruct the related convex body; prove the convergence; and design an algorithm based on the computational procedure, and then implement it and test it.
Aleksandrov’s projection theorem and its stability versions belong to the first topic (I). The spherical harmonic is the powerful tool therein. This technique can also solve problem (II) for brightness function (see, e.g., Lemma 4.4). Gardner and Milanfar [Reference Gardner and Milanfar13, Reference Gardner and Milanfar14] first studied the type (III) problem for brightness function. In fact, they provided algorithms to reconstruct the origin-symmetric convex body from its given brightness function. Before that, they gave an algorithm to reconstruct the polytope from the pre-given discrete measure. Gardner, Kiderlen, Milanfar, and Peyman [Reference Gardner, Kiderlen and Milanfar12] studied the convergence of related algorithms, which belongs to the topic (III). Bianchi, Gardner, and Kiderlen [Reference Bianchi, Gardner and Kiderlen2] studied problem (III) for covariograms of convex bodies. For studies about the reconstruction problems for other geometric functions (or measures), see, e.g., [Reference Kaasalainen and Lamberg20, Reference Kaasalainen and Lamberg21, Reference Little23, Reference Poonawala, Milanfar and Gardner26].
The Minkowski problems for geometric measures always involve type (I) and type (II). Except for the surface area measure, Aleksandrov, Fenchel, and Jessen also introduced the mixed area measures, and related problems are called the Christoffel–Minkowski problem. For these geometric measures, type (I) problems are completely solved, but type (II) problems are not (see [Reference Guan and Ma17] for a regular case). While the curvature measures were introduced by Federer about half a century ago, its dual counterpart was introduced quite recently by Huang, Lutwak, Yang, and Zhang [Reference Huang, Lutwak, Yang and Zhang19, Reference Lutwak, Yang and Zhang24]. The studies of the reconstruction problems of type (II) for the new geometric measures have caught great attention, in both convex geometry and partial differential equations (see, e.g., [Reference Böröczky, Henk and Pollehn4, Reference Böröczky, Lutwak, Yang, Zhang and Zhao3, Reference Chen, Huang and Zhao10, Reference Li, Sheng and Wang22, Reference Zhao30]). However, the problem of type (I) for the dual curvature measure is only solved in one special case on the plane [Reference Böröczky, Lutwak, Yang and Zhang5]. It seems that the problem of type (III) is only considered for surface area measure. For instance, see [Reference Gardner and Milanfar14, Reference Little23] for discrete case in $\mathbb {R}^n$ , and see [Reference Kaasalainen and Lamberg20] for smooth case in $\mathbb {R}^3$ . In [Reference Gardner and Milanfar14], the authors reconstructed polytopes from pre-given discrete measures, and they proved the convergence of their algorithm; in [Reference Kaasalainen and Lamberg20], by using the spherical harmonic, the authors gave an effective numerical algorithm for smooth curvature functions in $\mathbb {R}^3$ , but they did not rigorously prove the convergence.
Given a convex body K in $\mathbb {R}^n$ , a lightness function [Reference Campi, Gardner, Gronchi and Kiderlen9] is a binary function that gives the total reflected light resulting from illumination by a light source at infinity in the direction w that is visible when looking in the direction v. An interesting model appeared in [Reference Campi, Gardner, Gronchi and Kiderlen9] is the Lambertian lightness function, which means that it follows Lambert’s cosine law, and it is formally defined by
where $(t)_+$ and $(t)_-$ denote, respectively, the positive part and negative part of a real number t. In view of the general description of lightness function, the authors [Reference Campi, Gardner, Gronchi and Kiderlen9] also introduced a more general lightness function
where $f_w:\mathbb {S}^{n-1} \to [0,\infty )$ is positive and continuous in the open half-space $w^+ = \{x: x\cdot w>0\}$ and vanishes in its complement. This assumption is to fulfill that a watcher from the viewing direction w can only see a visible part of $\partial K$ . The authors [Reference Campi, Gardner, Gronchi and Kiderlen9] gave exhaustive studies about the first topic (I), for these lightness functions.
The aim of this is to study (the theoretical part of) problem (III) for surface area measure, and problems (II) and (III) for lightness functions. We explicitly restate them as follows:
-
• (III) Given the surface area measure $S(K,\cdot )$ , try to build a computational procedure to reconstruct the convex body, and prove the convergence.
-
• (II) What are the sufficient and necessary conditions of a binary function $G:\mathbb {S}^{n-1}\times \mathbb {S}^{n-1}\to \mathbb {R}$ so that there is a convex body K whose lightness function $Q_K$ is the pre-given G?
-
• (III) Given a lightness function $Q_K$ , try to build a computational procedure to reconstruct the convex body, and prove the convergence.
The spherical harmonics will be an important tool in this paper. They are the eigenvectors of spherical Laplace–Beltrami operator $\Delta _S:C^\infty (\mathbb {S}^{n-1})\to C^\infty (\mathbb {S}^{n-1})$ . For more detailed introduction, see Section 2.5. We also refer to [Reference Bourgain and Lindenstrauss6, Reference Groemer15, Reference Rubin27–Reference Seeley29] as good general references.
In Section 3, we give a computational procedure to reconstruct the convex body from a pre-given surface area measure $\mu $ on $\mathbb {S}^{n-1}$ . The main computations in the procedures will be the integrals of spherical harmonics and the mixed discriminants (which will be introduced in Section 2.3) of their Hessians. First, we consider a sequence of convex optimization problems ( $P_N$ ), and we prove in Theorem 3.5 that the solutions of ( $P_N$ ) will converge to a dilation of the convex body whose surface area measure is the pre-given $\mu $ . Noticing that ( $P_N$ ) has infinite many constraints, it is very difficult to solve ( $P_N$ ) directly. Therefore, we construct the subproblem ( $P_{N,k}$ ), which is an approximation problem of ( $P_N$ ). Nontrivially, we prove in Theorem 3.3 that ( $P_{N,k}$ ) are also convex optimization problems when k is sufficiently large, but with only finitely many constraints. Since problem ( $P_{N,k}$ ) are convex optimization problems with only finite many constraints, they can be efficiently solved by many existing optimization package, and we do not build the explicit algorithm. Finally, by combining Theorem 3.3 with Theorem 3.5, we conclude that the solution of ( $P_{N,k}$ ) (after normalization) converges (in Hausdorff metric) to the convex body K whose surface area measure is $\mu $ .
In Section 4.2, using the language of spherical harmonics, we give both the necessary and sufficient conditions of a binary function to be a lightness function, under the Sobolev setting. We use $H_k(\mathbb {S}^{n-1})$ to denote the Sobolev space of order k, and the detailed introduction can be seen in Section 2.4. We take the Lambertian lightness function $L_K(v,w)$ as an example, since it has more physical explanations. A trivial attempt is that if K is smooth and symmetric,
which is the $L^2$ cosine transform of $S(K,\cdot )$ . However, it is well known that the $L^p$ cosine transform is not surjective from $C_e^\infty (\mathbb {S}^{n-1})$ (even functions in $C^\infty (\mathbb {S}^{n-1})$ ) into $C_e^\infty (\mathbb {S}^{n-1})$ when p is even.
In view of this, the interesting but also reasonable thing in Section 4 is that the convex body K can be determined by the even expansion of the binary function $L_K$ without symmetric assumption, which is different from the brightness function. Moreover, the representations in the main theorems naturally lead to an algorithm to compute the curvature function of a convex body K. Our proof relies on some results that we slightly develop in Section 4.1, including the eigenvalues and the reverse of the partial cosine transform
in the Sobolev setting. We note that the Sobolev setting and the tool of spherical harmonics are very suitable and convenient to build algorithms to reconstruct convex bodies from known surface area measures or lightness functions.
In Section 5, we use the main results in Section 4.2 to construct a computational procedure to compute the surface area measure from a pre-given lightness function. Then, using the Minkowski procedure (Mink-Pro) in Section 3, we are able to reconstruct the convex body K (without the symmetric assumption). The convergence of the procedure is also shown.
The main computations in the two procedures of Sections 3 and 5 are merely about the spherical harmonics $Y_{m,j}$ ’s, including the mixed discriminants of their Hessians, and the scalar products of them with given function (or measure). These computations are independent of the indices $N,k$ of the problem ( $P_{N,k}$ ), and hence can be computed in advance and stored. In other words, the two procedures are robust and easily implementable.
2 Preliminaries
In this section, we collect some definitions, notations, and basic facts about convex bodies, the Minkowski problem, Sobolev space on $\mathbb {S}^{n-1}$ , and the spherical harmonics. We refer to [Reference Gardner11, Reference Gruber16, Reference Hebey18, Reference Schneider28, Reference Seeley29] as good general references.
2.1 Convex bodies
Let K be a nonempty convex compact set. Its support function $h_K:\mathbb {R}^n \to \mathbb {R}$ is defined by
Sometimes, we also view it as a spherical function by restricting it to $\mathbb {S}^{n-1}$ . The support function is always sublinear. Conversely, a sublinear function uniquely determines a nonempty convex compact set.
The Hausdorff distance between two nonempty compact convex sets $K_1,K_2$ is defined by
Given a function $f\in C(\mathbb {S}^{n-1})$ , denote $W[f]$ to be the Wullf shape determined by f:
If $W[f]$ has interiors, and $f_i\in C(\mathbb {S}^{n-1}), i=1,\ldots ,$ is a sequence of continuous functions that converges to $ f$ uniformly, then $W[f_i]\to W[f]$ in the Hausdorff metric.
We say the surface area measure $S(K,\cdot )$ has a density, if it is absolutely continuous with respect to the spherical Lebesgue measure. In this case, we denote the Radon–Nikodym derivative by
and we call $F_K$ the curvature function.
Given a function $h\in C^2(\mathbb {S}^{n-1}),$ denote $D^2h$ by
where $\nabla ^2 h$ is the Hessian of h with respect to the standard metric on $\mathbb {S}^{n-1}$ , and I denotes the identity map. If one extend the h to a positive-homogeneous function on $\mathbb {R}^n$ , then $D^2h(u)$ is exactly the restriction of its Hessian in $\mathbb {R}^n$ to the tangent space of $\mathbb {S}^{n-1}$ at u.
In fact, if $\partial K$ is $C^2$ and has positive Gauss curvature everywhere, then $S(K,\cdot )$ has a density,
which is the reciprocal Gauss curvature at $\nu _K^{-1}(u)\in \partial K.$
2.2 Minkowski problem and the Monge–Ampère equation
Solution to the classical Minkowski problem: Suppose $\mu$ is a finite Borel measure on satisfying:
(1) $\mu $ is not concentrated in a half-sphere;
(2) the centroid of $\mu $ is the origin:
Then, there is a unique (up to translations) convex body K in $\mathbb R^n$ so that $\mu (\cdot ) =S(K, \cdot )$ . Moreover, the convex body K is homothetic to the solution of the maximization problem
If the measure $\mu $ has a density function $g :\mathbb S^{n-1}\rightarrow \mathbb R$ with respect to the Lebesgue measure of the unit sphere $\mathbb {S}^{n-1}$ , the Minkowski problem is equivalent to the study of solution to the following Monge–Ampère equation on $\mathbb S^{n-1}$ :
Let $f\in L^1(\mathbb {S}^{n-1})$ be nonnegative, satisfying $\int _{\mathbb {S}^{n-1}} f>0$ and
Then, by the solution of the Minkowski problem, there is always a convex body K such that
2.3 Mixed discriminant and mixed volume
Let $A_1,\ldots ,A_m$ be self-adjoint linear transforms from a k-dimensional Hilbert V to itself. Then, the following is a k-homogeneous polynomial of $t_1,\ldots ,t_m \in \mathbb {R}$ :
Its coefficients $\mathrm {D} (A_{i_1},\ldots ,A_{i_{k}})$ , which are symmetric functions of $A_{i_1},\ldots ,A_{i_{k}}$ , are called the mixed discriminants of the linear transforms.
If we choose an orthonormal basis in V, these $A_i$ can be represented by symmetric $k\times k$ matrices, and we still denote them by $A_i$ . Denote the matrices $A_{r}=\left (a_{i j}^{r}\right )_{i, j=1}^{k}, r=1, \ldots , k$ , and the mixed discriminant $\mathrm {D} (A_1,\ldots ,A_k)$ can be computed by
where $P_k$ is the group of all permutations of the set $\{1,\ldots ,k\}.$ Clearly, $\mathrm {D} (A_1,\ldots ,A_k)$ is linear in each argument.
Let $K_1,\ldots ,K_m$ be convex bodies, and let $\lambda _1,\ldots ,\lambda _m\ge 0$ . Then, the volume of $\lambda _1K_1+\cdots \lambda _m K_m$ is an nth homogeneous polynomial of $\lambda _1,\ldots ,\lambda _m$ ,
The coefficients $\mathrm {V} \left (K_{i_{1}}, \ldots , K_{i_{n}}\right )$ , which will be called the mixed volumes, are nonnegative, symmetric in the indices, and dependent only on $K_{i_{1}}, \ldots , K_{i_{n}}$ .
If $\partial K$ is $C^2$ and has positive Gauss curvature everywhere, then
In addition, if $h_K=t_1h_1+\cdots t_mh_m$ , for some $h_i\in C^2$ and $t_i\in \mathbb {R}$ , we have
Here, $\mathrm {D} (D^2h_{i_1},\ldots ,D^2h_{i_{n-1}})(u)$ should be understood as the mixed discriminant of linear transforms restricted on $u^\perp $ .
Especially, notice that the spherical harmonics $(Y_{m,j})$ belong to $C^\infty (\mathbb {S}^{n-1})$ . When we are considering the numerical algorithm in Section 3, we represent the support function in finite-dimensional function space generated by the spherical harmonics in each step N, and (2.3) will be useful. Moreover, (2.3) also suggests the relation between mixed discriminant and mixed volume, which goes back to Hilbert–Minkowski–Aleksandrov (see [Reference Alexandrov1] for reference).
2.4 Sobolev space $H_k$
Let $k\ge 0$ be an integer. We denote $H_k(\mathbb S^{n-1})$ to be the Sobolev space, which is the completion of $C^\infty (\mathbb {S}^{n-1})$ with respect to the norm
where $\nabla ^m f$ denotes the mth covariant derivative of f. For a function $g\in H_k(\mathbb {S}^{n-1})$ and $m\le k$ , $\nabla ^m g\in L^2(\mathbb {S}^{n-1})$ denotes its weak covariant derivative, which can be defined as the limit of $\nabla ^m f_k$ in $L^2$ , where $f_k\in C^\infty (\mathbb {S}^{n-1})$ ( $k\in \mathbb {N}$ ) is an arbitrary Cauchy sequence with respect to the norm $\|\cdot \|_{H_k}$ . See [Reference Hebey18] for reference.
Let $\Delta _S: C^\infty (\mathbb {S}^{n-1})\to C^\infty (\mathbb {S}^{n-1})$ be the Laplace operator on $\mathbb {S}^{n-1}$ , which is self-adjoint. Extend it to $\Delta _S: H_{k+2}(\mathbb {S}^{n-1})\to H_k(\mathbb {S}^{n-1}).$ Then, for $f,g\in H_{k+2}(\mathbb {S}^{n-1})$ , it satisfies
Note that $(\cdot ,\cdot )$ always denotes the scalar product in $L^2(\mathbb {S}^{n-1})$ in this paper.
2.5 Spherical harmonics
The spherical harmonics are the eigenvectors (eigenfunctions) of spherical Laplace-Beltrami operator $\Delta _S:C^\infty (\mathbb {S}^{n-1})\to C^\infty (\mathbb {S}^{n-1})$ . It is well known that these eigenvectors can be chosen to be the restriction onto $\mathbb {S}^{n-1}$ of the harmonic polynomials (with respect to an orthonormal basis $e_1,\ldots ,e_n$ of $\mathbb {R}^n$ ) on $\mathbb {R}^n\backslash \{0\}.$ The eigenvalues of $\Delta _S$ are
With regard to $\lambda _m$ , the related eigenspace $E_m$ consists of all the harmonic polynomials on $\mathbb {R}^n\backslash \{0\}$ of degree m. The dimension of $E_m$ is
The compact operator theory in functional analysis tells that one can choose the eigenvectors to be a Hilbert basis of $L^2(\mathbb {S}^{n-1}).$ More precisely, we can choose $(Y_{m,j})$ such that the eigenspace of $\lambda _m$ is
and
For distinct $m,m'$ , we automatically have $\big (Y_{m,j},Y_{m',j'}\big ) = 0.$ Especially, when $m=1$ , the eigenfunctions $\{Y_{1,j}\}_{j=1}^n$ are chosen to be
so that $(Y_{1,k},Y_{1,j})=\delta _j^k.$ From this, it is easy to obtain the following formula that describes the barycenter of an integrable function g:
Since $(Y_{m,j})$ is a Hilbert basis, for each $g\in L^2(\mathbb S^{n-1})$ , we have
3 Reconstruction from a pre-given surface area measure
For a given finite Borel measure $\mu $ that is not concentrated in any hemisphere and satisfies
we shall give a computational procedure to reconstruct the convex body whose surface area measure is $\mu $ , and its convergence will also be proved.
For our convenience, only in this section, we relabel the spherical harmonics $\{Y_{m,j}: m=0,1,\ldots , j=1,\ldots ,N_m\}\setminus \{Y_{1,1},\ldots ,Y_{1,n}\}$ by $\{Y_m\}_{m=0}^\infty $ . That is to say, in $\{Y_i\}_{i=0}^\infty $ , there are no spherical harmonics of degree $1$ . As a result, if $h_Q = \sum _{i=0}^N x_i Y_i$ is a support function of a convex body Q. Then, by (2.4),
So, one convenience is to assign a specific translation of Q, since the surface area measure $S(Q,\cdot )$ is translation invariant.
Denote
and
Here, $\mathrm {D}$ is the mixed discriminant that can be computed by (2.2).
For a point $(x_0,\ldots ,x_N)\in \mathbb {R}^{N+1},$ let V be the volume function
3.1 Optimization problem ( $\overline P_N$ )
Originally, our aim in Step N is to solve the following optimization problem ( $P_N$ ):
where $\mathcal {C}_N$ is a convex set (will be proved later) defined by
Solve the optimization problem ( $\overline P_N$ ) and normalize the solution.
However, the problem ( $P_N$ ) cannot be solved directly by computer, since the requirement $(x_0,\ldots ,x_N)\in \overline {\mathcal {C}}_N $ has actually infinitely many inequality constraints, although it can be shown to be a convex optimization problem. Therefore, we consider the following optimization problem ( $P_{N,k}$ ), for $N,k\in \mathbb {N}$ . We will show in Theorem 3.3 that the subproblem ( $P_{N,k}$ ) is also a convex optimization problem with only finitely many constraints when k is large, and it follows from Theorems 3.3 and 3.5 that its solution will converge to the solution of the Minkowski problem as $N,k\to \infty $ .
Now, we put our computational procedure Mink-Pro and the subproblem ( $P_{N,k}$ ) as follows.
For each $k\in \mathbb {N}$ , we assign a kth partition $T_k = \{O_1^k,\ldots ,O_{2^k}^k\}$ of $\mathbb {S}^{n-1}$ such that
Extract one point from each $O_i^k$ to obtain a set $U_k\subset \mathbb {S}^{n-1}$ such that $U_k$ contains $2^k$ elements and
Define
and
Clearly, $\mathcal {C}_{N,k+1}\subset \mathcal {C}_{N,k}$ , and
Consider the modified optimization problem ( $ P_N^*$ ):
Denote its solution by $(x_0^*,\ldots ,x_N^*)$ , and define
3.2 Minkowski Procedure
-
• Input: Natural numbers $n,N,k$ , the general mixed volumes $\mathrm {V}_{i_1,\ldots ,i_n}$ and $ \mathrm {V}_i$ for the spherical harmonics $Y_0,\ldots ,Y_N$ defined by (3.1) and (3.2), the finite set $U_k$ satisfying (3.4) and (3.3).
-
• Task: Compute the function $h_{N,k}$ .
-
• Action:
-
1. Solve the following convex optimization problem ( $P_{N,k}$ ):
$$\begin{align*}\sup ~ V (x_0,\ldots,x_N), \end{align*}$$$$\begin{align*}\mathrm{subject~to} \quad(x_0,\ldots,x_N)\in \mathcal{C}_{N,k}.\end{align*}$$Denote the solution of ( $P_{N,k}$ ) by $(x_0^k,\ldots ,x_N^k).$
-
2. Compute
$$\begin{align*}h_{N,k} = \left(\sum_{i=0}^N x_i^kY_i\right) V(x_0^k,\ldots,x_N^k)^{- 1/({n-1})}.\end{align*}$$
-
Remark 3.1 1. We will show in Theorem 3.3 that ( $P_{N,k}$ ) and ( $P_{N}$ ) are convex optimization problems, for sufficiently large k. Thus, the maximizers $(x_0^k,\ldots ,x_0^k)$ and $(x_0^*,\ldots ,x_0^*)$ (theoretically) exist and are unique.
Especially, the approach to find the solution of ( $P_{N,k}$ ) is computationally tractable, since ( $P_{N,k}$ ) has only finitely many constraints. To write the detailed algorithm for the approach, see [Reference Boyd and Vandenberghe7, Reference Nocedal and Wright25] for references.
2. We will prove in Theorem 3.3(6) that the convex body, determined by $h_{N,k}$ ,
converges to $K_N$ as $k\to \infty $ . In addition, we will prove in Theorem 3.5 that $K_N\rightarrow K$ , when $N\to \infty $ , where $K_N$ is the convex body determined by $h_N$ , and $S(K,\cdot ) = \mu $ .
3. Recalling equation (2.3), $V(x_0,\ldots ,x_N)$ is actually the volume of the convex body Q whose support function is
for $(x_0,\ldots ,x_N)\in \mathcal {C}_{N}$ . The quantity $\sum _{i=0}^N \mathrm {V}_i\cdot x_i$ is $\int h_Q d\mu .$ This explains that Action 1 coincides with the optimization problem in solving the Minkowski problem. However, they are not completely the same, and we have to prove the convergence.
4. Although our procedure can reconstruct the body from a general measure, in practice, the smooth functions in the form
are more tractable. Since the $1$ -degree spherical harmonics $Y_{1,1},\ldots ,Y_{1,n}$ are removed from $Y_i$ ’s, such a function F always satisfies $\int _{\mathbb {S}^{n-1}} uF(u) du =0$ .
Lemma 3.1 Suppose K is a convex body such that
and
Then, there is a constant $c_1$ that only depends on $\mu $ , such that $\max \{|x|: x\in K\} \le c_1.$
Proof Since $\mu $ is not concentrated in any hemisphere, the function
is strictly positive on $S^{n-1}$ , and since the function is continuous on $S^{n-1}$ , which is compact, there exists $c_0>0$ such that
Denote $R_0 = \max \{|x| : x\in \bar K_N \}.$ Then, there exists $v_0\in \mathbb {S}^{n-1}$ such that $R_0\cdot v_0\in K.$ By (3.7), we must have $o\in K.$ Thus, we obtain
and hence $R_0\leq 1/c_0$ . Denote $c_1=1/c_0$ , and we complete the proof.
Lemma 3.2 (1) The sets $\overline {\mathcal {C}}_N, \mathcal {C}_N, \mathcal {C}_{N,k}\subset \mathbb {R}^{N+1}$ are nonempty, closed, and convex.
(2) $\overline {\mathcal {C}}_N$ is bounded.
Proof (1) Since $Y_0$ is a constant, $(c,0,\ldots ,0)\in \mathcal {C}_N$ for some sufficiently large constant c. Since $\mathcal {C}_N \subset \overline {\mathcal {C}}_N$ and $\mathcal {C}_N\subset \mathcal {C}_{N,k}$ , they are all nonempty. Since $ D^2\Big (\sum _{i=0}^N x_iY_i\Big ) (u)$ is continuous in $(x_0,\ldots ,x_N)$ , $\overline {\mathcal {C}}_N, \mathcal {C}_N$ , and $\mathcal {C}_{N,k}$ are closed.
If, at some $u\in \mathbb {S}^{n-1}$ ,
are positive semidefinite, then
is also positive semidefinite at this u, $\forall \lambda \in [0,1]$ . It follows immediately that $\overline {\mathcal {C}}_N, \mathcal {C}_N$ , and $\mathcal {C}_{N,k}$ are convex.
(2) By Lemma 3.1, for any $(x_0,\ldots ,x_N)\in \mathcal {C}_N$ , we have
where $c_1$ is a constant that only depends on $\mu $ . Then,
Since $Y_m\in C^\infty (\mathbb {S}^{n-1})$ are fixed functions, we deduce that $\overline {\mathcal {C}}_N$ is bounded.
Theorem 3.3 (1) $V(x_0,\ldots ,x_N)$ is strictly positive on $\mathcal {C}_N$ .
(2) The function $V^{1/n}$ is strictly concave on $\mathcal {C}_{N}$ .
(3) $\mathcal {C}_{N,k}\to \mathcal {C}_N$ in the Hausdorff metric; and $\mathcal {C}_N \subset \mathrm {int}~\overline {\mathcal {C}}_N$ .
(4) $(V)^{1/n}$ is strictly concave on $\mathcal {C}_{N,k}$ , for sufficiently large k.
(5) $(x_0^k,\ldots ,x_N^k)\to (x_0^*,\ldots ,x_N^*)$ as $k\to \infty $ .
(6) The convex body $K_{N,k}$ determined by $h_{N,k}$ converges to the convex body $K_N$ determined by $h_N$ (defined in (3.6)) as $k\to \infty $ .
Proof (1) We will use the fact that a smooth support function is either the support function of a convex body or the support function of a fixed point. By our construction, $\sum _{i=0}^N x_i Y_i$ is always smooth, and $\{Y_i\}_{i=0}^N$ does not contain the spherical harmonics of degree $1$ . Thus, for $(x_0,\ldots ,x_N)\in \mathcal {C}_N$ , it cannot be the support function of a fixed point. Thus, V is always positive.
(2) Since $(x_0,\ldots ,x_N)\in \mathcal {C}_N$ ,
is a support function of some convex compact set, which we denote by $K(x_0,\ldots ,x_N)$ . Now,
is its volume. Notice also that, for $(y_0,\ldots ,y_N)\in \mathcal {C}_N$ and $\lambda \in [0,1],$
Thus, statement (2) follows immediately from the Brunn–Minkowski inequality, statement (1), and the fact that
(3) Since $\mathcal {C}_{N,k}$ is convex and closed for each $k\in \mathbb {N}$ ,
and $\mathcal {C}_N$ is compact, statement (3) follows immediately.
(4) Since $\mathcal {C}_{N,k} \to \mathcal {C}_N\subset \mathrm {int} \overline {\mathcal {C}}_N$ , and $ \overline {\mathcal {C}}_N$ is compact, V is still positive on $\mathcal {C}_{N,k}$ when k is sufficiently large. Moreover, $\mathcal {C}_{N,k}$ will also be compact. Recall that
is a polynomial of degree n, where $\mathrm {V}_{i_1,\ldots ,i_n}$ are fixed constants. Thus, $V(\cdot ) C^\infty (\mathbb {R}^{N+1})$ . By the construction of $\mathcal {C}_{N,k}$ , when k is sufficiently large, there is a constant $\theta>0,$ such that the Hessian
everywhere on $\mathcal {C}_{N,k}$ .
(5) As a consequence of (3), for sufficiently large k, $\mathcal {C}_{N,k}$ is also compact. Thus, any subsequence of $(x_0^k,\ldots ,x_N^k)$ has a subsubsequence, which is also denoted by $(x_0^k,\ldots ,x_N^k)$ , converging to a point $(\bar x_0,\ldots ,\bar x_N) \in \bigcap _k \mathcal {C}_{N,k} =\mathcal {C}_N$ . It suffices to show that $(\bar x_0,\ldots ,\bar x_N)=(x_0^*,\ldots ,x_N^*)$ .
Since $(x_0^*,\ldots ,x_N^*)\in \mathcal {C}_{N,k}$ , for each $k\in \mathbb {N}$ , and $(x_0^k,\ldots ,x_N^k)$ is the maximizer of ( $P_{N,k}$ ) for sufficiently large k, we have
Taking limits in (3.8), we obtain
Since $(\bar x_0,\ldots ,\bar x_N)\in \mathcal {C}_N$ , and $(x_0^*,\ldots ,x_N^*)$ is the unique maximizer of ( $P_N^*$ ), we must have $(\bar x_0,\ldots ,\bar x_N)=(x_0^*,\ldots ,x_N^*)$ . This completes the proof.
(6) By (5), we have
uniformly on $\mathbb {S}^{n-1}$ . Recall the fact collected in Section 2.1: $f_i\to f$ uniformly and that $W[f]$ has interior imply that $W[f_i]\to W[f]$ in the Hausdorff metric. Now, (6) follows immediately from this property, the continuity of volume, and the definitions of $h_N$ and $h_{N,k}$ .
Lemma 3.4 [Reference Seeley29, Theorems 4 and 5]
Suppose $F\in C^\infty (\mathbb {S}^{n-1})$ satisfying $\int F(u)udu =0.$ Denote
Then:
-
(1) $F_N$ converges to F uniformly on $\mathbb {S}^{n-1}$ ;
-
(2) $D^2F_N$ converges to $D^2F$ uniformly on $\mathbb {S}^{n-1}$ .
Remark 3.2 This Lemma follows immediately from the estimates in [Reference Seeley29, Theorems 4 and 5], which suggest the relationship between the max norm and the $L^2$ -norm of Fourier expansion into spherical harmonics. The assumption $\int F(u)udu =0$ is because our $\{Y_i\}_{i=0}^\infty $ does not contain the spherical harmonics $Y_{1,1},\ldots ,Y_{1,n}$ . Except for this, the only difference is that the $D^\alpha $ in [Reference Seeley29, Theorem 4(b)] is slightly different from ours.
Let $f\in C^\infty (\mathbb {S}^{n-1})$ . Denote $\bar D^2 f(x/|x|)$ to be the Euclidean Hessian, and $D^2 f$ to be our notation in Section 2.1. Then, restricting to $x^\perp $ ,
and $\bar D^2 f(x/|x|) =0$ on the line $\mathbb {R} x$ .
Theorem 3.5 Denote the solution of ( $P_N$ ) by $(x_0^*,\ldots ,x_N^*),$ and
Let $K_N$ be the convex body determined by $h_N$ ,
Then, $K_N$ converges to K in the Hausdorff metric, as $N\to \infty $ . Here, K is the unique (theoretical) solution of the Minkowski problem
such that
Proof Denote
and denote $\bar K=\lambda K$ ( $\lambda>0$ ) to be the dilation of K such that $\int _{\mathbb {S}^{n-1}} h_{\bar K} d\mu = 1.$ To prove this theorem, it suffices to prove that $\bar K_N\to \bar K$ .
Step 1. Since $\{Y_m\}_{m=0}^\infty $ does not contain the $\{Y_{1,j}\}$ ’s, by equation (2.4), we always have
Step 2. It follows immediately from Lemma 3.1 that the sequence $\{\bar K_N\}_{N=1}^\infty $ is bounded.
Step 3. By Step 2 and Blaschke’s selection theorem, any subsequence of $\{\bar K_N\}_{N=1}^\infty $ has a subsubsequence $\bar K_{N_i}$ converging to a convex body $K^*$ in the Hausdorff metric. In this step, we will show that $K^*$ satisfies
and
The second equation follows immediately from the convergence. Let us consider (4.10).
-
• Step 3.1. Suppose Q is a smooth convex body that has positive Gauss curvature, and $\int _{\mathbb {S}^{n-1}}h_Qd\mu =1$ . We will show that $V(K^*)\ge V(Q).$ By the assumption, $h_Q\in C^\infty (\mathbb {S}^{n-1})$ , and
$$\begin{align*}\det (D^2h_Q)(u)> 0 , \quad \forall u\in\mathbb{S}^{n-1}.\end{align*}$$Denote$$\begin{align*}h_{Q_N} = \sum_{m=0}^N (h_Q,Y_{m})Y_m.\end{align*}$$Then, by Lemma 3.4(1), $h_{Q_N}$ converges to $h_Q$ uniformly. In addition, by Lemma 3.4(2), for sufficiently large N, we still have$$\begin{align*}\det (D^2h_{Q_N})(u)> 0, \quad \forall u\in\mathbb{S}^{n-1},\end{align*}$$and hence $h_{Q_N}$ is also support function. Therefore, if we denote$$\begin{align*}\bar Q_N =\frac1{\int_{\mathbb{S}^{n-1}}h_{Q_N}d\mu} Q_N, \end{align*}$$by the construction of $\bar K_N$ , we have $V(\bar K_N)\ge V(\bar Q_N)$ . Since $\bar K_N\to K^*$ , $\bar Q_N\to Q$ , we get$$\begin{align*}V(K^*)\ge V(Q).\end{align*}$$ -
• Step 3.2. For any convex body L, we can choose a sequence of smooth convex bodies $Q_m$ that have positive curvatures, and such that $Q_m\to L$ . By the approximation argument, we complete the proof of (4.10).
By the solution of the Minkowski problem, the convex body satisfying (4.10) and (3.11) is unique up to translations. Recalling the fact that
we have $\bar K= K^*$ . Since we have proved that any subsequence of $\{\bar K_N\}$ has a convergent subsequence converging to $\bar K$ , we have actually completed the proof of this theorem.
4 Reconstruction problems (II) for the lightness function
In this section, we consider the general lightness function $Q_K,$ which was defined by
where $f_w$ is positive and continuous in the open set $w^+\cap \mathbb {S}^{n-1}$ and vanishes in its complement. Before this, we need to develop some applications about the Funk–Hecke theorem.
4.1 Further applications of the Funk–Hecke theorem
The following is known as the Funk–Hecke theorem.
Lemma 4.1 [Reference Schneider28]
Let $f \in L_2[-1,1]$ be a function satisfying
If $Y_m \in \mathscr {F}^m$ is a spherical harmonic of degree m, then
for $v\in \mathbb {S}^{n-1}$ , where
where
Lemma 4.2 [Reference Schneider28]
For the function $f(t)=|t|$ , the eigenvalue $\lambda _m(f)$ in the Funk–Hecke theorem is
for even m, and $\lambda _{m}(|t|)=0$ for odd m. Furthermore, $|\lambda _m(|t|)^{-1}|=O\left (m^{(n+2)/2}\right )$ for even $m\rightarrow \infty .$
Since the lightness functions are defined in the form of partial cosine transform
we require the following Lemma to compute its eigenvalues.
Lemma 4.3 Let $f(t)=t_-$ , and denote $\beta _m $ to be the $\lambda _m(f)$ in the Funk–Hecke theorem. Then,
for even m, $\beta _{m}=0$ if m is odd and $m\geq 3,$ and
Applying Lemma 4.1, one can directly compute the coefficients $\beta _m$ , but it will be complicated when $m\ge 3$ is odd. However, we would like to provide the following simpler method, to understand the geometric fact that the partial cosine transform $\mathcal {C}_- g$ also defines a support function of a zonoid (which may not be symmetric with respect to the origin).
Proof Since $t_-= 1/2(|t| - t)$ , we have
Denote $a_{m,j} = (g,Y_{m,j})$ . By the above equation and Lemmas 4.1 and 4.2, we have
Using Lemma 4.1 again, for the function $(v\cdot u)_-,$ we know that
By (2.4), we have
By this and Lemma 4.2, comparing (4.2) with (4.3), we get the desired result.
Lemma 4.4 (Theorem 3.5.4 [Reference Schneider28])
If $\rho $ is an even signed measure on $\mathbb {S}^{n-1}$ with
then $\rho =0$ .
If G is an even real function on $\mathbb {S}^{n-1}$ of differentiability class $C^k$ , where $k\geq n+2$ is even, then there exist an even continuous function g on $\mathbb {S}^{n-1}$ such that
for $v\in \mathbb {S}^{n-1}$ .
The above Lemma and its proof can be found in [Reference Schneider28]. Since we want to work with the nonsmooth convex bodies, we slightly improve this result to the functions in Sobolev spaces. Precisely, we consider the functions in ${H}_{k}(\mathbb {S}^{n-1})$ . Although the approach is standard, we did not find the result as follows. So, for completeness, we give a proof here.
Lemma 4.5 Let $k\ge n+1$ be an even integer. If $G \in H_k(\mathbb S^{n-1})$ is an even function, then there exists an even function $g\in L^2(\mathbb {S}^{n-1})$ such that
for all $v\in \mathbb {S}^{n-1}$ . Especially, g can be formulated by the convergent function series
where $\lambda _m^{-1}(|t|)$ is provided by the Funk–Hecke theorem.
Proof Denote $k=2l$ . Since $G\in H_k(\mathbb {S}^{n-1}),$ we have
We aim to prove that the series
converges in $ L^2(\mathbb {S}^{n-1})$ .
Recall that
Since $\Delta _{S}$ is self-adjoint, we have
By the Cauchy–Schwarz inequality, we have
As a result, the right-hand side of (4.7) satisfies
Since $G\in H_k$ , it suffices to prove that
Here,
If $l>n/2+1/4$ , the inequality (4.8) holds for
Therefore, when $k=2l\ge n+1$ , we have the desired result.
Lemmas 4.3 and 4.5 immediately lead to the following.
Corollary 4.6 Let $k\ge n+1$ be an even integer, and let $G\in H_k(\mathbb {S}^{n-1})$ . Suppose $g\in L^2(\mathbb {S}^{n-1})$ satisfies
Then,
and
Moreover, $(G,Y_{m,j}) = 0$ , for each odd integer $m\geq 3$ .
Proof Denote $\tilde g(u) = g(u) + g(-u), a_{m,j} = (g,Y_{m,j}),$ and $b_{m,j} = (\tilde g,Y_{m,j}).$ When m is even, $Y_{m,j}$ is even, and thus $b_{m,j} = 2a_{m,j}$ . On one hand, by Lemmas 4.1 and 4.2, we have
On the other hand,
Comparing the coefficients, we have the following. For each odd integer $m\geq 3$ , we have $(G,Y_{m,j}) = 0$ . For $m=1$ , we have $(G,Y_{1,j}) =- ({\omega _n}/2) a_{1,j}.$ Applying (2.4) for g and G, respectively, we get
For even m, we have
Since $\tilde g$ is an even function, we get
The convergence of the above function series follows from Lemma 4.5 and the fact that $2\beta _m = \lambda _m(|t|)$ .
4.2 The necessary and sufficient conditions of being a lightness function
Given a binary function $G: \mathbb {S}^{n-1}\times \mathbb {S}^{n-1} \rightarrow \mathbb {R}$ so that $G(\cdot ,w)$ is Sobolev for a fixed $w\in \mathbb {S}^{n-1}$ , we shall use the notation
Lemma 4.7 Let K be a convex body, and let $k\ge n+1$ be an even integer. If, for each $i\in \{1,\ldots ,n\}$ , its lightness functions $Q_K(\cdot ,e_i), Q_K(\cdot ,-e_i)\in H_k(\mathbb {S}^{n-1})$ , then the surface area measure $S(K,\cdot )$ has a density in $L^2(\mathbb {S}^{n-1})$ .
Proof Denote
which is an even function in $H_k(\mathbb {S}^{n-1})$ . On one hand,
On the other hand, by Lemma 4.5, there is an even function $F_{K,i}\in L^2(\mathbb {S}^{n-1})$ such that
Since the cosine transform is injective (Lemma 4.4), we have
Note that, for each $u\in \mathbb {S}^{n-1},$ there exists an $i\in \{1,\ldots ,n\}$ such that $|u_i|\ge 1/\sqrt n.$ Since $|\tilde f_{e_i}|$ is positive and continuous in $\{u: |u_i|\geq 1/\sqrt n\}$ , it has a strict positive lower bound. Thus,
has a density in $L^2(\mathbb {S}^{n-1})$ , where the function $F_{K,i}/(f_{e_i}+ f_{-e_i})$ must be independent of i.
Theorem 4.8 Let K be a convex body, and let $k\ge n+1$ be an even integer. Suppose $Q_K(\cdot ,w ) \in H_k(\mathbb {S}^{n-1})$ , for any $w\in \mathbb {S}^{n-1}$ . Then:
-
(1) The curvature function $F_K$ satisfies
$$\begin{align*}\frac12 f_w(u)F_K(u) = 1_{w^+}(u)\cdot\sum_{m \mathrm{~even}}\sum_{j=1}^{N_m} \beta_m^{-1}\big(Q_K(\cdot,w),Y_{m,j}\big) Y_{m,j}(u).\end{align*}$$ -
(2) For each odd integer $m\geq 3$ , and for each fixed $w\in \mathbb {S}^{n-1}$ , the function $Q_K(\cdot ,w ): \mathbb {S}^{n-1}\rightarrow \mathbb {R}$ satisfies
$$\begin{align*}\big(Q_K(\cdot,w ),Y_{m,j}\big)=0, \quad \forall j=1,\ldots,N_{m}. \end{align*}$$ -
(3) Moreover, $F_K$ satisfies
$$\begin{align*}\int_{\mathbb{S}^{n-1}}vQ_K(v,w) dv = -\frac{\omega_n}2 \int_{\mathbb{S}^{n-1}} u f_w F_K(u)du , \quad \forall i=1,\ldots,n.\end{align*}$$
Proof By Lemma 4.7, $S(K,\cdot )$ has an $L^2(\mathbb {S}^{n-1})$ density $F_K:\mathbb {S}^{n-1}\rightarrow [0,\infty )$ . Substitute the function $f_w F_K$ into Corollary 4.6, and observe that
Then, we immediately obtain (2) and (3), and
Since $f_w(u)$ is concentrated on $w^+, f_w(-u) F_K(-u)$ is just a symmetric copy of $f_w(u) F_K(u)$ in $w^-$ . Thus,
This explains (1).
Theorem 4.9 Let $G: \mathbb {S}^{n-1}\times \mathbb {S}^{n-1} \rightarrow \mathbb {R}$ be a binary function, and let $k\ge n+1$ be an even integer. Suppose $G(\cdot ,w)\in H_k(\mathbb {S}^{n-1})$ , for any fixed $w\in \mathbb {S}^{n-1}$ . In addition, suppose G satisfies the following conditions:
-
(1) The following function
(4.9) $$ \begin{align} F(u) = \begin{cases} \sum\limits_{m~\mathrm{even}} \sum\limits_{k=1}^{N_m} \beta_m^{-1}\big(G(\cdot,w),Y_{m,k}\big) \frac{Y_{m,k}(u)}{f_w(u)} , \ \ & u\cdot w>0 \\ \phantom{\cdot} &\phantom{\cdot} \\ \sum\limits_{m~\mathrm{even}} \sum\limits_{k=1}^{N_m} \beta_m^{-1}\big(G(\cdot,-w),Y_{m,k}\big) \frac{Y_{m,k}(u)}{f_{-w}(u)} , \ \ & u\cdot w<0 \end{cases} \end{align} $$is nonnegative and independent of $w\in \mathbb {S}^{n-1}$ . -
(2) For each odd integer $m\geq 3$ , and for each fixed $w\in \mathbb {S}^{n-1}$ , the function $G(\cdot ,w): \mathbb {S}^{n-1}\rightarrow \mathbb {R}$ satisfies
$$\begin{align*}\big(G(\cdot,w),Y_{m,j}\big)=0, \quad \forall k=1,\ldots,N_{m}. \end{align*}$$ -
(3) The function F in (1) satisfies
$$\begin{align*}\int_{\mathbb{S}^{n-1}} uF(u) du =0.\end{align*}$$ -
(4) G and the function F in (1) satisfy
$$\begin{align*}\frac 2{\omega_n}\int_{\mathbb{S}^{n-1}} vG(v,w)dv = \int_{\mathbb{S}^{n-1}} uf_w(u)F(u)du.\end{align*}$$
Then, there is a convex body K such that
and K is unique up to translations.
Condition (1) means that for any $w_1,w_2\in \mathbb {S}^{n-1},$ the related two function series are equivalent almost everywhere.
Proof By Lemma 4.5, for each fixed $w\in \mathbb {S}^{n-1}$ , the function series
Let $e_1,\ldots ,e_n$ be an orthonormal basis of $\mathbb {R}^n$ , and denote
By (4.9) and (4.10), we infer that $F \in L^2(\Omega _i)$ , $i=1,\ldots ,n.$ Since $\{\Omega _i\}_{i=1}^n$ covers $\mathbb {S}^{n-1}$ , we see that $F \in L^2(\mathbb {S}^{n-1}).$ Since $\mathbb {S}^{n-1}$ is compact, F is also in $L^1(\mathbb {S}^{n-1})$ . By conditions (1) and (3), F is independent of w and is of centroid $0$ . Then, there is a unique (up to translations) convex body K that solves the classical Minkowski problem
Notice that
For this K, by Corollary 4.6, Lemma 4.3, condition (4), equation (2.4), and finally by condition (2), we have
That the convex body K is unique follows from Theorem 4.8 and the uniqueness of the Minkowski problem.
5 Reconstruction from the pre-given lightness function
In this section, for a pre-given lightness function $Q_K(v,w)$ together with a function $f_w$ as in (1.2), we give the following numerical procedure to reconstruct K.
Light-Pro
-
• Input: Natural numbers $n,M,N,k$ , the spherical harmonics $Y_{m,j}$ , where $m=1,\ldots ,M$ , and $j=1,\ldots ,N_m$ .
-
• Task: Compute the function $h_{M,N,k}$ .
-
• Action:
-
1. Fix a $w\in \mathbb {S}^{n-1}$ . Compute
$$\begin{align*}A_{m,j} = \beta_m^{-1} (Q(\cdot,w), Y_{m,j})\end{align*}$$and$$\begin{align*}B_{m,j} = \beta_m^{-1} (Q(\cdot,-w), Y_{m,j}).\end{align*}$$ -
2. Compute
(5.1) $$ \begin{align} F_M(u) = \begin{cases} \sum\limits_{\substack{m=0\\m~\mathrm{even}}}^M \sum\limits_{k=1}^{N_m} A_{m,j} \frac{Y_{m,k}(u)}{f_w(u)} , \ \ & u\cdot w>0, \\ \phantom{\cdot} &\phantom{\cdot} \\ \sum\limits_{\substack{m=0\\m~\mathrm{even}}}^M \sum\limits_{k=1}^{N_m} B_{m,j}\frac{Y_{m,k}(u)}{f_{-w}(u)} , \ \ & u\cdot w<0. \end{cases} \end{align} $$ -
3. Compute
(5.2) $$ \begin{align}\tilde F_M(u) = F_M(u) - \frac1{\omega_n}\int_{\mathbb{S}^{n-1}} (u\cdot v) F_M(v) dv.\end{align} $$ -
4. If $\tilde F_M\ge 0,$ then apply the Mink-Pro to find the solution of the ( $P_{N,k}$ ) problem for
$$\begin{align*}d\mu =\tilde F_M du,\end{align*}$$and denote the resulting function by $h_{M,N,k}$ .
-
Remark 5.1 1. On one hand, it follows from a direct computation that the $\tilde F_M$ defined in (5.2) satisfies
On the other hand, since $f_w$ is positive and continuous on $w^+\cap \mathbb {S}^{n-1},$ by Lemma 3.4 (1) and (5.1), we see that $ F_M\rightarrow F_K$ uniformly. By this and
we also have $\tilde F_M\rightarrow F_K$ uniformly, too.
2. If $\partial K$ is smooth and has positive curvature, we can prove in Theorem 5.1 that $K_{M,N,k}\to K$ as $M,N,k\to \infty $ , where $K_{M,N,k}$ is the convex body determined by $h_{M,N,k}$ .
Theorem 5.1 If $\partial K$ is smooth and has positive curvature, then $K_{M,N,k}\to K$ in Hausdorff metric as $M,N,k\to \infty $ .
Proof By Remark 5.1(1), $\tilde F_M\rightarrow F_K$ uniformly. Since $F_K$ is smooth and has positive curvature, when M is sufficiently large, we also have $\tilde F_M>0$ . In addition, for any $g\in C(\mathbb {S}^{n-1}),$
That is to say, $\tilde F_M(u)du\to F_K(u)du$ weakly. Let $K_M$ be the convex body whose curvature function is $\tilde F_M$ . Since $K_M,K$ are convex bodies such that
$K_M\rightarrow K$ in the Hausdorff metric. This is a classical result that can be proved by using Blaschke’s selection theorem (see [Reference Schneider28] for reference). The assertion follows from this and Theorem 3.3.
When $f_w=1_{w^+\cap \mathbb {S}^{n-1}}$ , the characteristic function of the open hemisphere, $Q_K$ becomes the partial brightness function (which is a binary function rather than $b_K $ )
Recalling that the brightness function $b_K$ can only determine an origin-symmetric convex body, for the partial brightness function, we are able to construct the convex body without symmetric assumption. It is reasonable because the change of viewing and illumination directions provides more information.
Acknowledgment
We thank Professor Boris Rubin for his clarifications on the cosine transforms. We are grateful to Professor Richard J. Gardner for his evaluations and suggestions on an earlier manuscript. We appreciate the anonymous referee’s insightful remarks, which helped us to make our paper better.