座標系の偏倍なんて嫌いだあ（T-T)

質問

はじめまして。kodemarimanと言います。
皆さんよろしくお願いします。

早速ですが、質問です。多分高校数学レベルの問題です。

右上の座標が(x1,y1)=(200,200)、左下の座標が(x2,y2)=(0,0)の正方形を考えます。

次の手順で座標を変換して、この正方形を描画すると、斜めになった長方形になります。
１：座標系を反時計周りに30度回転します。
２：座標系をx方向に1.2倍、y方向に1.5倍偏倍します。
３：正方形を書きます。

ですが、次の手順だと斜めになった長方形でなく、菱形になってしまいます。
１：座標系をx方向に1.2倍、y方向に1.5倍編倍します。
２：座標系を反時計周りに30度回転します。
３：正方形を書きます。

これって、なぜなんでなんでしょう？

目的としては、このような座標変換を繰り返し行っても
（例えば座標を回転して偏倍して平行移動して座標を回転して偏倍して平行移動して
座標を回転して偏倍して平行移動して・・・）
いつでも長方形（正方形）が描かれるような座標変換を行いたいのですが・・・

ちなみにテストに使ってるソフトはGhostScriptといいます。
これはポストスクリプトのプログラムを書くと画像にしてくれる
便利なソフトです。
↓URL
http://www.cs.wisc.edu/~ghost/index.htm

回答No.4

guiter · ベストアンサー

少し不親切でした。
n は座標変換の回数、1.2^n は siegmund さんの補足にあるようにｎ乗のことです。

以下、No.2の回答での行列の導き方です。
まず、初期状態として座標 O(0,0)、A(200,0)、B(200,200)、C(0,200)
を頂点とする正方形OABCを考えます。

ここで、n回座標変換を行なった後OAベクトルは
長さが 1.2 の n乗倍に、また 30*n[度]＝ nπ/6[rad]回転します。
すなわち、OA'ベクトル
　( 200*1.2^n*cos(nπ/6),200*1.2^n*sin(nπ/6) )
に写されます。
したがって、n回の座標変換を行なう行列をＡnとして
　┌ 200*1.2^n*cos(nπ/6) ┐　　　　┌200┐
　│　　　　　　　　　　　 │ ＝　Ａn│　　│　　　
　└ 200*1.2^n*sin(nπ/6) ┘　　　　└ 0 ┘
という式が成り立ちます。

同様にしてOCベクトルがOC'ベクトル
　( -200*1.5^n*sin(nπ/6),200*1.5^n*cos(nπ/6) )
に写されることから
　┌ -200*1.5^n*sin(nπ/6)┐　　　　┌ 0 ┐
　│　　　　　　　　　　　 │＝　Ａn│　　│　　　
　└ 200*1.5^n*cos(nπ/6)┘　　　　└200┘
という式が成り立ちます。
これらの式から行列ＡnがNo.2のように決まります。

また、変換後のOA'とOC'の内積をとると０になることから
変換後の２つのベクトルは垂直になっていることがわかりますね。
さらに、変換前
　OB = OA + OC
ですから変換後のOB'ベクトルは
　OB'= OA'+ OC'
というベクトルに写されます。
よって、四角形OA'B'C'は長方形になっています。
平行移動については行列で操作できないので最後にn回分まとめて動かすしかなさそうですね。

ところで、同じ辺を1.2倍し続けるのでいいのでしょうか？

補足

いえ、目的としては同じ辺を１．２倍しつづけるのではありません。

やりたいのは、ある座標系に描いた図形を、回転と拡大縮小の座標変換を行う事で
任意の大きさに拡大縮小して、かつ任意の角度に回転させることです。
また、そうして変換した座標系の中でさらに拡大や縮小、回転を行いたいのです。

で、どーもそういう作業（ポストスクリプトでプログラムを組んでいるんですが）をしていると、途中で図形がゆがんできてしまうので
困ってしまったんですよ。（＾＾；
(２回目の回転と拡大縮小で、なんかゆがんでしまう）

ところで、教えていただいたこの変換行列で試したところ、やぱり長方形
じゃなくて平行四辺形になってしまいました。理屈は納得できるのですが・・・
謎です。Ｇｈｏｓｔのバグなのかしら。（でもこれも長生きのソフトだからそれはないかも）
ほかに実験できそうなツール探して試してみようかな。