Mam shapefile o nazwie „ind_adm” i SpatialPointsDataFrame o nazwie „pnts”. „Pnts” zawiera punkty generowane losowo, a niektóre punkty pokrywają się z wielokątem. Zobacz obrazek poniżej.

Teraz chcę zrobić punkt w analizie wielokątów, tj. Chcę dowiedzieć się, które punkty leżą wewnątrz szarego wielokąta reprezentującego granicę Indii. Do tego używam funkcji over () w bibliotece sp.

pt.in.poly <- sp::over(ind_adm, pnts, fn = mean) #do the join

Jednak wyjście, które otrzymuję, jest

    >pt.in.poly
    values
    0 6.019467

Powinienem uzyskać indeks punktów, które są „w” wielokącie.

Gdzie się mylę?