Modifying fieldIndex to add more VIs

[pjjagwe]

Hi @filipematias23, I added more VIs to the Indices.txt file used to run the fieldIndex function.

When i run this code; X.Indices2 <- fieldIndex2(X.RemSoil$newMosaic,Red=1,Green=2,Blue=3,RedEdge=4,NIR=5, index = c("NDVI","NDRE", "BI", "SCI", "GLI", "HI","NGRDI", "SI", "VARI", "HUE", "BGI", "PSRI", "GNDVI", "RVI", "TVI", "CVI", "CIG", "CIRE", "DVI", "EVI", "GDVI","GLI","GOSAVI","GRVI", "GEMI","GARI" ,"GSAVI","IPVI","LAI", "MNLI","MSAVI2","MSR","NLI","OSAVI","RDVI","SAVI","TDVI","VARI","WDRVI","SR"))

I get an error: Error in eval(parse(text = as.character(Ind$eq[as.character(Ind$index) == : object 'EVI' not found In addition: Warning messages: 1: In for (f in 1:length(filename)) { : closing unused connection 4 (C:/Users/pjjagwe/AppData/Local/Temp/RtmpaOFv1e/raster/r_tmp_2023-08-25_181651.200041_16852_88055.gri) 2: In for (f in 1:length(filename)) { : closing unused connection 3 (C:/Users/pjjagwe/AppData/Local/Temp/RtmpaOFv1e/raster/r_tmp_2023-08-25_180154.128608_16852_62923.gri)

[pjjagwe]

@filipematias23 Could you kindly share some insight.

Thanks

[filipematias23]

Hi there, try using the parameter "myIndex" in fieldIndex(). You can write the new indices direct in the coding. Here is one example:

Test.Indices<- fieldIndex(mosaic = Test, Red = 1, Green = 2, Blue = 3, RedEdge=4, NIR=5, index = c("NGRDI","NDVI"), myIndex = c("Write your Index", "Write your Index", "2*Green/Blue"))

Hope it helps...!

All the best, Filipe

[pjjagwe]

Thanks @filipematias23