def test_BlocksCount(self):
sw = mm.sw_high(k=5, gdf=self.df_tessellation, ids="uID")
count = mm.BlocksCount(self.df_tessellation, "bID", sw, "uID").series
count2 = mm.BlocksCount(self.df_tessellation, self.df_tessellation.bID,
sw, "uID").series
unweigthed = mm.BlocksCount(self.df_tessellation,
"bID",
sw,
"uID",
weighted=False).series
check = 3.142437439120778e-05
check2 = 5.222222222222222
assert count.mean() == check
assert count2.mean() == check
assert unweigthed.mean() == check2
with pytest.raises(ValueError):
count = mm.BlocksCount(self.df_tessellation,
"bID",
sw,
"uID",
weighted="yes")
def time_BlocksCount(self):
mm.BlocksCount(self.df_tessellation, "bID", self.sw5, "uID")
tess = gpd.read_file('files/elements.gpkg', layer='tessellation')
blg = gpd.read_file('files/elements.gpkg', layer='buildings')
blocks = gpd.read_file('files/elements.gpkg', layer='blocks')
streets = gpd.read_file('files/elements.gpkg', layer='streets')
queen1 = mm.sw_high(k=1, gdf=tess, ids='uID')
queen3 = mm.sw_high(k=3, weights=queen1)
blg_queen = mm.sw_high(k=1, gdf=blg, ids='uID')
blg['ltbIBD'] = mm.MeanInterbuildingDistance(blg, queen1, 'uID', queen3).series
blg['ltcBuA'] = mm.BuildingAdjacency(blg, queen3, 'uID', blg_queen).series
# blg['temp_fa'] = mm.floor_area(blg, 'sdbHei', 'sdbAre')
# tess = tess.merge(blg[['temp_fa', 'uID']], on='uID', how='left')
# tess['licGDe'] = mm.density(tess, 'temp_fa', queen3, 'uID', 'sdcAre')
tess['ltcWRB'] = mm.BlocksCount(tess, 'bID', queen3, 'uID').series
tess.to_file('files/elements.gpkg', layer='tessellation', driver='GPKG')
blg.to_file('files/elements.gpkg', layer='buildings', driver='GPKG')
fo = libpysal.io.open('files/GRqueen1.gal', 'w')
fo.write(queen1)
fo.close()
fo = libpysal.io.open('files/GRqueen3.gal', 'w')
fo.write(queen3)
fo.close()
fo = libpysal.io.open('files/GRblg_queen.gal', 'w')
fo.write(blg_queen)
fo.close()