無印のアクリルフレームが捗る件と注意点

twitterでフォローしてる人が無印のアクリルフレームにアイカツ劇場版のスリーピングオーロラコーデの4枚を入れて飾っていて、いいなと思って真似したんだけど、うまくいかなかったこともあったのでそれも含めて、アイカツカードを飾るのって楽しいということを書こうと思う。

最初は4枚を横に並べるって構想で無印のアクリルフレームのパノラマを
買ってきた。無印のサイトから引っ張ってきたそのパノラマフレームがこれ。

これで4枚を横に並べて、データカードダスの筐体にあるサンプルみたいにして、パソコンのディスプレイの前に飾れば、作業中に目線を下に向けるだけでカードが目に入って机の上がすごい楽しくなるなって思って買った。ただ、買って入れようとして気づいたんだけど、これ端でネジを締めることでフレームの間に固定するから、端の2枚のせいでフレームが浮いちゃって真ん中の2枚は固定できない。ずっと入らない入らないって思ってたけど端にカード入れているからフレームが弧を描くようになっちゃってた。無印に電話したらやっぱりこの商品は1枚のパノラマサイズの写真とか用なのでそういうカード4枚的なのは無理ということで返品した。

で代わりに買ったのがA5サイズ。

横にしてちょうど4枚がきれいに収まった。参考にしたフォローしてる人ははがきサイズのフレーム使っていたんだけど、それだとスペースの問題で4枚を一部重ねてディスプレイしてて、重ねるってのが嫌で全部きれいに飾れるA5サイズにした。

やっぱりこうやってトップス・ボトムス・シューズ・アクセサリーの4枚をきれいに飾るとすごい楽しいし、きれい。アクリルフレームだから裏の説明も裏返せばすぐに読むことができる。映画は3回見に行って1セットは友人からもらって計4セットもあってどうしようと思ってたけど、あかりちゃん大好きだしこうやって飾れるってので有効活用できてよかった。PRも結構ダブったりするからトレード用とかじゃなければこうやって飾ったりシャドーアート（シャドーボックス）とかで数をつかって工作したりとかで有効活用したい。

月食でのターコイズフリンジの原因

天文学会に行って、そこであった高校生が発表するジュニアセッションにも行った。そのジュニアセッションのポスターで月食の本影の縁が青くなる、ターコイズフリンジ（Turquoise fringe）というのを測って検出しようとしたというところがあって結構面白かった。この現象自体知らなかったからそもそもリアルなのかってところからだったが、星ナビにも載っているしリアルらしい。そこでの彼らの説明だと成層圏を通っているから青くなるという説明だった。いやただやっぱりそこで思うのは月食が赤くなるのって青い光が大気中で吸収というか散乱されて赤い光が屈折することなんだからそこで青い光だけが直進して外側で縁になるの？って感じだった。まあそこでオゾン層が発光してるんじゃないかとかまた適当なこと言っちゃったんですが気になったので調べました。

結論から言うと高校生の彼らの言っていることが正しくて、成層圏中のオゾンによる吸収で青くなる。NASAのページでまさにそのページがあった。
http://science.nasa.gov/science-news/science-at-nasa/2008/13feb_lunareclipse/

The source of the turquoise is ozone. Eclipse researcher Dr. Richard Keen of the University of Colorado explains: "During a lunar eclipse, most of the light illuminating the moon passes through the stratosphere where it is reddened by scattering. However, light passing through the upper stratosphere penetrates the ozone layer, which absorbs red light and actually makes the passing light ray bluer." This can be seen, he says, as a soft blue fringe around the red core of Earth's shadow.

 英語版Wikipediaの成層圏にあったこの写真のまさにこの白い部分、が成層圏。でオレンジの部分が対流圏。青いところはさらに上の中間圏だそう。

オゾン層は10-50kmぐらいだそうで、ターコイズフリンジがそれに対応していると考えるとじゃあ赤い光はどこを通っているのか。上の写真の赤い部分がそれだと考えられるけど定量的には？

月食中の月のスペクトルの分光をした研究があってそれで月食中に観測した月のある点に地球大気のどこを通ってきたのかを計算した結果だと下の図になる。
Yan et al. 2014
High resolution transmission spectrum of the Earth's atmosphere -- Seeing Earth as an exoplanet using a lunar eclipse
http://arxiv.org/abs/1405.4780

