2008年8月の日記

2008年8月の日記	2008-08-30(Sat) 12:58

打ち上げ日	射場	ロケット	軌道	名前	開発	運用	重量など	参照
8月2日 17時30分成功	日本鹿児島県内之浦宇宙空間観測所	観測ロケット S-520-24号機			JAXA 帝京科学大学ロケット飛行中の微小重力環境を利用した結晶成長実験	JAXA ISAS
8月3日 12時34分失敗	マーシャル諸島共和国クェゼリン環礁レーガン試験場	SpaceX Falcon 1-3		Jumpstart PreSat NanoSail-D			Spaceflight Now
8月15日 5時44分成功	仏領ギアナクールーギアナ宇宙センター ELA-3	アリアンスペース Ariane 5 ECA	静止東経144度	通信衛星スーパーバード7号機 Superbird 7 Superbird C2 33274 2008-038A	日本三菱電機	日本三菱電機宇宙通信株式会社	5,000kg以下 DS-2000 28×Ku 15年	宇宙通信三菱電機三菱電機 Arianespace Gunter
静止西経125度	通信衛星 AMC-21 33277 2008-038D	Alcatel Alenia Space (prime) Orbital Sciences Corporation (OSC) (bus)	アメリカ SES Americom	2,500kg STar-2 24×Ku 16年	Arianespace Gunter
8月20日 7時43分成功	カザフスタンチュラタムバイコヌール宇宙基地	ILS Proton M/Breeze M	静止	通信衛星 Inmarsat-4 F3	EADS Astrium	Inmarsat	Eurostar 3000 6,000kg	ILS ILS Spaceflight Now
8月29日 16時15分成功	カザフスタンチュラタムバイコヌール宇宙基地	Dnepr	600km 98度	地球観測衛星 RapidEye×5	イギリス Surrey Satellite Technology Ltd.	ドイツ	150kg×5 7年以上	Spaceflight Now
← 7月の打ち上げ　／　9月の打ち上げ →

火星探査機フェニックスがロボットアームを使って火星の地表の深さ約5cmのところから採取した土壌をフェニックスが搭載しているTEGA(Thermal and Evolved-Gas Analyzer;熱発生気体分析装置)で分析したところ水蒸気が検出され、NASAは正式に火星に水が存在することを確認したと発表しました。6月1日に書いたように火星探査機フェニックスは火星に到着して間もなく自分の足元に氷と思しき滑らかな物体があるのを発見し、その物質が一部蒸発したことから水が含まれていることはほぼ確実とみられていました。火星周回軌道からの間接的な観測でも極地域に大量の氷があると予測されていましたが、地球以外の天体で採取した物質を直接分析装置にかけることで水の存在が確認されたのは今回が初めてです。5cm掘っただけで氷が出てくるとなると近い将来から遠い将来まで様々な火星探査・開発計画に弾みがつくのではないでしょうか。火星の水をめぐる利権争いなどが起こりませんように。

2008年8月2日（土）

日本時間の昨夜、カナダからグリーンランド、北極海、ロシア、モンゴル、中国にかけて皆既日食がありました。美しいダイヤモンドリングが観測されています。

5月26日に書いたように日本国内で見られる皆既日食が2009年7月22日にあります。

17時30分、JAXAは帝京科学大学と共同でロケット飛行中の微小重力環境を利用した結晶成長実験を行うための観測ロケットS-520-24号機を鹿児島県の内之浦宇宙空間観測所から打ち上げました。ロケットは最高高度293kmに達し、実験を終えて予定された海域に落下しました。

2008年8月4日（月）

日本時間3日12時34分、3度目の試みとなるFalcon 1ロケットの打ち上げがマーシャル諸島共和国クェゼリン環礁のレーガン試験場で行われましたが、ステージを分離することができず、打ち上げは失敗しました。

Googleの天気検索。日本語モードで都道府県名、市区町村名、駅名などに「天気」(またはweather)というキーワードを添えて検索すると日本国内のローカルな天気予報が表示される。例えば天気東京都。ちなみに英語モードを使うと海外の天気予報が表示される。例えばweather Beijing。日本語モードと英語モードでは天気予報のマークが違うので比べてみると面白い。晴れのマークの太陽の色は日本語モードでは赤くて英語モードでは黄色いのだ。マークなのだからわかりやすければよいのだけれど、日本人が太陽を赤く塗るのを変だと思っている外国人は少なくないかも知れない。

