Characterisation of structure across the nanometre scale is key to bridging the gap between the local atomic environment and macro-scale and can be achieved by means of scanning electron nanobeam diffraction (SEND). As a technique, SEND allows for a broad range of samples, due to being relatively tolerant of specimen thickness with low electron dosage. This, coupled with the capacity for automation of data collection over wide areas, allows for statistically representative probing of the microstructure. This paper outlines a versatile, data-driven approach for producing domain maps, and a statistical approach for assessing their applicability. The workflow utilises a Variational AutoEncoder to identify the sources of variance in the diffraction signal, and this, in combination with clustering techniques, is used to produce domain maps. This approach is agnostic to domain crystallinity, requires no prior knowledge of crystal structure, and does not require simulation of a library of expected diffraction patterns.