Lorentzian Iterative Hard Thresholding: Robust Compressed Sensing with Prior Information

Commonly employed reconstruction algorithms in compressed sensing (CS) use the $L_2$ norm as the metric for the residual error. However, it is well-known that least squares (LS) based estimators are highly sensitive to outliers present in the measurement vector leading to a poor performance when the noise no longer follows the Gaussian assumption but, instead, is better characterized by heavier-than-Gaussian tailed distributions. In this paper, we propose a robust iterative hard Thresholding (IHT) algorithm for reconstructing sparse signals in the presence of impulsive noise. To address this problem, we use a Lorentzian cost function instead of the $L_2$ cost function employed by the traditional IHT algorithm. We also modify the algorithm to incorporate prior signal information in the recovery process. Specifically, we study the case of CS with partially known support. The proposed algorithm is a fast method with computational load comparable to the LS based IHT, whilst having the advantage of robustness against heavy-tailed impulsive noise. Sufficient conditions for stability are studied and a reconstruction error bound is derived. We also derive sufficient conditions for stable sparse signal recovery with partially known support. Theoretical analysis shows that including prior support information relaxes the conditions for successful reconstruction. Simulation results demonstrate that the Lorentzian-based IHT algorithm significantly outperform commonly employed sparse reconstruction techniques in impulsive environments, while providing comparable performance in less demanding, light-tailed environments. Numerical results also demonstrate that the partially known support inclusion improves the performance of the proposed algorithm, thereby requiring fewer samples to yield an approximate reconstruction.Comment: 28 pages, 9 figures, accepted in IEEE Transactions on Signal Processin

PURIFY: a new algorithmic framework for next-generation radio-interferometric imaging

In recent works, compressed sensing (CS) and convex opti- mization techniques have been applied to radio-interferometric imaging showing the potential to outperform state-of-the-art imaging algorithms in the field. We review our latest contributions [1, 2, 3], which leverage the versatility of convex optimization to both handle realistic continuous visibilities and offer a highly parallelizable structure paving the way to significant acceleration of the reconstruction and high-dimensional data scalability. The new algorithmic structure promoted in a new software PURIFY (beta version) relies on the simultaneous-direction method of multipliers (SDMM). The performance of various sparsity priors is evaluated through simulations in the continuous visibility setting, confirming the superiority of our recent average sparsity approach SARA

On sparsity averaging

Recent developments in Carrillo et al. (2012) and Carrillo et al. (2013) introduced a novel regularization method for compressive imaging in the context of compressed sensing with coherent redundant dictionaries. The approach relies on the observation that natural images exhibit strong average sparsity over multiple coherent frames. The associated reconstruction algorithm, based on an analysis prior and a reweighted $\ell_1$ scheme, is dubbed Sparsity Averaging Reweighted Analysis (SARA). We review these advances and extend associated simulations establishing the superiority of SARA to regularization methods based on sparsity in a single frame, for a generic spread spectrum acquisition and for a Fourier acquisition of particular interest in radio astronomy.Comment: 4 pages, 3 figures, Proceedings of 10th International Conference on Sampling Theory and Applications (SampTA), Code available at https://github.com/basp-group/sopt, Full journal letter available at http://arxiv.org/abs/arXiv:1208.233

PURIFY: a new approach to radio-interferometric imaging

