F.12. cube
Этот модуль реализует тип данных cube для представления многомерных кубов.
F.12.1. Синтаксис
В Таблице F.7 показаны внешние представления типа cube. Буквы x, y и т. д. обозначают числа с плавающей точкой.
Таблица F.7. Внешние представления кубов
| Внешний синтаксис | Значение |
|---|---|
| Одномерная точка (или одномерный интервал нулевой длины) |
( | То же, что и выше |
| Точка в n-мерном пространстве, представленная внутри как куб нулевого объёма |
( | То же, что и выше |
( | Одномерный интервал, начинающийся в точке x и заканчивающийся в y, либо наоборот; порядок значения не имеет |
[( | То же, что и выше |
( | N-мерный куб, представленный парой диагонально противоположных углов |
[( | То же, что и выше |
В каком порядке вводятся противоположные углы куба, не имеет значения. Функции, принимающие тип cube, автоматически меняют углы местами, чтобы получить единое внутреннее представление «левый нижний — правый верхний». Когда эти углы совмещаются, в cube для экономии пространства хранится только один угол с флагом «является точкой».
Пробельные символы игнорируются, так что [( не отличается от x),(y)][ ( .x ), ( y ) ]
F.12.2. Точность
Значения хранятся внутри как 64-битные числа с плавающей точкой. Это значит, что числа с более чем 16 значащими цифрами будут усекаться.
F.12.3. Использование
В Таблице F.8 показаны операторы, предназначенные для работы с типом cube.
Таблица F.8. Операторы для кубов
| Оператор | Результат | Описание |
|---|---|---|
a = b | boolean | Кубы a и b идентичны. |
a && b | boolean | Кубы a и b пересекаются. |
a @> b | boolean | Куб a включает куб b. |
a <@ b | boolean | Куб a включён в куб b. |
a < b | boolean | Куб a меньше куба b. |
a <= b | boolean | Куб a меньше или равен кубу b. |
a > b | boolean | Куб a больше куба b. |
a >= b | boolean | Куб a больше или равен кубу b. |
a <> b | boolean | Куб a не равен кубу b. |
a -> n | float8 | Выдаёт n-ную координату куба (считая с 1). |
a ~> n | float8 | Выдаёт n-ную координату куба следующим образом: n = 2 * k - 1 обозначает нижнюю границу k-ой размерности, n = 2 * k обозначает верхнюю границу k-ой размерности. Этот оператор предназначен для поддержки KNN-GiST. |
a <-> b | float8 | Евклидово расстояние между a и b. |
a <#> b | float8 | Расстояние городских кварталов (метрика L-1) между a и b. |
a <=> b | float8 | Расстояние Чебышева (метрика L-inf) между a и b. |
(До версии PostgreSQL 8.2 операторы включения @> и <@ обозначались соответственно как @ и ~. Эти имена по-прежнему действуют, но считаются устаревшими и в конце концов будут упразднены. Заметьте, что старые имена произошли из соглашения, которому раньше следовали ключевые геометрические типы данных!)
Скалярные операторы упорядочивания (<, >= и т. д.) не имеют большого смысла ни для каких практических целей, кроме сортировки. Эти операторы сначала сравнивают первые координаты и если они равны, сравнивают вторые и т. д. Они предназначены в основном для поддержки класса операторов индекса-B-дерева для типа cube, который может быть полезен, например, если вы хотите создать ограничение UNIQUE для столбца типа cube.
Модуль cube также предоставляет класс операторов индекса GiST для значений cube. Индекс GiST для cube может применяться для поиска значений в выражениях с операторами =, &&, @> и <@ в предложениях WHERE.
GiST-индекс для cube может быть полезен и для поиска ближайших соседей с использованием операторов метрики <->, <#> и <=> в предложениях ORDER BY. Например, ближайшего соседа точки в трёхмерном пространстве (0.5, 0.5, 0.5) можно эффективно найти так:
SELECT c FROM test ORDER BY c <-> cube(array[0.5,0.5,0.5]) LIMIT 1;
Оператор ~> может также использоваться таким образом, чтобы эффективно выдавать первые несколько значений, отсортированных по выбранной координате. Например, чтобы получить первые несколько кубов, упорядоченных по возрастанию первой координаты (левого нижнего угла), можно использовать следующий запрос:
SELECT c FROM test ORDER BY c ~> 1 LIMIT 5;
А чтобы получить двумерные кубы, отсортированные по убыванию первой координаты правого верхнего угла:
SELECT c FROM test ORDER BY c ~> 3 DESC LIMIT 5;
В Таблице F.9 перечислены все доступные функции.
Таблица F.9. Функции для работы с кубами
| Функция | Результат | Описание | Пример |
|---|---|---|---|
cube(float8) | cube | Создаёт одномерный куб, у которого обе координаты равны. | cube(1) == '(1)' |
cube(float8, float8) | cube | Создаёт одномерный куб. | cube(1,2) == '(1),(2)' |
cube(float8[]) | cube | Создаёт куб нулевого объёма по координатам, определяемым массивом. | cube(ARRAY[1,2]) == '(1,2)' |
cube(float8[], float8[]) | cube | Создаёт куб с координатами правого верхнего и левого нижнего углов, определяемыми двумя массивами, которые должны быть одинаковой длины. | cube(ARRAY[1,2], ARRAY[3,4]) == '(1,2),(3,4)' |
cube(cube, float8) | cube | Создаёт новый куб, добавляя размерность к существующему кубу с одинаковым значением новой координаты для обеих углов. Это бывает полезно, когда нужно построить кубы поэтапно из вычисляемых значений. | cube('(1,2),(3,4)'::cube, 5) == '(1,2,5),(3,4,5)' |
cube(cube, float8, float8) | cube | Создаёт новый куб, добавляя размерность к существующему кубу. Это бывает полезно, когда нужно построить кубы поэтапно из вычисляемых значений. | cube('(1,2),(3,4)'::cube, 5, 6) == '(1,2,5),(3,4,6)' |
cube_dim(cube) | integer | Возвращает число размерностей куба. | cube_dim('(1,2),(3,4)') == '2' |
cube_ll_coord(cube, integer) | float8 | Возвращает значение n-ной координаты левого нижнего угла куба. | cube_ll_coord('(1,2),(3,4)', 2) == '2' |
cube_ur_coord(cube, integer) | float8 | Возвращает значение n-ной координаты правого верхнего угла куба. | cube_ur_coord('(1,2),(3,4)', 2) == '4' |
cube_is_point(cube) | boolean | Возвращает true, если куб является точкой, то есть если два определяющих его угла совпадают. | |
cube_distance(cube, cube) | float8 | Возвращает расстояние между двумя кубами. Если оба куба являются точками, вычисляется обычная функция расстояния. | |
cube_subset(cube, integer[]) | cube | Создаёт новый куб из существующего, используя список размерностей из массива. Может применяться для получения координат углов в одном измерении, для удаления измерений и изменения их порядка. | cube_subset(cube('(1,3,5),(6,7,8)'), ARRAY[2]) == '(3),(7)' cube_subset(cube('(1,3,5),(6,7,8)'), ARRAY[3,2,1,1]) == '(5,3,1,1),(8,7,6,6)' |
cube_union(cube, cube) | cube | Создаёт объединение двух кубов. | |
cube_inter(cube, cube) | cube | Создаёт пересечение двух кубов. | |
cube_enlarge(c cube, r double, n integer) | cube | Увеличивает размер куба на заданный радиус r как минимум в n измерениях. Если радиус отрицательный, куб, наоборот, уменьшается. Все определённые измерения изменяются на величину радиуса r. Координаты левого нижнего угла уменьшаются на r, а координаты правого верхнего увеличиваются на r. Если координата левого нижнего угла становится больше соответствующей координаты правого верхнего (это возможно, только когда r < 0), обоим координатам присваивается их среднее значение. Если n превышает число определённых измерений и куб увеличивается (r > 0), добавляются дополнительные размерности, недостающие до n; начальным значением для дополнительных координат считается ноль. Эта функция полезна для создания окружающих точку прямоугольников для поиска ближайших точек. | cube_enlarge('(1,2),(3,4)', 0.5, 3) == '(0.5,1.5,-0.5),(3.5,4.5,0.5)' |
F.12.4. Поведение по умолчанию
Я полагаю, что это объединение:
select cube_union('(0,5,2),(2,3,1)', '0');
cube_union
-------------------
(0, 0, 0),(2, 5, 2)
(1 row)не противоречит здравому смыслу, как и это пересечение
select cube_inter('(0,-1),(1,1)', '(-2),(2)');
cube_inter
-------------
(0, 0),(1, 0)
(1 row)Во всех бинарных операциях с кубами разных размерностей, я полагаю, что куб с меньшей размерностью является декартовой проекцией; то есть в опущенных в строковом представлении координатах предполагаются нули. Таким образом, показанные выше вызовы равнозначны следующим:
cube_union('(0,5,2),(2,3,1)','(0,0,0),(0,0,0)');
cube_inter('(0,-1),(1,1)','(-2,0),(2,0)');В следующем предикате включения применяется синтаксис точек, хотя фактически второй аргумент представляется внутри кубом. Этот синтаксис избавляет от необходимости определять отдельный тип точек и функции для предикатов (cube,point).
select cube_contains('(0,0),(1,1)', '0.5,0.5');
cube_contains
--------------
t
(1 row)F.12.5. Замечания
Примеры использования можно увидеть в регрессионном тесте sql/cube.sql.
Во избежание некорректного применения этого типа, число размерностей кубов искусственно ограничено значением 100. Если это ограничение вас не устраивает, его можно изменить в cubedata.h.
F.12.6. Благодарности
Первый автор: Джин Селков мл. <selkovjr@mcs.anl.gov>, Аргоннская национальная лаборатория, Отдел математики и компьютерных наук
Я очень благодарен в первую очередь профессору Джо Геллерштейну (https://dsf.berkeley.edu/jmh/) за пояснение сути GiST (http://gist.cs.berkeley.edu/) и его бывшему студенту, Энди Донгу, за пример, написанный для Illustra. Я также признателен всем разработчикам Postgres в настоящем и прошлом за возможность создать свой собственный мир и спокойно жить в нём. Ещё я хотел бы выразить признательность Аргоннской лаборатории и Министерству энергетики США за годы постоянной поддержки моих исследований в области баз данных.
Небольшие изменения в этот пакет внёс Бруно Вольф III <bruno@wolff.to> в августе/сентябре 2002 г. В том числе он перешёл от одинарной к двойной точности и добавил несколько новых функций.
Дополнительные изменения внёс Джошуа Рейх <josh@root.net> в июле 2006 г. В частности, он добавил cube(float8[], float8[]), подчистил код и перевёл его на протокол вызовов версии V1 с устаревшего протокола V0.
F.12. cube
This module implements a data type cube for representing multidimensional cubes.
F.12.1. Syntax
Table F.7 shows the valid external representations for the cube type. x, y, etc. denote floating-point numbers.
Table F.7. Cube External Representations
| External Syntax | Meaning |
|---|---|
| A one-dimensional point (or, zero-length one-dimensional interval) |
( | Same as above |
| A point in n-dimensional space, represented internally as a zero-volume cube |
( | Same as above |
( | A one-dimensional interval starting at x and ending at y or vice versa; the order does not matter |
[( | Same as above |
( | An n-dimensional cube represented by a pair of its diagonally opposite corners |
[( | Same as above |
It does not matter which order the opposite corners of a cube are entered in. The cube functions automatically swap values if needed to create a uniform “lower left — upper right” internal representation. When the corners coincide, cube stores only one corner along with an “is point” flag to avoid wasting space.
White space is ignored on input, so [( is the same as x),(y)][ ( . x ), ( y ) ]
F.12.2. Precision
Values are stored internally as 64-bit floating point numbers. This means that numbers with more than about 16 significant digits will be truncated.
F.12.3. Usage
Table F.8 shows the operators provided for type cube.
Table F.8. Cube Operators
| Operator | Result | Description |
|---|---|---|
a = b | boolean | The cubes a and b are identical. |
a && b | boolean | The cubes a and b overlap. |
a @> b | boolean | The cube a contains the cube b. |
a <@ b | boolean | The cube a is contained in the cube b. |
a < b | boolean | The cube a is less than the cube b. |
a <= b | boolean | The cube a is less than or equal to the cube b. |
a > b | boolean | The cube a is greater than the cube b. |
a >= b | boolean | The cube a is greater than or equal to the cube b. |
a <> b | boolean | The cube a is not equal to the cube b. |
a -> n | float8 | Get n-th coordinate of cube (counting from 1). |
a ~> n | float8 | Get n-th coordinate of cube in following way: n = 2 * k - 1 means lower bound of k-th dimension, n = 2 * k means upper bound of k-th dimension. This operator is designed for KNN-GiST support. |
a <-> b | float8 | Euclidean distance between a and b. |
a <#> b | float8 | Taxicab (L-1 metric) distance between a and b. |
a <=> b | float8 | Chebyshev (L-inf metric) distance between a and b. |
(Before PostgreSQL 8.2, the containment operators @> and <@ were respectively called @ and ~. These names are still available, but are deprecated and will eventually be retired. Notice that the old names are reversed from the convention formerly followed by the core geometric data types!)
The scalar ordering operators (<, >=, etc) do not make a lot of sense for any practical purpose but sorting. These operators first compare the first coordinates, and if those are equal, compare the second coordinates, etc. They exist mainly to support the b-tree index operator class for cube, which can be useful for example if you would like a UNIQUE constraint on a cube column.
The cube module also provides a GiST index operator class for cube values. A cube GiST index can be used to search for values using the =, &&, @>, and <@ operators in WHERE clauses.
In addition, a cube GiST index can be used to find nearest neighbors using the metric operators <->, <#>, and <=> in ORDER BY clauses. For example, the nearest neighbor of the 3-D point (0.5, 0.5, 0.5) could be found efficiently with:
SELECT c FROM test ORDER BY c <-> cube(array[0.5,0.5,0.5]) LIMIT 1;
The ~> operator can also be used in this way to efficiently retrieve the first few values sorted by a selected coordinate. For example, to get the first few cubes ordered by the first coordinate (lower left corner) ascending one could use the following query:
SELECT c FROM test ORDER BY c ~> 1 LIMIT 5;
And to get 2-D cubes ordered by the first coordinate of the upper right corner descending:
SELECT c FROM test ORDER BY c ~> 3 DESC LIMIT 5;
Table F.9 shows the available functions.
Table F.9. Cube Functions
| Function | Result | Description | Example |
|---|---|---|---|
cube(float8) | cube | Makes a one dimensional cube with both coordinates the same. | cube(1) == '(1)' |
cube(float8, float8) | cube | Makes a one dimensional cube. | cube(1,2) == '(1),(2)' |
cube(float8[]) | cube | Makes a zero-volume cube using the coordinates defined by the array. | cube(ARRAY[1,2]) == '(1,2)' |
cube(float8[], float8[]) | cube | Makes a cube with upper right and lower left coordinates as defined by the two arrays, which must be of the same length. | cube(ARRAY[1,2], ARRAY[3,4]) == '(1,2),(3,4)' |
cube(cube, float8) | cube | Makes a new cube by adding a dimension on to an existing cube, with the same values for both endpoints of the new coordinate. This is useful for building cubes piece by piece from calculated values. | cube('(1,2),(3,4)'::cube, 5) == '(1,2,5),(3,4,5)' |
cube(cube, float8, float8) | cube | Makes a new cube by adding a dimension on to an existing cube. This is useful for building cubes piece by piece from calculated values. | cube('(1,2),(3,4)'::cube, 5, 6) == '(1,2,5),(3,4,6)' |
cube_dim(cube) | integer | Returns the number of dimensions of the cube. | cube_dim('(1,2),(3,4)') == '2' |
cube_ll_coord(cube, integer) | float8 | Returns the n-th coordinate value for the lower left corner of the cube. | cube_ll_coord('(1,2),(3,4)', 2) == '2' |
cube_ur_coord(cube, integer) | float8 | Returns the n-th coordinate value for the upper right corner of the cube. | cube_ur_coord('(1,2),(3,4)', 2) == '4' |
cube_is_point(cube) | boolean | Returns true if the cube is a point, that is, the two defining corners are the same. | |
cube_distance(cube, cube) | float8 | Returns the distance between two cubes. If both cubes are points, this is the normal distance function. | |
cube_subset(cube, integer[]) | cube | Makes a new cube from an existing cube, using a list of dimension indexes from an array. Can be used to extract the endpoints of a single dimension, or to drop dimensions, or to reorder them as desired. | cube_subset(cube('(1,3,5),(6,7,8)'), ARRAY[2]) == '(3),(7)' cube_subset(cube('(1,3,5),(6,7,8)'), ARRAY[3,2,1,1]) == '(5,3,1,1),(8,7,6,6)' |
cube_union(cube, cube) | cube | Produces the union of two cubes. | |
cube_inter(cube, cube) | cube | Produces the intersection of two cubes. | |
cube_enlarge(c cube, r double, n integer) | cube | Increases the size of the cube by the specified radius r in at least n dimensions. If the radius is negative the cube is shrunk instead. All defined dimensions are changed by the radius r. Lower-left coordinates are decreased by r and upper-right coordinates are increased by r. If a lower-left coordinate is increased to more than the corresponding upper-right coordinate (this can only happen when r < 0) than both coordinates are set to their average. If n is greater than the number of defined dimensions and the cube is being enlarged (r > 0), then extra dimensions are added to make n altogether; 0 is used as the initial value for the extra coordinates. This function is useful for creating bounding boxes around a point for searching for nearby points. | cube_enlarge('(1,2),(3,4)', 0.5, 3) == '(0.5,1.5,-0.5),(3.5,4.5,0.5)' |
F.12.4. Defaults
I believe this union:
select cube_union('(0,5,2),(2,3,1)', '0');
cube_union
-------------------
(0, 0, 0),(2, 5, 2)
(1 row)
does not contradict common sense, neither does the intersection
select cube_inter('(0,-1),(1,1)', '(-2),(2)');
cube_inter
-------------
(0, 0),(1, 0)
(1 row)
In all binary operations on differently-dimensioned cubes, I assume the lower-dimensional one to be a Cartesian projection, i. e., having zeroes in place of coordinates omitted in the string representation. The above examples are equivalent to:
cube_union('(0,5,2),(2,3,1)','(0,0,0),(0,0,0)');
cube_inter('(0,-1),(1,1)','(-2,0),(2,0)');
The following containment predicate uses the point syntax, while in fact the second argument is internally represented by a box. This syntax makes it unnecessary to define a separate point type and functions for (box,point) predicates.
select cube_contains('(0,0),(1,1)', '0.5,0.5');
cube_contains
--------------
t
(1 row)
F.12.5. Notes
For examples of usage, see the regression test sql/cube.sql.
To make it harder for people to break things, there is a limit of 100 on the number of dimensions of cubes. This is set in cubedata.h if you need something bigger.
F.12.6. Credits
Original author: Gene Selkov, Jr. <selkovjr@mcs.anl.gov>, Mathematics and Computer Science Division, Argonne National Laboratory.
My thanks are primarily to Prof. Joe Hellerstein (https://dsf.berkeley.edu/jmh/) for elucidating the gist of the GiST (http://gist.cs.berkeley.edu/), and to his former student Andy Dong for his example written for Illustra. I am also grateful to all Postgres developers, present and past, for enabling myself to create my own world and live undisturbed in it. And I would like to acknowledge my gratitude to Argonne Lab and to the U.S. Department of Energy for the years of faithful support of my database research.
Minor updates to this package were made by Bruno Wolff III <bruno@wolff.to> in August/September of 2002. These include changing the precision from single precision to double precision and adding some new functions.
Additional updates were made by Joshua Reich <josh@root.net> in July 2006. These include cube(float8[], float8[]) and cleaning up the code to use the V1 call protocol instead of the deprecated V0 protocol.