To further reveal the interplay between chemical substitution patterns and the resulting crystal packing types for azapentacenes in an automated way, a kernel-based machine learning technique was applied in conjunction with clustering and dimensionality reduction analysis to produce a sketch-map for the individual crystal packing landscapes for pentacene and two azapentacene molecules. These methods provided valuable insight into the profound changes in the overall crystal energy landscapes with changes in molecular substitution pattern.