tb 0.9 → 1.0

data/lib/tb/json.rb

@@ -1,4 +1,4 @@
-# Copyright (C) 2012 Tanaka Akira  <akr@fsij.org>
+# Copyright (C) 2012-2014 Tanaka Akira  <akr@fsij.org>
 #
 # Redistribution and use in source and binary forms, with or without
 # modification, are permitted provided that the following conditions
@@ -28,23 +28,25 @@
 
 require 'json'
 
-class Tb
-  class JSONReader
-    include Tb::Enumerable
-
-    def initialize(string)
-      @ary = JSON.parse(string)
+module Tb
+  class JSONReader < Tb::HashReader
+    def initialize(io)
+      ary = JSON.parse(io.read)
+      super lambda { ary.shift }
     end
+  end
 
-    def header_and_each(header_proc)
-      header_proc.call(nil) if header_proc
-      @ary.each {|obj|
-        yield obj
+  class JSONWriter < Tb::HashWriter
+    def initialize(io)
+      io << "[\n"
+      sep = ""
+      super lambda {|hash|
+        io << sep << JSON.pretty_generate(hash)
+        sep = ",\n"
+      },
+      lambda {
+        io << "\n]\n"
       }
     end
-
-    def each(&block)
-      header_and_each(nil, &block)
-    end
   end
 end

data/lib/tb/ltsv.rb

@@ -1,6 +1,6 @@
 # lib/tb/ltsv.rb - LTSV related fetures for table library
 #
-# Copyright (C) 2013 Tanaka Akira  <akr@fsij.org>
+# Copyright (C) 2013-2014 Tanaka Akira  <akr@fsij.org>
 #
 # Redistribution and use in source and binary forms, with or without
 # modification, are permitted provided that the following conditions
@@ -30,38 +30,7 @@
 
 require 'stringio'
 
-class Tb
-  def Tb.load_ltsv(filename, &block)
-    Tb.parse_ltsv(File.read(filename), &block)
-  end
-
-  def Tb.parse_ltsv(ltsv)
-    assoc_list = []
-    ltsv_stream_input(ltsv) {|assoc|
-      assoc_list << assoc
-    }
-    fields_hash = {}
-    assoc_list.each {|assoc|
-      assoc.each {|key, val|
-        fields_hash[key] ||= fields_hash.size
-      }
-    }
-    header = fields_hash.keys
-    t = Tb.new(header)
-    assoc_list.each {|assoc|
-      t.insert Hash[assoc]
-    }
-    t
-  end
-
-  def Tb.ltsv_stream_input(ltsv)
-    ltsvreader = LTSVReader.new(ltsv)
-    while assoc = ltsvreader.shift
-      yield assoc
-    end
-    nil
-  end
-
+module Tb
   def Tb.ltsv_escape_key(str)
     if /[\0-\x1f":\\\x7f]/ =~ str
       '"' +
@@ -138,70 +107,19 @@ class Tb
     end
   end
 
-  class LTSVReader
-    include Tb::Enumerable
-
-    def initialize(input)
-      if input.respond_to? :to_str
-        @input = StringIO.new(input)
-      else
-        @input = input
-      end
-    end
-
-    def shift
-      line = @input.gets
-      return nil if !line
-      line = line.chomp("\n")
-      line = line.chomp("\r")
-      ary = line.split(/\t/, -1)
-      assoc = ary.map {|str|
-        /:/ =~ str
-        key = $`
-        val = $'
-        key = Tb.ltsv_unescape_string(key)
-        val = Tb.ltsv_unescape_string(val)
-        [key, val]
-      }
-      assoc
-    end
-
-    def header_and_each(header_proc)
-      header_proc.call(nil) if header_proc
-      while assoc = self.shift
-        yield Hash[assoc]
-      end
-    end
-
-    def each(&block)
-      header_and_each(nil, &block)
-    end
-  end
-
-  def Tb.ltsv_stream_output(out)
-    gen = Object.new
-    gen.instance_variable_set(:@out, out)
-    def gen.<<(assoc)
-      @out << Tb.ltsv_assoc_join(assoc) << "\n"
-    end
-    yield gen
-  end
-
-  # :call-seq:
-  #   generate_ltsv(out='') {|recordids| modified_recordids }
-  #   generate_ltsv(out='')
-  #
-  def generate_ltsv(out='', fields=nil, &block)
-    recordids = list_recordids
-    if block_given?
-      recordids = yield(recordids)
-    end
-    Tb.ltsv_stream_output(out) {|gen|
-      recordids.each {|recordid|
-        gen << get_record(recordid)
-      }
+  def Tb.ltsv_split_line(line)
+    line = line.chomp("\n")
+    line = line.chomp("\r")
+    ary = line.split(/\t/, -1)
+    assoc = ary.map {|str|
+      /:/ =~ str
+      key = $`
+      val = $'
+      key = Tb.ltsv_unescape_string(key)
+      val = Tb.ltsv_unescape_string(val)
+      [key, val]
     }
-    out
+    assoc
   end
 
   def Tb.ltsv_assoc_join(assoc)
@@ -209,4 +127,25 @@ class Tb
       Tb.ltsv_escape_key(key) + ':' + Tb.ltsv_escape_value(val)
     }.join("\t")
   end
+
+  class LTSVReader < Tb::HashReader
+    def initialize(io)
+      super lambda {
+        line = io.gets
+        if line
+          Hash[Tb.ltsv_split_line(line)]
+        else
+          nil
+        end
+      }
+    end
+  end
+
+  class LTSVWriter < Tb::HashWriter
+    def initialize(io)
+      super lambda {|hash|
+        io << (Tb.ltsv_assoc_join(hash) + "\n")
+      }
+    end
+  end
 end

data/lib/tb/ndjson.rb

@@ -0,0 +1,63 @@
+# Copyright (C) 2014 Tanaka Akira  <akr@fsij.org>
+#
+# Redistribution and use in source and binary forms, with or without
+# modification, are permitted provided that the following conditions
+# are met:
+#
+#  1. Redistributions of source code must retain the above copyright
+#     notice, this list of conditions and the following disclaimer.
+#  2. Redistributions in binary form must reproduce the above
+#     copyright notice, this list of conditions and the following
+#     disclaimer in the documentation and/or other materials provided
+#     with the distribution.
+#  3. The name of the author may not be used to endorse or promote
+#     products derived from this software without specific prior
+#     written permission.
+#
+# THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS
+# OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
+# WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+# ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY
+# DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+# DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE
+# GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+# INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
+# WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE
+# OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE,
+# EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+require 'json'
+
+# NDJSON is Newline Delimited JSON
+# http://ndjson.org/
+#
+# Tb::NDJSONReader accepts empty lines.
+
+module Tb
+  class NDJSONReader < Tb::HashReader
+    # io.gets should returns a string.
+    def initialize(io)
+      super lambda {
+        while true
+          line = io.gets
+          if line.nil? || /\S/ =~ line
+            break
+          end
+        end
+        if line
+          JSON.parse(line)
+        else
+          nil
+        end
+      }
+    end
+  end
+
+  class NDJSONWriter < Tb::HashWriter
+    def initialize(io)
+      super lambda {|hash|
+        io << (JSON.generate(hash) + "\n")
+      }
+    end
+  end
+end

data/lib/tb/numericreader.rb

@@ -0,0 +1,66 @@
+# lib/tb/numericreaderm.rb - reader mixin for table without header
+#
+# Copyright (C) 2014 Tanaka Akira  <akr@fsij.org>
+#
+# Redistribution and use in source and binary forms, with or without
+# modification, are permitted provided that the following conditions
+# are met:
+#
+#  1. Redistributions of source code must retain the above copyright
+#     notice, this list of conditions and the following disclaimer.
+#  2. Redistributions in binary form must reproduce the above
+#     copyright notice, this list of conditions and the following
+#     disclaimer in the documentation and/or other materials provided
+#     with the distribution.
+#  3. The name of the author may not be used to endorse or promote
+#     products derived from this software without specific prior
+#     written permission.
+#
+# THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS
+# OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
+# WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+# ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY
+# DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+# DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE
+# GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+# INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
+# WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE
+# OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE,
+# EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+class Tb::NumericReader
+  include Tb::EnumerableWithEach
+
+  def initialize(get_array)
+    @get_array = get_array
+  end
+
+  def header_known?
+    false
+  end
+
+  def get_named_header
+    []
+  end
+
+  def get_hash
+    ary = @get_array.call
+    if !ary
+      return nil
+    end
+    hash = {}
+    ary.each_with_index {|v, i|
+      field = (i+1).to_s
+      hash[field] = v
+    }
+    hash
+  end
+
+  def header_and_each(header_proc)
+    header_proc.call(get_named_header) if header_proc
+    while hash = get_hash
+      yield hash
+    end
+    nil
+  end
+end

data/lib/tb/numericwriter.rb

@@ -0,0 +1,61 @@
+# lib/tb/numericwriterm.rb - writer mixin for table without header
+#
+# Copyright (C) 2014 Tanaka Akira  <akr@fsij.org>
+#
+# Redistribution and use in source and binary forms, with or without
+# modification, are permitted provided that the following conditions
+# are met:
+#
+#  1. Redistributions of source code must retain the above copyright
+#     notice, this list of conditions and the following disclaimer.
+#  2. Redistributions in binary form must reproduce the above
+#     copyright notice, this list of conditions and the following
+#     disclaimer in the documentation and/or other materials provided
+#     with the distribution.
+#  3. The name of the author may not be used to endorse or promote
+#     products derived from this software without specific prior
+#     written permission.
+#
+# THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS
+# OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
+# WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+# ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY
+# DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+# DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE
+# GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+# INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
+# WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE
+# OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE,
+# EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+require 'tempfile'
+
+class Tb::NumericWriter
+  def initialize(put_array, put_finish=nil)
+    @put_array = put_array
+    @put_finish = put_finish
+  end
+
+  def header_required?
+    false
+  end
+
+  def header_generator=(gen)
+  end
+
+  def put_hash(hash)
+    ary = []
+    hash.each {|k, v|
+      if /\A[1-9][0-9]*\z/ !~ k
+        raise ArgumentError, "numeric field name expected: #{k.inspect}"
+      end
+      ary[k.to_i-1] = v
+    }
+    @put_array.call ary
+    nil
+  end
+
+  def finish
+    @put_finish.call if @put_finish
+  end
+end

data/lib/tb/pnm.rb

@@ -1,6 +1,6 @@
 # lib/tb/pnm.rb - tools for (very small) PNM images.
 #
-# Copyright (C) 2010-2012 Tanaka Akira  <akr@fsij.org>
+# Copyright (C) 2010-2014 Tanaka Akira  <akr@fsij.org>
 #
 # Redistribution and use in source and binary forms, with or without
 # modification, are permitted provided that the following conditions
@@ -28,43 +28,47 @@
 # OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE,
 # EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 
-class Tb
-  # :call-seq:
-  #   Tb.load_pnm(pnm_filename) -> tb
-  #
-  def Tb.load_pnm(pnm_filename)
-    pnm_content = File.open(pnm_filename, "rb") {|f| f.read }
-    Tb.parse_pnm(pnm_content)
-  end
+module Tb
+  # practical only for (very) small images.
+  class PNMReader < HashReader
+    def initialize(pnm_io)
+      pnm_io.binmode
+      content = pnm_io.read
+      initialize0(content)
+      header = @ary.shift
+      super lambda {
+        a = @ary.shift
+        if a
+          h = {}
+          a.each_with_index {|v, i|
+            h[header[i]] = v if !v.nil?
+          }
+          h
+        else
+          nil
+        end
+      }
+    end
 
-  # :call-seq:
-  #   Tb.parse_pnm(pnm_content) -> tb
-  #
-  def Tb.parse_pnm(pnm_content)
-    reader = PNMReader.new(pnm_content)
-    header = reader.shift
-    t = Tb.new(header)
-    reader.each {|ary|
-      t.insert_values header, ary
-    }
-    t
-  end
+    def header_known?
+      true
+    end
 
-  # :call-seq:
-  #   Tb.pnm_stream_input(pnm_io) {|ary| ... }
-  #
-  def Tb.pnm_stream_input(pnm_io)
-    pnm_io.binmode
-    content = pnm_io.read
-    PNMReader.new(content)
-  end
+    def get_named_header
+      ['type', 'x', 'y', 'component', 'value']
+    end
 
-  # practical only for (very) small images.
-  class PNMReader
     WSP = /(?:[ \t\r\n]|\#[^\r\n]*[\r\n])+/
 
-    def initialize(pnm_content)
-      pnm_content.force_encoding("ASCII-8BIT") if pnm_content.respond_to? :force_encoding
+    def initialize0(pnm_content)
+      if pnm_content.respond_to? :to_str
+        pnm_content.force_encoding("ASCII-8BIT") if pnm_content.respond_to? :force_encoding
+      else
+        # IO
+        pnm_content.binmode
+        pnm_content = pnm_content.read
+      end
+
       if /\A(P[63])(#{WSP})(\d+)(#{WSP})(\d+)(#{WSP})(\d+)[ \t\r\n]/on =~ pnm_content
         magic, wsp1, w, wsp2, h, wsp3, max, raster = $1, $2, $3.to_i, $4, $5.to_i, $6, $7.to_i, $'
         pixel_component = %w[R G B]
@@ -162,197 +166,199 @@ class Tb
       end
       iter
     end
+  end
 
-    def shift
-      @ary.shift
+  class PNMWriter
+    def initialize(io)
+      @io = io
+      @ary = []
+      @fields = {}
     end
 
-    def each
-      while ary = self.shift
-        yield ary
-      end
-      nil
+    def header_required?
+      false
     end
 
-    def to_a
-      result= []
-      each {|ary| result << ary }
-      result
+    def header_generator=(gen)
     end
-  end
 
-  # :call-seq:
-  #   generate_pnm(out='')
-  #
-  def generate_pnm(out='')
-    undefined_field = ['x', 'y', 'component', 'value'] - self.list_fields
-    if !undefined_field.empty?
-      raise ArgumentError, "field not defined: #{undefined_field.inspect[1...-1]}"
+    def put_hash(pairs)
+      @ary << pairs
+      pairs.each {|k, v|
+        @fields[k] = true
+      }
     end
-    pnm_type = nil
-    width = height = nil
-    comments = []
-    max_value = nil
-    max_x = max_y = 0
-    values = { 0.0 => true, 1.0 => true }
-    components = {}
-    self.each {|rec|
-      case rec['component']
-      when 'pnm_type'
-        pnm_type = rec['value']
-      when 'width'
-        width = rec['value'].to_i
-      when 'height'
-        height = rec['value'].to_i
-      when 'max'
-        max_value = rec['value'].to_i
-      when 'comment'
-        comments << rec['value']
-      when 'R', 'G', 'B', 'V'
-        components[rec['component']] = true
-        x = rec['x'].to_i
-        y = rec['y'].to_i
-        max_x = x if max_x < x
-        max_y = y if max_y < y
-        values[rec['value'].to_f] = true
-      end
 
-    }
-    if !(components.keys - %w[V]).empty? &&
-       !(components.keys - %w[R G B]).empty?
-      raise ArgumentError, "inconsistent color component: #{components.keys.sort.inspect[1...-1]}"
-    end
-    case pnm_type
-    when 'P1', 'P4' then raise ArgumentError, "unexpected compoenent for PBM: #{components.keys.sort.inspect[1...-1]}" if !(components.keys - %w[V]).empty?
-    when 'P2', 'P5' then raise ArgumentError, "unexpected compoenent for PGM: #{components.keys.sort.inspect[1...-1]}" if !(components.keys - %w[V]).empty?
-    when 'P3', 'P6' then raise ArgumentError, "unexpected compoenent for PPM: #{components.keys.sort.inspect[1...-1]}" if !(components.keys - %w[R G B]).empty?
-    end
-    comments.each {|c|
-      if /[\r\n]/ =~ c
-        raise ArgumentError, "comment cannot contain a newline: #{c.inspect}"
+    def finish
+      undefined_field = ['x', 'y', 'component', 'value'] - @fields.keys
+      if !undefined_field.empty?
+        raise ArgumentError, "field not defined: #{undefined_field.inspect[1...-1]}"
       end
-    }
-    if !width
-      width = max_x + 1
-    end
-    if !height
-      height = max_y + 1
-    end
-    if !max_value
-      min_interval = 1.0
-      values.keys.sort.each_cons(2) {|v1, v2|
-        d = v2-v1
-        min_interval = d if d < min_interval
+      pnm_type = nil
+      width = height = nil
+      comments = []
+      max_value = nil
+      max_x = max_y = 0
+      values = { 0.0 => true, 1.0 => true }
+      components = {}
+      @ary.each {|rec|
+        case rec['component']
+        when 'pnm_type'
+          pnm_type = rec['value']
+        when 'width'
+          width = rec['value'].to_i
+        when 'height'
+          height = rec['value'].to_i
+        when 'max'
+          max_value = rec['value'].to_i
+        when 'comment'
+          comments << rec['value']
+        when 'R', 'G', 'B', 'V'
+          components[rec['component']] = true
+          x = rec['x'].to_i
+          y = rec['y'].to_i
+          max_x = x if max_x < x
+          max_y = y if max_y < y
+          values[rec['value'].to_f] = true
+        end
       }
-      if min_interval < 0.0039 # 1/255 = 0.00392156862745098...
-        max_value = 0xffff
-      elsif min_interval < 1.0 || !(components.keys & %w[R G B]).empty?
-        max_value = 0xff
-      else
-        max_value = 1
+      if !(components.keys - %w[V]).empty? &&
+         !(components.keys - %w[R G B]).empty?
+        raise ArgumentError, "inconsistent color component: #{components.keys.sort.inspect[1...-1]}"
       end
-    end
-    if pnm_type
-      if !pnm_type.kind_of?(String) || /\AP[123456]\z/ !~ pnm_type
-        raise ArgumentError, "unexpected PNM type: #{pnm_type.inspect}"
+      case pnm_type
+      when 'P1', 'P4' then raise ArgumentError, "unexpected compoenent for PBM: #{components.keys.sort.inspect[1...-1]}" if !(components.keys - %w[V]).empty?
+      when 'P2', 'P5' then raise ArgumentError, "unexpected compoenent for PGM: #{components.keys.sort.inspect[1...-1]}" if !(components.keys - %w[V]).empty?
+      when 'P3', 'P6' then raise ArgumentError, "unexpected compoenent for PPM: #{components.keys.sort.inspect[1...-1]}" if !(components.keys - %w[R G B]).empty?
       end
-    else
-      if (components.keys - ['V']).empty?
-        if max_value == 1
-          pnm_type = 'P4' # PBM
-        else
-          pnm_type = 'P5' # PGM
+      comments.each {|c|
+        if /[\r\n]/ =~ c
+          raise ArgumentError, "comment cannot contain a newline: #{c.inspect}"
         end
-      else
-        pnm_type = 'P6' # PPM
+      }
+      if !width
+        width = max_x + 1
       end
-    end
-    header = "#{pnm_type}\n"
-    comments.each {|c| header << '#' << c << "\n" }
-    header << "#{width} #{height}\n"
-    header << "#{max_value}\n" if /P[2536]/ =~ pnm_type
-    if /P[14]/ =~ pnm_type # PBM
-      max_value = 1
-    end
-    bytes_per_component = bytes_per_line = component_fmt = component_template = nil
-    case pnm_type
-    when 'P1' then bytes_per_component = 1; raster = '1' * (width * height)
-    when 'P4' then bytes_per_line = (width + 7) / 8; raster = ["1"*width].pack("B*") * height
-    when 'P2' then bytes_per_component = max_value.to_s.length+1; component_fmt = "%#{bytes_per_component}d"; raster = (component_fmt % 0) * (width * height)
-    when 'P5' then bytes_per_component, component_template = max_value < 0x100 ? [1, 'C'] : [2, 'n']; raster = "\0" * (bytes_per_component * width * height)
-    when 'P3' then bytes_per_component = max_value.to_s.length+1; component_fmt = "%#{bytes_per_component}d"; raster = (component_fmt % 0) * (3 * width * height)
-    when 'P6' then bytes_per_component, component_template = max_value < 0x100 ? [1, 'C'] : [2, 'n']; raster = "\0" * (bytes_per_component * 3 * width * height)
-    else
-      raise
-    end
-    raster.force_encoding("ASCII-8BIT") if raster.respond_to? :force_encoding
-    self.each {|rec|
-      c = rec['component']
-      next if /\A[RGBV]\z/ !~ c
-      x = rec['x'].to_i
-      y = rec['y'].to_i
-      next if x < 0 || width <= x
-      next if y < 0 || height <= y
-      v = rec['value'].to_f
-      if v < 0
-        v = 0
-      elsif 1 < v
-        v = 1
+      if !height
+        height = max_y + 1
       end
-      case pnm_type
-      when 'P1'
-        v = v < 0.5 ? '1' : '0'
-        raster[y * width + x] = v
-      when 'P4'
-        xhi, xlo = x.divmod(8)
-        i = y * bytes_per_line + xhi
-        byte = raster[i].ord
-        if v < 0.5
-          byte |= 0x80 >> xlo
+      if !max_value
+        min_interval = 1.0
+        values.keys.sort.each_cons(2) {|v1, v2|
+          d = v2-v1
+          min_interval = d if d < min_interval
+        }
+        if min_interval < 0.0039 # 1/255 = 0.00392156862745098...
+          max_value = 0xffff
+        elsif min_interval < 1.0 || !(components.keys & %w[R G B]).empty?
+          max_value = 0xff
         else
-          byte &= 0xff7f >> xlo
+          max_value = 1
         end
-        raster[i] = [byte].pack("C")
-      when 'P2'
-        v = (v * max_value).round
-        raster[(y * width + x) * bytes_per_component, bytes_per_component] = component_fmt % v
-      when 'P5'
-        v = (v * max_value).round
-        raster[(y * width + x) * bytes_per_component, bytes_per_component] = [v].pack(component_template)
-      when 'P3'
-        v = (v * max_value).round
-        i = (y * width + x) * 3
-        if c == 'G' then i += 1
-        elsif c == 'B' then i += 2
+      end
+      if pnm_type
+        if !pnm_type.kind_of?(String) || /\AP[123456]\z/ !~ pnm_type
+          raise ArgumentError, "unexpected PNM type: #{pnm_type.inspect}"
         end
-        raster[i * bytes_per_component, bytes_per_component] = component_fmt % v
-      when 'P6'
-        v = (v * max_value).round
-        i = (y * width + x) * 3
-        if c == 'G' then i += 1
-        elsif c == 'B' then i += 2
+      else
+        if (components.keys - ['V']).empty?
+          if max_value == 1
+            pnm_type = 'P4' # PBM
+          else
+            pnm_type = 'P5' # PGM
+          end
+        else
+          pnm_type = 'P6' # PPM
         end
-        raster[i * bytes_per_component, bytes_per_component] = [v].pack(component_template)
+      end
+      header = "#{pnm_type}\n"
+      comments.each {|c| header << '#' << c << "\n" }
+      header << "#{width} #{height}\n"
+      header << "#{max_value}\n" if /P[2536]/ =~ pnm_type
+      if /P[14]/ =~ pnm_type # PBM
+        max_value = 1
+      end
+      bytes_per_component = bytes_per_line = component_fmt = component_template = nil
+      case pnm_type
+      when 'P1' then bytes_per_component = 1; raster = '1' * (width * height)
+      when 'P4' then bytes_per_line = (width + 7) / 8; raster = ["1"*width].pack("B*") * height
+      when 'P2' then bytes_per_component = max_value.to_s.length+1; component_fmt = "%#{bytes_per_component}d"; raster = (component_fmt % 0) * (width * height)
+      when 'P5' then bytes_per_component, component_template = max_value < 0x100 ? [1, 'C'] : [2, 'n']; raster = "\0" * (bytes_per_component * width * height)
+      when 'P3' then bytes_per_component = max_value.to_s.length+1; component_fmt = "%#{bytes_per_component}d"; raster = (component_fmt % 0) * (3 * width * height)
+      when 'P6' then bytes_per_component, component_template = max_value < 0x100 ? [1, 'C'] : [2, 'n']; raster = "\0" * (bytes_per_component * 3 * width * height)
       else
         raise
       end
-    }
-    if pnm_type == 'P1'
-      raster.gsub!(/[01]{#{width}}/, "\\&\n")
-      if 70 < width
-        raster.gsub!(/[01]{70}/, "\\&\n")
-      end
-      raster << "\n" if /\n\z/ !~ raster
-    elsif /P[23]/ =~ pnm_type
-      components_per_line = /P2/ =~ pnm_type ? width : 3 * width
-      raster.gsub!(/  +/, ' ')
-      raster.gsub!(/( \d+){#{components_per_line}}/, "\\&\n")
-      raster.gsub!(/(\A|\n) +/, '\1')
-      raster.gsub!(/.{71,}\n/) {
-        $&.gsub(/(.{1,69})[ \n]/, "\\1\n")
+      raster.force_encoding("ASCII-8BIT") if raster.respond_to? :force_encoding
+      @ary.each {|rec|
+        c = rec['component']
+        next if /\A[RGBV]\z/ !~ c
+        x = rec['x'].to_i
+        y = rec['y'].to_i
+        next if x < 0 || width <= x
+        next if y < 0 || height <= y
+        v = rec['value'].to_f
+        if v < 0
+          v = 0
+        elsif 1 < v
+          v = 1
+        end
+        case pnm_type
+        when 'P1'
+          v = v < 0.5 ? '1' : '0'
+          raster[y * width + x] = v
+        when 'P4'
+          xhi, xlo = x.divmod(8)
+          i = y * bytes_per_line + xhi
+          byte = raster[i].ord
+          if v < 0.5
+            byte |= 0x80 >> xlo
+          else
+            byte &= 0xff7f >> xlo
+          end
+          raster[i] = [byte].pack("C")
+        when 'P2'
+          v = (v * max_value).round
+          raster[(y * width + x) * bytes_per_component, bytes_per_component] = component_fmt % v
+        when 'P5'
+          v = (v * max_value).round
+          raster[(y * width + x) * bytes_per_component, bytes_per_component] = [v].pack(component_template)
+        when 'P3'
+          v = (v * max_value).round
+          i = (y * width + x) * 3
+          if c == 'G' then i += 1
+          elsif c == 'B' then i += 2
+          end
+          raster[i * bytes_per_component, bytes_per_component] = component_fmt % v
+        when 'P6'
+          v = (v * max_value).round
+          i = (y * width + x) * 3
+          if c == 'G' then i += 1
+          elsif c == 'B' then i += 2
+          end
+          raster[i * bytes_per_component, bytes_per_component] = [v].pack(component_template)
+        else
+          raise
+        end
       }
-      raster << "\n" if /\n\z/ !~ raster
+      if pnm_type == 'P1'
+        raster.gsub!(/[01]{#{width}}/, "\\&\n")
+        if 70 < width
+          raster.gsub!(/[01]{70}/, "\\&\n")
+        end
+        raster << "\n" if /\n\z/ !~ raster
+      elsif /P[23]/ =~ pnm_type
+        components_per_line = /P2/ =~ pnm_type ? width : 3 * width
+        raster.gsub!(/  +/, ' ')
+        raster.gsub!(/( \d+){#{components_per_line}}/, "\\&\n")
+        raster.gsub!(/(\A|\n) +/, '\1')
+        raster.gsub!(/.{71,}\n/) {
+          $&.gsub(/(.{1,69})[ \n]/, "\\1\n")
+        }
+        raster << "\n" if /\n\z/ !~ raster
+      end
+      @io << (header+raster)
     end
-    out << (header+raster)
   end
+
 end