Convex Optimization
- ISSN: 10556788
- ISBN: 9780521833783
- DOI: 10.1080/10556781003625177
- PubMed: 20876008
Abstract
We are developing a dual panel breast-dedicated PET system using LSO scintillators coupled to position sensitive avalanche photodiodes (PSAPD). The charge output is amplified and read using NOVA RENA-3 ASICs. This paper shows that the coincidence timing resolution of the RENA-3 ASIC can be improved using certain list-mode calibrations. We treat the calibration problem as a convex optimization problem and use the RENA-3s analog-based timing system to correct the measured data for time dispersion effects from correlated noise, PSAPD signal delays and varying signal amplitudes. The direct solution to the optimization problem involves a matrix inversion that grows order (n3) with the number of parameters. An iterative method using single-coordinate descent to approximate the inversion grows order (n). The inversion does not need to run to convergence, since any gains at high iteration number will be low compared to noise amplification. The system calibration method is demonstrated with measured pulser data as well as with two LSO-PSAPD detectors in electronic coincidence. After applying the algorithm, the 511keV photopeak paired coincidence time resolution from the LSO-PSAPD detectors under study improved by 57%, from the raw value of 16.30.07 ns FWHM to 6.920.02 ns FWHM (11.520.05 ns to 4.890.02 ns for unpaired photons).
Convex Optimization
objective function. We denote a general vector optimization problem as
minimize (with respect to K) f0(x)
subject to fi(x) ≤ 0, i = 1, . . . ,m
hi(x) = 0, i = 1, . . . , p.
(4.56)
Here x ∈ Rn is the optimization variable, K ⊆ Rq is a proper cone, f0 : Rn → Rq
is the objective function, fi : Rn → R are the inequality constraint functions, and
hi : Rn → R are the equality constraint functions. The only difference between this
problem and the standard optimization problem (4.1) is that here, the objective
function takes values in Rq, and the problem specification includes a proper cone
K, which is used to compare objective values. In the context of vector optimization,
the standard optimization problem (4.1) is sometimes called a scalar optimization
problem.
We say the vector optimization problem (4.56) is a convex vector optimization
problem if the objective function f0 isK-convex, the inequality constraint functions
f1, . . . , fm are convex, and the equality constraint functions h1, . . . , hp are affine.
(As in the scalar case, we usually express the equality constraints as Ax = b, where
A ∈ Rp×n.)
What meaning can we give to the vector optimization problem (4.56)? Suppose
x and y are two feasible points (i.e., they satisfy the constraints). Their associated
objective values, f0(x) and f0(y), are to be compared using the generalized inequal-
ity K . We interpret f0(x) K f0(y) as meaning that x is ‘better than or equal’ in
value to y (as judged by the objective f0, with respect to K). The confusing aspect
of vector optimization is that the two objective values f0(x) and f0(y) need not be
comparable; we can have neither f0(x) K f0(y) nor f0(y) K f0(x), i.e., neither
is better than the other. This cannot happen in a scalar objective optimization
problem.
4.7.2 Optimal points and values
We first consider a special case, in which the meaning of the vector optimization
problem is clear. Consider the set of objective values of feasible points,
O = {f0(x) | ∃x ∈ D, fi(x) ≤ 0, i = 1, . . . ,m, hi(x) = 0, i = 1, . . . , p} ⊆ Rq,
which is called the set of achievable objective values. If this set has a minimum
element (see §2.4.2), i.e., there is a feasible x such that f0(x) K f0(y) for all
feasible y, then we say x is optimal for the problem (4.56), and refer to f0(x) as
the optimal value of the problem. (When a vector optimization problem has an
optimal value, it is unique.) If x? is an optimal point, then f0(x?), the objective
at x?, can be compared to the objective at every other feasible point, and is better
than or equal to it. Roughly speaking, x? is unambiguously a best choice for x,
among feasible points.
A point x? is optimal if and only if it is feasible and
O ⊆ f0(x?) +K (4.57)
Applications
Sign up today - FREE
Mendeley saves you time finding and organizing research. Learn more
- All your research in one place
- Add and import papers easily
- Access it anywhere, anytime





