def test_SharedWallsRatio(self): self.df_buildings["swr"] = mm.SharedWallsRatio(self.df_buildings).series self.df_buildings["swr_array"] = mm.SharedWallsRatio( self.df_buildings, perimeters=self.df_buildings.geometry.length ).series nonconsecutive = self.df_buildings.drop(2) result = mm.SharedWallsRatio(nonconsecutive).series check = 0.3424804411228673 assert self.df_buildings["swr"][10] == check assert self.df_buildings["swr_array"][10] == check assert result[10] == check
def test_SharedWallsRatio(self): self.df_buildings["swr"] = mm.SharedWallsRatio(self.df_buildings, "uID").series self.df_buildings["swr_uid"] = mm.SharedWallsRatio( self.df_buildings, range(len(self.df_buildings))).series self.df_buildings["swr_array"] = mm.SharedWallsRatio( self.df_buildings, "uID", self.df_buildings.geometry.length).series nonconsecutive = self.df_buildings.drop(2) result = mm.SharedWallsRatio(nonconsecutive, "uID").series check = 0.3424804411228673 assert self.df_buildings["swr"][10] == check assert self.df_buildings["swr_uid"][10] == check assert self.df_buildings["swr_array"][10] == check assert result[10] == check
def time_SharedWallsRatio(self): mm.SharedWallsRatio(self.df_buildings)
import momepy as mm import geopandas as gpd import libpysal 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") blg["mtbSWR"] = mm.SharedWallsRatio(blg, "uID", "sdbPer").series queen_1 = libpysal.weights.contiguity.Queen.from_dataframe(tess, ids="uID") blg["mtbAli"] = mm.Alignment(blg, queen_1, "uID", "stbOri").series blg["mtbNDi"] = mm.NeighbourDistance(blg, queen_1, "uID").series tess["mtcWNe"] = mm.Neighbours(tess, queen_1, "uID", weighted=True).series tess["mdcAre"] = mm.CoveredArea(tess, queen_1, "uID").series str_q1 = libpysal.weights.contiguity.Queen.from_dataframe(streets) streets["misRea"] = mm.Reached(streets, tess, "nID", "nID", spatial_weights=str_q1, mode="count").series streets["mdsAre"] = mm.Reached(streets, tess, "nID", "nID",