Multireference Alignment Is Easier With an Aperiodic Translation Distribution

被引:30
作者
Abbe, Emmanuel [1 ,2 ]
Bendory, Tamir [1 ]
Leeb, William [1 ,3 ]
Pereira, Joao M. [1 ]
Sharon, Nir [1 ,4 ]
Singer, Amit [1 ,5 ]
机构
[1] Princeton Univ, Program Appl & Computat Math, Princeton, NJ 08544 USA
[2] Princeton Univ, Dept Elect Engn, Princeton, NJ 08544 USA
[3] Univ Minnesota Twin Cities, Sch Math, Minneapolis, MN 55455 USA
[4] Tel Aviv Univ, Sch Math Sci, IL-6997801 Tel Aviv, Israel
[5] Princeton Univ, Dept Math, Princeton, NJ 08544 USA
关键词
Multireference alignment; spectral algorithm; method of moments; spiked covariance model; non-convex optimization; expectation-maximization; cryo-EM; MAXIMUM-LIKELIHOOD; SYNCHRONIZATION; RECONSTRUCTION; ALGORITHMS;
D O I
10.1109/TIT.2018.2889674
中图分类号
TP [自动化技术、计算机技术];
学科分类号
0812 ;
摘要
In the multireference alignment model, a signal is observed by the action of a random circular translation and the addition of Gaussian noise. The goal is to recover the signal's orbit by accessing multiple independent observations. Of particular interest is the sample complexity, i.e., the number of observations/samples needed in terms of the signal-to-noise ratio (SNR) (the signal energy divided by the noise variance) in order to drive the mean-square error to zero. Previous work showed that if the translations are drawn from the uniform distribution, then, in the low SNR regime, the sample complexity of the problem scales as omega(1/SNR3). In this paper, using a generalization of the Chapman-Robbins bound for orbits and expansions of the chi(2) divergence at low SNR, we show that in the same regime the sample complexity for any aperiodic translation distribution scales as omega(1/SNR2). This rate is achieved by a simple spectral algorithm. We propose two additional algorithms based on non-convex optimization and expectation-maximization. We also draw a connection between the multireference alignment problem and the spiked covariance model.
引用
收藏
页码:3565 / 3584
页数:20
相关论文
共 51 条
[1]  
Abbe E, 2018, IEEE INT SYMP INFO, P561, DOI 10.1109/ISIT.2018.8437646
[2]  
Abbe E, 2017, IEEE INT SYMP INFO, P1316, DOI 10.1109/ISIT.2017.8006742
[3]   Fundamental Limits in Multi-Image Alignment [J].
Aguerrebere, Cecilia ;
Delbracio, Mauricio ;
Bartesaghi, Alberto ;
Sapiro, Guillermo .
IEEE TRANSACTIONS ON SIGNAL PROCESSING, 2016, 64 (21) :5707-5722
[4]  
[Anonymous], 1998, STAT SHAPE ANAL
[5]  
[Anonymous], [No title captured]
[6]  
[Anonymous], 2015, Non-unique games over compact groups and orientation estimation in cryo-em
[7]  
BANDEIRA A. S., 2016, P C LEARN THEOR, V49, P361
[8]  
Bandeira A. S., 2017, OPTIMAL RATES ESTIMA
[9]  
Bandeira Afonso S, 2014, P 5 C INN THEOR COMP, P459
[10]   2.2 Å resolution cryo-EM structure of β-galactosidase in complex with a cell-permeant inhibitor [J].
Bartesaghi, Alberto ;
Merk, Alan ;
Banerjee, Soojay ;
Matthies, Doreen ;
Wu, Xiongwu ;
Milne, Jacqueline L. S. ;
Subramaniam, Sriram .
SCIENCE, 2015, 348 (6239) :1147-1151