In this study, a temporal shape–based fusion method using a spatially and temporally moving window is proposed to incorporate time lag of fine and coarse resolution observations, and to fully utilize target fine resolution pixel and similar coarse resolution pixels in the process. This method provides high accuracy fused images with Pearson’s r of ~0.95, root mean square error of ~0.04, and bias of ~0.01 for commonly used fine spatial resolution satellites, including Landsat 7 and 8, Sentinel 2, and Gaofen 1, over different heterogeneous regions, such as urban, mountain, forest, and savanna regions. The fused fine resolution Enhanced Vegetation Index (EVI) time series using different fine spatial resolution satellites data as input are all highly correlated with the PhenoCam monitored green chromatic coordinate, with no temporal lag. Compared with commonly used data fusion method, this method provides equivalent and slightly higher accuracy because both neighboring similar pixels and the annual temporal variation are fully considered. This temporal shape–based fusion method does not require each input fine resolution image to be cloud-free; therefore, it can be used at a large spatial scale without further preprocessing and generates continuous datasets over a long-time range with only one input preparation process. The factors that could affect the method accuracy are the cloud detection accuracy of fine resolution data and the temporal continuity of the coarse resolution data. The method may also be used to produce spatially and temporally continuous surface reflectance and other surface reflectance derived indices.