I am using Landsat 5, 7, 8 images and RStudio. My goal is to compute the disturbance in a given area using the Vegetation Change Tracker (VCT) algorithm. I would like to ask you if you could provide me a general guidance on how to construct it using R Studio. I computed the IFZ of each individual image, I created the LTSS, I have the NDVI and the NBR of each year. My problem is on how to create the rules required form the VCT (the rules can be seen in the image below).

The acronyms CLIV, MNCLIV and MNCHIV means: - the pixel was classified as a forest at least in 2 sequentially images, - the minimum years of sequentially images needed for a pixel to be classified as a forest (2 years) and - minimum number of consecutive high IFZ values for a pixel to be classified as a non-forest pixel. The paper named An automated approach for reconstructing recent forest disturbance history using dense Landsat time series stack.