Geographically Constrained Transformer for Spatiotemporal Reconstruction of 2 m NDVI in Complex Coastal Landscapes

High-resolution Normalized Difference Vegetation Index (NDVI) data are essential for monitoring fine-scale coastal environmental dynamics, yet persistent cloud cover, rapid geomorphic change, and strong spatial heterogeneity limit the availability of temporally continuous observations. Existing spatiotemporal fusion approaches can partially address these limitations, but many rely primarily on data-driven feature learning and do not explicitly incorporate geographic information, leading to boundary blurring, structural inconsistency, and sensitivity to background noise in complex coastal environments. This study presents a geographically constrained Transformer-based framework for 2 m NDVI spatiotemporal reconstruction in coastal landscapes named Coastal-Prior-Embedded Global–Local Fusion Transformer (Coastal-GLFT). The approach integrates high-resolution Gaofen-6 panchromatic and multispectral imagery with high-frequency wide-field-view observations and auxiliary geographic datasets describing elevation, coastline proximity, and land use/land cover. Geographic priors were incorporated as explicit spatial constraints, while a spatiotemporal gating mechanism and global–local fusion architecture were used to improve the representation of temporal variation and multi-scale spatial structure. The method was evaluated using a multi-temporal dataset for the Yellow River Delta comprising 49 high-resolution scenes and 137 coarse-resolution scenes acquired between 2020 and 2025. Compared with representative physics-based, convolutional neural network, generative adversarial network, and Transformer-based fusion methods, the proposed approach reduced reconstruction error by approximately 5–72%, increased signal fidelity by approximately 1–12%, and improved structural similarity by approximately 2–52%. Compared with the strongest Transformer-based baseline, SwinSTFM, Coastal-GLFT reduced RMSE from 0.0896 to 0.0855, increased PSNR from 36.19 dB to 37.09 dB, and improved SSIM from 0.8551 to 0.8742. Qualitative analysis further demonstrated improved preservation of boundary structure, spatial continuity, and heterogeneous coastal features, including aquaculture ponds, tidal creeks, and fragmented wetlands. These results indicate that integrating geographic constraints with multi-scale Transformer-based reconstruction can improve the fidelity and structural consistency of high-resolution NDVI reconstruction in complex coastal environments. The framework provides a basis for fine-scale coastal vegetation monitoring and land-cover analysis, while future work should assess transferability across diverse coastal systems and improve computational scalability.