Postgresql + PostGIS: найти запись по столбцам геометрии

Question

У меня есть таблица:

CREATE TABLE main.address_list 
(
    id serial NOT NULL, 
    address_string text NOT NULL, 
    address_position geometry(Point,4326) 
    CONSTRAINT pk_address_list PRIMARY KEY (id) 
) 

Я хочу найти некоторые записи в этой таблице, используя address_position колонки:

select * 
     , public.ST_Y(al.address_position) latitude 
     , public.ST_X(al.address_position) longitude 
from address_list al 
where al.address_position = public.ST_SetSRID(public.ST_MakePoint(56.187339 --longitude$n 
                    , 57.964295 --latitude$n 
                   ), 4326) 

Но я получил следующее ошибка:

ERROR: operator is not unique: public.geometry = public.geometry 
LINE 5: where al.address_position = public.ST_SetSRID(public.ST_Make... 
           ^
HINT: Could not choose a best candidate operator. You might need to add explicit type casts. 

Как я могу получить явный тип приведения?

Answer 1

Вы пытаетесь сравнить два geometries с оператором equals-to =. Это не поддерживается; вместо этого вы должны использовать некоторую функцию пространственных отношений. Вероятно, лучше всего использовать минимальное расстояние:

select * 
     , public.ST_Y(al.address_position) latitude 
     , public.ST_X(al.address_position) longitude 
from address_list al 
where public.ST_DWithin(al.address_position::geography 
         , public.ST_SetSRID(public.ST_MakePoint(56.187339 --longitude$n 
                  , 57.964295 --latitude$n 
             ), 4326)::geography 
         , 1.); -- maximum separation of 1m 

В качестве альтернативы, вы можете сравнить значения широты и долготы для равенства:

select * 
     , public.ST_Y(al.address_position) latitude 
     , public.ST_X(al.address_position) longitude 
from address_list al 
where public.ST_X(al.address_position) = 56.187339 --longitude$n 
    and public.ST_Y(al.address_position) = 57.964295; --latitude$n 

Но вы должны иметь в виду, что сравнение значений с плавающей точкой для равенство обычно не является безопасным вариантом из-за представления в двоичной форме.

Answer 2

Я преобразовал столбец геометрии в WKT (хорошо известный текст) и сравните эти значения.

select * 
     , public.ST_Y(al.address_position) latitude 
     , public.ST_X(al.address_position) longitude 
     , public.ST_AsText(al.address_position) 
from address_list al 
where public.ST_AsText(al.address_position) = public.ST_AsText(public.ST_SetSRID(public.ST_MakePoint(56.187339 --longitude$n 
                            , 57.964295 --latitude$n 
                            ), 4326) 
                   ) 

И добавил индекс функции на колонке address_position

CREATE INDEX i_address_position_as_text 
    ON address_list 
    USING btree 
    (public.ST_AsText(address_position));