I am trying to calculate the habitat fragmentation of a forested region in Africa in order to protect the areas less fragmented. I have used the tool "Landscape Fragmentation Tool V2 requiring Spatial Analyst" but I have an error during its execution. Could you recommended me another tool to calculate fragmentation of habitats using arcgis 10.1?
The input that I used was a land cover map with 1 and 2 data, and the edge width was 50.
This is the error:
Executing: Landscape_Fragmentation_Tool Bool_fragm 50 C:\CW2\project_2 fragmentation
Start Time: Mon May 12 22:51:18 2014
Running script Landscape_Fragmentation_Tool...
OUTPUT DATASET NAME IS:
C:\CW2\project_2\fragmentation.img
DATA EXTRACTED FROM LAND COVER
EDGE FRAGMENTATION IDENTIFIED
INTERIOR, PATCH, AND PERFORATED FRAGMENTATION IDENTIFIED
CATEGORIZING CORE FOREST PATCHES...
File "C:\arcgis_extras\Landscape_Fragmentation_SA\LFT_SA.py", line 254, in <module>
gp.RegionGroup_sa("core1_", core_grp, "EIGHT", "WITHIN")
<class 'arcgisscripting.ExecuteError'>: ERROR 010213: Error in reading raster c:\fragdata\core1_.
ERROR 010067: Error in executing grid expression.
Failed to execute (RegionGroup).
Completed script Landscape_Fragmentation_Tool...
Failed to execute (Landscape_Fragmentation_Tool).
Failed at Mon May 12 23:33:04 2014 (Elapsed Time: 41 minutes 46 seconds)
Thanks!
UPDATE:
Here you can find more lines in case you have some idea:
#-------------------------------------------------------------------
# CLASSIFY CORE FOREST BASED ON PATCH SIZE
gp.workspace = TempWS
gp.AddMessage("CATEGORIZING CORE FOREST PATCHES...")
print "CATEGORIZING CORE FOREST PATCHES..."
gp.Reclassify_sa (fragMap2_f, "Value", "4 1", "core1_", "NODATA")
# region group core class...
core_grp = TempWS + "\\core_grp.img"
gp.RegionGroup_sa("core1_", core_grp, "EIGHT", "WITHIN") <-- The error is here
# convert core groups to txt file...
coreGrpFile = "%s\\coreGrpFile.txt" % TempWS
gp.RasterToASCII_conversion(core_grp, coreGrpFile)
# convert 4 class frag map to txt file...
fragMap_4c = "%s\\fragMap_4c.txt" % TempWS
gp.RasterToASCII_conversion(fragMap2_f, fragMap_4c)
cur = gp.SearchCursor(core_grp)
row = cur.next()
remap_dct = {}