Python相当于R的ConnCompLabel

时间:2017-06-10 16:04:50

标签: python r

我有一个用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)
  }
}

0 个答案:

没有答案