arcpy.gp.RegionGroup_sa("binary","C:/FELIX_DATA/SCRARCH/GROUPS","EIGHT","WITHIN","NO_LINK","#")
arcpy.gp.SingleOutputMapAlgebra_sa("$$XMAP","C:/FELIX_DATA/SCRARCH/X","#")
arcpy.gp.SingleOutputMapAlgebra_sa("$$YMAP","C:/FELIX_DATA/SCRARCH/Y","#")
arcpy.gp.ZonalStatistics_sa("GROUPS","VALUE","X","C:/FELIX_DATA/SCRARCH/xMean","MEAN","DATA")
arcpy.gp.ZonalStatistics_sa("GROUPS","VALUE","Y","C:/FELIX_DATA/SCRARCH/yMean","MEAN","DATA")
arcpy.gp.RasterCalculator_sa("""Power(("X" - "xMean")*("X" - "xMean")+("Y" - "yMean")*("Y" - "yMean"),0.5)""","C:/FELIX_DATA/SCRARCH/distance")
arcpy.gp.ZonalStatisticsAsTable_sa("GROUPS","VALUE","distance","C:/FELIX_DATA/SCRARCH/stats.dbf","DATA","MEAN")
arcpy.RasterToPolygon_conversion("GROUPS","C:/FELIX_DATA/SCRARCH/polygons.shp","SIMPLIFY","VALUE")
arcpy.AddField_management("POLYGONS","R_GYR","FLOAT","#","#","#","#","NULLABLE","NON_REQUIRED","#")
arcpy.AddJoin_management("POLYGONS","GRIDCODE","stats","VALUE","KEEP_ALL")
arcpy.CalculateField_management("POLYGONS","polygons.R_GYR","[stats.MEAN]","VB","#")

outRegionGrp = RegionGroup(binary, "EIGHT")