Accurate delineation of fields is difficult in fragmented landscapes where single-date images provide no seasonal cues and supervised models require labels. We propose a method that explicitly represents phenology to improve zero-shot delineation. Using 22 cloud-free PlanetScope scenes over a 5 × 5 km area, a single harmonic model is fitted to the NDVI per pixel to obtain the phase, amplitude and mean. These values are then mapped into cylindrical colour spaces (Hue–Saturation–Value, Hue–Whiteness–Blackness, Luminance-Chroma-Hue). The resulting recoloured composites are segmented using the Segment Anything Model (SAM), without fine-tuning. The results are evaluated object-wise, object-wise grouped by area size, and pixel-wise. Pixel-wise evaluation achieved up to F1 = 0.898, and a mean Intersection-over-Union (mIoU) of 0.815, while object-wise performance reached F1 = 0.610. HSV achieved the strongest area match, while HWB produced the fewest fragments. The ordinal time-of-day basis provided better parcel separability than the annual radian adjustment. The main errors were over-segmentation and fragmentation. As the parcel size increased, the IoU increased, but the precision decreased. It is concluded that recolouring using harmonic NDVI time series is a simple, scalable, and interpretable basis for field delineation that can be easily improved.

Operational agricultural monitoring in the Central European lowlands requires timely parcel boundaries; however, unmarked field edges produce minimal spectral contrast in single-date imagery. Previous works demonstrated that harmonic NDVI encoding enables zero-shot field delineation using foundational models, but the influence of the spectral index choice on temporal boundaries remained unquantified. This study systematically evaluates eleven vegetation indices—NDVI, GNDVI, NDRE, EVI, EVI2, SAVI, MSAVI, NDWI, CIg, CIre, and NDYVI—within a fixed harmonic phenology encoding pipeline. A one-year PlanetScope time series (15 × 15 km, Slavonija, Croatia) was decomposed via annual sinusoidal regression to extract per-pixel phase, amplitude, and mean parameters. These harmonic descriptors were mapped to HSV colour channels and segmented using the Segment Anything Model without fine-tuning. Official agricultural parcels (PAAFRD, 2025) provided ground truth for pixel-wise, object-wise, and size-stratified evaluation. Performance stratified into three tiers based on object-wise metrics. Soil-adjusted and enhanced-greenness indices (MSAVI, EVI, EVI2, and SAVI) achieved F1 = 0.51–0.52, and mIoU = 0.70–0.71, statistically outperforming standard ratio formulations (NDVI: F1 = 0.49) and chlorophyll indices (CIg, CIre: F1 = 0.45–0.47). Pixel-wise scores remained compressed (F1 > 0.88 across all indices), indicating consistent interior coverage but index-dependent boundary precision. Error analysis revealed scale-dependent patterns: merging dominated small parcels (<10,000 m²), while fragmentation increased with parcel size. Results demonstrate that spectral formulation is a systematic design factor in phenology-based delineation, with soil background correction and dynamic range compression improving seasonal trajectory separability. The harmonic parameters generated by this framework provide feature-ready input for crop classification, suggesting that integrated boundary extraction and crop mapping workflows merit further investigation.