３次元の回転 (原点を通る任意方向の回転軸まわり)
― 座標系に依存しないベクトル表現と回転行列 ―

注意

• 対象読者：(３次元) ベクトルの内積・外積，２次元の回転行列を理解している人． (大学１年の物理数学レベル)
  (テンソル解析の初歩 (具体的にはエディントンのイプシロンを用いた３次元外積の表現) を知っていればなお良い．)
• 「オイラー角」や「ロール ピッチ ヨー (カルダン角)」等で検索して来た人は↓必読！
  ７０秒で分る、使える、四元数・４元数・クォータニオン・ Quaternionで回転 (産業技術総合研究所 中田 亨)
  ⇒ オイラー角から四元数への変換やその逆変換という問題 (背景説明篇)

  この問題はしばしば相談を受けるのですが、
  数学的理由と、数学の外の事情で難しいのです。

• ワールド座標 ⇔ カメラ座標 (視野座標) の変換式を求めようとしてここに迷い込んでくる人がたまにいるみたいだけど， ややこしいカメラの回転を考える必要はありません． そういう人はこっち．
• 画像表示を ON にしないと，図と数式が表示されません．

1. 問題
2. 解答
   1. 記号の定義
   2. ベクトルによる表現 (座標系非依存)
   3. 回転行列による表現 (正規直交座標系限定)
   4. 四元数 (クォータニオン) による表現 (研究中)
3. 参考・関連図書
4. サイト内関連ページ
5. 外部へのリンク
6. 更新履歴

３次元空間内において，原点を通る任意方向の回転軸のまわりに 点を角度θだけ回転させたときの位置を求めよ．

	O：原点 OA：回転軸 P：回転前の点の位置 P'：回転後の点の位置 Q：点 P を含む回転面と回転軸の交点 a：回転軸の方向・向きを表す単位ベクトル θ：回転角 回転方向は，座標系が右手系ならば a に対して右ネジの向き，左手系ならば左ネジの向きとする． p ≡ O→P p' ≡ O→P' u ≡ Q→P v ≡ a × u (u を +90°回転させたベクトル) 注意：右手系と左手系では外積の向きの定義が逆になる．

O→Q は単位ベクトル a に対する p の平行成分，u は垂直成分なので，

……………………………………… (2.2.1)
………………… (2.2.2)

v は定義より，

…………………(2.2.3)

a⊥u かつ |a|＝1 なので，|u|＝|v| である点に注意．

以上と u⊥v であることを考慮すると，回転後の位置 p' は２次元の回転の式を使って次のように書けることは明らか．

…………………(2.2.4)

ねっ，簡単でしょ？

なお，以上の計算式は座標系の取り方には依存しないので， 直交座標 (デカルト座標)でなくても (斜交座標や極座標 (円柱座標，球座標など) でも) 使用できる． ただし図にも描いたとおり，座標系の掌性 (chirality， キラリティ， カイラリティ：右手系と左手系の区別) により回転方向の定義が逆になる点に注意．
(そうなる理由は，右手系と左手系では外積の向き (図の v の向き) の定義が逆になるから．)

今までの計算式は座標系の取り方に依存しなかったが，ここから先は正規直交座標系限定．
(よく知られている外積の成分の計算式が正規直交座標系を前提としているため．)

ベクトルは縦ベクトル，^t を転置とする．例えば，

また行列 [a×] を，任意のベクトル a，p について次式が成り立つように定義する．
("[a×]" はここで勝手に定義する記号なので，よそを探してもないよ．)
(2008/08/17(日) 追記：外積行列 (cross-product matrix) というらしい．)

…………………(2.3.1)

つまり，

…………………(2.3.2)

これを用いると，(2.2.4) は次のように書ける．

…………………(2.3.3)
(I₃ は３次の単位行列)

したがって回転行列 R(a, θ) は次のようになる．

…………………(2.3.4)

成分で書けば次のとおり．

…………………(2.3.5)
(δ_ij は クロネッカーのデルタ，ε_ikj は エディントンのイプシロン．)

(余談：英語版 Wikipedia には，ε_ijk を３次元行列 (３次行列ではない) で表現した図がある．
また，ε_ijk と δ_ij の面白い関係式も載っている．)