2008年8月5日（火）

Google マップとGoogle Earthのストリートビューが日本とオーストラリアに対応。主要都市の周辺だけとはいえ、対応地域を示す青い線(Google Earthではカメラのマーク)が予想以上に郊外まで広がっていて自宅のある団地の中まで辿れたのでちょっと焦った。道の名前は「都道○号線」ではなくて「○○街道」で表示して欲しいかも。

2008年8月7日（木）

2008年8月8日（金）

12時57分頃、東京都と神奈川県で震度4の地震。震源は東京都多摩東部、深さ40km、M4.5。津波の心配はなし。

2008年8月9日（土）

Tentative Name.さんから。3月17日に書いた欧州合同原子核研究機構（CERN）の大型ハドロン衝突型加速器（LHC）が初めてビームを循環させる日が9月10日に決定したそうです。2008年の目標は1ビームあたり5TeV(5兆電子ボルト)で当日の入射エネルギーは450GeV(4,500億電子ボルト)。人類史上最強の加速器ははたして何を見せてくれるのでしょうか。

2008年8月12日（火）

北島も谷本もスゴイ。口パクはよくあることだと思うけれど花火がCGだったのはがっかり。

2008年8月14日（木）

X68トラブル110番より。推測ですが、JPEGED.R V1.26で作成したJPEGファイルを表示できるブラウザと表示できないブラウザがあるのは、直流成分のハフマン符号の最長符号のひとつとしてすべてのビットが1の符号を許容するものと許容しないものがあるからだと思います。JPEGED.Rの元になったJPEG.X v1.00はそのような符号は出力しないようですが、それとは別に符号長の問題があってJPEG.Xで作成したJPEGファイルについても表示できるブラウザと表示できないブラウザがあるようです。また、JPEG.Xが横長のアスペクト比としてX密度2:Y密度3を出力していたところがJPEGED.Rでは縦長のX密度3:Y密度2に変更されており、アスペクト比を再現するSusie32などでJPEGED.Rで作成したJPEGファイルを表示すると異常に縦に伸びた画像が表示されてしまいます。

2008年8月15日（金）

日本時間5時44分、三菱電機と宇宙通信の通信衛星スーパーバード7号機(Superbird-7)とSES Americomの通信衛星AMC-21が仏領ギアナクールーのギアナ宇宙センター ELA-3からアリアンスペースのAriane 5 ECAロケットで打ち上げられました。2機の衛星は近地点高度250km、遠地点高度35,931kmの静止トランスファ軌道でロケットから切り離され、打ち上げは成功しました。スーパーバード7号機は技術試験衛星VIII型「きく8号」(KIKU-8 (ETS-VIII))や運輸多目的衛星「ひまわり7号」(MTSAT-2)で実証されている三菱電機のDS2000衛星バスにKuバンドの中継器28本(27MHz×20、36MHz×8)を搭載した日本初の国産商業衛星です。打ち上げ重量は4,820kg。スーパーバードC号機の後継機スーパーバードC2号機として東経144度の静止軌道に配備されます。AMC-21は打ち上げ重量が約2,500kgで静止位置は西経125度。

2008年8月16日（土）

14日の続き。X680x0で作ったJPEGファイルがIEやFirefoxなどで表示されないことがあるのはJPEGファイルの構造が微妙に間違っているからです。下のPerlスクリプトはそのようなJPEGファイルを修正して表示できるようにします。JPEGファイルは普通は非可逆圧縮なので画像を一旦ビットマップに展開して元と同じサイズに再圧縮すると画質が劣化してしまうことがありますが、このスクリプトを使った場合は画質は変化しません。アスペクト比の情報も修正していますがウェブブラウザ側が対応していないので普通は正方形で表示されます。電脳倶楽部の画廊で表示できなかったJPEGファイルを10個くらい試しただけなのでバグが残っているかも知れません。

correctjpeg.pl

