recent observations of the cosmic microwave background ( @xmath0 ) and of the large scale structure of the universe ( @xmath3 ) have led to the definition of a concordance cosmological model .
the analysis of the @xmath0 temperature data over the whole celestial sphere from the wilkinson microwave anisotropy probe ( @xmath4 ) satellite experiment has been playing a dominant role in designing this precise picture of the universe @xcite . according to this model , the cosmic structures and the @xmath0 originate from gaussian adiabatic perturbations seeded in the early phase of inflation of the universe .
however , cosmological scenarios motivated by theories for the unification of the fundamental interactions predict the existence of topological defects , resulting from phase transitions at the end of inflation @xcite .
these defects would have participated to the formation of the cosmic structures , also imprinting the @xmath0 .
cosmic strings are a line - like version of such defects which are also predicted in the framework of fundamental string theory @xcite . as a consequence ,
the issue of their existence is a central question in cosmology today .
cosmic strings are parametrized by a string tension @xmath5 , i.e. a mass per unit length , which sets the overall strength of a string network .
their main signature in the @xmath0 is characterized by temperature steps along the string positions .
this localized effect , known as the kaiser - stebbins ( @xmath1 ) effect @xcite , hence implies a non - gaussian imprint of the string network in the @xmath0 .
the most numerous strings appear at an angular size around @xmath6 @xmath7 on the celestial sphere .
@xmath0 experiments with an angular resolution much below @xmath6 @xmath7 are thus required in order to resolve the width of cosmic strings .
experimental constraints have been obtained on a possible string contribution in terms of upper bounds on the string tension @xmath5 @xcite . in this context , even though observations largely fit with an origin of the cosmic structures in terms of adiabatic perturbations , room is still available for the existence of cosmic strings .
the observed @xmath0 signal can be modeled as a linear superposition of a statistically isotropic but non - gaussian string signal @xmath8 proportional to an unknown string tension , with statistically isotropic gaussian noise @xmath9 comprising the standard components of the @xmath0 induced by adiabatic perturbations , as well as instrumental noise . in a first recent work @xcite we developed an effective method for mapping the string network potentially imprinted in @xmath0 temperature data , in the perspective of forthcoming arcminute resolution experiments performing a mapping of planar patches of the sky in real space , with standard radiometers .
we took a bayesian approach to this denoising problem , based on statistical models for both the string signal and the noise .
the denoising is done in the wavelet domain , using a steerable wavelet transform well adapted for representing the strongly oriented features present in the string signal .
we showed that the string signal coefficients are well modeled by generalized gaussian distributions ( @xmath10 s ) , which are fitted at each wavelet scale using a training simulation borrowed from the set of realistic string simulations recently produced by @xcite .    in the present work we aim at extending the approach described above to the reconstruction of a string signal from radio - interferometric observations of the @xmath0 probing small patches of the celestial sphere at the same arcminute resolution .
under the most standard assumptions , radio interferometers provide a noisy and incomplete fourier coverage of the planar signal . in these terms
, the signal reconstruction amounts to solving an ill - posed inverse problem encompassing both denoising and deconvolution . relying on the @xmath1 effect in the @xmath0 ,
the string network itself can be mapped by the magnitude of the gradient @xmath11 of the string signal , which is very sparse or compressible in nature , with significant amplitudes only along the strings themselves .
the theory of compressed sensing offers a new framework for solving ill - posed inverse problem when the signals are sparse or compressible @xcite .
a band - limited signal may be expressed as the @xmath12-dimensional vector of its values sampled at the nyquist - shannon rate . by definition ,
a signal is sparse in some basis if its expansion contains only a small number @xmath13 of non - zero coefficients .
more generally it is compressible if its expansion only contains a small number of significant coefficients , i.e. if a large number of its coefficients bear a negligible value .
the theory of compressed sensing demonstrates that , for sparse or compressible signals , a small number @xmath14 of random measurements , in a sensing basis incoherent with the sparsity basis , will suffice for an accurate and stable reconstruction of the signals .
the mutual coherence between two bases may be defined as the maximum complex modulus of the scalar product between unit - norm vectors of the two bases .
a random selection of fourier measurements of a signal sparse in real space or in its gradient are particular examples of good sensing procedures . in this framework
, the reconstruction of the signals may be approached in different ways , but all of them encompass a regularization of the ill - posed inverse problem by the introduction of a sparsity or compressibility constraint . in a second recent work @xcite , we presented results showing that compressed sensing offers powerful image reconstruction techniques for radio - interferometric data , which are stable relative to noise , as well as stable relative to non - exact sparsity , i.e. compressibility of the signals .
these techniques are based on global minimization problems , which are solved by convex optimization algorithms @xcite .
we particularly illustrated the versatility of the scheme relative to the inclusion of specific prior information on the signal in the minimization problems .
we here wish to apply these results for the reconstruction of a string signal in the @xmath0 . as for the denoising algorithm proposed in @xcite , the reconstruction algorithm that we present
may be considered as a modular component of a larger data analysis .
it may notably find use as a pre - processing step for other methods for cosmic string detection based on explicit edge detection @xcite .
our reconstruction algorithm is tested under different conditions with astrophysical noise components including various contributions to the standard components of the @xmath0 , i.e. primary and secondary anisotropies .
our analyses rely on the set of realistic string simulations recently produced by @xcite .
the remainder of this paper is organized as follows . in section [ sec : string signal and noise ] , we discuss the string signal and noise in the @xmath0 . in section [ sec : interferometric inverse problem ] , we describe the inverse problem posed in the context of radio interferometry . in section [ sec : compressed sensing reconstruction ] , we define the minimization problems to be solved for signal reconstruction in the context of compressed sensing , notably relying on a prior statistical model of the string signal fitted from simulations . in section [ sec : analyses and results ] , we describe our numerical analyses and results . we finally conclude in section [ sec : conclusion ] .
in this section , we describe the string signal in the @xmath0 , as well as the astrophysical noise made up of the standard primary and secondary anisotropies .
we also describe the numerical simulations used for our subsequent analyses .      in an inflationary cosmological model ,
the phase transitions responsible for the formation of a cosmic string network occur after the end of inflation , so as to produce observable defects . from the epoch of last scattering until today
, the cosmic string network continuously imprints the @xmath0 .
the so - called scaling solution for the string network implies that the most numerous strings are imprinted just after last scattering and have a typical angular size around @xmath6 @xmath7 , of the order of the horizon size at that time @xcite .
longer strings are also imprinted in the later stages of the universe evolution , but in smaller number , according to the number of corresponding horizon volumes required to fill the sky .
the main signature of a cosmic string in the @xmath0 is described by the @xmath1 effect according to which a temperature step is induced along the string position .
the relative amplitude of this step is given by @xmath15 where @xmath16 and @xmath17 with @xmath18 standing for the string velocity transverse to the line of sight and @xmath19 for the speed of light , and where @xmath20 is a dimensionless parameter uniquely associated with the string tension @xmath5 through @xmath21 with @xmath22 standing for the gravitational constant . in the following we call @xmath20 the string tension .
on small angular scales , realistic simulations can be produced by stacking @xmath0 maps induced in different redshift ranges between last scattering and today .
the simulations that we use in this work have been produced by this technique @xcite .
the string signal is understood as a realization of a statistically isotropic but non - gaussian process on the celestial sphere with an overall amplitude rescaled by the string tension @xmath20 , and characterized by a nearly scale - free angular power spectrum : @xmath23 , where the positive integer index @xmath24 stands for the angular frequency index on the sphere .
an analytical expression of this spectrum was provided for @xmath24 larger than a few hundreds by @xcite , on the basis of their simulations .
we consider radio - interferometric @xmath0 observations with a small field of view corresponding to an angular opening @xmath25 on the celestial sphere . in this context
, the small portion of the celestial sphere accessible is identified to a planar patch of size @xmath26 , and we may consider planar signals functions of a two - dimensional position vector @xmath27 .
the corresponding spatial frequencies may be denoted as two - dimensional vectors @xmath28 with a radial component given by the norm @xmath29 of the vector . in this euclidean limit
, the radial component identifies with the angular frequency on the celestial sphere , below some band limit @xmath30 set by the resolution of the interferometer under consideration : @xmath31 . for @xmath29 larger than a few hundreds , the nearly scale - free planar power spectrum of the string signal @xmath32 reads as : @xmath33 with @xmath34 for @xmath35 .
the observed @xmath0 signal can be understood as a linear superposition of the string signal and statistically isotropic noise of astrophysical and instrumental origin @xmath36 with some angular power spectrum @xmath37 . in the euclidean limit
, the corresponding planar power spectrum for the noise @xmath36 may be written as @xmath38 for @xmath35 .
let us note that we consider zero mean signals , identifying perturbations around statistical means .
angular power spectra of the string signal and noise ( borrowed from @xcite ) as functions of the angular frequency in the range @xmath39 $ ] in @xmath40 axes scaling .
the spectrum of the string signal is represented for a string tension @xmath41 in terms of its analytical expression valid at high angular frequencies ( red straight line ) .
the noise spectra ( ordered by decreasing amplitude at low angular frequencies ) are : the primary @xmath0 anisotropies ( black solid line ) and the gravitational lensing correction ( black dashed line ) , the thermal @xmath42 effect in the rayleigh - jeans limit ( blue solid line ) , the non - linear kinetic @xmath42 effect ( green solid line ) , and the ostriker - vishniac effect ( magenta solid line).,width=302 ]    current @xmath0 experiments achieve an angular resolution on the celestial sphere of the order of @xmath43 arcminutes , corresponding to a limit angular frequency not far above @xmath44 . at such resolutions ,
the standard components of the @xmath0 primarily contain the statistically isotropic and gaussian primary anisotropies induced by adiabatic perturbations at last scattering . in this context
, any possible string signal is confined to amplitudes largely dominated by these primary anisotropies .
the constraints mainly come from a best fit analysis of the angular power spectrum of the overall @xmath0 signal in the @xmath4 temperature data @xcite .
the tightest of these constraints @xcite gives the following upper bound at @xmath45 per cent confidence level : @xmath46    forthcoming experiments will provide access to higher angular resolution . among other instruments , radio interferometers such as the arcminute microkelvin imager ( ami ) will probe the @xmath0 at angular resolutions below @xmath6 @xmath47 @xcite . at these resolutions , the secondary @xmath0 anisotropies induced by interaction of @xmath0 photons with the evolving universe after last scattering will dominate the primary anisotropies , and must be accounted for
.    the primary anisotropies exhibit exponential damping at high angular frequencies .
this contrasts with the slow decay of the nearly scale - free angular power spectrum of the string signal , which thus dominates over the primary anisotropies at high enough angular frequencies .
the secondary anisotropies include gravitational effects such as the integrated sachs - wolfe ( @xmath48 ) effect , the rees - sciama ( @xmath49 ) effect , and gravitational lensing , as well as re - scattering effects such as the thermal and kinetic sunyaev - zeldovich ( @xmath42 ) effects .
the @xmath42 effects dominate these secondary anisotropies @xcite . the @xmath48 and @xmath49 effects associated with the time evolution of the standard gravitational potentials can be neglected at these angular frequencies .
one may thus restrict the secondary anisotropies considered in the noise to the linear ( ostriker - vishniac ) and non - linear kinetic @xmath42 effects , as well as the thermal @xmath42 effect .
the @xmath42 effects are actually non - gaussian , spatially dependent , and the kinetic and thermal effects are correlated . as a simplifying assumption , we treat these two effects as two independent statistically isotropic gaussian noise components .
the effect of gravitational lensing is very small relative to the @xmath42 effects , but we still take it into account as a correction to the angular power spectrum of the primary anisotropies .    at arcminute resolution ,
the thermal and kinetic @xmath42 effects have standard deviations around @xmath43 and @xmath50 @xmath51 respectively .
they also have a slow decay at high angular frequencies and will dominate the string signal for string tension values below the current experimental upper bound .
arcminute @xmath0 experiments are in fact primarily dedicated to the detection of these secondary anisotropies . unlike the other effects considered , which have the same black body spectrum as the primary anisotropies , the thermal @xmath42 effect on the @xmath0 temperature depends on the frequency of observation .
its amplitude decreases between the rayleigh - jeans limit ( null frequency ) and @xmath52 @xmath53 where it is expected to vanish , before increasing again at higher frequencies .
figure [ fig : spectra ] represents the angular power spectra as functions of the angular frequency @xmath24 , for a string signal with string tension @xmath41 , for the primary @xmath0 anisotropies and the correction due to gravitational lensing , for the ostriker - vishniac and the non - linear kinetic @xmath42 effect , and for the thermal @xmath42 effect in the rayleigh - jeans limit .
these spectra are explicitly borrowed from @xcite .    in this context , the performance of the reconstruction algorithm to be defined will be studied in the following limits . in a first noise condition , we will consider the secondary anisotropies as a statistically isotropic gaussian noise , with power spectrum given by the rayleigh - jeans limit , that is added to the primary anisotropies . in a second noise condition , we will also assume an observation frequency around @xmath52 @xmath53 taking advantage of the frequency dependence of the thermal @xmath42 effect , and include in the noise secondary anisotropies in absence of this effect .
this is equivalent to including only the kinetic @xmath42 effect and gravitational lensing in the secondary anisotropies .
instrumental noise is considered to be negligible and simply discarded .
these two different noise conditions are respectively denoted as @xmath54 ( secondary anisotropies with thermal @xmath42 effect ) and @xmath55 ( secondary anisotropies without thermal @xmath42 effect ) in the following .
analyzing these limits can reveal to what extent the kinetic and thermal @xmath42 effects hamper the reconstruction of the string signal , as a function of the string tension . in a third noise condition
, the reconstruction performance will also be examined in the limit where the noise only includes primary anisotropies , neglecting even instrumental noise , and assuming secondary anisotropies have been correctly separated .
this case is denoted as @xmath56 ( primary anisotropies without instrumental noise ) and will be studied in order to understand the behaviour of the reconstruction algorithm in ideal conditions .    for the sake of our analyses , foreground emissions such as galactic dust or point
sources @xcite are disregarded .
simulated temperature maps ( in @xmath57 ) of the string signal in the @xmath0 and of its magnitude of gradient for a tension @xmath58 are represented in the top and bottom left panels respectively .
temperature maps of the primary @xmath0 anisotropies , i.e. the noise in the @xmath56 condition , and of its magnitude of gradient are represented in the top and bottom right panels respectively.,title="fig:",width=151 ]   simulated temperature maps ( in @xmath57 ) of the string signal in the @xmath0 and of its magnitude of gradient for a tension @xmath58 are represented in the top and bottom left panels respectively .
temperature maps of the primary @xmath0 anisotropies , i.e. the noise in the @xmath56 condition , and of its magnitude of gradient are represented in the top and bottom right panels respectively.,title="fig:",width=151 ] +    simulated temperature maps ( in @xmath57 ) of the string signal in the @xmath0 and of its magnitude of gradient for a tension @xmath58 are represented in the top and bottom left panels respectively .
temperature maps of the primary @xmath0 anisotropies , i.e. the noise in the @xmath56 condition , and of its magnitude of gradient are represented in the top and bottom right panels respectively.,title="fig:",width=151 ]   simulated temperature maps ( in @xmath57 ) of the string signal in the @xmath0 and of its magnitude of gradient for a tension @xmath58 are represented in the top and bottom left panels respectively .
temperature maps of the primary @xmath0 anisotropies , i.e. the noise in the @xmath56 condition , and of its magnitude of gradient are represented in the top and bottom right panels respectively.,title="fig:",width=151 ]    we use simulations of the string signal borrowed from the set of realistic simulations produced by @xcite .
the specific simulations considered are defined on planar patches of size @xmath26 for a field of view defined by an angular opening @xmath59 @xmath7 .
the finite size of the patch induces a discretization of the spatial frequencies which may be defined on a grid with @xmath60 uniformly sampled points @xmath61 with @xmath62 , below some band limit @xmath30 .
the original maps are sampled at nyquist - shannon rate on grids with @xmath60 uniformly sampled points @xmath63 with @xmath62 .
we consider here samples with @xmath64 , so that the corresponding pixels have an angular size around @xmath65 , with a band limit @xmath66 on each each component of the spatial frequencies .    for each noise component , a simulation
may easily be produced by taking the fourier transform of gaussian white noise , renormalizing each fourier coefficient by the square root of the corresponding power spectrum , and inverting the fourier transform @xcite . in each noise condition
considered , an overall noise simulation is obtained by simple superposition of the required independent components simulated .
the power spectrum of the noise @xmath67 is the sum of the individual spectra .
our approach for signal reconstruction notably relies on the introduction of specific statistical prior information on the signal in the regularization of the ill - posed inverse problem .
a number of @xmath68 simulations of the string signal @xmath8 are used as training data for fitting a prior statistical model for the gradient of the signal .
only @xmath6 additional simulation is reserved for testing the algorithm . in all noise conditions considered and for each string tension ,
this test string signal simulation is combined with @xmath69 independent realizations of the noise in order to produce multiple test simulations . for illustration ,
figure [ fig : maps - gradients ] represents the test simulation of the string signal and one simulation of the noise , as well as corresponding maps of the magnitude of gradient .
note that due to the finite size of the simulations , different string signal simulations exhibit slightly different estimated power spectra .
for simplicity of our argument though both training and test simulations have been normalized so as to exhibit the same total power for the same string tension ( see @xcite for more details ) .
in this section , we describe the ill - posed inverse problem posed for reconstruction of the string signal in the context of radio - interferometric observations of the @xmath0 . we also describe a modified inverse problem including a whitening operation on the measured visibilities and discuss basic ideas for signal reconstruction .
a map of the primary beam considered is represented in the left panel , with darkness increasing with the value of the function @xmath70 in the range [ @xmath71 .
a map of the mask made up of a random selection of @xmath72 per cent of the spatial frequencies up to the band limit in the fourier plane is represented in the right panel , with a unit value at the position of the selected frequencies identified by black points , and zero otherwise.,title="fig:",width=151 ]   a map of the primary beam considered is represented in the left panel , with darkness increasing with the value of the function @xmath70 in the range [ @xmath71 .
a map of the mask made up of a random selection of @xmath72 per cent of the spatial frequencies up to the band limit in the fourier plane is represented in the right panel , with a unit value at the position of the selected frequencies identified by black points , and zero otherwise.,title="fig:",width=151 ]    aperture synthesis in radio interferometry is a powerful technique in radio astronomy dating back to more than sixty years ago @xcite . in this context
, the small portion of the celestial sphere around the pointing direction tracked by a radio telescope array during observation defines the original real planar signal or image probed @xmath8 .
the field of view observed is limited by a primary beam represented by a function @xmath70 with values in the range [ @xmath71 .
each telescope pair at one instant of observation identifies a baseline defined as the relative position between the two telescopes . to each baseline
is associated one measurement called visibility . in the simplest setting one also considers baselines with negligible component in the pointing direction of the instrument . under this additional assumption ,
if the signal is made up of incoherent sources , each visibility corresponds to the value of the fourier transform of the signal multiplied by the primary beam at a spatial frequency @xmath28 identified by the components of the baseline in the plane of the signal .
we assume that the spatial frequencies @xmath28 probed by all telescope pairs during the observation belong to the discrete grid of points @xmath73 ( see @xcite for more details ) .
the fourier coverage provided by the @xmath74 spatial frequencies probed @xmath75 , with @xmath76 , can simply be identified by a binary mask in the fourier plane equal to @xmath6 for each spatial frequency probed and @xmath77 otherwise .
the visibilities measured may be denoted by a vector of @xmath74 complex fourier coefficients @xmath78 .
formally , the visibilities may equivalently be denoted by a vector of @xmath79 real measures @xmath80 consisting of the real and imaginary parts of the complex measures . in this
discrete setting , the fourier coverage is in general incomplete in the sense that the number of real constraints @xmath79 is smaller than the number of unknowns @xmath12 : @xmath81 .
an ill - posed inverse problem is thus defined for the reconstruction of the signal , seen as a vector of real sampled values @xmath82 , from the measured visibilities @xmath83 . in the context of our analysis ,
the noise affecting the signal consists in the vector of real sampled values for the standard @xmath0 components @xmath84 .
the inverse problem therefore reads as : @xmath85 where the measurement matrix @xmath86 identifies the complete linear relation between the signal and the visibilities . the matrix @xmath87 is the diagonal matrix implementing the primary beam .
the unitary matrix @xmath88 implements the discrete fourier transform providing the fourier coefficients .
the matrix @xmath89 is the rectangular binary matrix implementing the mask characterizing the interferometer .
it contains only one non - zero value on each line , at the index of the fourier coefficient corresponding to each of the spatial frequencies probed @xmath75 . neglecting small correlations introduced between spatial frequencies by the primary beam ,
the power spectrum of the noise @xmath90 is nonetheless modified as @xmath91 .
the noise @xmath92 affecting the visibilities is gaussian noise with different variances @xmath93 for different components @xmath94 .
let us recall that @xmath90 is real so that the components of @xmath95 are independent only when assuming that the mask @xmath96 contains no pair of antipodal spatial frequencies @xmath97 .
for each such pair , the two corresponding noise components are complex conjugate of one another . in this work
@xmath74 denotes the number of independent complex visibilities , while @xmath98 , with @xmath99 , will stand for the total number of complex visibilities including possible pairs of antipodal spatial frequencies .
the measurement matrix may be augmented by a whitening operation which will be useful in a reconstruction perspective .
it consists in dividing each measured visibility @xmath100 by the standard deviation of the real and imaginary parts of the corresponding noise component @xmath94 .
this operation is simply implemented in terms of the diagonal matrix @xmath101^{1/2}\delta_{bb'}\}_{1\leq b , b'\leq ( m/2)}$ ] .
this gives rise to the modified inverse problem @xmath102 where the whitened visibilities @xmath103 are now affected by identically distributed independent ( i.i.d . )
gaussian noise with unit variance of the real and imaginary part of each component @xmath94 , i.e. white gaussian noise .    for each test simulation ,
the visibilities are computed in accordance with relation ( [ eq : iip1 ] ) .
the primary beam @xmath70 considered for our observational set up is a gaussian function with a full width at half maximum ( @xmath104 ) equal to the angular opening of the field of view @xmath105 .
the function has a maximum value equal to unity .
two fourier coverages are considered independently , defined by masks made up of a random selection of @xmath106 spatial frequencies ( see @xcite for more details ) with either @xmath72 per cent or @xmath107 per cent of the fourier plane probed up to the band limit @xmath30 .
note that , after discarding redundant visibilities associated with pairs of antipodal spatial frequencies in the two masks considered , these two percentages respectively correspond to selections of @xmath79 spatial frequencies with @xmath108 and @xmath109 per cent .
the primary beam as well as the mask for a fourier coverage of @xmath72 per cent are illustrated in figure [ fig : obssetup ] .
the visibilities are also whitened in accordance with relation ( [ eq : iip1 ] ) by the introduction of the suitable matrix @xmath110 in each noise condition .      in the perspective of the reconstruction of the signal @xmath111 , relation ( [ eq : iip1 ] ) represents the measurement constraint . considering a candidate reconstruction @xmath112 ,
the residual noise reads as @xmath113 .
the residual noise level estimator , defined as twice the negative logarithm of the likelihood associated with @xmath112 , simply reads as @xmath114 the notation @xmath115 for a scalar stands for the complex modulus when applied to a complex number and for the absolute value when applied to a real number .
the @xmath116 norm of the residual noise is the standard norm of the corresponding vector : @xmath117 .
this noise level estimator follows a chi - square distribution with @xmath79 degrees of freedom .
typically , this estimator should be minimized by the good candidate reconstruction .
the measurement constraint on the reconstruction may be defined as a bound @xmath118 , with @xmath119 corresponding to some suitable large percentile of the @xmath120 distribution in order to avoid noise over - fitting .    in this context , many signals may formally satisfy the measurement constraint . a regularization scheme that encompasses
enough prior information on the original signal is needed in order to find a unique solution .
all image reconstruction algorithms will differ through the kind of regularization considered .
firstly , the simplest possible reconstruction procedure would consist in assuming that the non - probed spatial frequencies have a null value , and simply inverting the fourier transform after dividing each visibility by the variance of the noise in order to enhance the signal - to - noise ratio .
this reconstruction can by written as @xmath121 , where the notation @xmath122 stands for the conjugate transpose operation .
this map simply identifies with the convolution of the original signal and noise @xmath123 by a beam @xmath124 , where @xmath125 stands for a spike at the center of the field of view : @xmath126 , with @xmath127 identifying the central pixel . in the standard vocabulary of radio interferometry , the beam
is called the dirty beam , and the convolved map is called the dirty map in recognition of its visual appearance contrasting with the original signal .
a normalization constant needs to be chosen so as to get a unit central value of the beam .
in other words , no denoising or deconvolution is performed by this simple procedure that we call @xmath128 .
the only important operation in this reconstruction is the enhancement of the signal - to - noise ratio by the matrix @xmath129 .
secondly , the standard @xmath2 algorithm is essentially a matching pursuit ( @xmath130 ) algorithm @xcite which performs the deconvolution by iterative removal of the dirty beam , thus leading from the dirty map to a clean map @xcite .
formally , the algorithm should be stopped when the @xmath120 constraint above is satisfied .
in this section , we describe the standard global minimization problem @xmath131 for the reconstruction of the string signal in the context of compressed sensing .
we model the prior statistical distribution of the gradient of the string signal with a @xmath10 , fitted on the basis of the training simulations .
we also propose to enhance the @xmath131 problem on the basis of this statistical prior information , hence defining a @xmath132 minimization problem .
we finally propose a power spectral model ( psm ) used for independent estimation of the string tension , which may notably be useful for accelerating the convergence of the iterative algorithms solving the minimization problems .      for the reconstruction of a signal of sparsity @xmath133 , the theory of compressed sensing proves that the number @xmath79 of random measurements required , in a sensing basis incoherent with the sparsity basis , roughly scales with the sparsity , well below the number of samples @xmath12 of the signal at nyquist - shannon rate : @xmath134 @xcite .
in the case of radio - interferometric data , the measurements are assumed to be random fourier samples which represent a good sensing procedure for a string signal @xmath8 exhibiting sparse magnitude of the gradient @xmath11 . in this case , in the presence of i.i.d .
gaussian noise , the so - called total variation minimization problem ( @xmath131 ) applies @xcite , which corresponds to the minimization of the @xmath135 norm of the real signal under a @xmath136 constraint on the candidate reconstruction @xmath112 : @xmath137 the vector of sample values of the magnitude of the gradient associated with the vector @xmath111 is denoted as @xmath138 . by definition , the @xmath135 norm of a signal is the @xmath139 norm of the magnitude of its gradient , i.e. @xmath140 .
we set @xmath119 to be the @xmath141 percentile of the @xmath120 distribution .
the iterative algorithm considered for solving this convex optimization problem is based on the douglas - rachford splitting method @xcite in the framework of proximal operator theory @xcite .
[ c ] [ c ]   positive tail of modeled prior @xmath10 distribution @xmath142 for a string tension @xmath143 .
the @xmath10 is superimposed on the histogram of the aggregate of the magnitude of the gradient of the @xmath68 test simulations.,title="fig:",width=226 ]    in a bayesian perspective , the minimization of the @xmath135 norm of the signal amounts to assuming that the signal gradient follows a laplacian distribution with no spatial correlations , and to minimizing the negative logarithm of this prior distribution , i.e. to maximizing the distribution , subject to the @xmath120 constraint .
this promotes compressibility in a generic sense .
a laplacian distribution is indeed more peaked and more heavy - tailed than a gaussian prior for example .
however , with realistic simulations of the string signal at hand one can study more precisely the signal distribution . in that perspective
, the signal gradient can be modeled at each sample point by the positive tail of a generalized gaussian distribution ( @xmath10 ) : @xmath144 where @xmath145 is the gamma function , and where @xmath146 and @xmath147 are respectively called shape and scale parameters .
the shape parameter @xmath146 can be considered as a continuous measure of the compressibility of the underlying distribution . setting @xmath148 recovers a gaussian distribution , and @xmath149 recovers a laplacian distribution
. letting @xmath146 approach @xmath77 yields compressible distributions .
this parameter determines the kurtosis @xmath150 , i.e. the ratio of the fourth central moment to the square of the variance ( second central moment ) , by @xmath151^{2}}.\label{eq : csr3}\ ] ] by nature this kurtosis is independent of the string tension @xmath20 .
the scale parameter @xmath147 is linearly proportional to the standard deviation @xmath152 of the distribution .
the corresponding variance reads as : @xmath153^{2}= \rho^{2}\frac{\gamma\left(3q^{-1}\right)}{\gamma\left(q^{-1}\right)}s^{2}.\label{eq : csr4}\ ] ] the shape and scale parameters are estimated by a moment method on the aggregate of the magnitude of the gradient of the @xmath68 test simulations .
the values with two significant figures are @xmath154 and @xmath155 .
the shape parameter value shows that the signal gradient exhibits a much more compressible gradient than a generic laplacian distribution .
figure [ fig : histogram ] shows the positive tail of the modeled prior @xmath10 distribution @xmath142 for a string tension @xmath143 .
the @xmath10 is superimposed on the histogram of the aggregate of the magnitude of the gradient of the @xmath68 test simulations .
qualitatively , we see that the prior distribution is well modeled by a @xmath10 .
accounting for the specific statistical prior information in the regularization of the inverse problem should enhance the reconstruction performance , beyond the standard @xmath131 minimization . in the bayesian perspective
discussed , we define a statistical total variation @xmath132 minimization problem that consists in minimizing the logarithm of the @xmath10 prior distribution modeled , subject to the @xmath120 constraint : @xmath156 again , we set @xmath119 to be the @xmath141 percentile of the @xmath120 distribution . because @xmath157 the problem is non - convex in nature , but a re - weighted scheme allows convergence of the algorithm through iterations of simple @xmath131 minimizations @xcite .    note that no prior knowledge of the string tension @xmath20 is formally required for the reconstruction .
however we must acknowledge the fact that convergence appears to be faster when some prior estimation of @xmath20 is used to renormalize @xmath112 in the minimization term of relations ( [ eq : csr1 ] ) and ( [ eq : csr5 ] ) , as expected from relation ( [ eq : csr2 ] ) . as discussed in section [ sec : analyses and results ] , such an estimation of @xmath20 will also define a threshold below which one should not expect good reconstruction of the string network .
our estimation procedure is described below .      by bayes theorem ,
the posterior probability distribution for @xmath20 given the measured visibilities @xmath158 is simply obtained from the likelihood @xmath159 and the prior probability distribution @xmath160 on @xmath20 . for complete consistency
, the likelihood should be calculated using the model established for the signal gradient and noise . however , while this model by construction accounts for the non - gaussianity , i.e. sparsity of the string signal , it ignores the spatial correlations .    as in @xcite , we have observed that a likelihood yielding a precise localization of the string tension value can actually be obtained using a power spectral model .
such a model assumes that both the string signal and the noise arise from statistically isotropic gaussian random processes , such that their fourier coefficients are independent gaussian variables . under this model ,
as the string signal and noise are independent , the observed signal @xmath123 has a power spectrum , as modified by the primary beam , of the form : @xmath161 with @xmath162 defined in section [ sub : reconstruction basics ] and similarly @xmath163 . in this setting the likelihood can be computed easily in terms of the measured visibilities @xmath83 as : @xmath164 the posterior probability distribution for @xmath20 given the measured visibilities thus reads as @xmath165 with normalization @xmath166 .
we take the prior @xmath160 to be flat in an interval @xmath167 $ ] , with an upper bound @xmath168 large relative to the upper bound ( [ eq : ssn4 ] ) .
we define an estimation @xmath169 of the string tension as the expectation value of this posterior probability distribution : @xmath170.\label{eq : csr9}\ ] ]
in this section , we firstly describe our numerical analyses , including the performance criteria of reconstruction , and discuss their results .
we also compare the corresponding eye visibility thresholds on the string tension with the detectability threshold obtained on the basis of the psm .
we finally propose a small discussion of our results .
the overall reconstruction is effective for mapping the string network if the magnitude of gradient of the reconstructed signal closely resembles the magnitude of gradient of the true string signal . a simple qualitative measure of the performance is given by whether the string network is visible in the magnitude of gradient of the reconstructed signal .
we define the eye visibility threshold as the minimum string tension around which the overall denoising and mapping by the magnitude of the gradient begins to exhibit string features visible by eye .
we will augment this qualitative assessment of the reconstruction performance with three quantitative measures , namely the signal - to - noise ratio , the correlation coefficient and the kurtosis of the magnitude of the gradient of the reconstructed string signal re - multiplied by the primary beam .
the kurtosis is known to be a good statistic for discriminating between models with and without cosmic strings @xcite .
the signal - to - noise ratio is defined in terms of the original and reconstructed string signal gradients as : @xmath171 where @xmath172 and @xmath173 respectively stand for the standard deviations of the discrepancy signal @xmath174 and of the original signal @xmath175 .
the standard deviations are estimated from the sample variances on the basis of the signal realizations concerned . with this definition
, the @xmath176 is measured in decibels ( @xmath177 ) .
we will consider that the reconstruction is meaningful in terms of signal - to - noise ratio for the values of @xmath20 where this statistic is positive .
the correlation coefficient is also defined in terms of the original and reconstructed string signal gradients as : @xmath178 where @xmath179 stands for the covariance between @xmath175 and @xmath180 , and @xmath173 and @xmath181 respectively stand for the standard deviations of the signals @xmath175 and @xmath180 .
the covariance and variances are again estimated from the sample covariance and variances on the basis of the signal realizations concerned .
we will consider that the reconstruction is meaningful in terms of correlation coefficient for the values of @xmath20 where this statistic is also positive .
note that the variances and covariances discussed obviously depend on the string tension , but not their dimensionless ratios . for this reason
we have omitted this dependence in the definitions ( [ eq : r1 ] ) and ( [ eq : r2 ] ) .
analogously , kurtoses of the magnitude of the gradient of the reconstructed string signal re - multiplied by the primary beam are estimated from the sample kurtoses on the basis of the signal realizations concerned .
the kurtosis of the magnitude of gradient of a string signal as expected on the basis of our test string signal simulation takes the value @xmath182 with two significant numbers .
we compare the reconstruction results of the simple @xmath128 procedure and of the standard @xmath183 algorithm with those of the @xmath131 and @xmath132 algorithms .
reconstructions are performed for string tensions equi - spaced in logarithmic scaling in the range @xmath184 $ ] , corresponding to ratio values for @xmath20 of @xmath185 , @xmath186 , @xmath187 , @xmath188 , and @xmath189 in each order of magnitude . for each string tension , each algorithm , each noise condition ( @xmath56 , @xmath55 , or @xmath54 ) and each coverage condition ( @xmath72 or @xmath107 per cent ) analysed , we consider that the quantitative measures described above indicate effective performance when this performance is significant over the @xmath69 corresponding simulations .
.[tab : detectionthresholds ] psm detectability thresholds and eye visibility thresholds on the string tension , for each of the noise and coverage conditions considered .
all values are given with two significant figures . [ cols="^,^,^,^",options="header " , ]      magnitude of the gradient of reconstructed string signals re - multiplied by the primary beam , for a tension @xmath190 , in the noise condition @xmath56 and for a coverage at @xmath72 per cent .
the reconstructions considered are those associated with the simple @xmath128 procedure ( top left panel ) , @xmath183 ( top right panel ) , @xmath131 ( bottom left panel ) , and @xmath132 ( bottom right panel).,title="fig:",width=151 ]   magnitude of the gradient of reconstructed string signals re - multiplied by the primary beam , for a tension @xmath190 , in the noise condition @xmath56 and for a coverage at @xmath72 per cent .
the reconstructions considered are those associated with the simple @xmath128 procedure ( top left panel ) , @xmath183 ( top right panel ) , @xmath131 ( bottom left panel ) , and @xmath132 ( bottom right panel).,title="fig:",width=151 ] +   magnitude of the gradient of reconstructed string signals re - multiplied by the primary beam , for a tension @xmath190 , in the noise condition @xmath56 and for a coverage at @xmath72 per cent .
the reconstructions considered are those associated with the simple @xmath128 procedure ( top left panel ) , @xmath183 ( top right panel ) , @xmath131 ( bottom left panel ) , and @xmath132 ( bottom right panel).,title="fig:",width=151 ]   magnitude of the gradient of reconstructed string signals re - multiplied by the primary beam , for a tension @xmath190 , in the noise condition @xmath56 and for a coverage at @xmath72 per cent .
the reconstructions considered are those associated with the simple @xmath128 procedure ( top left panel ) , @xmath183 ( top right panel ) , @xmath131 ( bottom left panel ) , and @xmath132 ( bottom right panel).,title="fig:",width=151 ]    [ b ] [ b ] [ b ] [ ] [ b ] [ b ] [ b ]      +      +      +    in all noise and coverage conditions , the simple @xmath128 procedure is already effective at enhancing the visibility of strings by eye thanks to the multiplication by @xmath129 in fourier space , but the resulting maps are dirty .
the string tensions representing the eye visibility thresholds are listed in table [ tab : detectionthresholds ] .
these thresholds are well below the experimental upper bound in the noise condition @xmath56 , while for the noise conditions @xmath55 and @xmath54 they are respectively of the same order and larger .
this evolution is naturally due to the larger amount of noise at high spatial frequencies , where the signal gradient lives , in the presence of secondary anisotropies .
the eye visibility thresholds for @xmath183 , @xmath131 , and @xmath132 are the same as for @xmath128 .
however it clearly appears that these more evolved algorithms offer much better reconstructions .
this can already be acknowledged by simple eye inspection , as illustrated in figure [ fig : pa - in - rec ] for the magnitude of the gradient of the reconstructed signals re - multiplied by the primary beam in the noise condition @xmath56 at @xmath72 per cent and for a string tension @xmath191 . this qualitative assessment is confirmed by our quantitative measures for the reconstruction performance in all noise and coverage conditions and for all string tensions .
figure [ fig : stat ] represents the @xmath176 , correlation coefficient , and kurtosis of the magnitude of the gradient of reconstructed string signals re - multiplied by the primary beam as functions of the string tension in all noise conditions and for a coverage at @xmath72 per cent , for @xmath183 , @xmath131 , and @xmath132 . the corresponding curves for @xmath128 are actually not shown for readability reasons as they are very far from optimal values . in this context
, the major interest of the analysis actually resides in the comparison of the standard @xmath183 algorithm with the @xmath131 and @xmath132 algorithms .
firstly , the reconstructions illustrated in figure [ fig : pa - in - rec ] seem very similar for @xmath183 , @xmath131 , and @xmath132 .
note however that the magnitude of the gradient of the original signal ( not shown on the figure ) re - multiplied by the primary beam lies roughly in the range [ @xmath192 .
the range of amplitudes reported on the maps therefore indicate that the dynamic range of the reconstructions is better for @xmath131 than for @xmath183 , and even better for @xmath132 .
this is a generic behaviour that can be observed in all noise and coverage conditions and at all string tensions .
secondly , as shown in figure [ fig : stat ] in all noise conditions and for a coverage at @xmath72 per cent , the eye visibility thresholds identify with the string tension where the @xmath176 and correlation coefficient get down to very small values , and the kurtosis exhibits a large variability over the @xmath69 simulations around the expected value .
this allows us to identify rigorously the eye visibility threshold as the string tension below which the algorithms are not effective anymore at recovering the string network .
thirdly , as also shown in figure [ fig : stat ] in all noise conditions and for a coverage at @xmath72 per cent , the @xmath176 and correlation coefficient are always significantly larger for @xmath131 than for @xmath183 .
they are always at least as good for @xmath132 as for @xmath131 and often significantly larger .
however , it appears that the kurtosis is less optimally recovered with @xmath131 than with @xmath183 .
but the @xmath132 approach regularizes the kurtosis significantly better than @xmath183 in most cases .
note that these last two conclusions also hold for a @xmath107 per cent coverage , hence suggesting their stability and robustness relative to the coverage .
the reconstruction performance is always significantly better at @xmath107 per cent than at @xmath72 per cent , but the improvement is not drastic .
this suggests that a small number of visibilities can already lead to good reconstructions of the string signal .
in the perspective of the evaluation of the performance of the reconstruction algorithms , it is natural to ask whether the requirement for a complete reconstruction of the signal strongly limits the lowest string tension accessible in comparison with the lowest detectable string tension on the basis of the psm .    in this regard , for any possible string tension @xmath20
, the probability distribution for the estimate @xmath169 on the basis of the psm may be obtained from simulations .
one can identify the critical value @xmath193 such that , for null string tension , the probability that the estimate @xmath169 is larger than @xmath193 is small : @xmath194 , for some suitable positive value @xmath195 much smaller than unity .
the test for the hypothesis of null string tension is then defined as follows .
for estimated values @xmath196 , the hypothesis of null string tension may be rejected with a significance level @xmath195 . on the contrary for estimated values @xmath197 , the hypothesis of null string tension may not be rejected .
we define the psm detectability threshold @xmath198 such that , for the string tension @xmath198 , the probability that the estimate @xmath169 is larger than @xmath193 is large : @xmath199 , for some suitable positive value @xmath200 much smaller than unity .
consequently , for string tensions larger than @xmath198 , the probability of rejecting a null string tension on the basis of the hypothesis test defined is larger than @xmath201 , called the power of the test .
the value @xmath198 is the smallest string tension that can be discriminated from the hypothesis of null string tension for given values of significance and power of the test .
it may thus be understood as a detectability threshold determined on the basis of the psm .
the psm detectability thresholds in the various noise and coverage conditions considered are reported in table [ tab : detectionthresholds ] for @xmath202 . in the noise condition @xmath56 , detectability thresholds are well below the experimental upper bound , while for the noise conditions @xmath55 and @xmath54 they are respectively of the same order and larger . in all cases ,
the psm detectability threshold @xmath203 is only slightly below the eye visibility threshold associated with the reconstruction algorithms discussed , by a factor smaller than @xmath204 .
this discrepancy between these two thresholds confirms that the detection problem alone can be solved with the psm only down to slightly lower string tensions than the more difficult reconstruction problem .
as already emphasized in @xcite , it is one thing to estimate a single global parameter such as the string tension on the basis of a psm , but quite another to explicitly reconstruct the string network itself .
these comparison results are thus already very positive .
our results strongly support the idea of the superiority of the compressed sensing approaches for a string signal reconstruction from radio - interferometric data , in particular when specific statistical prior information is introduced in the global minimization problem defined for reconstruction .
however , at the high angular resolution considered and in the noise condition @xmath54 , the eye visibility thresholds , as well as the psm detectability thresholds , are above the current best experimental upper bound on the string tension .
one first approach to circumvent this difficulty in the perspective of the reconstruction of a possible string network from forthcoming radio - interferometric observations is to consider an observation frequency around @xmath52 @xmath53 so as to eliminate the thermal @xmath42 effect and to fall in the noise condition @xmath55 .
but even in this case the eye visibility thresholds remain slightly higher than the experimental upper bound .    in this context , recent work by the authors might help to drastically enhance the signal reconstruction performance in the context of compressed sensing approaches @xcite .
radio interferometers with small field of view and baselines with non - negligible and constant component in the pointing direction should be considered . in this context
, the visibilities measured essentially identify with a noisy and incomplete fourier coverage of the product of the planar signal multiplied by the primary beam with a linear chirp modulation .
theoretical and numerical results show the universality of the corresponding spread spectrum phenomenon relative to the sparsity basis , in terms of the achievable quality of reconstruction through global minimization problems such as @xmath131 and @xmath132 .
preliminary numerical results regarding the application of such a framework to the reconstruction of a string signal in the @xmath0 look very promising .    in another line of thoughts ,
global component separation techniques might be envisaged in order to simultaneously extract all non - gaussian components of the @xmath0 temperature data , including the @xmath42 effects and the string signal .
finally , let us recall that the string signal simulations used for modelling the prior statistical distribution of the signal gradient as @xmath10 s and for computing its power spectrum are realistic .
nevertheless , further analysis should be performed in order to study the stability of the @xmath132 reconstruction quality relative to a discrepancy between the statistical distribution of the true signal and our model . analyzing the difference between the @xmath131 and @xmath132 reconstruction qualities
already provides useful information in that respect , as the @xmath131 minimization problem captures the essence of the compressiblity of the signal gradient , but with a non - exact @xmath10 shape parameter @xmath149 . in this context
, one can argue in favour of a relative stability of the reconstruction quality in terms of the @xmath176 and correlation coefficient as these measures are indeed smaller for @xmath131 than for @xmath132 , but still larger than for @xmath183 .
this stability is lost in terms of the kurtosis though .
the @xmath132 algorithm proposed for the reconstruction of a string signal from arcminute resolution radio - interferometric data of the @xmath0 is designed in the framework of compressed sensing .
it notably relies on a model of the prior statistical distribution of the signal , fitted on the basis of realistic simulations .
the algorithm shows superior performance relative to the standard @xmath183 algorithm .
secondary @xmath0 anisotropies strongly hamper the signal reconstruction though , and further work is still required in order to obtain good performance at string tensions below the current experimental upper bound .
the authors wish to thank m. j. fadili for private communication of results on optimization by proximal methods .
the authors also thank a. a. fraisse , c. ringeval , d. n. spergel , and f. r. bouchet for kindly providing simulations of a string signal .
the authors also thank the reviewer for his valuable comments .
y. w. is postdoctoral researcher of the belgian national science foundation ( f.r.s .- fnrs ) .
jones m. e. , 2002 , in chen l .- w .
, ma c .- p . , ng k .- w . , pen u .-
, eds , asp conf . ser .
257 , amiba 2001 : high - z clusters , missing baryons , and @xmath0 polarization .
, san francisco , p. 35