The ability to study cancer-immune cell communication across the whole tumor section without tissue dissociation is needed, especially for cancer immunotherapy development, which requires understanding of molecular mechanisms and discovery of more druggable targets. In this work, we assembled and evaluated an integrated experimental framework and analytical process to enable genome-wide scale discovery of ligand-receptors potentially used for cellular crosstalks, followed by targeted validation. We assessed the complementarity of four different technologies: single-cell RNA sequencing and Spatial transcriptomic (measuring over >20,000 genes), RNA In Situ Hybridization (RNAscope, measuring 4-12 genes) and Opal Polaris multiplex protein staining (4-9 proteins). To utilize the multimodal data, we implemented existing methods and also developed STRISH (Spatial TRanscriptomic In Situ Hybridization), a computational method that can automatically scan across the whole tissue section for local expression of gene (e.g. RNAscope data) and/or protein markers (e.g. Polaris data) to recapitulate an interaction landscape across the whole tissue. We evaluated the approach to discover and validate cell-cell interaction in situ through in-depth analysis of two types of cancer, basal cell carcinoma and squamous cell carcinoma, which account for over 70% of cancer cases. We showed that inference of cell-cell interactions using scRNA-seq data can misdetect or detect false positive interactions. Spatial transcriptomics still suffers from misdetecting lowly expressed ligand-receptor interactions, but reduces false discovery. RNAscope and Polaris are sensitive methods for defining the location of potential ligand receptor interactions, and the STRISH program can determine the probability that local gene co-expression reflects true cell-cell interaction. We expect that the approach described here will be widely applied to discover and validate ligand receptor interaction in different types of solid cancer tumors.