#Correct JPEG files created by X680x0 JPEG.X, JPEGED.R
#  2008-08-16 Makoto Kamada
#    Pixel density should not be 3:2(portrait) but 2:3(landscape) (JPEGED.R v1.26)
#    Huffman code whose all bits are one is not allowed (JPEGED.R v1.26)
#    Too long symbols are not allowed (JPEG.X v1.00)
#    DHT segment should have proper size (JPEGED.R v1.00+0.16)
#    JPEG file should end at EOI
use strict;
$| = 1;
@ARGV or die "Correct JPEG files created by X680x0 JPEG.X, JPEGED.R\n" .
"usage: perl $0 JPEG-file-name ...\n";
my $ncorrected = 0;
for my $name (@ARGV) {
  my $namebak = "$name.bak";  #backup file name
  my $nametmp = "$name.tmp.$$";  #temporary file name
  my $in = 0;  #is input file opened?
  my $out = 0;  #is output file opened?
  my $corrected = 0;  #is file corrected?
  my $eoi = 0;  #is EOI appeared?
  my $message = '';  #error message
thisfile:
  {
    -d $name and $message = 'directory', last thisfile;
    -f $name or $message = 'not found', last thisfile;
    open IN, $name or $message = 'cannot open', last thisfile;
    $in = 1;
    open OUT, "> $nametmp" or $message = "cannot create temporary file $nametmp", last thisfile;
    $out = 1;
    my $nmcu;  #number of mcus
    my $ncomponent;  #number of components
    my @nblock;  #number of blocks
    my @ht = ('', '', '', '');  #huffman table DY,DC,AY,AC
    my $b;  #input buffer
    #segment loop
    while (2 == read IN, $b, 2) {
      print OUT $b;
      0xff == ord $b or $message = sprintf ('missing segment marker, 0x%02x appeared', ord $b), last thisfile;
      my $m = vec $b, 1, 8;  #segment marker
      $m == 0xd8 and next;  #SOI
      if ($m == 0xd9) {  #EOI
        $eoi = 1;
        #JPEG file should end at EOI
        read IN, $b, 1 and $corrected = 1;
        last;
      }
      2 == read IN, $b, 2 or $message = 'unexpected EOF', last thisfile;
      my $l = vec $b, 0, 16;  #length of segment, inclusive
      $l -= 2;  #length of segment, exclusive
      $l >= 0 or $message = 'unknown format', last thisfile;
      $l == read IN, $b, $l or $message = 'unexpected EOF', last thisfile;  #segment body
      if ($m == 0xe0) {  #APP0
        "JFIF\x00" eq substr $b, 0, 5 or $message = 'not JFIF', last thisfile;
        #Pixel density should not be 3:2(portrait) but 2:3(landscape)
        my ($xdensity, $ydensity) = unpack 'n2', substr $b, 8, 4;  #pixel density
        if ($xdensity == 3 && $ydensity == 2) {
          vec ($b, 9, 8) = 2;
          vec ($b, 11, 8) = 3;
          $corrected = 1;
        }
      } elsif ($m == 0xc0) {  #SOF0 (baseline)
        8 == ord $b or $message = 'unknown format', last thisfile;
        my ($height, $width) = unpack 'n2', substr $b, 1, 4;  #image size (height,width)
        $ncomponent = vec $b, 5, 8;  #number of components (1=Y,3=YCbCr)
        my $maxhs = 0;
        my $maxvs = 0;
        @nblock = ();
        for (my $k = 0; $k < $ncomponent; $k++) {
          my $vs = vec $b, 7 + 3 * $k, 8;
          my $hs = $vs >> 4;
          $vs &= 0x0f;
          $maxhs >= $hs or $maxhs = $hs;
          $maxvs >= $vs or $maxvs = $vs;
          push @nblock, $hs * $vs;
        }
        my $mcuwidth = 8 * $maxhs;
        my $mcuheight = 8 * $maxvs;
        my $hmcu = int (($width + $mcuwidth - 1) / $mcuwidth);  #number of mcu, horizontal
        my $vmcu = int (($height + $mcuheight - 1) / $mcuheight);  #number of mcu, vertical
        $nmcu = $hmcu * $vmcu;  #number of mcu
      } elsif ($m == 0xc2) {  #SOF2 (progressive)
        $message = 'progressive mode, not supported';
        last thisfile;
      } elsif ($m == 0xc4) {  #DHT
        my $j = 0;
        while ($j + 17 < $l) {
          my $n = vec $b, $j++, 8;
          $n = ($n & 0x10) >> 3 | ($n & 0x01);  #0:DY,1:DC,2:AY,3:AC
          my @s = ();  #length of huffman code
          for (my $k = 0; $k < 16; $k++) {
            push @s, ($k + 1) x vec $b, $j + $k, 8;
          }
          $j += 16;
          my @c = ();  #huffman code
          for (my ($c, $s) = (0, 1); @c < @s; $c <<= 1, $s++) {
            while ($s[@c] == $s) {
              $c < 1 << $s or $message = 'huffman code collides', last thisfile;
              push @c, $c++;
            }
          }
          for (my $k = 0; $k < @s; $k++) {
            my $s = $s[$k];  #length of huffman code
            my $c = $c[$k];  #huffman code
            my $w = vec $b, $j + $k, 8;
            my $z = $w >> 4;  #zero run-length
            $w &= 0x0f;  #length of symbol
            my $e = { size => $s, code => $c, width => $w, zero => $z, zrl => 0, eob => 0,
              newsize => $s, newcode => $c, error => 0 };
            #Huffman code whose all bits are one is not allowed
            if ($c + 1 == 1 << $s) {
              $e->{newsize} = $s + 1;
              $e->{newcode} = $c << 1;
              vec ($b, $j - 17 + $s, 8)--;
              vec ($b, $j - 16 + $s, 8)++;
              $corrected = 1;
            }
            if ($n < 2) {  #D
              #Too long symbols are not allowed
              if ($z) {
                $e->{width} = $w = 0;
                $e->{zero} = $z = 0;
                $e->{error} = 1;
                vec ($b, $j + $k, 8) = 0;
                $corrected = 1;
              }
            } else {  #A
              if ($w == 0) {
                if ($z == 15) { #ZRL
                  $e->{zrl} = 1;
                } elsif ($z == 0) { #EOB
                  $e->{eob} = 1;
                } else {
                  $message = sprintf 'unknown symbol code 0x%x%x', $z, $w;
                  last thisfile;
                }
              }
            }
            $ht[$n] or $ht[$n] = [('') x 256];
            my $ht = $ht[$n];
            if ($s > 8) {
              $s -= 8;
              my $ch = $c >> $s;
              $c &= (1 << $s) - 1;
              'HASH' eq ref $ht->[$ch] and $message = 'huffman code collides', last thisfile;
              'ARRAY' eq ref $ht->[$ch] or $ht->[$ch] = [('') x 256];
              $ht = $ht->[$ch];
            }
            $c <<= 8 - $s;
            my $c2 = $c + (256 >> $s);
            while ($c < $c2) {
              'HASH' eq ref $ht->[$c] and $message = 'huffman code collides', last thisfile;
              'ARRAY' eq ref $ht->[$c] and $message = 'huffman code collides', last thisfile;
              $ht->[$c++] = $e;
            }
          }
          $j += @s;
        }
        $j <= $l or $message = 'broken huffman table', last thisfile;
        #DHT segment should have proper size
        if ($j < $l) {
          $b = substr $b, 0, $j;
          $corrected = 1;
        }
      } elsif ($m == 0xda) {  #SOS
        print OUT pack 'n', 2 + length $b;  #segment length
        print OUT $b;  #segment body
        my $q = 0;  #current bits (input)
        my $f = 0;  #number of filled bits (input)
        my $p = 0;  #previous byte (input)
        my $r = 0;  #current bits (output)
        my $g = 0;  #number of filled bits (output)
        for (my $mcu = 0; $mcu < $nmcu; $mcu++) {
          for (my $component = 0; $component < $ncomponent; $component++) {  #0:Y,1:Cb,2:Cr
            my $nblock = $nblock[$component];
            for (my $block = 0; $block < $nblock; $block++) {
              for (my $scan = 0; $scan < 64; $scan++) {
                my $n = ($scan == 0 ? 0 : 2) | ($component == 0 ? 0 : 1);  #0:DY,1:DC,2:AY,3:AC
                #input huffman code
                if ($f < 8) {
                  if ($p == 0xff) {
                    read IN, $b, 1 or $message = 'unexpected EOF', last thisfile;
                    0x00 == ord $b or $message = 'unknown format', last thisfile;
                  }
                  read IN, $b, 1 or $message = 'unexpected EOF', last thisfile;
                  $q |= ($p = ord $b) << (32 - ($f += 8));
                }
                my $e = $ht[$n]->[($q >> 24) & 0xff];
                if ('ARRAY' eq ref $e) {
                  if ($f < 16) {
                    if ($p == 0xff) {
                      read IN, $b, 1 or $message = 'unexpected EOF', last thisfile;
                      0x00 == ord $b or $message = 'unknown format', last thisfile;
                    }
                    read IN, $b, 1 or $message = 'unexpected EOF', last thisfile;
                    $q |= ($p = ord $b) << (32 - ($f += 8));
                  }
                  $e = $e->[($q >> 16) & 0xff];
                }
                $e or $message = 'unknown format', last thisfile;
                $e->{error} and $message = 'too long symbol used', last thisfile;
                $q = ($q << $e->{size}) & 0xffffffff;
                $f -= $e->{size};
                #output huffman code
                $r |= $e->{newcode} << (32 - ($g += $e->{newsize}));
                while ($g >= 8) {
                  my $t = ($r >> 24) & 0xff;
                  print OUT chr $t;
                  $t == 0xff and print OUT chr 0x00;
                  $r = ($r << 8) & 0xffffffff;
                  $g -= 8;
                }
                $e->{eob} and last;
                if ($e->{zrl}) {
                  $scan += 15;
                } else {
                  #input symbol
                  while ($f < $e->{width}) {
                    if ($p == 0xff) {
                      read IN, $b, 1 or $message = 'unexpected EOF', last thisfile;
                      0x00 == ord $b or $message = 'unknown format', last thisfile;
                    }
                    read IN, $b, 1 or $message = 'unexpected EOF', last thisfile;
                    $q |= ($p = ord $b) << (32 - ($f += 8));
                  }
                  my $d = ($q >> (32 - $e->{width})) & ((1 << $e->{width}) - 1);  #symbol
                  $q = ($q << $e->{width}) & 0xffffffff;
                  $f -= $e->{width};
                  #output symbol
                  $r |= $d << (32 - ($g += $e->{width}));
                  while ($g >= 8) {
                    my $t = ($r >> 24) & 0xff;
                    print OUT chr $t;
                    $t == 0xff and print OUT chr 0x00;
                    $r = ($r << 8) & 0xffffffff;
                    $g -= 8;
                  }
                  $scan += $e->{zero};
                }
              }
            }
          }
        }
        $f >= 8 and seek IN, -($f >> 3), 1;
        if ($g > 0) {
          my $t = ($r >> 24) & (0x100 - (0x100 >> $g));
          print OUT chr $t;
          $t == 0xff and print OUT chr 0x00;
        }
        next;
      }
      print OUT pack 'n', 2 + length $b;  #segment length
      print OUT $b;  #segment body
    }
  }
  if ($out) {
    #JPEG file should end at EOI
    if (!$message && !$eoi) {
      print OUT "\xff\xd9";  #EOI
      $corrected = 1;
    }
    close OUT;
  }
  $in and close IN;
  $message or ($corrected or $message = 'not corrected');
  if ($message) {
    $out and (unlink $nametmp or die "$0: $name: cannot delete $nametmp\n");
    warn "$0: $name: $message\n";
  } else {
    $ncorrected++;
    rename $name, $namebak or die "$0: $name: cannot rename $name to $namebak\n";
    rename $nametmp, $name or die "$0: $name: cannot rename $nametmp to $name\n";
    warn "$0: $name: corrected\n";
  }
}
@ARGV > 1 and warn "$ncorrected file(s) corrected\n";

