4

私は geopandas をテストして、非常に単純なものを作成しています。差分メソッドを使用して、円の内側にある GeoDataFrame のいくつかのポイントを削除します。

これが私のスクリプトの始まりです：

%matplotlib inline
# previous line is because I used ipynb
import pandas as pd
import geopandas as gp
from shapely.geometry import Point
[...]
points_df = gp.GeoDataFrame(csv_file, crs=None, geometry=geometry)

points_df の最初の行は次のとおりです。

    Name        Adress      geometry
0   place1      street1     POINT (6.182674 48.694416)
1   place2      street2     POINT (6.177306 48.689889)
2   place3      street3     POINT (6.18 48.69600000000001)
3   place4      street4     POINT (6.1819 48.6938)
4   place5      street5     POINT (6.175694 48.690833)

次に、最初の GeoDF のいくつかのポイントを含むポイントを追加します。

base = points_df.plot(marker='o', color='red', markersize=5)

center_coord = [Point(6.18, 48.689900)]
center = gp.GeoDataFrame(crs=None, geometry=center_coord)
center.plot(ax=base, color = 'blue',markersize=5)

circle = center.buffer(0.015)
circle.plot(ax=base, color = 'green')

iPython ノートブックによって表示される結果は次のとおりです。

ポリゴンとポイント

ここでの目標は、緑の円の内側にある赤い点を削除することです。そのためには差分法で十分だと思いました。しかし、私が書くとき：

selection = points_df['geometry'].difference(circle)
selection.plot(color = 'green', markersize=5)

結果は... points_df で何も変わっていません:

変更なし

difference() メソッドはポリゴン GeoDataFrames でのみ機能し、ポイントとポリゴンの混合は不可能だと思います。しかし、多分私は何かを逃した!

この場合、円内の点の存在をテストする関数は差分法よりも優れていますか?

4

1 に答える 1

5

difference() メソッドはポリゴン GeoDataFrames でのみ機能し、ポイントとポリゴンの混合は不可能だと思います。

それが問題のようです。ポイントでオーバーレイを使用することはできません。

また、その種の空間操作については、単純な空間結合が最も簡単なソリューションのようです。

最後の例から始めます ;):

%matplotlib inline
import pandas as pd
import geopandas as gp
import numpy as np
import matplotlib.pyplot as plt
from shapely.geometry import Point

# Create Fake Data
df = pd.DataFrame(np.random.randint(10,20,size=(35, 3)), columns=['Longitude','Latitude','data'])

# create Geometry series with lat / longitude
geometry = [Point(xy) for xy in zip(df.Longitude, df.Latitude)]

df = df.drop(['Longitude', 'Latitude'], axis = 1)

# Create GeoDataFrame
points = gp.GeoDataFrame(df, crs=None, geometry=geometry)

# Create Matplotlib figure
fig, ax = plt.subplots()

# Set Axes to equal (otherwise plot looks weird)
ax.set_aspect('equal')

# Plot GeoDataFrame on Axis ax
points.plot(ax=ax,marker='o', color='red', markersize=5)

# Create new point
center_coord = [Point(15, 13)]
center = gp.GeoDataFrame(crs=None, geometry=center_coord)

# Plot new point
center.plot(ax=ax,color = 'blue',markersize=5)
# Buffer point and plot it
circle = gp.GeoDataFrame(crs=None, geometry=center.buffer(2.5))

circle.plot(color = 'white',ax=ax)

問題

ポイントがポリゴンの内側にあるか外側にあるかを判断する方法に関する問題が残ります...それを達成する1つの方法は、ポリゴン内のすべてのポイントを結合し、すべてのポイントと内のポイントの違いでDataFrameを作成することですサークル：

# Calculate the points inside the circle 

pointsinside = gp.sjoin(points,circle,how="inner")

# Now the points outside the circle is just the difference 
# between  points and points inside (see the ~)

pointsoutside = points[~points.index.isin(pointsinside.index)]


# Create a nice plot 
fig, ax = plt.subplots()
ax.set_aspect('equal')
circle.plot(color = 'white',ax=ax)
center.plot(ax=ax,color = 'blue',markersize=5)
pointsinside.plot(ax=ax,marker='o', color='green', markersize=5)

pointsoutside.plot(ax=ax,marker='o', color='yellow', markersize=5)

print('Total points:' ,len(points))
print('Points inside circle:' ,len(pointsinside))
print('Points outside circle:' ,len(pointsoutside))

合計点: 35

円内のポイント: 10

円の外側のポイント: 25

問題が解決しました ;）

于 2016-11-27T12:59:33.993 に答える