オウムアムアはエイリアンの探査機か？

　 １１月１８日付でオウムアムアに関するブログを書いた後、新たな新聞記事も出たので、現段階でオウムアムアとは何者なのかについて整理してみることにした（オウマウマウは私の読み違いらしいので、訂正します）。

　すでに記載した通り、２０１７年１０月に、太陽系外からやってきたとみられる天体が観測され、オウムアムアと名付けられて、世界中の天文学関係者の興味を引くことになった。

　この天体が何者なのかについて、二つの見解が対立している。

　一つはＮＡＳＡが主張するもので、天体は彗星であり、太陽光によって温められ、ガスを噴き出して加速したと結論づけた。

　もう一つは、米ハーバード・スミソニアン天体物理学センターの研究者らが唱えるもので、天体は、太陽光の圧力を受けて進む宇宙帆船であるという。

　両見解から推察できることは、この天体は、通常の小惑星などとは異なり、太陽系内で加速されているという観測データが得られているということである。この天体が太陽系に侵入するときの初速度は２６ｋｍ／ｓ程度で、特に注目するような速度ではないようだ。

　ハーバードの研究者は、この天体が彗星とする説は、観測データと整合するような説明が難しいとしている。

　研究者は、この天体は、太陽光を受けて推進できることから、薄さ０．３～０．９ｍｍの曲面をもつシートでできた帆船であると推定している。

　この天体が星間を航行している間に、微細な星間物質が付着して小惑星のような外観になったという。

　研究者によれば、オウムアムアは、薄いシート曲面をもつ星間物質を代表するものであり、まだ未知のプロセスによって星間あるいは原始惑星円盤内で自然に生成されたものか、宇宙人によって構築された宇宙帆船の可能性があるという。

　もしそれが宇宙人によってつくられたものであるならば、オウムアムアは、宇宙帆船の一つである可能性がある。宇宙帆船は、恒星系内に存在する惑星の間の物品輸送に利用されるであろう。そうであれば、不要となった宇宙帆船が廃船として宇宙空間に放出されても不思議ではないだろう。

　ＪＡＸＡが打ち上げたＩＫＡＲＯＳは、そのような宇宙帆船が可能であることを実証するテスト機であり、燃料なしで宇宙空間を移動できる。

　ハーバードの研究者は、ＩＫＡＲＯＳから多くの示唆を受けたことは疑いない。

　しかしながら、オウムアムアがエイリアンの探査機と考えるのは早計と考える。オウムアムアから発信された電波らしいものが検出されていないので、探査機の可能性は低い。また、オウムアムアは宇宙人由来の廃船と考えるのもＳＦめいていて信じ難いものである。

　ハーバードの宇宙帆船説は、仮説を提唱しただけで、誰もこの仮説を実証することができない状態にあるので、まだ科学の領域には入っていないものと考える。

　少なくとも、コンピュータ・シミュレーションによって、理論上のオウムアムア・モデルの加速性能が観測データと整合することを示して欲しいと思うのである。

　また、研究者らは、宇宙空間にはオウムアムアのような宇宙帆船またはその廃船が多数放浪していると考えているので、オウムアムアに続く第二の太陽系外天体の観測が期待される。

　参考文献
　２０１８年１１月２８日付の朝日新聞夕刊の記事のほかに、
　Could Solar Radiation Pressure Explain ‘Oumuamua’s Ｐｅｃｕｌｉａｒ Acceleration ?, The Astrophysical Journal Letters, 868 No.1, 2018 Nov.20

画像データに関する主成分分析

　多数の顔写真が登録してあるとする。ある人の顔の画像を得たとき、この人の顔写真が登録されているか否かをチェックすることによって本人の認証を行うシステムがある。

　たとえば多数の顔写真｛ｆi｝，ｉ＝１，．．．，Ｎがあって、顔の大きさがほぼ一定になるように調節してあるとする。画像ｆiはｎ×ｎ画像のｎ＾２個の画素に１からｎ＾２までの通し番号をつけ、その値を縦に並べたｎ＾２次元列ベクトルとする。

　｛ｆi｝の平均をとれば、「平均顔」の画像ｇ0が得られる。各画像からｇ0を引いた
ｆi’＝ｆi－ｇ0
は、平均顔からの差異を表す画像である。そのような代表的差画像の第１次近似をｇ1とする。さらに差画像の第２次近似をｇ2とする。これを続けて差画像の第３次以下の近似ｇ3，ｇ4，．．．を作る。