なお，この行列を同次座標用に拡張したものは， OpenGL の glRotate (glRotatef, glRotatex) 関数でも用いられている． また DirectX の XMMatrixRotationAxis / XMMatrixRotationNormal (DirectXMath / XNA Math) や D3DXMatrixRotationAxis (D3DX) もこの行列を返すはず．

OpenCV では cv::Rodrigues()．(2015/12/26(土) 追記)

2008/06/18(水) 追記

ロドリゲスの回転公式 (Rodrigues' Rotation Formula) というらしい．

2008/08/17(日) 追記

エディントンのイプシロンを用いた式は， ３次元外積の成分を統一的に表現することができるので代数的に扱うには都合がよい． しかし数値計算 (CG を含む) では， 無駄に積和演算と多次元配列のアクセスを増やすので遅い． クロネッカーのデルタについても同様． sinθ [a×] と cosθ I₃ の非零要素だけを個別に加算するようハードコードする方がよい．

2010/02/27(土) 追記

ベクトルの回転 (物理のかぎしっぽ) より引用．

注)有限回転の式を，ロドリグの公式と呼びましたが， この式を導いたのはロドリグではないという説が有力です． この式自体は，実はもっと昔から知られていたようなのですが， ベクトルという概念がまだ無かったので，違った説明のされかたをされていました． ベクトルという概念を前面に押し出して，この公式を初めて導いたのはギブスのようです． こういった事情により，ロドリグの公式という言い方を避けて， ベクトルの回転公式などと呼ぶ人もいるようですが， 未だに定着している名前はありません．

(研究中)

楽天で買う 価格：4,830円（税込、送料別）

Androidゲームプログラミング A to Z posted with amazlet at 12.05.20 Mario Zechner インプレスジャパン 売り上げランキング: 26792 Amazon.co.jp で詳細を見る

---

---

---

---

楽天で買う 価格：3,900円（税込、送料別）

微分幾何学と接続―技術者のための (One Pointテキストシリーズ 11) posted with amazlet at 10.11.29 棚橋 隆彦 三恵社 売り上げランキング: 545133 Amazon.co.jp で詳細を見る

---

楽天で買う 価格：2,730円（税込、送料別）

C++による簡単実習 3次元CG入門 posted with amazlet at 11.09.03 小笠原 祐治 森北出版 売り上げランキング: 285345 Amazon.co.jp で詳細を見る

理論的な説明 ＋ サンプルプログラム．

---

楽天で買う 価格：2,940円（税込、送料別）

初めてのOpenGL ES posted with amazlet at 12.12.09 山下 武志 オライリージャパン 売り上げランキング: 85713 Amazon.co.jp で詳細を見る

---

---

---

---

---

---

---

---

楽天で買う 価格：1,575円（税込、送料別）

内積・外積・空間図形を通して ベクトルを深く理解しよう (数学のかんどころ　1) posted with amazlet at 11.11.03 飯高 茂 共立出版 売り上げランキング: 190055 Amazon.co.jp で詳細を見る

---

楽天で買う 価格：2,835円（税込、送料別）

なっとくするベクトル (なっとくシリーズ) posted with amazlet at 11.11.03 小野寺 嘉孝 講談社 売り上げランキング: 122943 Amazon.co.jp で詳細を見る

---

楽天で買う 価格：2,835円（税込、送料別）

なっとくする行列・ベクトル (なっとくシリーズ) posted with amazlet at 11.11.03 川久保 勝夫 講談社 売り上げランキング: 172378 Amazon.co.jp で詳細を見る

---

ベクトル・行列がビジュアルにわかる線形代数と幾何 posted with amazlet at 11.11.03 江見 圭司 江見 善一 共立出版 売り上げランキング: 158714 Amazon.co.jp で詳細を見る

---

楽天で買う 価格：2,100円（税込、送料別）

行列とベクトルのはなし―線形代数の基礎 posted with amazlet at 11.11.03 大村 平 日科技連出版社 売り上げランキング: 248023 Amazon.co.jp で詳細を見る

---

楽天で買う 価格：1,155円（税込、送料別）

