我有一个用R
编写的代码,我需要在Python
中重写它。
现在,我被ConnCompLabel
SDMTools
个问津卡在R
上
这是我需要重写的代码:
indice<-function(NDVI,threshold,theo=NULL)
{
nCell=sum(values(NDVI)>threshold,na.rm=T)
nPotCell=sum(!is.na(values(NDVI)))
if(nCell<=1){
return(1)
}else{
require(SDMTools)
if(is.null(theo)){
theo=load("DATA/B2-fragTheo.RData")
theo=get(theo)
}
ccl.mat = ConnCompLabel(NDVI>threshold)
nFrag=max(values(ccl.mat),na.rm=T)
nFragTheo=(theo$nfragtheo[(theo$ncell)/10000==round(nCell/nPotCell,4)])*nPotCell/10000
return(nFrag/nFragTheo)
}
}