У мене є форм-формат полігонів та ще один файл CSV, який містить список точок як пари (Lat, Lng).
Я хочу перевірити для кожної пари (лат, 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