Продолжаем соединять точки в контур

Article
09/11/2011

В предыдущем посте мы рассмотрели задачу превращения последовательности точек в LineString, а замкнутой ломаной в Polygon на основе WKT/WKB-представлений геометрических величин. Здесь мы рассмотрим решение той же задачи с использованием еще одного способа представления - GML.

Пусть имеем прежнюю таблицу точек #points (см. Скрипт 1 предыдущего поста)

select id, p.ToString() from #points

-------------------------------------------------------------------------------------

1 POINT (0 0)

2 POINT (1 1)

3 POINT (1 -1)

4 POINT (-1 -1)

5 POINT (-1 1)

Скрипт 1

из которых требуетсяполучить LINESTRING(0 0, 1 1, 1 -1, -1 -1, -1 1, 0 0). Наряду с текстовым (WKT) и бинарным (WKB)-представлениями для выражения геопространственных значений существует XML-вариант под названием GML (Geography Markup Language). Подобно методам STAsText() или AsTextZM(), он же ToString() в диалекте SQL Server, STGeomFromText() c его частностями STPointFromText(), STLineFromText() и др.; STAsBinary(), STGeomFromWKB() и его частностям; существуют методы для конвертации в/из GML: AsGml() / GeomFromGML().

Построение ломаной из таблицы точек с помощью GML

Рассмотрим GML-представление вышепомянутой ломаной:

declare @l geometry = geometry::STGeomFromText('LINESTRING(0 0, 1 1, 1 -1, -1 -1, -1 1)', 0)

select @l.AsGml()

-------------------------------------------------------------------------------------

</LineString>

Скрипт 2

Совершенно очевидно, как из таблицы Скрипт 1 получить XML Скрипт 2. Превращаем в каждой записи координаты точки в строку:

select id, p.ToString(), Format(p.STX, 'N8') + ' ' + Format(p.STY, 'N8') from #points

Рис.1

Складываем строки позаписьно

declare @s nvarchar(max) = (select ',' + Format(p.STX, 'N8') + ' ' + Format(p.STY, 'N8') from #points order by id for xml path(''))

select @s

Рис.2

И делаем требуемый XML по образу Скрипт 2:

declare @s nvarchar(max) = (select ' ' + Format(p.STX, 'N8') + ' ' + Format(p.STY, 'N8') from #points order by id for xml path(''))

declare @x xml = '<LineString xmlns="http://www.opengis.net/gml">

<posList>' + stuff(@s, 1, 1, '') + '</posList>

</LineString>'

select @x

Рис.3

который легким движением руки превращается в геометрию:

select geometry::GeomFromGml(@x, 0).ToString()

Рис.4

Построение многоугольника из таблицы точек с помощью GML

Совершенно аналогично вместо ломаной получить многоугольник, ограниченный создаваемым ею контуром. Посмотрим, как он выглядит в GML-представлении:

declare @l geometry = geometry::STGeomFromText('POLYGON((0 0, 1 1, 1 -1, -1 -1, -1 1, 0 0))', 0)

select @l.AsGml()

-------------------------------------------------------------------------------------

</LinearRing>

</exterior>

</Polygon>

Скрипт 3

и создадим из таблицы Скрипт 1 такой же. Кроме тэгов необходимо иметь в виду еще одно изменение. Контур многоугольника должен быть замкнут, т.е. заканчиваться той же точкой, что и начинался. Отсюда set @s +=… Остальное понятно.

declare @s nvarchar(max) = (select Format(p.STX, 'N8') + ' ' + Format(p.STY, 'N8') + ' ' from #points order by id for xml path(''))

set @s += (select top 1 Format(p.STX, 'N8') + ' ' + Format(p.STY, 'N8') from #points order by id)

declare @x xml = '<Polygon xmlns="http://www.opengis.net/gml">

</LinearRing>

</exterior>

</Polygon>'

select geometry::GeomFromGml(@x, 0).ToString()

Рис.5

Вариации способа

Получить начальную точку не запросом, а подстрокой

Можно взять начальную точку в конец не запросом из таблицы точек #points, а подстрокой из почти сформированной строки @s. Начальной точкой в ней будет все до второго пробела:

declare @s nvarchar(max) = (select Format(p.STX, 'N8') + ' ' + Format(p.STY, 'N8') + ' ' from #points order by id for xml path(''))

set @s += left(@s, charindex(' ', @s, charindex(' ', @s) + 1) - 1) --первая точка = левая часть до 2-го пробела

select @s

Рис.6

Получение многоугольника из замкнутой ломаной

Можно также взять готовый LineString (Рис.4) и превратить ее в Polygon, чтобы GML стал выглядеть, как в Скрипт 3. Жалко, в T-SQL нет возможности переименовывать узлы, только replace value of. Из-за этого XML проще построить заново, чем видоизменять структуру существующего.

declare @l geometry = geometry::STGeomFromText('LINESTRING(0 0, 1 1, 1 -1, -1 -1, -1 1)', 0)

declare @s varchar(max) = @l.AsGml().value('declare default element namespace "http://www.opengis.net/gml";

(LineString/posList)[1]', 'varchar(max)') + ' ' +

@l.STPointN(1).AsGml().value('declare default element namespace "http://www.opengis.net/gml";

(Point/pos)[1]', 'varchar(50)')

select geometry::GeomFromGml('<Polygon xmlns="http://www.opengis.net/gml">

</LinearRing>

</exterior>

</Polygon>', 0).ToString()

Рис.7

Здесь во втором операторе скрипта из LineString XML-запросом выбирается строчка с координатами и через пробел к ней добавляются координаты первой точки (STPointN(1)), которые тоже получаются XML-запросом (value()) из ее GML-представления (AsGml()). Затем строка с координатами вставляется в GML-описание многоугольника аналогично Рис.5, из которого методом GeomFromGml создается, собственно, геометрический объект.

Продолжение следует.

Алексей Шуленин

Продолжаем соединять точки в контур

Построение ломаной из таблицы точек с помощью GML

Построение многоугольника из таблицы точек с помощью GML

Вариации способа

Получить начальную точку не запросом, а подстрокой

Получение многоугольника из замкнутой ломаной

Additional resources