会社の人から「何やってんの、Rawtherapee使えばそんなの楽勝っすよ」と言われたので試してみた。rawtherapee hdriでぐぐると、日本語ではその人のページと2chのログくらいしか検索されないんですけど… さっそくインストールして使ってみた。Windows版にはgtk+のライブラリやら何やらも詰め合わされて入っている。Linuxアプリをgtk+ライブラリ詰め合わせにしてWindows版、Windowsアプリをwine詰め合わせにしてLinux版、で両方配布するのが今後は主流になるのかもね。言語ファイルのenglishをjapaneseという名前にコピーして、自分で翻訳していけばUIも日本語が出るようになるので便利だ(？)。私はこういう技術に疎いので妥当な訳語が分からないですが。 rawのデコードにはdcrawを使っているような雰囲気。機能がありすぎてどれを使えばいいのかよくわからないが、「Color Boost」で最高にブーストすると(?)、不自然に色が強調された、それっぽい絵になった。jpegでもColor Boostすれば、それっぽくなるなぁ。 今までいじっていたufraw(こちらはdcrawのフロントエンドか？)でも似たような感じにまではできる。機能はRawTherapeeのほうが多そうです。 それで、RawTherapeeを使って、jpeg1枚とraw1枚でどういう差が出るのか、対決してみようかと思いました。 このシリーズで初めて画像を出してみよう。元の絵はさきほど撮った室内の様子です。自宅の写真はちょっと恥ずかしいですけど。キャプションをそれっぽくするために、無駄にポラロイド風にして傾けてしまった(笑)。RawTherapeeの設定は、以下のようなもの。操作としては普通に読み込んで、ColorBoostでめいっぱい効果が出るようにしてみただけです。

VMWare経由で使っているFedora 8をx86_64版にしてみた。ホストOSが32bitでもゲストOSを64bitにできるようです。えらいです。

VMを生成するときに警告が出てIntel VTが有効になってないとか言われたので、BIOSを眺めて有効にした。今まで有効だと思って「さすがに速いな…」とか呟いていたのだけど、間抜けでした。VMWareはコードのロード時にバイナリを書き換えて実行するときのオーバーヘッドを減らしているので、無効でもそこそこ速いはずです、というのは言い訳です。

確かに、こころなしか速くなってる気がしますね。気のせいかもしれません。

聞いていたとおり、ポインタとlongが64ビットになってますね。あと特に何かせずに普通にoff_tを使っても64ビットになっているので、面倒だったFILE_OFFSET_BITS=64も必要なくなります。あとは変わらず、かな。

# for i in char short int long "long long" quad_t "void *" off_t; do sh sizeof.sh "$i" sys/types.h ; done
sizeof(char)=1
sizeof(short)=2
sizeof(int)=4
sizeof(long)=8
sizeof(long long)=8
sizeof(quad_t)=8
sizeof(void *)=8
sizeof(off_t)=8

スクリプトは以下の通り。

理論はわかりませんが、少しずつためていたサンプル画像を観察すると、だいたいEV補正値が0.5違う場合、JPEGで見えるRGBの値(0～255)で20程度違いがあります。20という値が正確なのかどうかはわかりませんが、色の平均をとってみると、そんな感じで出る画像が多い。EV補正±0.3の場合は10程度の差に見える。 そのため、オートブラケットでEV補正±0.5で3枚撮った場合は、RGBの範囲が0～255だったのが-20～275くらいに広がるという感じにまではできそうです。 明るいほうと暗いほうに、だいたい7～8%ずつ増えて、256階調→約300階調くらいになる程度ですね。ちょっと物足りない感じもしますが、あくまでも「お手軽」だからこのくらいでいっか、という気もするところ。 0～99を暗め明るめの画像から、100～199(元画像の値で80～179)は中くらいの画像から、200～300(160～255)を明るめ暗めの画像から採用すれば、そこそこ色の範囲が広くなるのではないかと思います。 (追記) 2008/02/20(水) 19:44 「明るめ」と「暗め」が逆でしたね。

今回はメタデータについて。Raw(DNG)にもjpegにもEXIF情報が入っていて、撮影したときの情報を得ることができる。実際に出ている色のみで処理するのもいいけど、こういうメタデータをもとに明るさ情報を取り出してどうのこうのする、というのもアリだと思います。 EXIF情報を眺めると、ブラケット撮影で変更しながら撮影された明るさの情報が入っているのは、ExposureBiasValueという部分です。EV値(意味はわかりませんが高いほど明るい画像になるようです)の差分が入っている。これは…分数型かなぁ。 もともとのEV値は直接は入っていなくて、計算しなくてはならないようだ。BrightnessValue+log2(ISOSpeedRatings/3.125)+ExposureBiasValueの値がEV値、という情報をWebのどこか(URL忘れました)で仕入れて計算。 私のカメラでは、DNGでもJPEGでもこのスクリプトで食えています。