2008年8月17日（日）

JPEGED.Rで作られたJPEGファイルで、交流成分のハフマン符号表に1が16個並んでいる符号が含まれているものを見つけました(132号のF14B_JR.JPGや135号のBUS_STOP.JPGなど)。そのようなファイルに昨日のスクリプトを適用するとファイルが壊れてしまうので、とりあえずエラーにして回避するようにしました。後で修正を試みます。変換されて新しいファイルができたのに表示できるようにならなかったときは、変換後のファイルを捨てて.bakファイルをリネームして元に戻すようにしてください。

2008年8月18日（月）

現地時間16日、イランが国産初の衛星打ち上げロケットSafir Omid(Safir 1, Messenger of Hope, Envoy of Hope, Ambassador of Peace)の打ち上げに成功したそうです。ロケットは長さ22m、直径1.25m、重さ26トンで2段式とのこと。なお、一部報道で人工衛星の打ち上げに成功したと報じられましたが、今回は人工衛星(として機能するもの)は積んでいなかったようです。

X68000 LIBRARYにjpegtran.xを収録しました。JPEGファイルをロスレスで変換するツールで、主にプログレッシブ形式のJPEGファイルをベースライン形式のJPEGファイルに変換したいときに使うことが想定されています。また、jpeged.rが出力した微妙に間違っていてウェブブラウザで表示できないJPEGファイルをjpegtran.xで空変換すると表示できるようになる場合があります(X68トラブル110番のbigslopeさんの書き込みより)。

