Non-small cell lung cancer (NSCLC) is the most widely distributed tumor in the world, and its immunotherapy is not practical. Neutrophil is one of a tumor's most abundant immune cell groups. This research aimed to investigate the complex communication network in the immune microenvironment (TIME) of NSCLC tumors to clarify the interaction between immune cells and tumors and establish a prognostic risk model that can predict immune response and prognosis of patients by analyzing the characteristics of Neutrophil differentiation. Integrated single-cell RNA sequencing (scRNA-seq) data from NSCLC samples and Bulk RNA-seq were used for analysis. Twenty-eight main cell clusters were identified, and their interactions were clarified. Subsequently, four subsets of Neutrophils with different differentiation states were found, closely related to immune regulation and metabolic pathways. Based on the ratio of four housekeeping genes (ACTB, GAPDH, TFRC, TUBB), six Neutrophil differentiation-related genes (NDRGs) prognostic risk models, including MS4A7, CXCR2, CSRNP1, RETN, CD177, and LUCAT1, were constructed by Elastic Net and Multivariate Cox regression, and patients' total survival time and immunotherapy response were successfully predicted and validated in three large cohorts. Finally, the causes of the worse prognosis of NSCLC caused by six prognostic genes were explored, and the small molecular compounds targeted at the anti-tumor effect of prognostic genes were screened. This study clarifies the TIME regulation network in NSCLC and emphasizes the critical role of NDRGs in predicting the prognosis of patients with NSCLC and their potential response to immunotherapy, thus providing a promising therapeutic target for NSCLC.