細野真宏のベクトル〈平面図形〉が本当によくわかる本―数B (1週間集中講義シリーズ) posted with amazlet at 11.11.03 細野 真宏 小学館 売り上げランキング: 32875 Amazon.co.jp で詳細を見る

---

楽天で買う 価格：1,155円（税込、送料別）

細野真宏の ベクトル〈空間図形〉が本当によくわかる本 1週間集中講義シリーズ posted with amazlet at 11.11.03 細野 真宏 小学館 売り上げランキング: 8392 Amazon.co.jp で詳細を見る

---

楽天で買う 価格：1,365円（税込、送料別）

決定版 志田晶の ベクトルが面白いほどわかる本 posted with amazlet at 11.11.03 志田 晶 中経出版 売り上げランキング: 43835 Amazon.co.jp で詳細を見る

---

坂田アキラの三角関数・指数・対数が面白いほどわかる本 新装版―数学2対応 (数学が面白いほどわかるシリーズ) posted with amazlet at 11.11.03 坂田 アキラ 中経出版 売り上げランキング: 17311 Amazon.co.jp で詳細を見る

---

山本俊郎のベクトル原則編が面白いほどわかる本 (数学が面白いほどわかるシリーズ) posted with amazlet at 11.11.03 山本 俊郎 中経出版 売り上げランキング: 142524 Amazon.co.jp で詳細を見る

---

山本俊郎のベクトル実戦編が面白いほどわかる本 (数学が面白いほどわかるシリーズ) posted with amazlet at 11.11.03 山本 俊郎 中経出版 売り上げランキング: 331107 Amazon.co.jp で詳細を見る

---

イラスト・図解 はじめての行列とベクトル posted with amazlet at 11.11.03 長谷川 勝也 技術評論社 売り上げランキング: 220444 Amazon.co.jp で詳細を見る

---

楽天で買う 価格：788円（税込、送料別）

なるほど高校数学　ベクトルの物語 (ブルーバックス) posted with amazlet at 11.11.03 原岡 喜重 講談社 売り上げランキング: 229182 Amazon.co.jp で詳細を見る

---

楽天で買う 価格：788円（税込、送料別）

なるほど高校数学 三角関数の物語 (ブルーバックス) posted with amazlet at 11.11.03 原岡 喜重 講談社 売り上げランキング: 158574 Amazon.co.jp で詳細を見る

---

楽天で買う 価格：1,575円（税込、送料別）

いかにして問題をとくか posted with amazlet at 11.09.03 G. ポリア 丸善 売り上げランキング: 41 Amazon.co.jp で詳細を見る

• 幾何学・CG のアルゴリズム集
• ３点の座標から簡単に角度と回転方向を求める．(2・3・N次元，外積を用いる方法)
• 視野変換・射影変換 (ワールド座標系 ⇒ カメラ座標系 ⇒ 画面座標系)
• DirectX Graphics に関するメモ (Direct3D，DXGI ほか)
• OpenCV に関するメモ

