Lacosmicx

I have finished the optimized LA Cosmic implementation. It uses Cython for the main part of the algorithm and uses C for the individual routines. Everything is multithreaded using OpenMP (the Cython standard). For a 4K x 2K image, this implementation is ~35 times faster than cosmics.py written by Malte Tewes (which is already significantly faster than the original IRAF version). 

Given the complexity of this code, should this be included in imageutils or should it be its own package? 

Any other comments are welcome.
Sure. I'm happy to do that. I will put in the full reference in the
documentation of the main function as well.

On Mon, Nov 24, 2014 at 11:49 AM, Erik Tollerud <notifications@github.com>
wrote:

> In imageutils/lacosmicx/__init__.py:
>
> > +Author : Curtis McCully
> > +Date : October 2014
> > +
> > +Lacosmicx is designed to detect cosmic rays in images (numpy arrays),
> > +based on Pieter van Dokkum's L.A.Cosmic algorithm.
> > +
> > +Much of this was originally adapted from cosmics.py written by Malte Tewes.
> > +I have ported all of the slow functions to Cython/C, and optimized
> > +where I can. This is designed to be as fast as possible so some of the
> > +readability has been sacrificed, specifically in the C code.
> > +
> > +L.A.Cosmic = LAplacian Cosmic ray detection
> > +
> > +U{http://www.astro.yale.edu/dokkum/lacosmic/}
> > +
> > +(article : U{http://arxiv.org/abs/astro-ph/0108003})
>
> Perhaps you should also link to the journal article (
> http://www.jstor.org/stable/10.1086/323894)? Or, if it's open access,
> maybe *instead* of the arXiv version? (I can't tell if it's open right
> now, because I'm on a network that has paid-for access...)
>
> —
> Reply to this email directly or view it on GitHub
> <https://github.com/astropy/imageutils/pull/35/files#r20816434>.
>



-- 
~Dr. Curtis McCully
Department of Physics and Astronomy
Rutgers, The State University of New Jersey