Эта формула никак не относится к вегетационному индексу, да и к самому растру она не привязана.
Описание всех присутствующих в ней переменных есть
тут.
ESUN - константа, есть в справочных таблицах на сайтах операторов съемочной системы (в данном случае USGS), если специфичный сенсор или новый, можно и рассчитывать на основе одной из моделей.
d - есть в справочных таблицах для дней юлианского календаря, есть упрощенный вариант расчета по формуле, есть более продвинутые варианты расчета (пример по ссылке выше)
Ѳ - формула по ссылке выше, sun elevation есть в метаданных к каждому снимку, также sun elevation можно рассчитать зная координаты, дату и время съемки.