• ３次元の回転
  • Rodrigues' rotation formula (Wikipedia)
  • Rodrigues' Rotation Formula (Wolfram MathWorld)
  • ベクトルの回転 (物理のかぎしっぽ)
  • 旧・加納裕のBLOGです
    • search : 任意軸周りの回転
    • 任意軸周りの回転
    • 任意軸周りの回転(2)
    • 任意軸周りの回転(3)
    • 任意軸周りの回転(4)
    • 任意軸周りの回転(5)
    • 任意軸周りの回転(6)
    • 任意軸周りの回転(7)
  • ７０秒で分る、使える、四元数・４元数・クォータニオン・ Quaternionで回転
    (産業技術総合研究所 中田 亨)
    • 回転篇
      • 四元数を使った回転の取り扱い手順
      • 四元数の長所と欠点
      • 四元数回転のソフトウエアのソースコード
    • 一般変位篇：回転以外の変形である拡大縮小（と併進？）の手順
    • 所望の四元数を作ろう篇：ある回転を実施する四元数の作り方
    • オイラー角から四元数への変換やその逆変換という問題（背景説明篇）

      この問題はしばしば相談を受けるのですが、
      数学的理由と、数学の外の事情で難しいのです。

      「オイラー角」や「ロール ピッチ ヨー (カルダン角)」等で検索して来た人は必読！

    • 四元数ムダ知識
  • 四元数を用いた三次元回転計算の幾何学的意味 (HMMNRST)
  • 北海道大学 講義資料
    ３次元の回転及び角速度： 数学的準備，シンプルローテーション，回転行列，オイラーパラメータ， ロドリゲスパラメータ，オイラー角，連続する回転操作，角速度の表現，演習ほか．
    • 宇宙機ダイナミクス
      • 第３回 (2007/06/12) 資料 (p.12〜)
      • 第４回 (2007/06/19) 資料
      • 第５回 (2007/06/26) 資料
    • ロボット情報学
      順運動学，逆運動学．
  • 人間−ロボット情報学2007 講義資料 (東北大学)
    ３次元回転姿勢と角速度に関する補足： ロドリゲスパラメータ，オイラーパラメータ，２×２複素行列ほか．
  • Mini Tips 1：任意軸周りの回転行列を求める (伊藤ソフトウェア研究所)
    回転行列の導出と，それを作成する Java のサンプルプログラム．
  • TMPSwiki
    • 3D空間における回転の表現形式 (TMPSwiki)
    • クォータニオン逆運動学
      
      • 対数クォータニオンを用いた数値的逆運動学法 (インバースキネマティクス，IK)
      • ３次元回転表現の相互変換 (回転行列，クォータニオン，オイラー角，etc.)
  • DOSEI日記
    • 回転行列と回転ベクトル
    • 回転行列の表現
    • 外積が行列で書けること
  • u_1rohのカタチ
    • クォータニオンの探検
    • クォータニオンとオイラーの公式に関する疑問
    • クォータニオンによる回転を美しく理解する
  • 四元数と3次元空間中の回転 (++C++;// 未確認飛行 C)
  • コラム/QuaternionPowers
    

    GameDev?.orgで2003年2月に投稿された“Quaternion Powers”（Sobeit Void氏著）の無許可な和訳です。

    元記事：http://www.gamedev.net/reference/articles/article1095.asp
  • 回転行列の表現方法
    3次元図形変換の行列表現：平行移動，拡大・縮小・反転，回転 (X/Y/Z 軸まわり，任意ベクトルまわり，オイラー角， ロール・ピッチ・ヨー，四元数)
  • 座標変換
    ２次元と３次元の座標変換 (ロボット用なので回転と平行移動のみ)， 同次変換，オイラー角．
  • 回転行列からオイラー角のパラメータ抽出を行う (It_lives_vainlyの日記)
  • 剛体の運動には必ず回転軸があることの証明をば (::Rukeのページ:: )
  • 回転を与える変換行列の解析 (島田 静雄，CAD・CGのための基礎数学)
    
    • 回転行列から回転軸の方向と回転角を求める．
    • 回転行列から座標軸まわりの回転角を求める．
    • 回転中心を求める．

    注意

    • FireFox では数式 (画像) が見えない．IE では見える．
      (画像のパス名の区切りに '/' ではなく '\' を使っているため．)
    • 一部しか読んでないけど，式 (4.41) は間違ってる．
  • 有顔ベクトルを用いた３次元の回転の表現
    • Inverse Kinematicsを用いたヒューマンフィギュアの動作決定に関する研究
      (筑波大学システム工学研究科 南城康之 (修論))
  • miscellaneous
    CG の各種アルゴリズムを ActionScript で実装したサンプルがある．
    • 回転行列から回転軸と回転角を求める
    • 回転行列の分解
  • 球座標の回転
  • 3 次元のアフィン行列を求める問題について (新潟工科大学 情報電子工学科 竹野茂治)
    「1 はじめに」から引用．

    以前機械系の方と、3 次元の合同アフィン変換 (回転 ＋ 平行移動からなる合同変換) の、回転と平行移動を求める方法に関して話を聞いたことがあり、 その方法が多少面白い方法であったので、ここにまとめておく。

    PDF ファイル

  • SO(3)の生成子 (Quod Erat Demonstrandum?)
• 右手系と左手系
  • OpenGL / DirectXにおける３Ｄ基礎概念の対比
    (The Fifth DIMENSION：実践 ３Ｄプログラミングの集い)
    OpenGL と DirectX の相違点．右/左手系，回転方向，行/列ベクトルなど．
    (Web魚拓：元のページは Infoseek isweb のサービス終了に伴い消滅．)
  • 座標系 (Direct3D 9) (MSDN)
    右手系の CG データを DirectX の (左手) 座標系に変換する方法．
    • ポリゴンの各三角形の頂点順序が (表側から見て) 右回りになるように並べ替える．
    • ビュー行列でＺ成分の符号を反転させる．
  • 左手座標系と右手座標系間の変換 (篤のメモ帳)
  • 右手座標系と左手座標系（クォータニオン）
  • OpenGL を Direct3D 互換で使う (ホイール欲しい ハンドル欲しい)
    • OpenGL を左手座標系で使用する方法
    • OpenGL と Direct3D の座標系の違い
    • OpenGL と Direct3D のシャドウマップ
  • 軸性ベクトル (擬ベクトル) に関する詳しい説明
    • ベクトルの外積　軸性ベクトルについて (OKWave)
    • NO.1825　線形な（正直な？）エッセイ（？）其の１　2010.1.4.　DDT
      (Weekend Mathematics コロキウム室 NO.251)
• Wikipedia
  • ドット積
  • クロス積
  • クロネッカーのデルタ
  • エディントンのイプシロン (Levi-Civita symbol)
    英語版には ε_ijk を３次元行列 (３次行列ではない) で表現した図がある．
    また，ε_ijk と δ_ij の面白い関係式も載っている．
  • テンソル
  • 四元数
  • Axis angle
  • オイラー角
  • 方位角 (Azimuth)
  • 仰俯角 (仰角・俯角，Altura)
  • 天球座標系
  • 経緯度
    • 緯度
    • 経度
    • ISO 6709
      

      ISO 6709は経緯度や高度など座標値による地点位置の表現についての国際標準規格である。 表題は Standard representation of geographic point location by coordinates， 国際規格分類コードICS 35.040 (文字集合と情報符号化)。

  • Euler-Rodrigues parameters
  • スピノール
  • Orientation entanglement (方位のもつれ)
• ベクトル解析 (物理のかぎしっぽ)
  • 四元数
  • 軸性ベクトルと極性ベクトル
  • 擬テンソル
• ベクトル・複素数・クォータニオン (PDF)
• Quaternion 〜いや、もう皆使えるだろ〜 (t-pot)
  • Quaternion(４元数)とは
  • 逆Quaternion
  • Quaternion と回転の関係
  • 任意軸回りの回転
  • 球面線型補間
  • Quaternion による回転の行列表記
  • プログラム (Direct3D)
• Quaternion (Wolfram MathWorld)
• 外積について (高校生のための微分幾何)
• ２次元ベクトルの外積の効用 (線形代数学の教科内容の改善に向けて)
  ２次元外積のわかりやすい解説．
• 固有ベクトル、行列の幾何学的解釈 (アクションゲーム大好き)
• 四次元への扉 (岡田好一ホームページ)
  四次元の超立体について．
• 座標変換とスピノール
  Ｎ次元の回転について．
• 高次元回転行列の補間とその応用
  (名古屋大学大学院情報科学研究科 高橋友和ほか， Visual Computing / グラフィクスと CAD 合同シンポジウム 2007)

このページの主な更新は Blog でお知らせします．

1. 2008/06/08(日) 作成開始
2. 2008/06/10(火) 公開
3. 2008/06/15(日) タイトル変更
4. 2008/06/18(水) ロドリゲスの回転公式について追記．
5. 2008/08/17(日) 数値計算に関する注意と外積行列に関する追記．
6. 2008/09/05(金) 参考・関連図書を追加．
7. 2009/02/05(木) 左手系の場合の図を追加．
8. 2010/02/27(土) 回転公式の名称について追記．
9. 2010/06/12(土) 参考・関連図書を追加．
10. 2011/11/03(木) 参考・関連図書を追加．
11. 2012/12/09(日)
    • 回転の計算式の座標系依存性に関する記述を追加．
    • XMMatrixRotation{Axis,Normal}() に関する記述を追加．
12. 2013/01/01(火) (本当の) ４次元の回転についてちょっと追記．