การดำเนินการ "spatial join"

การทำงานกับข้อมูลเชิงพื้นที่ใน Python

Dani Arribas-Bel

Geographic Data Science Lab (University of Liverpool)

ความสัมพันธ์เชิงพื้นที่ I

ความสัมพันธ์เชิงพื้นที่ II

เมืองใดบ้างที่อยู่ในบราซิล?

brazil = countries.loc[22, 'geometry']
cities[cities.within(brazil)]

               name                                       geometry
169        Brasília  POINT (-47.91799814700306 -15.78139437287899)
238  Rio de Janeiro  POINT (-43.22696665284366 -22.92307731561596)
239       São Paulo  POINT (-46.62696583905523 -23.55673372837896)

แล้วถ้าต้องการทราบว่าแต่ละเมืองอยู่ในประเทศใดล่ะ?

Spatial Join

SPATIAL JOIN = การถ่ายโอนแอตทริบิวต์จากเลเยอร์หนึ่งไปยังอีกเลเยอร์ตามความสัมพันธ์เชิงพื้นที่

Spatial Join ด้วย GeoPandas

joined = geopandas.sjoin(cities,
                         countries[['name', 'geometry']],
                         op="within")

joined.head()

        name_left                                     geometry name_right
0    Vatican City  POINT (12.45338654497177 41.90328217996012)      Italy
1      San Marino    POINT (12.44177015780014 43.936095834768)      Italy
226          Rome    POINT (12.481312562874 41.89790148509894)      Italy
2           Vaduz  POINT (9.516669472907267 47.13372377429357)    Austria
212        Vienna  POINT (16.36469309674374 48.20196113681686)    Austria

มาฝึกกันเถอะ!

การทำงานกับข้อมูลเชิงพื้นที่ใน Python