Image recovery in optical interferometry is an ill-posed nonlinear inverse problem arising from incomplete power spectrum and bi-spectrum measurements. We formulate a linear version of the problem for the order-3 tensor formed by the tensor product of the signal with itself. This linear problem is regularized by standard convex \ell_1 relaxations of sparsity and low rank constraints and solved using the most advanced algorithms in convex optimization. We show preliminary results on small size synthetic images as a proof of concept.