3D electron diffraction (3D ED) has shown great potential in crystal structure determination in materials, small organic molecules, and macromolecules. In this work, an automated, low-dose and low-bias 3D ED protocol has been implemented to identify six phases from a multiple-phase melt-crystallisation product of an active pharmaceutical ingredient, griseofulvin (GSF). Batch data collection under low-dose conditions using a widely available commercial software was combined with automated data analysis to collect and process over 230 datasets in three days. Accurate unit cell parameters obtained from 3D ED data allowed direct phase identification of GSF Forms III, I and the known GSF inclusion complex (IC) with polyethylene glycol (PEG) (GSF-PEG IC-I), as well as three minor phases, namely GSF Forms II, V and an elusive new phase, GSF-PEG IC-II. Their structures were then directly determined by 3D ED. Furthermore, we reveal how the stabilities of the two GSF-PEG IC polymorphs are closely related to their crystal structures. These results demonstrate the power of automated 3D ED for accurate phase identification and direct structure determination of complex, beam-sensitive crystallisation products, which is significant for drug development where solid form screening is crucial for the overall efficacy of the drug product.