The problem of segmenting an image into several modalities representing different textures can be modelled using Gaussian mixtures. Moreover, texture image patches when translated, rotated or scaled lie in low dimensional subspaces of the high-dimensional space spanned by the grey values. These two aspects make the mixture of local subspace models worth consideration for segmenting this type of images. In recent years a number of mixtures of local PCA models have been proposed. Most of these models require the user to set the number of subspaces and subspace dimensionalities. To make the model autonomous, we propose a greedy EM algorithm to find a suboptimal number of subspaces, besides using a global retained variance ratio to estimate for each subspace the dimensionality that retains the given variability ratio. We provide experimental results for testing the proposed method on texture segmentation.