2008年8月19日（火）

日本時間7時43分、Inmarsatの通信衛星Inmarsat-4 F3がカザフスタンチュラタムのバイコヌール宇宙基地からILSのProton M/Breeze Mロケットで打ち上げられました。9時間3分後の16時46分に衛星が切り離され、打ち上げは成功しました。

2008年8月22日（金）

15日に打ち上げられた宇宙通信の通信衛星スーパーバード7号機がもうすぐ静止化します。静止位置は東経144度の赤道上空。翼を広げた大きさは31.6mあります。

2008年8月24日（日）

定番のFTPクライアントFFFTPのVer 1.96cが公開されました。5月に指摘されていたディレクトリトラバーサルの脆弱性を修正したそうです。

2008年8月25日（月）

ついに1,000万桁を超える素数が発見されたのかも知れないし、間違いかも知れない。

2006年9月4日に知られている44番目のメルセンヌ素数M32,582,657=2^32,582,657-1(980万8,358桁)が発見されてから2年近く経っているので今度こそ大台に乗ることを期待したいところですが、必ずしも小さい順に見つかるとは限らないので知られている44番目よりも小さい可能性もあります。

2008年8月27日（水）

GIMPSのトップページに情報が出ました。素数候補の発見は23日でしたがGeorge Woltman氏(GIMPSの創設者)が休暇中だったため検証は26日に始まりました。検証には異なるプログラムが用いられており、2週間くらいかかるそうです。

