1月 就職しました．まぁぶっちゃけ転職という感じですが．まったく知らない研究分野のまったくよくわからないプロジェクトになんだかよくわからないうちに流れ着いたという感じです．着任する前は出来るだけ早く来いということだったのに着任してみると居室…

一昨日は共同研究者の息子と会いました．どうやら私達のプロジェクトを彼の観点から手伝ってくれるそうです．共同研究者が「すごくできるんだから」と自慢していたのでかなり期待していったのですが，大分色々なことを知らないみたいでした． でも，それはた…

やっぱりいろいろと書き残しておかないともったいないなと思って久しぶりに blog を再開しようかと思います．いつまで続くかしらん．

昼過ぎに来たメールに「（ぼくのなまえ）さんのメールを見ていますと、数十年前○○さんと見た夢が実現出来そうな気がしてきます。」と書いてあったのを見て，なんだかすごくやる気が出て来ました．今の仕事を彼らのためにやってる訳ではないので目的はそれぞ…

日本にいる頃から4-5時くらいに起きる生活だったのですが，未だに早朝起きるクセが抜けません．夜まで起きてると色々と不安になることが多いので夜は早く寝る方です．でも北京の朝はとても寒いのでさすがにつらいです． 今日は朝すごく悲しい夢を見ました．…

昨日は一昨日もらったデータを2日間かけて解析しました．二種の時系列データを適当に解析してみていろいろ眺めてみると結構素性がいいもので安心しました．自然界の話なのにここまで素性がいいと気持ちがいいものです．

http://d.hatena.ne.jp/seaborgium/20111227/1324993156

週末に EOS 6D の実験しようと思ったんですけど，火曜日くらいに満充電にして電源切っておいておいていたのに電池切れしていて使えませんでした．GPS logging モードは電池を猛烈に食うみたいです．勘弁してよね．というわけでまた来週．

EOS 6D で撮られたデータを，手っ取り早く python で処理したいと思います．fits にさえなっていればいつも通りの処理で出来るのでとりあえず fits にします． git clone git://github.com/eaydin/cr2fits.git このソフトはシェルでしかえ使えないようになっ…

去年，EF8-15mm F4L フィッシュアイ USM を買いました．このレンズ，APS-C でも使えるのですが，フルサイズじゃないと全周魚眼にならないので，手持ちの EOS Kiss X3 じゃ我慢できず，結局フルサイズセンサーの載っている EOS 6D を買ってしまいました．高か…

さて，何をするにしてもコンピュータから制御できないと話にならないので，早速 Mac から制御するようにしたいと思います． $ port install gphoto2昔 gphoto2 で 5D MkII を操作する実験をするときに作ったメモを参考にいろいろやりたいと思います． 実験室…

iPhone5 を買いました．iPhone 4 からの進化とそれから，softbank から au への MNP をして au LTE が使えるようになったおかげで外出先ではとても快適なモバイルライフがおくられるようになりました．特にまな板と化していた wifi 版 iPad がテザリングを使…

この歳になるまできちんと理解していなかったので一応まとめておきたいと思います．平均値に並ぶ基本的な統計量として，分散が上げられると思います．この分散の計算は で得られるわけです．ところが，この手の計算を習うときに，よく N で割るのではなく， …