A点が4.2kmでB点が17.8km。つまりオゾン層の下の対流圏。

なのでやっぱり月食中が赤いのは地球大気下層、高度数kmから十数kmのところを通ってきているから（だから火山とかの影響も受ける）。赤い光は屈折して影の中にも光が入る。影のごく外側は高度の高いところを通ってきていてオゾン層で青くなる。高くなると空気の密度が小さくなって屈折率も小さくなるからほとんど直進して縁になるってことかな。逆にいうと影のどこを見るかで地球大気のどの高度を見るかに対応しているってわけだな。

あと補足だけど、なんで月食中の月を分光して観測しているかというと将来的に地球のような惑星が見つかってその大気を観測する時は月食のように、惑星をかすめてくる大気を調べることになる。その時にその地球のような惑星が本当に地球と同じようなのかって確かめるために月食で同じように観測して地球のテンプレを作っているってわけです。

院試対策投稿に対する今のコメントとか

なんかGoogleで「院試　対策」とかググると
http://tgrr-putan.blogspot.jp/2012/09/blog-post_13.html
が結構上に来て、見ている人も多いらしい。
この間の天文台のスプリングスクールでも見ましたって言ってたのがいたし、というか後輩でもいた。これから院試に向けてたぶんプレビュー数また増えるし、院試勉強どうしたらいいですかってのよく聞くんで一応、またコメントしときます。

あの投稿での方法はあくまでも付け焼刃的方法で、あれが参考になるとか言ってるやつはどうしようもないので、ちゃんと勉強しましょうということです。というかこのことについても記事の中で書いてるはずだけど。普通にこれまでの講義や演習のノートを見返して問題を解き直す、演習書の問題を解くということをしてください。

アイカツ！LIVE☆イリュージョンに行ってきた

日曜日に行ってきたけど本当によかった。もう一度見たい。やって欲しい。もうこれに尽きるんだけど、だらだらと自分が書いて思い出してまた読み返した時に思い返したいというのも考えて書いておく。

最初の芸能人はカードが命というらしい、アイカツ！システムで変身する曲が流れてカウントダウンが始まったあたりですでに泣いた。なんかアイカツ世界へ行けるんだなっていうワクワク感がこの上なくあって体が震えてたと思う。で、ステージ真ん中の扉が開いてアイドルが出ていて踊ってて本当に感動した。なんかありえないというか、あのテレビとかで映ってるあのステージが目の前にあってテレビで映ってるあの観客がここだと思うと本当に夢のようだった。

あと、色が無線制御で変わるペンライト！　すごい。自前でサイリウム持ち込んでた人いたけど正直邪魔だった。あれ本当に揃った会場見てみたい。たぶんファミリータイムだとそうなってたのかな。サイリウムを扇状にしてジュリアナみたいに振ってたオタク死すべし。曲に合わせてペンライト変わってたし、しかも手拍子とかそういうレベルの早さで。

だから一番楽しかったのはあかりちゃんのペンライトの使い方講座。あかりちゃんが「右、左、右、左って振ってね」ってので右左に振るとペンライトの色が右と左で色が変わったってのが楽しい。あかりちゃんめっちゃかわいかった。2015シーズン本当に楽しみです。

後ろの大学生っぽい二人組とかは「あおい姐さん」とか叫んでて、自分も「あかりちゃん応援してるよー」ってあかりちゃん出てくるときに言おうかと思ったら、もう感極まってて声が出なかった。

アンコールのSHINING LINE*、フルバージョンだったし、あかりちゃんも出てきたしで本当によかった。泣いてた。本当にあかりちゃん応援してます。

なんか本当に夢のようで、曲名とかもライブブックみないと思い出せないんだよね。ただ最後のいちごちゃんとあかりちゃんが挨拶するところで、ああバトンが受け継がれていくんだなあって感じはあった。twitterのフォローしてる人とかはSHINING LINE*の最後のサビであかりちゃんが出てきてバトンが渡されたって書いてあったけど、思い出せない。たぶん泣いてたから。SHINING LINE*好きな曲だったから。

