A fast version of the k-means classification algorithm for astronomical applications

Context. K-means is a clustering algorithm that has been used to classify large datasets in astronomical databases. It is an unsupervised method, able to cope very different types of problems. Aims. We check whether a variant of the algorithm called single-pass k-means can be used as a fast alternative to the traditional k-means. Methods. The execution time of the two algorithms are compared when classifying subsets drawn from the SDSS-DR7 catalog of galaxy spectra. Results. Single-pass k-means turn out to be between 20 % and 40 % faster than k-means and provide statistically equivalent classifications. This conclusion can be scaled up to other larger databases because the execution time of both algorithms increases linearly with the number of objects. Conclusions. Single-pass k-means can be safely used as a fast alternative to k-means

Interpretation of HINODE SOT/SP asymmetric Stokes profiles observed in quiet Sun network and internetwork

We present the first interpretation of the Stokes profile asymmetries measured in the FeI 630 nm lines by SOT/SP, in both quiet Sun internetwork (IN) and network regions. The inversion is carried out under the hypothesis of MISMA, where the unresolved structure is assumed to be optically thin. We analyze a 29.52"x31.70" subfield carefully selected to be representative of the properties of a 302"x162" quiet Sun field-of-view at disk center. The inversion code is able to reproduce the observed asymmetries in a very satisfactory way. The inversion code interprets 25% of inverted profiles as emerging from pixels in which both positive and negative polarities coexist. kG field strengths are found at the base of the photosphere in both network and IN; in the case of the latter, both kG fields and hG fields are admixed. When considering the magnetic properties at the mid photosphere most kG fields are gone, and the statistics is dominated by hG fields. We constrain the magnetic field of only 4.5% of the analyzed photosphere (and this percentage reduces to 1.3% when referred to all pixels, including those with low polarization not analyzed). The rest of the plasma is consistent with the presence of weak fields not contributing to the detected polarization signals. The average flux densities derived in the full subfield and in IN regions are higher than the ones derived from the same dataset by Milne-Eddington inversion. The existence of large asymmetries in SOT/SP polarization profiles is uncovered. These are not negligible in quiet Sun data. The MISMA inversion code reproduces them in a satisfactory way, and provides a statistical description of the magnetized IN and network which partly differs and complements the results obtained so far. From this it follows the importance of having a complete interpretation of the line profile shapes.Comment: 11 pages, 9 figures, 1 table - Accepted for publication on A&

Magnetic flux in the inter-network quiet Sun from comparison with numerical simulations

Khomenko et al. estimate the mean magnetic field strength of the quiet Sun to be 20 G. The figure is smaller than several existing estimates, and it comes from the comparison between observed Zeeman polarization signals and synthetic signals from numerical simulations of magneto-convection. The numerical simulations require an artificially large magnetic diffusivity, which smears out magnetic structures smaller than the grid scale. Assuming a turbulent cascade for the unresolved artificially smeared magnetic fields, we find that their unsigned magnetic flux is at least as important as that explicitly shown in the simulation. The unresolved fields do not produce Zeeman polarization but contribute to the unsigned flux.Since they are not considered by Khomenko et al., their mean magnetic field strength has to be regarded as a lower limit. This kind of bias is not specific of a particular numerical simulation or a spectral line. It is to be expected when observed quiet Sun Zeeman signals are compared with synthetic signals from simulations.Comment: Accepted A&A. 4 pages, 0 figure

Asymmetries of the Stokes V profiles observed by HINODE SOT/SP in the quiet Sun

We present the first classification of SOT/SP circular polarization measurements with the aim of highlighting exhaustively the whole variety of Stokes V shapes emerging from the quiet Sun. k-means is used to classify HINODE SOT/SP Stokes V profiles observed in the quiet Sun network and internetwork (IN). We analyze a 302 x 162 square arcsec field-of-view (FOV) which can be considered a complete sample of quiet Sun measurements performed at at the disk center with 0.32 arcsec angular resolution and 0.001 polarimetric sensitivity. Such a classification allows us to divide the whole dataset in classes, with each class represented by a cluster profile, i.e., the average of the profiles in the class. The set of 35 cluster profiles derived from the analysis completely characterizes SOT/SP quiet Sun measurements. The separation between network and IN profile shapes is evident - classes in the network are not present in the IN, and vice versa. Asymmetric profiles are approximatively 93 % of the total number of profiles. Among these, approximatively 34 % of the profiles are strongly asymmetric profiles, and they can be divided in three families: blue-lobe, red-lobe, and Q-like profiles. The blue-lobe profiles tend to be associated with upflows (granules), whereas the red-lobe and Q-like ones appear in downflows (intergranular lanes). Such profiles need to be interpreted considering model atmospheres different from a uniformly magnetized Milne-Eddington (ME) atmosphere, i.e., characterized by gradients and/or discontinuities in the magnetic field and velocity along the line-of-sight (LOS).Comment: 11 pages, 4 figures, accepted for publication in Astronomy and Astrophysic

