Я работаю с набором данных (~ 40000 строк), который содержит информацию xy для каждой строки. Поскольку этот набор данных был составлен за ~ 20 лет, произошел сдвиг в данных, которые использовали люди, отправляющие образцы. Но не было очень хорошего контроля качества для этого набора данных в течение примерно 5 лет go, поэтому существует около 30000 записей, которые имеют координаты UTM xy (и информацию о зоне), но не имеют данных о базовой точке.
I Я хотел бы решить эту проблему, но я не совсем уверен, как подойти к этому. Учитывая, что исходные данные WGS84 изначально были сопоставлены с исходными данными NAD83 (хотя с тех пор они немного разошлись), я подумал, что будет проще взять неизвестные значения UTM базовых данных, спроецировать их с помощью spTransform в rgdal, построить исходный образец координаты и прогнозируемые координаты, и если они были на расстоянии более нескольких метров, добавьте значение в столбец в кадре данных, который отмечает, что исходный UTM, скорее всего, был NAD27. Данные выглядят так, хотя есть еще около 40 столбцов:
ID Easting Northing Zone
1 721234 3811234 12
2 721334 3811334 12
Однако я не могу понять, как заставить работать часть сравнения. У кого-нибудь есть предложения или кто-нибудь сделал что-то подобное?
Спасибо.