Почему MySQL находит точку внутри полигона, хотя она находится не в нем?

Question

[Глеб Николаев]

Собственно есть полигон и точка:

Делаю такой запрос и мне возвращается "yes" (якобы точка внутри полигона). Что я делаю не так?

set @r = st_GeomFromText('POLYGON((
37.22216838998539 55.84054891160275, 
37.21838662328132 55.83319286655621, 
37.241389247793016 55.825028405079834, 
37.23314950169925 55.81497747433474, 
37.2333211630762 55.81246433423204, 
37.23928639592532 55.78352732389286, 
37.33145779202875 55.769560201254166, 
37.3683649880737 55.7713020105184, 
37.39943569730221 55.78020337060204, 
37.44464477187915 55.80740526357401, 
37.454772793119375 55.82122732160703, 
37.465587459867436 55.83697652282129, 
37.459922634427976 55.83581729079571, 
37.45631774551197 55.83523766176832, 
37.43589004165455 55.836590116002256, 
37.42865282177775 55.83428693095061, 
37.42607790112345 55.83090550567636, 
37.39346223950235 55.84790625380606, 
37.387677439690705 55.85716977857572, 
37.36999631786453 55.868271531747475, 
37.33513310159018 55.880326311160644, 
37.31625035012533 55.885536584352174, 
37.28723957742026 55.876562899285055, 
37.22539627677173 55.85184165763322, 
37.22216838998539 55.84054891160275))');
set @p = st_GeomFromText('POINT(37.441008 55.872223)');
select if(mbrcontains(@r, @p), 'yes', 'no') as contain;

Answer 1

[Глеб Николаев]

Всем спасибо! Нашел решение. Это ST_Within (с версии mysql 5.6) Обращаю внимание что для ее правильной работы нужно поменять местами переменные из примера в моем вопросе т.е

select if(ST_Within(@p, @r), 'yes', 'no') as contain;

Comments

[xmoonlight]

а это разве не на плоскости?

[Глеб Николаев]

xmoonlight: полигон на плоскости, да. В данном конкретном случае это регион на карте Москвы

[xmoonlight]

Глеб Николаев: а полигон составлен на поверхности гео-сферы и затем преобразован или же просто как секущая сферу плоскость, проходящая через полярные координаты (lat,lon)? Т.е. вопрос у меня - в расхождении.

[Глеб Николаев]

xmoonlight: Полигон составлен на плоскости. Прям на самой плоской :) бумажной карте Москвы , по gps координатам. Это северо-западная часть Москвы+Красногорск

[xmoonlight]

Глеб Николаев: а-а... а я думал что переводили как-то так. :)

Answer 2

[xmoonlight]

MBR - как я понимаю, это лишь прямоугольная ограничивающая область для вписанного в неё полигона.
Т.е., иными словами, Вам не подходит такой вариант детектирования принадлежности точки полигону.

Comments

[Глеб Николаев]

Не подскажете какие-нибудь альтернативы?

[xmoonlight]

Глеб Николаев: да я сам обрадовался - думал ща скопирастю)))
Пока - не сталкивался, но вещь нужная.

[longclaps]

Глеб Николаев: xmoonlight: Когда-то писал такое ручками. Разбивал полигон на треугольнки (для невыпуклого это гемор), принадлежость точки треугольнику - просто.

[Глеб Николаев]

Нашел здесь функцию MBRWithin, говорят это поиск непосредственно в полигоне. Читаю сейчас мануал по ней

[xmoonlight]

longclaps: Глеб Николаев: ну раз на треугольники - то и concave (не выпуклый) уж сразу надо)
Знаю, что concave-фигуры "пилят" на треугольники через "обход" произвольного полигона, а потом их конкатенируют через рёбра в одну группу. (но там была 2D-графика, а здесь - GEO-координаты)

[longclaps]

xmoonlight: Да я для невыпуклого и решал. Отдельный гемор - когда разрез (при разделе на треугольники) проходит через точку (задачка была олимпиадная, кейсы всякие гнусные).