Какой тип данных для широты и долготы?

В PostGIS для точек с широтой и долготой существует тип данных geography.

Чтобы добавить столбец:

alter table your_table add column geog geography;

Чтобы вставить данные:

insert into your_table (geog) values ('SRID=4326;POINT(longitude latitude)');

4326 - это идентификатор пространственной привязки, в котором указаны данные в градусах долготы и широты, как и в GPS. Подробнее об этом: http://epsg.io/4326.

Порядок - это долгота, широта - поэтому, если вы построите его как карту, это будет (x, y).

Чтобы найти ближайшую точку, вам нужно сначала создать пространственный индекс:

create index on your_table using gist (geog);

а затем запросите, скажем, 5 ближайших к заданной точке:

select * 
from your_table 
order by geog <-> 'SRID=4326;POINT(lon lat)' 
limit 5;

Я настоятельно рекомендую PostGis. Он специфичен для такого типа данных и имеет готовые методы для расчета расстояния между точками, среди других операций ГИС, которые вы можете найти полезными в будущем.

Если вам не нужны все функции, предлагаемые PostGIS, Postgres (в настоящее время) предлагает модуль расширения под названием earthdistance. Он использует точку или cube в зависимости от вашей точности, необходимой для расчета расстояния.

Теперь вы можете использовать функцию earth_box, например, для запроса точек на определенном расстоянии от местоположения.

В PostGIS геометрия предпочтительнее географии (круглая модель земли), потому что вычисления намного проще, а значит, быстрее. Он также имеет МНОЖЕСТВО доступных функций, но менее точен на очень больших расстояниях.

Импортируйте поля длинного и широтного формата CSV в DECIMAL(10,6) столбец. 6 цифр - это точность 10 см, этого должно хватить для большинства случаев использования. Затем передайте импортированные данные в правильный SRID.

Неправильный путь!

/* try what seems the obvious solution */
DROP TABLE IF EXISTS public.test_geom_bad;
-- Big Ben, London
SELECT ST_SetSRID(ST_MakePoint(-0.116773, 51.510357),4326) AS geom
INTO public.test_geom_bad;

ПРАВИЛЬНЫЙ способ

/* add the necessary CAST to make it work */
DROP TABLE IF EXISTS public.test_geom_correct;
SELECT ST_SetSRID(ST_MakePoint(-0.116773, 51.510357),4326)::geometry(Geometry, 4326) AS geom
INTO public.test_geom_correct;

Убедитесь, что SRID не равен нулю!

/* now observe the incorrect SRID 0 */
SELECT * FROM public.geometry_columns
WHERE f_table_name IN ('test_geom_bad','test_geom_correct');

Проверьте порядок вашего параметра долгой широты с помощью средства просмотра WKT и

SELECT ST_AsEWKT(geom) FROM public.test_geom_correct

Затем проиндексируйте его для лучшей производительности

CREATE INDEX idx_target_table_geom_gist
    ON target_table USING gist(geom);

Используйте тип данных Точка, чтобы сохранить долготу и широту в одном столбце:

CREATE TABLE table_name (
    id integer NOT NULL,
    name text NOT NULL,
    location point NOT NULL,
    created_on timestamp with time zone NOT NULL DEFAULT CURRENT_TIMESTAMP,
    CONSTRAINT table_name_pkey PRIMARY KEY (id)
)

Создайте индексирование для столбца "местоположение":

CREATE INDEX ON table_name USING GIST(location);

Индекс GiST позволяет оптимизировать поиск «ближайшего соседа»:

SELECT * FROM table_name ORDER BY location <-> point '(-74.013, 40.711)' LIMIT 10;

Примечание. Первый элемент точки - это долгота, а второй - широта.

Для получения дополнительной информации см. Операторы запроса < / а>.