Как используя PostGIS вычислить пересечение 3-х окружностей?

[aquam]

Начальные данные: есть три окружности заданные точкой центра (широта и долгота) и радиусом в километрах.
Вопрос №1. Как вычислить полигоны окружностей и занести их в PostGIS?
Пока что нашёл функцию make_ellipse(cx float8, cy float8, a float8, b float8, slices float8),
но не знаю как вычислить a и b (у меня в километрах).

Вопрос №2. Как на основе сделанных полигонов окружностей вычислить полигон пересечения?

Спасибо!

[KolesovDmitry]

Начальные данные: есть три окружности заданные точкой центра (широта и долгота) и радиусом в километрах.
Вопрос №1. Как вычислить полигоны окружностей и занести их в PostGIS?
Пока что нашёл функцию make_ellipse(cx float8, cy float8, a float8, b float8, slices float8),
но не знаю как вычислить a и b (у меня в километрах).

Если у вас данные в широта/долгота, а радиус в километрах, то можно сначала перепроецировать ваши данные в прямоугольную систему координат (или как это там называется?) при помощи ST_Transform, а потом считать в километрах.

Вопрос №2. Как на основе сделанных полигонов окружностей вычислить полигон пересечения?

Попробуйте ST_Intersection.

[aquam]

Спасибо! Буду пробовать.

[aquam]

Можно поподробней c первым вопросом где перепроецирование в прямоугольную систему координат?
Получается мне нужно использовать что-то вида make_ellipse(37, 53, ST_Transform(радиус), ST_Transform(радиус), 18) для получения полигона-окружности с центром в точке (37, 53)?

[KolesovDmitry]

Можно поподробней c первым вопросом где перепроецирование в прямоугольную систему координат?
Получается мне нужно использовать что-то вида make_ellipse(37, 53, ST_Transform(радиус), ST_Transform(радиус), 18) для получения полигона-окружности с центром в точке (37, 53)?

Не, не так. Я бы попробовал что-то в этом роде:
* сначала создаете точку в системе координат широта/долгота
* потом перепроецируете эту точку при помощи ST_Transform в нужную систему координат
* в новой системе координат делаете хоть make_ellips, хоть буферные зоны вокруг нее.

[Denis Rykov]

Пример:

Код: Выделить всё

SELECT AsText(ST_Buffer(ST_Transform(ST_SetSRID(ST_MakePoint(37,53),4326),3395),1000000));

[aquam]

Спасибо!
Отличные советы!