真値 0 を持ち，分散 1 で観測される x という変数を 5 回観測したとする．その5回の観測量から，元の分布関数の特長を調べたいので bootstrap する． import numpy import pylab n = 5 trial = 1000 a=numpy.random.normal(0,1,n) b=numpy.array([a[map(lam…

lambda で演算する方 import numpy import time f = lambda x: x**2 def test0(n): """空""" t1=time.time() t2=time.time() print t2-t1 def test1(n): """内包表現""" t1=time.time() a=[ f(x) for x in range(n) ] t2=time.time() print t2-t1 def test2(…

ちょっとしたテスト import numpy import time def test0(n): """空""" t1=time.time() t2=time.time() print t2-t1 def test1(n): """内包表現""" t1=time.time() a=[ x for x in range(n) ] t2=time.time() print t2-t1 def test2(n): """append""" t1=tim…

最近考えた TSV (tab separated values) 形式で書かれたデータファイルを読んで辞書のリストに登録する． import re tsvfile = "filename" f=open(tsv,"r") s=re.compile("[\t\s]") objects=map( lambda x: { "id" : int(x[0]), # ここら辺はそれぞれのデー…

たぶんどれも一緒だけど http://www.cv.nrao.edu/fits/documents/wcs/wcs.all.ps http://www.aoc.nrao.edu/~egreisen/AnAI.pdf 細かく描いてある．FITS の手引きもこれを参考にしたのではないだろうか．

python を使った天文用 plotter として大変便利な astLib ですが，scalebar の長さが違うという重要なバグがありましたのでメモしておきます． なにが間違っているかと言うと，scalebar を描く時の終点を決めるところで，ra 方向にも角度をそのまま足してし…

いつも忘れてしまうのでメモしておく matplotlib を使って図を作ってるときに，tex 形式のテキストをいれたいときに行う環境変数の設定． from matplotlib import rc rc('text', usetex=True) これで使いたいところで r"$\log\left( h^{-1} M/ M_{\odot}\rig…

これが分かりやすいかは別として，テキストファイルととかから読み込んだデータを配列に格納する． p=re.compile("[\t\n]") data = numpy.array( \ map( lambda x: map(lambda y: float(y), p.split(x)), \ [ \ filter( lambda x: x is not "", x) \ for x i…

Synthetic photometry refers to an algorithm which provides magnitudes and colors computed from a spectral energy distribution and the response functions of a photometric system. http://adsabs.harvard.edu/abs/1986HiA.....7..799B何のことかと…

はじめに iPhone や softbank の無線LANサービス加入者は JR 東海の N700 系で，東京新大阪間の新幹線内無線LANサービスを無料で使うことができます．しかし，新大阪までの間を iPhone 画面でインターネットやってると目が壊れるので，母艦の Mac でやりたい…

Basically, it does not support to use the deimos pipeline on the 64bit Linux.Some modification would enable us to use the pipeline: evilmake in idlutils add -m32 in X_CFLAGS, X_FFLAGS options add -melf_i386 in X_LD_FLAGS mrdfits of idlutil…

印象的な変更点 HDR (High Dynamic Range) photo support Bluetooth AVRCP profile HDR は室内で撮った写真なのに窓の外の風景も一緒に映るという機能． AVRCPはこれまで不可能だった Bluetooth デバイスからのリモートコントロール．これで先送り巻き戻しも…

http://www.ucolick.org/%7Ephillips/deimos_ref/guider_picts.html に書いてある．しかし，2番はそのままでは通らない．dss サーバが multi-extension fits をはく模様． ERROR: FXF: EOF encountered while reading FITS fileとか言います．というわけで，…

データのソートについてまとめてある： http://blog.livedoor.jp/yawamen/archives/51492356.htmlたとえば 3つの tuple をもつリストをソートする： for x, y, nu in sorted(peak, key=lambda x:x[2], reverse=True): print n, x, y, nu 相変わらず lambda …

のはみなさんもご承知でしょうけれども，パソコンも大変らしいです． 半年くらい前に買った MacBook Air がもう壊滅的な遅さになって， さすがにやっぱり前の MacBook Pro のスペックと同等なのでそんなもんかぁ・・・ふぅ，25MHz のマシンに比べたらずいぶ…

6/25くらいに iPhone 4 予約しました． 本当はまだ au の二年周期契約がまだ数ヶ月残ってたし，おサイフケータイもないから iPhone には移れないと思っていました． でも，柏の方の研究会に行くために秋葉原を通る用事があって，しかも，研究会が終わってか…

日記を付けることの重要性（要するにlogなんだけど）を思い知った． たとえば2009-07-29 - CGIBOYがサービスやめるからってそれくらいであきらめたりはしない。とか昔に書いたやり方が探せなくて困っていたけど，ここに来らわかった． ファイル更新日時とか…