Як відфільтрувати неправильні GPS-точки?

У мене база даних GPS-точок. Немає треків, лише окуляри. Мені потрібно обчислити деяке значення на кожні 100 метрів, але іноді GPS давав неправильні координати, які лежать далеко від реальних точок GPS, і замість обчислення значень для маленького квадрата я мушу обчислити його для дійсно великої прямокутної площі.

Який найкращий алгоритм фільтрації неправильних GPS-точок?

Я зробив скріншот, щоб допомогти зрозуміти:

Запустіть Ігор Анселіна Локального Морана проти очок і викиньте що-небудь із z-балом нижче -1,96. Це статистичний метод визначення місця просторів. Ви повинні переконатися, що для всіх точок є значення, пов'язане з їх просторовим положенням.

Але перевіряючи інструменти 10.1 після коментаря whuber, я розумію, що якщо ви використовуєте ArcGIS 10.1, доступний інструмент аналізу групування, що дійсно ви хочете зробити.

Я думаю, ви хочете зробити групуючий аналіз з просторовим обмеженням триангуляції Делоне. Проблема перешкод тут полягає в тому, що вам потрібно мати кількість груп розподілу, рівних або більших, ніж кількість відключених груп (якщо хтось із залишків є природними сусідами один для одного). Інакше люди, які не мають природних сусідів, не знайдуть групи з аналізу групування.

Виходячи з цього, я думаю, що триангуляція Делоні може бути джерелом алгоритму фільтрування, але я ще не впевнений.

Ще одне оновлення: Після копання у Partition.py, скрипт, який запускає інструмент аналізу групування, я думаю, що можна використовувати алгоритм для відключених груп у поєднанні з частиною NoNeighbors, хоча у мене виникають проблеми з викопанням цієї частини сценарій.

Це може допомогти скласти список людей, що вижили:

SELECT p1.point_id 
FROM p1 AS points, p2 AS points
WHERE p1.point_id <> p2.point_id AND
ST_Distance(p1.geom, p2.geom) > 10000

Тут point_id буде основним ключем у вашій таблиці очок. Функція відстані знайде точки, де найближча більша за 10000 метрів. (Ви, звичайно, можете поставити будь-яке значення, відповідне)

Якщо вищезазначене працює, то перейдіть до пункту DELETE, наприклад:

DELETE FROM points WHERE point_id IN (
-- SELECT as above
SELECT ....
);

Я спробую надати більш практичну відповідь, щоб допомогти вам виконати роботу. (вибачте, якщо шукаєте обговорення алгоритмів)

Сценарій 1: Ви згадуєте «GPS точки», тож якщо у вас є доступ до оригінальних точок GPS, робота стає набагато простішою. Ви можете викинути точки з високим HDOP / VDOP або кількістю супутників на огляді - що могло б спричинити помилку спочатку. Безкоштовний інструмент, як gpsbabel, має вбудовані такі фільтри. http://www.gpsbabel.org/htmldoc-development/Data_Filters.html

Сценарій 2: Ви просто маєте набір точок. Потім проблемою стає виявлення просторових видатків. У цій галузі є багато досліджень, і я бачу багато робіт на цю тему з веб-пошуку. Якщо ви хочете очистити свої дані, ви можете використовувати алгоритм vASS.lierlier GRASS, який повинен працювати у вашому випадку на основі екрана, яким ви поділилися. http://grass.osgeo.org/gdp/html_grass63/v.outlier.html

Я думаю, у вас є непотрібні дані. Реально, якщо ви дбаєте про те, що деякі дані неправильні, і ви не можете достовірно визначити кожну неправильну точку, використовуючи якийсь інший фактор, то у вашому аналізі з'являться погані дані.

Якщо це має значення, то, ймовірно, слід розглянути питання про викид все, з'ясувати першопричину (наприклад, погані точки GPS з багатостороннього шляху), вирішити цю першопричину (наприклад, додати антену задушення або краще тип GPS, або все, що найкраще виправити є), а потім повторно збирати дані.

Якщо погані дані не мають значення, просто використовуйте їх і ігноруйте помилки.