We have previously used cryo-electron tomography combined with sub-volume averaging and classification to obtain 3D structures of macromolecular assemblies in cases where a single dominant species was present, and applied these methods to the analysis of a variety of trimeric HIV-1 and SIV envelope glycoproteins (Env). Here, we extend these studies by demonstrating automated, iterative, missing wedge-corrected 3D image alignment and classification methods to distinguish multiple conformations that are present simultaneously. We present a method for measuring the spatial distribution of the vector elements representing distinct conformational states of Env. We identify data processing strategies that allow clear separation of the previously characterized closed and open conformations, as well as unliganded and antibody-liganded states of Env when they are present in mixtures. We show that identifying and removing spikes with the lowest signal-to-noise ratios improves the overall accuracy of alignment between individual Env sub-volumes, and that alignment accuracy, in turn, determines the success of image classification in assessing conformational heterogeneity in heterogeneous mixtures. We validate these procedures for computational separation by successfully separating and reconstructing distinct 3D structures for unliganded and antibody-liganded as well as open and closed conformations of Env present simultaneously in mixtures.