Breast cancer (BC) remains the deadliest cancer for women worldwide. Neoadjuvant immunotherapies have demonstrated improved responses for some patients. Unfortunately, no robust method exists for predicting which patients will respond to immunotherapy. Imaging of diagnostic BC biopsies has revealed that the spatial distribution of tumor infiltrating lymphocytes (TILs) and other immune cells within and around the tumor can help stratify BC patients into responders and non-responders. However, clinical microscopy cannot differentiate between subtypes of TILs; numerous markers are needed to capture the heterogeneity of cancer cells and immune cells in the TME. Highly multiplexed fluorescence microscopy, or high-plex IF, has emerged as a workhorse in data collection for spatial proteomics. We present a pilot study of the TME of BC patients treated with neoadjuvant immunotherapy. Specifically in this abstract, we discuss computer vision methods for analyzing the cellular constituents probed in these complex and rich images. We discuss image stitching and channel registration for high-plex modalities, deep learning algorithms for cell detection and segmentation, and pseudo-spectral angle mapping (pSAM) for cell classification. We present strategies for accurate quantification of these images, facilitating investigations into immune activity in breast tumors with high phenotypic accuracy.