In a recent article series, the authors have promoted convex optimization algorithms for radio-interferometric imaging in the framework of compressed sensing, which leverages sparsity regularization priors for the associated inverse problem and defines a minimization problem for image reconstruction. This approach was shown, in theory and through simulations in a simple discrete visibility setting, to have the potential to outperform significantly CLEAN and its evolutions. In this work, we leverage the versatility of convex optimization in solving minimization problems to both handle realistic continuous visibilities and offer a highly parallelizable structure paving the way to significant acceleration of the reconstruction and high-dimensional data scalability. The new algorithmic structure promoted relies on the simultaneous-direction method of multipliers (SDMM), and contrasts with the current major-minor cycle structure of CLEAN and its evolutions, which in particular cannot handle the state-of-the-art minimization problems under consideration where neither the regularization term nor the data term are differentiable functions. We release a beta version of an SDMM-based imaging software written in C and dubbed PURIFY (http://basp-group.github.io/purify/) that handles various sparsity priors, including our recent average sparsity approach SARA. We evaluate the performance of different priors through simulations in the continuous visibility setting, confirming the superiority of SARA

Actinomicetos termoalcalófilos del área subtropical de Jujuy, Argentina

El objetivo de este trabajo fue examinar los actinomicetos termoalcalófilos presentes en el área subtropical de Jujuy, Argentina, caracterizada por el cultivo de la caña de azúcar. Se aislaron en medio con novobiocina las especies Laceyella putida, Laceyella sacchari, Thermoactinomyces intermedius, Thermoactinomyces vulgaris y Thermoflavimicrobium dichotomicum a partir de la rizósfera de plantas y de renuevos de caña de azúcar, así como de suelos de pastura y de monte natural. El pH de los suelos era casi neutro a ligeramente alcalino, excepto en un solo caso en que el suelo estaba acidificado por licor láctico. El número de actinomicetos encontrados sobre los tejidos vivos y en el bagazo recién obtenido o almacenado según el método de Ritter fue pequeño en comparación con el observado sobre las hojas en descomposición. L. sacchari predominó respecto de T. vulgaris. Se aislaron especies termoalcalófilas de Laceyella, Thermoactinomyces, Thermoflavimicrobium, Saccharomonospora, Streptomyces y Thermononospora de los residuos compostados de caña de azúcar utilizando medio sin novobiocina. El aire capturado cerca de pilas de bagazo en compostaje contenía esporos de estos organismos.The objective of this study was to examine the alkalithermophilic actinomycete communities in the subtropical environment of Jujuy, Argentina, characterized by sugarcane crops. Laceyella putida, Laceyella sacchari, Thermoactinomyces intermedius, Thermoactinomyces vulgaris and Thermoflavimicrobium dichotomicum were isolated on the media with novobiocin, from sugar cane plants and renewal rhizospheres, and grass and wood soils. Soil pH was almost neutral or lightly alkaline, except for grass soil acidified by lactic liquor. A smaller number of actinomycetes was found on the living plants and bagasse (recently obtained or stored according to the Ritter method) with respect to decomposed leaves on the soil. Thermophilic species of Laceyella, Thermoactinomyces, Thermoflavimicrobium, Saccharomonospora, Streptomyces and Thermononospora were isolated on the media without novobiocin, from composted sugar cane residues. Air captured near composted bagasse piles, contained alkalithermophilic actinomycete spores.Fil: Carrillo, L.. Universidad Nacional de Jujuy. Facultad de Ciencias Agrarias; ArgentinaFil: Benitez Ahrendts, Marcelo Rafael. Universidad Nacional de Jujuy. Facultad de Ciencias Agrarias; ArgentinaFil: Maldonado, Marcos Javier. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tucumán. Planta Piloto de Procesos Industriales Microbiológicos; Argentina. Universidad Nacional de Jujuy. Facultad de Ciencias Agrarias; Argentin

A randomised primal-dual algorithm for distributed radio-interferometric imaging

Next generation radio telescopes, like the Square Kilometre Array, will acquire an unprecedented amount of data for radio astronomy. The development of fast, parallelisable or distributed algorithms for handling such large-scale data sets is of prime importance. Motivated by this, we investigate herein a convex optimisation algorithmic structure, based on primal-dual forward-backward iterations, for solving the radio interferometric imaging problem. It can encompass any convex prior of interest. It allows for the distributed processing of the measured data and introduces further flexibility by employing a probabilistic approach for the selection of the data blocks used at a given iteration. We study the reconstruction performance with respect to the data distribution and we propose the use of nonuniform probabilities for the randomised updates. Our simulations show the feasibility of the randomisation given a limited computing infrastructure as well as important computational advantages when compared to state-of-the-art algorithmic structures.Comment: 5 pages, 3 figures, Proceedings of the European Signal Processing Conference (EUSIPCO) 2016, Related journal publication available at https://arxiv.org/abs/1601.0402

The Acoustics of the Choir in Spanish Cathedrals

One of the most significant enclosures in worship spaces is that of the choir. Generally, from a historical point of view, the choir is a semi-enclosed and privileged area reserved for the clergy, whose position and configuration gives it a private character. Regarding the generation and transformation of ecclesial interior spaces, the choir commands a role of the first magnitude. Its shape and location produce, on occasions, major modifications that significantly affect the acoustics of these indoor spaces. In the case of Spanish cathedrals, whose design responds to the so-called “Spanish type”, the central position of the choir, enclosed by high stonework walls on three of its sides and with numerous wooden stalls inside, breaks up the space in the main nave, thereby generating other new spaces, such as the trascoro. The aim of this work was to analyse the acoustic evolution of the choir as one of the main elements that configure the sound space of Spanish cathedrals. By means of in situ measurements and simulation models, the main acoustic parameters were evaluated, both in their current state and in their original configurations that have since disappeared. This analysis enabled the various acoustic conditions existing between the choir itself and the area of the faithful to be verified, and the significant improvement of the acoustic quality in the choir space to become apparent. The effect on the acoustic parameters is highly significant, with slight differences in the choir, where the values are appropriate for Gregorian chants, and suitable intelligibility of sung text. High values are also obtained in the area of the faithful, which lacked specific acoustic requirements at the time of construction

Uso de tubos timpánicos auto elaborados en pacientes con otitis media serosa realizado en el servicio de otorrinolaringología del hospital Antonio Lenin Fonseca durante el período de Mayo a Diciembre del 2014

Se realizó el estudio de “Uso de tubos timpánicos autoelaborados en pacientes con otitis media serosa en el servicio de otorrinolaringología del Hospital Antonio Lenin Fonseca durante el período de mayo a diciembre del 2014.” El objetivo era conocer los resultados que presentaron los pacientes al colocárseles tubos timpánicos autoelaborados en el servicio de Otorrinolaringología del Hospital Antonio Lenin Fonseca en el período de mayo a diciembre del 2014. Es un estudio descriptivo, longitudinal, prospectivo y observacional, estudio de una cohorte, realizado en el servicio de otorrinolaringología del Hospital Antonio Lenín Fonseca de Managua en el 2014, cuyo universo está constituido por pacientes diagnosticados con Otitis Media serosa que requirieron de cirugía en el período de estudio, tomando de muestra a todos los pacientes a los que se les colocó tubos timpánicos durante el período en estudio(N = 10). Se aplicó un muestreo por conveniencia. En los pacientes predominó la edad de 5-10 años (50%), de escolaridad primaria (60%), sin diferencias en cuanto al sexo y la procedencia. En los 17 oídos estudiados se encontró membrana timpánica retraída (100%), la cual estaba levemente enrojecida (47%). El timpanograma nos graficó una curva tipo B en cada oído. La complicación más frecuente en los oídos derechos fue tejido de granulación (40%) así también en los oídos izquierdos (55.6%). Estos tubos fueron expulsados antes de los 3 meses (60%) Dentro de los datos y resultados obtenidos, podemos observar la producción de tejido de granulación en la membrana timpánica