2008年8月28日（木）

Jason Papadopoulosさんによる素因数分解プログラムMsieveのバージョン1.37が公開されました。

現在2つの検証プログラムが動いており、終わるのは9月12日頃と9月16日頃とのこと。

2008年8月29日（金）

国際宇宙ステーションISSは2006年6月25日に打ち上げられたロシアのスパイ衛星Cosmos 2421(Kosmos 2421)(の打ち上げロケット?)から発生したデブリ(33246 2006-026RU)との衝突を避けるために日本時間28日1時11分から5分02秒間ATV Jules Verneのスラスタを逆噴射して高度を1.77km下げました。ATVはISSの進行方向から見て最後尾にドッキングしているのであらかじめISS全体の向きを180度反転してから逆噴射しました。現在のISSの構成でこんな芸当ができるモジュールはATVだけです。ISSは放っておくと希薄大気との摩擦抵抗で次第に高度が下がってしまうのでスペースシャトルやATVのスラスタを用いて定期的に巨体を押し上げていますが、デブリを回避するために軌道を変更したのは2003年5月30日以来5年ぶりで、さらにISSにブレーキをかけて高度を下げた例は8年ぶりとのことです。Cosmos 2421から発生した別のデブリ(33248 2006-026RW)も明日10時09分頃ISSに接近すると予測されており、双方の軌道が確定してデブリ回避マヌーバ(DAM;Debris Avoidance Maneuver)が必要と判断されれば明日8時頃に再びATVのスラスタを噴射することになるそうです。

6月10日に「もうすぐ完成」と書いた上海環球金融中心(高さ492m)が完成し、明日30日に地面からの高さが世界一高い展望台(高さ472m)が一般にお披露目されるそうです。

2008年8月30日（土）

日本時間29日16時15分、ドイツの地球観測衛星RapidEye×5機がカザフスタンチュラタムのバイコヌール宇宙基地からDneprロケットで打ち上げられました。5機の衛星は高度600km、傾斜角98度の軌道に投入され、打ち上げは成功しました。150kgの衛星を多数用いることで地球上のどこでも1日以内に撮影することができ、農作物を状況を迅速に把握します。