Non-negative Matrix Factorization (NMF) is a robust approach to learning spatially localized parts-based subspace patterns in applications such as document analysis, image interpretation, and gene expression analysis. NMF-based decomposition capabilities are lacking in the present ITK toolkit. We provide a generic framework for such decompositions. A specific implementation using a Kulback-Liebler type divergence function is provided to illustrate a possible extension of the base class along with test images to illustrate usage. We have found this method to be robust to noisy image data and show examples from our on-going research using the Allan Brain Atlas data to illustrate its ability to analyze higher dimension data.