物販はライブブックとトートバッグとサインライブタオルとあかりちゃんの衣装のキーホルダー。ただ、twitter見てるとハコライブ買っておけばよかったかなとか思った。というか買おうか迷ってる感じ。

twitterにも書いたけど、やっぱりアイドル物でこういうライブやったのアイカツ！が初めてだけど、アイマスが先陣切って欲しかったなと思いました。自分も途中からアイマスにはまってた層だけどあのライブ動画の感動が元だから。twitterでもアイマス好きだった人がかなりアイカツ！やってるのもそういうライブCG動画の進化がアイカツ！のほうが上だったからだようなあと思う。

とりあえず2015シーズンのあかりちゃん見るまでは死ねないので研究とかいろいろと頑張ろうと思いました。BGMはSHINING LINE*で「かなえたい今日の前で　ココロが試される」んだから。

pythonで-infとか入っているデータを読み込んでnanに置き換え、その後nanに平均を代入する方法

pythonで文字と数字が混ざっていたり、空白があるデータ
#value1 value2 value3
1 2 3
4    5
6 -inf error
とかをpythonというかnumpyで読ませるにはgenfromtxtを使う。
空白や文字は自動的にnanになる。たぶん上の場合だと4 5 nanになるはず。「,」区切りで4,,5なら4 nan 5になる。

nanなら後はnumpyにあるnanmeanとかnanstdやscipyのscipy.stats.nanmedianとか使えば平均、標準偏差、中央値が出せるんだけど、-inf入ってるとこれは文字じゃなくてマイナス無限大として認識されるからやっかいって話。

genfromtxt自体は
data = numpy.genfromtxt('data.txt', comments="#", dtype="d")
みたいな感じ。dはdouble型のとき。dtypeは"f8"とかでfloat, "S5"でstring 10文字みたいにいじれる。さらにdtype=[('x','f8'), ('y','f8'), ('z', 'f8')]みたいにするとラベルつけられる。
http://retrofocus28.blogspot.com/2014/03/numpy.html
参照。
あとヘッダーがある場合はskip_header=5とかすると5行分飛ばす。

文字とか欠けているところがなかったらloadtxtでやったほうが楽かもしれん。
http://kaiseki-web.lhd.nifs.ac.jp/documents/Python/textfile.htm
http://pacocat.blogspot.com/2011/12/numpyloadtxt.html
参照。

で、置き換えするためにnumpy.whereで-infが入っているarrayを探す。
index = numpy.where(data, == float('-inf') )
で-infがある場所をindexする。たぶん上のdataの例なら
index=(array([2]), array([1]) )
か。

このindexを使ってfor文で-infをnanに置き換えしていく。
for i in range(len(index[0]))
  data[index[0][i]][index[1][i]] = float('nan')
という感じ。
index[0]がdata[0]に対応してindex[0][i]でindex[0]のarrayから1つずつ要素を取り出してくる。index[1]とdata[1]の関係も同様。で、＝で代入。こんな感じでたぶん他で代入するときにもできるはず。

で、平均を測る。ここで測る平均は列の平均を測りたい。
average = numpy.nanmean(data, axis=0)
axisをつけないと全部の要素の平均、axis=0で列、axis=1で行の平均。というかpythonのarrayでの行と列、どっちがどっちだかいまだによくわかってねえ。
http://docs.scipy.org/doc/numpy-dev/reference/generated/numpy.nanmean.html
を見ること。

でnanに置き換えしたようにfor文使って平均に置き換えしていく。
for i in range(len(index[0]))
  data[index[0][i]][index[1][i]] = average[index[1][i]]

列の平均なのでindex[1]なのがミソ。

これでnumpy.savetxtで書きだす。

numpy.savetxt("test.txt", data)

たぶんこれもヘッダーつけたりとかオプションある。

http://docs.scipy.org/doc/numpy/reference/generated/numpy.savetxt.html

http://airtoxin.wordpress.com/2013/02/02/numpy-savetxt%E3%81%A7csv%E5%87%BA%E5%8A%9B%E3%82%92%E3%81%99%E3%82%8B%E6%96%B9%E6%B3%95/

http://conta.hatenablog.com/entry/2013/12/02/183538

とか参照。

-infが入っててすげー疲れた。python雑魚が成長できたと思いました。