asciiの並びでは、「-」は数字やアルファベット、「/」や「.」よりも前、「_」は数字やアルファベット大文字より後ろでアルファベット小文字より前、という特性がある。

!"#$%&’()*+,-./0(略)9:;<=>?@A(略)Z[]^_`a(略)z{|}~

そのため、あるものに追加して属性を書く目的では、「-」よりも「_」が適している。何を言おうとしているのかというと、以下の2つのファイルがあるとすると、

• 01234.jpg
• 01234-1.jpg

名前でソートしたときの並び順が

• 01234-1.jpg
• 01234.jpg

である。「.(ドット)」よりも「-(ハイフン)」が前にあるためだ(拡張子を特別扱いするような面倒なソート方式を使っていれば別)。そのため、「_0」しかないときでも省略することができない。これが「-」ではなく「_」なら、

• 01234.jpg
• 01234_1.jpg

という順番になり、「_0」を省略しても構わないことになる。

今からasciiコードや拡張子の区切り文字を再定義できればなぁ…と思う人ってけっこういるんじゃないかな(「.」はもっと前にあるべき文字だよ…)。でも実現できるもんではないので、我々はあるものに合わせていくしかない。

xyzzyを作業環境にしようと思い立って、auto-save-buffers.lというのを見つけてきて(require ‘auto-save-buffers)してみたら、「名前が衝突するため、exportできません」と言われてしまう。見ると、(in-package “editor”)としているので、editorという名前空間に自前のシンボルをexportしようとしてエラーになっているようにも見えた。

調べてみると、xyzzyのlispの名前空間は勝手に作っていいものではなくて、いくつか決まったものがあってネストされた状態になっている(?)みたいな。そこで、(in-package “user”)に変更したら、うまいこと動いてくれたようだ。xyzzyのlispはEmacs-lispよりもcommon-lispに近いらしいですね。ちゃんと分数(実数)型になるらしい。Lisp系の言語は以前にちょっとかじってた程度なんで分からないけど。

実数型はなんで流行らないのかな。誤差を考えなくて済むので便利だと思うんですけどね。

あとはTortoiseHGを入れてmercurialで作業環境を管理するようにして、Linuxとの同期をとったりできるようにした。auto-save-buffersもこういうものを前提にしたものですしね。Windowsで作業するというのがイメージできていなかったが、そう悪くない環境になるような気がする。xyzzyのPythonのモードも悪くない出来だと思います。

HDRIに憧れる少年がいたとする。まあ実際は30代のオッサンですけど。自分のコンパクトデジカメでできるだけHDRIな人になるにはどうしたらいいのか。いくつか方法を考えてみたいなぁと。 HDRIっていうのはいろんな文脈で使われるようだが、私が言っているのはこっちに近いかな。写真っぽくない、絵画みたいな雰囲気。といっても非現実的になるわけじゃなくて、人間の目で見たものにより近い感じになる。写真だと逆光だのなんだので条件が悪いと猛烈に色が飛んで、目にははっきり見えるものがぜんぜん写らないじゃないですか。芸術を追求したいわけじゃなくて、あれがむかつくんですね。だからHDRIに憧れるといっても、「スクリプトを通せばある程度きれいになります」というレベルでいい。 そんなことができるのか。 このカメラにはオートブラケットの機能がついている。明るさ(用語はわからないが)を変えながら3回シャッターが切られて、３つの画像が保存される。これを組み合わせて、明るい部分を暗い画像から、暗い部分を明るい画像から採用すればそれっぽくなるような気がする。これならPythonのPILで簡単なスクリプトを書くか、Gimp(Script-Fu?)で頑張ればどうにかなるレベル。3枚のうち、暗すぎる部分と明るすぎる部分を捨てながら色を採用orブレンドしていけばいい。ただ、何の単位かわからないが、±0.5までしかないのであまりアグレッシブな(?)絵にはならないかもしれないし、三脚ではなく手で持って液晶を見ながらいい加減に撮るし、ずれると悲惨かもしれないね。ちょっとブレやすい程度で済めばいいけど。 ズレの解消法として、3つの画像をautopano-SIFTとhuginで重ねて、enblendの部分を自分で書くとか？ 自分で書かなくても、enblendのオプションでどうにかなるのかもしれない。ありものをつなぎ合わせるだけでスクリプトにもしやすい(huginの部分はnonaを直接呼べばいいだけ?)。調べてみる価値はあるだろう。 また、このカメラはRaw(AdobeのDNG形式＝tiffの拡張)でも保存できる。Rawでオートブラケットは効かないようだが、とにかくjpegでは飛んでしまう色まで記録できるように作られているはずだから、HDRIに近づけるだろう。OpenEXRに変換すれば、exrtoolsあたりでがんばれば、それなりになるかもしれないし、Pythonのモジュールもあるからどうにかなる可能性が高まる。 ざっと考えるとこの2つ。他にもやりようがあるのかもしれないが。 どっちがいいのかな。最初に気になるのは、Rawのデータのビット数。これが8ビットだと、何やってるんだという話になる。さすがにないか。調べる方法が分からないが手探りでやってみると、

