У меня есть шейп-файл полигонов и другой файл CSV, который содержит список точек в виде (Lat, Lng) пар.
Я хочу проверить для каждой пары (lat, lng) из файла CSV, в какой полигон он попадает внутрь ..
Шейп-файл проецируется, и файл proj выглядит так:
PROJCS["Transverse_Mercator",GEOGCS["GCS_OSGB 1936",
DATUM["D_OSGB_1936",SPHEROID["Airy_1830",6377563.396,299.3249646]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",49],PARAMETER["central_meridian",-2],PARAMETER["scale_factor",0.9996012717],PARAMETER["false_easting",400000],PARAMETER["false_northing",-100000],UNIT["Meter",1]]
Мой план следующий:
- Прочтите шейп- файл, используя
readShapePoly
функцию изMapTools
пакета R. - Считать координаты точек из CSV-файла в информационный кадр и преобразовать его в SpatialPointsDataFrame
- Используйте
over
функцию, чтобы определить, в какой полигон он попадает.
Для этого мне нужно указать proj4string
время загрузки шейп-файла на шаге 1, а также преобразовать координаты из файла CSV в ту же систему проекции, используя spTransform
функцию, прежде чем применять over
функцию на шаге 3, так как для этого необходимо, чтобы точки и многоугольники быть под той же системой проекции.
Любая идея о том, что должно правильное значение для содержимого файла proj, показанного выше?
readOGR
функции Я всегда получаю не удается открыть файл ошибки
summary
функции для SpatialPolygonDataFrame
объекта дало мне правильное значение дляproj4string