Основы машинной графики
Изменим теперь положение координатных осей левосторонней системы координат таким образом, чтобы ось OY была направлена в плоскость рисунка. Тогда положительная полуось OZ будет направлена горизонтально вправо, а положительная полуось OX - вертикально вверх (рис. 1.10).
Получить формулу вращения точки вокруг оси OY на угол можно заменив x на z, y на x в формуле двумерного поворота. При этом координата точки y при таком вращении не изменяется. В результате чего формула вращения точки (x,y,z) вокруг оси OY на угол
будет иметь следующий вид:
![\left \{ \begin{array}{ccc} x'=xcos\beta-ysin\beta \\ y'=y \\ z'=xsin\beta+zcos\beta \end{array}\right](/sites/default/files/tex_cache/aed07e2e2a5f7adc9a3fd03a725a773c.png)
![(x', y', z' ) = (x, y, z) \left( \begin{array}{ccc} cos\beta & 0 & sin\beta \\ 0 & 1 & 0 \\ -sin\beta & 0 & cos\beta \end{array} \right)](/sites/default/files/tex_cache/eb11c1c20281a23aa1b7f81e96aeac47.png)
Аналогично поступаем с осью вращения OX. Изменим положение координатных осей так, чтобы ось OX была направлена в плоскость рисунка, ось OY - горизонтально вправо, ось OZ - вертикально вверх (рис. 1.11).
Заменив в формуле двумерного поворота y на z, x на y, получим формулу вращения точки (x,y,z) вокруг оси OX на угол
![\left\{ \begin{array}{ccc} x'=x \\ y'=ycos\alpha-zsin\alpha \\ z'=ysin\alpha+zcos\alpha \end{array}\right](/sites/default/files/tex_cache/3856370bba05973a7affc4d118f3b27c.png)
![(x', y', z' ) = (x, y, z) \left( \begin{array}{ccc} 1 & 0 & 0 \\ 0 & cos\alpha & sin\alpha \\ 0 & -sin\alpha & cos\alpha \end{array} \right)](/sites/default/files/tex_cache/84307bb749670c792073660608d5de66.png)
Способ двумерного плоского вращения вокруг произвольной точки может быть обобщен на случай вращения вокруг произвольной оси трехмерного пространства. Пусть произвольная ось вращения задается вектором , причем
- точка, определяющая начало вектора, а
- конец вектора (рис. 1.12)
Вращение вокруг задаваемой оси (вектора ) на угол
выполняется в несколько этапов:
- Перенос вектора
так, чтобы начало вектора (точка
) совпала с началом системы координат. Это осуществляется с помощью операции сдвига T(-a,-b,-c) ;
- Поворот вокруг оси OY на угол
так, чтобы вектор (m,l,n) оказался в плоскости OYZ:
;
- Поворот вокруг оси OX на угол
так, чтобы вектор (m',l',n') совпал с осью OZ:
;
- Поворот вокруг оси OZ на заданный угол
;
- Выполнение преобразования, обратного, произведенному на шаге 3. Т.е. поворот вокруг оси OX на угол
;
- Выполнение преобразования, обратного, произведенному на шаге 2. Т.е. поворот вокруг оси OY на угол
;
- Выполнение преобразования, обратного, произведенному на шаге 1. Т.е. сдвиг на вектор (a,b,c): T(a,b,c)
Данный алгоритм вращения вокруг произвольной оси можно записать с помощью произведения серии элементарных матриц: , где V - исходная точка, V' - точка после поворота.
Остается определить чему равны углы поворотов и
(рис. 1.13).
Из простых тригонометрических соотношений можно получить следующие формулы:
![cos\beta=\frac{n}{\sqrt{m^{2}+n^{2}}}, cos\alpha=\frac{\sqrt{m^{2}+n^{2}}}{\sqrt{m^{2}+n^{2}+l^{2}}}](/sites/default/files/tex_cache/d54a8ee71c51c64d11ffe7c0c96761e0.png)
Как видно, операции трехмерного масштабирования и вращения могут быть реализованы с помощью умножения вектор-строки (точки) на матрицу преобразования. Операция же сдвига реализуется через сложение двух вектор-строк. Аналогично тому, как все двумерные преобразования (сдвиг, масштабирование и вращение) описываются матрицами размером 3х3 (через однородные координаты), трехмерные преобразования могут быть представлены в виде матриц размером 4х4. И тогда точка трехмерного пространства (x,y,z) записывается в однородных координатах как Wx,Wy,Wz,W, где . Если
, то для получения трехмерных декартовых координат точки (x,y,z) первые три однородные координаты нужно разделить на W. Отсюда следует, что две точки
и
в пространстве однородных координат описывают одну и ту же точку трехмерного пространства в том и только том случае, когда
для любой константы c не равной нулю. Таким образом, преобразование точки трехмерного пространства P=(x,y,z) в новую точку P'=(x',y',z') с использованием однородных координат можно записать как:
![(x', y', z', 1) = (x, y, z, 1) \left( \begin{array}{cccc} a_{1} & a_{2} & a_{3} & 0 \\ a_{4} & a_{5} & a_{6} & 0 \\ a_{7} & a_{8} & a_{9} & 0 \\a_{10} & a_{11} & a_{12} & 1 \\ \end{array} \right).](/sites/default/files/tex_cache/5606be65fd1b6e63a2d390946eddbf8f.png)
Уравнения трехмерного сдвига, масштабирования и вращения записываются в виде матриц преобразования однородных координат следующим образом:
![(x', y', z', 1) = (x, y, z, 1) \left( \begin{array}{cccc} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\Tx & Ty & Tz & 1 \\ \end{array} \right)=P*T(Tx,Ty,Tz)](/sites/default/files/tex_cache/bb87c407b703fd19427b97b9d308cd4f.png)
![(x', y', z', 1) = (x, y, z, 1) \left( \begin{array}{cccc} Sx & 0 & 0 & 0 \\ 0 & Sy & 0 & 0 \\ 0 & 0 & Sz & 0 \\0 & 0 & 0 & 1 \\ \end{array} \right)=P*S(Sx,Sy,Sz)](/sites/default/files/tex_cache/474b05d30cf430851293cf2f9e6838b2.png)
![(x', y', z', 1) = (x, y, z, 1) \left( \begin{array}{cccc} 1 & 0 & 0 & 0 \\ 0 & cos\alpha & sin\alpha & 0 \\ 0 & -sin\alpha & cos\alpha & 0 \\0 & 0 & 0 & 1 \\ \end{array} \right)=P*R_{x}(\alpha)](/sites/default/files/tex_cache/f0a808e87440e02a646a62d100aad360.png)
![(x', y', z', 1) = (x, y, z, 1) \left( \begin{array}{cccc} cos\beta & 0 & -sin\beta & 0 \\ 0 & 1 & 0 & 0 \\ sin\beta & 0 & cos\beta & 0 \\0 & 0 & 0
& 1 \\ \end{array} \right)=P*R_{y}(\beta)](/sites/default/files/tex_cache/d1c9cc04962a3418d7a72bb2f87d5d49.png)
![(x', y', z', 1) = (x, y, z, 1) \left( \begin{array}{cccc} cos\Upsilon & sin\Upsilon & 0 & 0 \\ -sin\Upsilon & cos\Upsilon & 0 & 0 \\ 0 & 0 & 1 & 0 \\0 & 0 & 0 & 1 \\ \end{array} \right)=P*R_{z}(\Upsilon),](/sites/default/files/tex_cache/b44f4f0067556b3b8eaae852a5fd33c1.png)
![R_{x}(\alpha),R_{y}(\beta),R_{z}(\Upsilon)](/sites/default/files/tex_cache/e5a76a9d896b1cb1f144398deade9950.png)
![\alpha, \beta;, \Upsilon;](/sites/default/files/tex_cache/c0ace837000fbd18119323c205167551.png)
Как и в двумерном случае, матричный подход позволяет совместить два или более элементарных преобразования в одно. Таким образом, последовательное применение двух преобразований и
может быть заменено применением одного преобразования T, причем матрица T будет равна произведению матриц преобразований
и
. Это легко можно увидеть на простом примере. Пусть точка (x,y,z) трансформируется в точку (x',y',z') с помощью преобразования
:
. Применяя затем преобразование
к точке (x',y',z'), получим точку
Теперь подставляя первое выражение во второе, получим:
. Причем порядок применения преобразований должен быть сохранен при перемножении соответствующих матриц.
Процесс вывода трехмерной графической информации по существу является более сложным, чем соответствующий двумерный процесс. Сложность, характерная для трехмерного случая, обуславливается тем, что поверхность вывода не имеет графического третьего измерения. Такое несоответствие между пространственными объектами и плоскими изображениями устраняется путем введения проекций, которые отображают трехмерные объекты на двумерной проекционной картинной плоскости. В процессе вывода трехмерной графической информации мы задаем видимый объем в мировом пространстве, проекцию на картинную плоскость и поле вывода на видовой поверхности. В общем случае объекты, определенные в трехмерном мировом пространстве, отсекаются по границам трехмерного видимого объема и после этого проецируются. То, что попадает в пределы окна, которое само является проекцией видимого объема на картинную плоскость, затем преобразуется в поле вывода и отображается на графическом устройстве. В общем случае операция проекции преобразует точки, заданные в системе координат размерности n, в точки системы координат размерности меньшей, чем n. В нашем случае точка трехмерного пространства отображается в двумерное пространство. Проекция трехмерного объекта строится при помощи прямых проецирующих лучей, которые называются проекторами и которые выходят из центра проекции, проходят через каждую точку объекта и, пересекая картинную плоскость, образуют проекцию. На рис. 1.14 представлены две различные проекции одного и того же отрезка и проекторы, проходящие через его конечные точки.
Определенный таким образом класс проекций известен под названием плоских геометрических проекций, т.к. проецирование осуществляется на плоскость, а не на искривленную поверхность и в качестве проекторов используют прямые линии. Плоские геометрические проекции можно подразделить на два основных класса: центральные (перспективные) и параллельные (ортогональные). Различие между ними определяется соотношением между центром проекции и проекционной плоскостью. Так, если расстояние между ними, конечно, то проекция будет центральной, если же оно бесконечно, то – параллельной. При описании центральной проекции мы явно задаем ее центр проекции, в то время как для параллельной проекции мы указываем лишь направление проецирования. Центр проекции порождает визуальный эффект, аналогичный тому, к которому приводят фотографические системы и используется в случаях, когда желательно достичь некоторой степени реализма. Следует заметить, что размер центральной проекции объекта изменяется обратно пропорционально расстоянию от центра проекции до объекта. Параллельная проекция порождает менее реалистичное изображение, т.к. отсутствует перспективное "укорачивание" объекта. Проекция фиксирует истинные размеры объекта, и параллельные линии остаются параллельными.
В общем случае задача получения центральной проекции заключается в том, чтобы определить проекцию точки объекта, расположенную в произвольном месте трехмерного пространства, на некоторую плоскость в этом же пространстве, называемую картинной. Нахождение центральной проекции является частным случаем задачи определения пересечения луча L с плоскостью в трехмерном пространстве (рис. 1.15)
В машинной графике задача вычисления центральной проекции, как правило, сильно упрощена. В данном случае центр проекции, который также называют точкой зрения, находится на одной из осей системы координат, картинная (проекционная) плоскость перпендикулярна оптической оси. Как правило, точку зрения (центр проекции) располагают на оси OZ, тогда картинная плоскость будет параллельна плоскости OXY системы координат (рис. 1.16).
В нашем случае точка C=(0,0,c) - центр проекции (положение наблюдателя), плоскость z=0 - картинная плоскость. Пусть точка P=(x,y,z) имеет проекцию P'=(x',y',0). Рассмотрим два подобных треугольника CPQ и CP'Q', и запишем отношение катетов: . Рассмотрим два других подобных треугольника CQ'O и CQB, и запишем отношения сторон для них:
. С другой стороны имеем:
. Так как OQ'=x', BQ=x, P'Q'=y', PQ=y имеем
![\left\{ \begin{array}{ccc} \frac{x'}{x}=\frac{c}{c-z} \\ \frac{y'}{y}=\frac{c}{c-z} \end{array}\right](/sites/default/files/tex_cache/c8d28626cd6023fe94705698bc381db3.png)
![\left\{ \begin{array}{ccc} x'=x\frac{c}{c-z}=x\frac{1}{1-\frac{z}{c}} \\ y'=y\frac{c}{c-z}=y\frac{1}{1-\frac{z}{c}} \end{array}\right](/sites/default/files/tex_cache/0923ab07d6afaefd5fa759512482dc34.png)
Если теперь c стремится к бесконечности, то получим формулу параллельной проекции: .
Следующим шагом необходимо спроецированное изображение перевести в координаты экрана. Это можно проделать следующим образом:
![\left\{ \begin{array}{ccc} X=X_{0}+x'l \\ Y=Y_{0}+y'l \end{array}\right,](/sites/default/files/tex_cache/fccae73f47cf0a13aa3bbb6d97f3d992.png)
![(X_{0},Y_{0})](/sites/default/files/tex_cache/a3d00ad0f1c2bb3a2213060fca4af88a.png)
Существует связь однородных координат с операцией центральной и параллельной проекциями, которая может быть выражена так: .
Для перехода от однородных координат к обычным, необходимо разделить все компоненты точки на четвертую координату: .
Для параллельной проекции матрица преобразования будет иметь вид: .
Таким образом, шаг проецирования можно описать в терминах матричной операции умножения. В результате этого мы можем объединить вместе операции преобразования объекта (сдвиг, масштабирование, вращение) и операцию проецирования в одну общую матрицу преобразования. Аналогично можно поступить с приведением спроецированных точек к экранным координатам:
![(X, Y, Z, H) = (x', y', z', 1) \left( \begin{array}{cccc} l & 0 & 0 & 0 \\ 0 & l & 0 & 0 \\ 0 & 0 & 0 & 0 \\X_{0} & Y_{0} & 0 & 1 \\ \end{array} \right)=(x'l+X_{0},y'l+Y_{0},0,1)](/sites/default/files/tex_cache/d042c72a2941243932c1ac4bb47ef007.png)
![\left\{ \begin{array}{ccc} X_{экр}=x'l+X_{0} \\ Y_{экр}=y'l+Y_{0} \end{array}\right](/sites/default/files/tex_cache/2770ed8eb7ef78c8a771c68f563411df.png)
Таким образом, все операции преобразования объекта трехмерного пространства на картинную плоскость (экран) можно описать в терминах матричных умножений.