カメラのスペックが「焦点距離 5.9mm（35mm換算値28mm）」だから、huginのcrop factorは28/5.9≒4.75のようだ。毎回計算するのがばからしい。

手すりの部分が大きな失敗

雪の降る中、サーバを移行しました。今までバックアップと兼任だったマシンをサーバに降格させたんですが、途中でバックアップがない状態でファイルを間違って消したり(アップした画像などは復活できませんでした)、いろいろ騒ぎながらもメールとWebサーバの移行ができました。今度のサーバはThinkPad X30(2002年発売の10周年記念の限定モデル)です。

メール関係ではbogofilterをやめてSpamAssassinに。ThunderbirdのSPAM除けがSpamAssassinのヘッダを簡単に(チェックボックスにチェックするだけで)認識してくれるのでそうした。ThunderbirdのSPAM除けと2重チェックになってくれて、良いかもしれない。

OS関係ではCentOS5.1にしておいた。仮想マシン上でサーバを動かして、次回のハードウェアリプレースのときにはマイグレーションでやろうとも思ったけど、Intel VTに対応したマシンでもないし、子供に邪魔されたり奥さんにコキ使われながら移行する作業というのも悪くないので、やめておいた。

私も仕事の関係でExcelなどを使うことがあるが、Excelの日付時刻情報の使いにくさは何とかならないものか。まあ、他の表計算のソフトの状況は知らないのだけど。 何を言いたいのかというと、Excelの日付時刻情報は浮動小数点数で表現されている。1.0が1日。これが曲者だ。1秒はどの値かというと、1/(24*60*60)=約0.00001157407407407407であり、10進の浮動小数点数では正確に表現することはできない。人類は1日の中の時間を表現する時、24進数と60進数を併用している。秒以下の部分(ms, us, ns…)は10進数だが、24進数と60進数を10進数(実際は2進数)の小数で表現するのは無理があると思う。 Excelで日付時刻情報を縦もしくは横軸にとってグラフ化するとき、10秒単位で目盛をつけたいと思ったとする。軸のプロパティから目盛単位を変更するのだけど、どうしたらいいのだろう？ 確か「00:00:10」が正解なのだが(ちょっと手元にないので確認できない)、なんだかなぁという気がする。目盛単位を修正するとき、得体の知れない小数に直されているので、ちょっと10秒を3秒にしたいんだ、というときに10を3にするような手軽さはなく、結局分かりにくいので、何だっけ、SECとかMINUTEとかの関数を組み合わせて秒単位の整数を作ってグラフ化するのが楽ということに(私の中では)なった。 コンピュータで日付時刻を表現するとき、少なくとも私の分野ではコンピュータにおける紀元(EPOCH)からの秒数で表現するという手法が一般的だ。EPOCHは1970年1月1日の00:00:00 GMT。秒以下の精度が必要な場合は、Cのstruct timevalのように秒とusの2つの整数を合わせて使うか、Pythonのように浮動小数点数にするか(usの精度を確保するには仮数部が31+20=51ビット以上あればいいので普通の倍精度浮動小数点数で大丈夫)。これがもっとも単純だ。60進数、24進数といえども整数ならば単純で、整数を整数(60や24)で割ったり剰余を取れば簡単に処理できる。 まあよくよく考えると、精度としては1.0が1日を示すか1秒を示すかに優劣はない。でもなんか気持ち悪いし、実際に処理しにくかったことがあったので書いてみました。