Do galaxies form a spectroscopic sequence?

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society ©: 2011 RAS © The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.We identify a spectroscopic sequence of galaxies, analogous to the Hubble sequence of morphological types, based on the Automatic Spectroscopic K-means-based (ASK) classification. Considering galaxy spectra as multidimensional vectors, the majority of the spectral classes are distributed along a well-defined curve going from the earliest to the latest types, suggesting that the optical spectra of normal galaxies can be described in terms of a single affine parameter. Optically bright active galaxies, however, appear as an independent, roughly orthogonal branch that intersects the main sequence exactly at the transition between early and late typesThis work has been funded by projects AYA2007-67965-C03-03, AYA 2007-67965-C03-01, AYA2007-67752-C03-01 and CSD 2006-00070 (Spanish Ministry of Science, Technology and Innovation

Strength distribution of solar magnetic fields in photospheric quiet Sun regions

The magnetic topology of the solar photosphere in its quietest regions is hidden by the difficulties to disentangle magnetic flux through the resolution element from the field strength of unresolved structures. The observation of spectral lines with strong coupling with hyperfine structure, like the observed MnI line at 553.7 nm, allows such differentiation. The main aim is to analyse the distribution of field strengths in the network and intranetwork of the solar photosphere through inversion of the MnI line at 553.7 nm. An inversion code for the magnetic field using the Principal Component Analysis (PCA) has been developed. Statistical tests are run on the code to validate it. The code has to draw information from the small-amplitude spectral feature oppearing in the core of the Stokes V profile of the observed line for field strengths below a certain threshold, coinciding with lower limit of the Paschen-Back effect in the fine structure of the involved atomic levels. The inversion of the observed profiles, using the circular polarization (V) and the intensity (I), shows the presence of magnetic fields strengths in a range from 0 to 2 kG, with predominant weak strength values. Mixed regions with mean strength field values of 1130 and 435 Gauss are found associated with the network and intranetwork respectively. The MnI line at 553 nm probes the field strength distribution in the quiet sun and shows the predominance of weak, hectoGauss fields in the intranetwork, and strong, kiloGauss fields in the network. It also shows that both network and intranetwork are to be understood at our present spatial resolutions as field distributions of which we hint the mean properties.Comment: 10 pages, 6 figure

Automated Detection and Tracking of Solar Magnetic Bright Points

Magnetic Bright Points (MBPs) in the internetwork are among the smallest objects in the solar photosphere and appear bright against the ambient environment. An algorithm is presented that can be used for the automated detection of the MBPs in the spatial and temporal domains. The algorithm works by mapping the lanes through intensity thresholding. A compass search, combined with a study of the intensity gradient across the detected objects, allows the disentanglement of MBPs from bright pixels within the granules. Object growing is implemented to account for any pixels that might have been removed when mapping the lanes. The images are stabilized by locating long-lived objects that may have been missed due to variable light levels and seeing quality. Tests of the algorithm employing data taken with the Swedish Solar Telescope (SST), reveal that ~90% of MBPs within a 75"x 75" field of view are detected

Relationship between Hubble type and spectroscopic class in local galaxies

We compare the Hubble type and the spectroscopic class of the galaxies with spectra in SDSS/DR7. As it is long known, elliptical galaxies tend to be red whereas spiral galaxies tend to be blue, however, this relationship presents a large scatter, which we measure and quantify in detail. We compare the Automatic Spectroscopic K-means based classification (ASK) with most of the commonly used morphological classifications. All of them provide consistent results. Given a spectral class, the morphological type wavers with a standard deviation between 2 and 3 T types, and the same large dispersion characterizes the variability of spectral classes fixed the morphological type. The distributions of Hubble types given an ASK class are very skewed -- they present long tails that go to the late morphological types for the red galaxies, and to the early morphological types for the blue spectroscopic classes. The scatter is not produced by problems in the classification, and it remains when particular subsets are considered. A considerable fraction of the red galaxies are spirals (40--60 %), but they never present very late Hubble types (Sd or later). Even though red spectra are not associated with ellipticals, most ellipticals do have red spectra: 97 % of the ellipticals in the morphological catalog by Nair & Abraham, used here for reference, belong to ASK 0, 2 or 3. It contains only a 3 % of blue ellipticals. The galaxies in the green valley class (ASK~5) are mostly spirals, and the AGN class (ASK 6) presents a large scatter of Hubble types from E to Sd. From redshift 0.25 to now the galaxies redden from ASK 2 to ASK 0, as expected from the passive evolution of their stellar populations. Two of the ASK classes (1 and 4) gather edge-on spirals, a property of interest in studies requiring knowing the intrinsic shape of a galaxy (e.g., weak lensing calibration).Comment: Accepted for publication in ApJ. 16 pages. 12 Figs. 2 summary table

Solar Intranetwork Magnetic Elements: bipolar flux appearance

The current study aims to quantify characteristic features of bipolar flux appearance of solar intranetwork (IN) magnetic elements. To attack such a problem, we use the Narrow-band Filter Imager (NFI) magnetograms from the Solar Optical Telescope (SOT) on board \emph{Hinode}; these data are from quiet and an enhanced network areas. Cluster emergence of mixed polarities and IN ephemeral regions (ERs) are the most conspicuous forms of bipolar flux appearance within the network. Each of the clusters is characterized by a few well-developed ERs that are partially or fully co-aligned in magnetic axis orientation. On average, the sampled IN ERs have total maximum unsigned flux of several 10^{17} Mx, separation of 3-4 arcsec, and a lifetime of 10-15 minutes. The smallest IN ERs have a maximum unsigned flux of several 10^{16} Mx, separations less than 1 arcsec, and lifetimes as short as 5 minutes. Most IN ERs exhibit a rotation of their magnetic axis of more than 10 degrees during flux emergence. Peculiar flux appearance, e.g., bipole shrinkage followed by growth or the reverse, is not unusual. A few examples show repeated shrinkage-growth or growth-shrinkage, like magnetic floats in the dynamic photosphere. The observed bipolar behavior seems to carry rich information on magneto-convection in the sub-photospheric layer.Comment: 26 pages, 14 figure