　このとき、ｇ1，ｇ2，ｇ3，．．．が正規直交系であれば、画像ｆi’をｇ1，ｇ2，ｇ3，．．．の一次結合で表現することができ、次のように表記できる。
　　　ｆi＝ｇ0＋ｃi1ｇ1＋ｃi2ｇ2＋ｃi3ｇ3＋．．．　　　（１）

　これは、個々の顔をまず平均顔で近似し、次に主たる特徴を補正し、さらに二次的な特徴を補正し、以下順に高次の特徴を補正しているとみなせる。

　正規直交基底ｇ1，ｇ2，ｇ3，．．．を求めるために、統計学で主成分分析とよばれる技法を用いる。

　ｆi’がｎ個の画素から成る列ベクトルとすると、｛ｆi｝についてｎ×ｎの分散・共分散行列を計算する。次にこの行列についてｎ個の固有値を計算し、各固有値に対する固有ベクトルを求める。このｎ個の固有ベクトルが正規直交基底ｇ1，ｇ2，．．．，ｇnとなる。

　簡単のために画像は白黒画像とし、画像はｎ＝４，８画素から成るものとし、画像数Ｎ＝３２とし、分散・共分散行列を計算し、固有値を求めた。ただし、ｇ0が０ベクトルとなるように画像標本を作成した。

　ｎ＝４の場合の計算結果は、分散・共分散行列が単位行列となるので、ただ一通りの固有値＝１しか得られない。

　標本数Ｎを非常に大きな数にすれば、データはランダム化されるので、ｇ0が０ベクトルとなり、各分散が１、すべての共分散が０になる結果、分散・共分散行列が単位行列になることが予想される。

　このように、画像一般を対象としたとき、ｎ個の画素は同等のチャネルとなり、同じカテゴリとなるので、特徴による差別化ができず、主成分分析の技法で扱うのは難しいことを知る。

　画素数が４，８程度であれば、画像ベクトルの総当たりで内積を計算することにより、各々４個と８個の固有ベクトルを見つけることができる。

　ｎ＝４の場合に、見つけた固有ベクトルの組を行ベクトルの形式で表現すると、次のようになる。
　　　（１１１１）；（１１－１－１）；（１－１１－１）；（１－１－１１）

　また、次の例も固有ベクトルの組であるので、一組のみとは限らない。
　　　（１１１－１）；（１１－１１）；（１－１１１）；（１－１－１１）

　ｎ＝８の場合には、ｎ＝４で見つけた最初の固有ベクトルの組を利用して、次のような固有ベクトルの組を作成できる。
　　　（１１１１１１１１）；（１１－１－１１１－１－１）；（１－１１－１１－１１－１）；（１－１－１１１－１－１１）；（１１１１－１－１－１－１）；（１１－１－１－１－１１１）；（１－１１－１－１１－１１）；（１－１－１１－１１１－１）

　ｎ＝４，８のいずれの場合についても、任意の画像ベクトルは、各固有ベクトルに重みを掛けたものの一次結合で表現できる。従って、すべての画像ベクトルは、原理的には（１）式のように表現できるものと予想される。

　何のことはない。結果的には１０月２８日付のブログで作成を試みた意味のない基底画像の段階に戻っていたのだ。

　冒頭の顔画像の識別問題について言及すれば、顔画像の中で識別のキーとなるような特徴領域のカテゴリを何か所か選び、各特徴領域を構成する画素パターンをまとめて一つの成分として分析するのではなかろうか。いずれにしても特徴領域の抽出には、事前の統計処理などそれなりのスキルとノウハウが必要になるのであろう。

　こうなると、単純な「平均顔」の画像をデザインし、いくつかの特徴領域のカテゴリを選択し、各特徴領域のデータをランダム化して、１０００個程度の変形顔の画像をつくってみたくなる。ここでは、主成分分析の手法を逆用し、単純かつ人為的な分散・共分散行列を作成することになるだろう。変形顔の作成に当たり、遺伝的アルゴリズムが使えるのか否かも検討対象となる。

　参考文献
　金谷健一著「インターネット時代の数学―重ね合わせの原理と応用」（共立出版）
　長谷川勝也著「ホントにわかる多変量解析」（共立出版）