rgeo 0.3.3 → 0.3.4

data/lib/rgeo/cartesian/feature_methods.rb

@@ -1,15 +1,15 @@
 # -----------------------------------------------------------------------------
-# 
+#
 # Cartesian common methods
-# 
+#
 # -----------------------------------------------------------------------------
-# Copyright 2010 Daniel Azuma
-# 
+# Copyright 2010-2012 Daniel Azuma
+#
 # All rights reserved.
-# 
+#
 # Redistribution and use in source and binary forms, with or without
 # modification, are permitted provided that the following conditions are met:
-# 
+#
 # * Redistributions of source code must retain the above copyright notice,
 #   this list of conditions and the following disclaimer.
 # * Redistributions in binary form must reproduce the above copyright notice,
@@ -18,7 +18,7 @@
 # * Neither the name of the copyright holder, nor the names of any other
 #   contributors to this software, may be used to endorse or promote products
 #   derived from this software without specific prior written permission.
-# 
+#
 # THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
 # AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
 # IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
@@ -35,29 +35,29 @@
 
 
 module RGeo
-  
+
   module Cartesian
-    
-    
+
+
     module GeometryMethods  # :nodoc:
-      
-      
+
+
       def srid
         factory.srid
       end
-      
-      
+
+
       def envelope
         BoundingBox.new(factory).add(self).to_geometry
       end
-      
-      
+
+
     end
-    
-    
+
+
     module PointMethods  # :nodoc:
-      
-      
+
+
       def distance(rhs_)
         rhs_ = ::RGeo::Feature.cast(rhs_, @factory)
         case rhs_
@@ -69,14 +69,14 @@ module RGeo
           super
         end
       end
-      
-      
+
+
     end
-    
-    
+
+
     module LineStringMethods  # :nodoc:
-      
-      
+
+
       def _segments
         unless @segments
           @segments = (0..num_points-2).map do |i_|
@@ -85,8 +85,8 @@ module RGeo
         end
         @segments
       end
-      
-      
+
+
       def is_simple?
         segs_ = _segments
         len_ = segs_.length
@@ -111,15 +111,15 @@ module RGeo
         end
         true
       end
-      
-      
+
+
       def length
         @segments.inject(0.0){ |sum_, seg_| sum_ + seg_.length }
       end
-      
-      
+
+
     end
-    
+
   end
-  
+
 end

data/lib/rgeo/cartesian/interface.rb

@@ -1,15 +1,15 @@
 # -----------------------------------------------------------------------------
-# 
+#
 # Cartesian toplevel interface
-# 
+#
 # -----------------------------------------------------------------------------
-# Copyright 2010 Daniel Azuma
-# 
+# Copyright 2010-2012 Daniel Azuma
+#
 # All rights reserved.
-# 
+#
 # Redistribution and use in source and binary forms, with or without
 # modification, are permitted provided that the following conditions are met:
-# 
+#
 # * Redistributions of source code must retain the above copyright notice,
 #   this list of conditions and the following disclaimer.
 # * Redistributions in binary form must reproduce the above copyright notice,
@@ -18,7 +18,7 @@
 # * Neither the name of the copyright holder, nor the names of any other
 #   contributors to this software, may be used to endorse or promote products
 #   derived from this software without specific prior written permission.
-# 
+#
 # THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
 # AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
 # IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
@@ -35,15 +35,15 @@
 
 
 module RGeo
-  
+
   module Cartesian
-    
+
     class << self
-      
-      
+
+
       # Creates and returns a cartesian factory of the preferred
       # Cartesian implementation.
-      # 
+      #
       # The actual implementation returned depends on which ruby
       # interpreter is running and what libraries are available.
       # RGeo will try to provide a fully-functional and performant
@@ -51,13 +51,13 @@ module RGeo
       # implementation will be returned.
       # In practice, this means it returns a Geos implementation if
       # available; otherwise it falls back to the simple implementation.
-      # 
+      #
       # The given options are passed to the factory's constructor.
       # What options are available depends on the particular
       # implementation. See RGeo::Geos.factory and
       # RGeo::Cartesian.simple_factory for details. Unsupported options
       # are ignored.
-      
+
       def preferred_factory(opts_={})
         if ::RGeo::Geos.supported?
           ::RGeo::Geos.factory(opts_)
@@ -66,14 +66,14 @@ module RGeo
         end
       end
       alias_method :factory, :preferred_factory
-      
-      
+
+
       # Returns a factory for the simple Cartesian implementation. This
       # implementation provides all SFS 1.1 types, and also allows Z and
       # M coordinates. It does not depend on external libraries, and is
       # thus always available, but it does not implement many of the more
       # advanced geometric operations. These limitations are:
-      # 
+      #
       # * Relational operators such as Feature::Geometry#intersects? are
       #   not implemented for most types.
       # * Relational constructors such as Feature::Geometry#union are
@@ -85,12 +85,12 @@ module RGeo
       # * Equality and simplicity evaluation are implemented for some but
       #   not all types.
       # * Assertions for polygons and multipolygons are not implemented.
-      # 
+      #
       # Unimplemented operations may raise Error::UnsupportedOperation
       # if invoked.
-      # 
+      #
       # Options include:
-      # 
+      #
       # [<tt>:srid</tt>]
       #   Set the SRID returned by geometries created by this factory.
       #   Default is 0.
@@ -130,41 +130,41 @@ module RGeo
       #   configuration parameters for WKRep::WKTGenerator.new.
       #   Default is the empty hash, indicating the default configuration
       #   for WKRep::WKBGenerator.
-      
+
       def simple_factory(opts_={})
         Cartesian::Factory.new(opts_)
       end
-      
-      
+
+
       # Returns a Feature::FactoryGenerator that creates preferred
       # factories. The given options are used as the default options.
-      # 
+      #
       # A common case for this is to provide the <tt>:srs_database</tt>
       # as a default. Then, the factory generator need only be passed
       # an SRID and it will automatically fetch the appropriate Proj4
       # and CoordSys objects.
-      
+
       def preferred_factory_generator(defaults_={})
         ::Proc.new{ |c_| preferred_factory(defaults_.merge(c_)) }
       end
       alias_method :factory_generator, :preferred_factory_generator
-      
-      
+
+
       # Returns a Feature::FactoryGenerator that creates simple factories.
       # The given options are used as the default options.
-      # 
+      #
       # A common case for this is to provide the <tt>:srs_database</tt>
       # as a default. Then, the factory generator need only be passed
       # an SRID and it will automatically fetch the appropriate Proj4
       # and CoordSys objects.
-      
+
       def simple_factory_generator(defaults_={})
         ::Proc.new{ |c_| simple_factory(defaults_.merge(c_)) }
       end
-      
-      
+
+
     end
-    
+
   end
-  
+
 end

data/lib/rgeo/coord_sys.rb

@@ -1,15 +1,15 @@
 # -----------------------------------------------------------------------------
-# 
+#
 # Coordinate systems for RGeo
-# 
+#
 # -----------------------------------------------------------------------------
-# Copyright 2010 Daniel Azuma
-# 
+# Copyright 2010-2012 Daniel Azuma
+#
 # All rights reserved.
-# 
+#
 # Redistribution and use in source and binary forms, with or without
 # modification, are permitted provided that the following conditions are met:
-# 
+#
 # * Redistributions of source code must retain the above copyright notice,
 #   this list of conditions and the following disclaimer.
 # * Redistributions in binary form must reproduce the above copyright notice,
@@ -18,7 +18,7 @@
 # * Neither the name of the copyright holder, nor the names of any other
 #   contributors to this software, may be used to endorse or promote products
 #   derived from this software without specific prior written permission.
-# 
+#
 # THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
 # AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
 # IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
@@ -35,32 +35,32 @@
 
 
 module RGeo
-  
-  
+
+
   # This module provides data structures and tools related to coordinate
   # systems and coordinate transforms. It comprises the following parts:
-  # 
+  #
   # RGeo::CoordSys::Proj4 is a wrapper around the proj4 library, which
   # defines a commonly-used syntax for specifying geographic and projected
   # coordinate systems, and performs coordinate transformations.
-  # 
+  #
   # The RGeo::CoordSys::CS module contains an implementation of the CS
   # (coordinate systems) package of the OGC Coordinate Transform spec.
   # This includes classes for representing ellipsoids, datums, coordinate
   # systems, and other related concepts, as well as a parser for the WKT
   # format for specifying coordinate systems.
-  # 
+  #
   # The RGeo::CoordSys::SRSDatabase module contains tools for accessing
   # spatial reference databases, from which you can look up coordinate
   # system specifications. You can access the <tt>spatial_ref_sys</tt>
   # table provided with OGC-compliant spatial databases such as PostGIS,
   # read the databases provided with the proj4 library, or access URLs
   # such as those provided by spatialreference.org.
-  
+
   module CoordSys
   end
-  
-  
+
+
 end
 
 

data/lib/rgeo/coord_sys/cs/entities.rb

@@ -1,16 +1,15 @@
-# -*- coding: utf-8 -*-
 # -----------------------------------------------------------------------------
-# 
+#
 # OGC CS objects for RGeo
-# 
+#
 # -----------------------------------------------------------------------------
-# Copyright 2010 Daniel Azuma
-# 
+# Copyright 2010-2012 Daniel Azuma
+#
 # All rights reserved.
-# 
+#
 # Redistribution and use in source and binary forms, with or without
 # modification, are permitted provided that the following conditions are met:
-# 
+#
 # * Redistributions of source code must retain the above copyright notice,
 #   this list of conditions and the following disclaimer.
 # * Redistributions in binary form must reproduce the above copyright notice,
@@ -19,7 +18,7 @@
 # * Neither the name of the copyright holder, nor the names of any other
 #   contributors to this software, may be used to endorse or promote products
 #   derived from this software without specific prior written permission.
-# 
+#
 # THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
 # AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
 # IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
@@ -36,57 +35,57 @@
 
 
 module RGeo
-  
+
   module CoordSys
-    
-    
+
+
     module CS
-      
-      
+
+
       # An axis orientation constant for AxisInfo.
       # Unknown or unspecified axis orientation. This can be used for
       # local or fitted coordinate systems.
       AO_OTHER = 0
-      
+
       # An axis orientation constant for AxisInfo.
       # Increasing ordinates values go North. This is usually used for
       # Grid Y coordinates and Latitude.
       AO_NORTH = 1
-      
+
       # An axis orientation constant for AxisInfo.
       # Increasing ordinates values go South. This is rarely used.
       AO_SOUTH = 2
-      
+
       # An axis orientation constant for AxisInfo.
       # Increasing ordinates values go East. This is rarely used.
       AO_EAST = 3
-      
+
       # An axis orientation constant for AxisInfo.
       # Increasing ordinates values go West. This is usually used for
       # Grid X coordinates and Longitude.
       AO_WEST = 4
-      
+
       # An axis orientation constant for AxisInfo.
       # Increasing ordinates values go up. This is used for vertical
       # coordinate systems.
       AO_UP = 5
-      
+
       # An axis orientation constant for AxisInfo.
       # Increasing ordinates values go down. This is used for vertical
       # coordinate systems.
       AO_DOWN = 6
-      
-      
+
+
       # A datum type constant for HorizontalDatum.
       # Lowest possible value for horizontal datum types.
       HD_MIN = 1000
-      
+
       # A datum type constant for HorizontalDatum.
       # Unspecified horizontal datum type. Horizontal datums with this
       # type should never supply a conversion to WGS84 using Bursa Wolf
       # parameters.
       HD_OTHER = 1000
-      
+
       # A datum type constant for HorizontalDatum.
       # These datums, such as ED50, NAD27 and NAD83, have been designed
       # to support horizontal positions on the ellipsoid as opposed to
@@ -94,7 +93,7 @@ module RGeo
       # support a horizontal component of a position in a domain of
       # limited extent, such as a country, a region or a continent.
       HD_CLASSIC = 1001
-      
+
       # A datum type constant for HorizontalDatum.
       # A geocentric datum is a "satellite age" modern geodetic datum
       # mainly of global extent, such as WGS84 (used in GPS), PZ90 (used
@@ -103,30 +102,30 @@ module RGeo
       # position (through ellipsoidal heights). The regional realizations
       # of ITRF, such as ETRF, are also included in this category.
       HD_GEOCENTRIC = 1002
-      
+
       # A datum type constant for HorizontalDatum.
       # Highest possible value for horizontal datum types.
       HD_MAX = 1999
-      
+
       # A datum type constant for VerticalDatum.
       # Lowest possible value for vertical datum types.
       VD_MIN = 2000
-      
+
       # A datum type constant for VerticalDatum.
       # Unspecified vertical datum type.
       VD_OTHER = 2000
-      
+
       # A datum type constant for VerticalDatum.
       # A vertical datum for orthometric heights that are measured along
       # the plumb line.
       VD_ORTHOMETRIC = 2001
-      
+
       # A datum type constant for VerticalDatum.
       # A vertical datum for ellipsoidal heights that are measured along
       # the normal to the ellipsoid used in the definition of horizontal
       # datum.
       VD_ELLIPSOIDAL = 2002
-      
+
       # A datum type constant for VerticalDatum.
       # The vertical datum of altitudes or heights in the atmosphere.
       # These are approximations of orthometric heights obtained with the
@@ -135,11 +134,11 @@ module RGeo
       # millibars (used to measure pressure levels), or theta value (units
       # used to measure geopotential height).
       VD_ALTITUDE_BAROMETRIC = 2003
-      
+
       # A datum type constant for VerticalDatum.
       # A normal height system.
       VD_NORMAL = 2004
-      
+
       # A datum type constant for VerticalDatum.
       # A vertical datum of geoid model derived heights, also called
       # GPS-derived heights. These heights are approximations of
@@ -147,7 +146,7 @@ module RGeo
       # (h) by the use of the given geoid undulation model (N) through
       # the equation: H=h-N.
       VD_GEOID_MODE_DERIVED = 2005
-      
+
       # A datum type constant for VerticalDatum.
       # This attribute is used to support the set of datums generated for
       # hydrographic engineering projects where depth measurements below
@@ -156,61 +155,61 @@ module RGeo
       # (approximately) to the actual equipotential surfaces of the
       # earth's gravity field, using such procedures as echo-sounding.
       VD_DEPTH = 2006
-      
+
       # A datum type constant for VerticalDatum.
       # Highest possible value for vertical datum types.
       VD_MAX = 2999
-      
+
       # A datum type constant for LocalDatum.
       # Lowest possible value for local datum types.
       LD_MIN = 10000
-      
+
       # A datum type constant for LocalDatum.
       # Highest possible value for local datum types.
       LD_MAX = 32767
-      
-      
+
+
       # This is a base class for all OGC coordinate system objects.
       # This includes both interfaces and data types from the OGC
       # Coordinate Transformation spec.
-      # 
+      #
       # This is a non-instantiable abstract class.
-      
+
       class Base
-        
-        
+
+
         def inspect  # :nodoc:
           "#<#{self.class}:0x#{object_id.to_s(16)} #{to_wkt}>"
         end
-        
-        
+
+
         # Tests for equality. Two objects are defined as equal if they
         # have the same type (class) and the same WKT representation.
-        
+
         def eql?(rhs_)
           rhs_.class == self.class && rhs_.to_wkt == self.to_wkt
         end
         alias_method :==, :eql?
-        
-        
+
+
         # Returns the default WKT representation.
-        
+
         def to_s
           to_wkt
         end
-        
-        
+
+
         # Computes the WKT representation. Options include:
-        # 
+        #
         # [<tt>:standard_brackets</tt>]
         #   If set to true, outputs parentheses rather than square
         #   brackets. Default is false.
-        
+
         def to_wkt(opts_={})
           opts_[:standard_brackets] ? _to_wkt('(', ')') : _to_wkt('[', ']')
         end
-        
-        
+
+
         def _to_wkt(open_, close_)  # :nodoc:
           content_ = _wkt_content(open_, close_).map{ |obj_| ",#{obj_}" }.join
           if defined?(@authority) && @authority
@@ -220,14 +219,14 @@ module RGeo
           end
           "#{_wkt_typename}#{open_}#{@name.inspect}#{content_}#{authority_}#{close_}"
         end
-        
-        
+
+
         # Marshal support
-        
+
         def marshal_dump  # :nodoc:
           to_wkt
         end
-        
+
         def marshal_load(data_)  # :nodoc:
           data_ = data_['wkt'] if data_.is_a?(::Hash)
           temp_ = CS.create_from_wkt(data_)
@@ -239,10 +238,10 @@ module RGeo
             raise ::TypeError, 'Bad Marshal data'
           end
         end
-        
-        
+
+
         # Psych support
-        
+
         def init_with(coder_)  # :nodoc:
           temp_ = CS.create_from_wkt(coder_.type == :scalar ? coder_.scalar : coder_['wkt'] )
           if temp_.class == self.class
@@ -253,33 +252,33 @@ module RGeo
             raise ::TypeError, 'Bad YAML data'
           end
         end
-        
+
         def encode_with(coder_)  # :nodoc:
           coder_['wkt'] = to_wkt
         end
-        
-        
+
+
         class << self
-          
+
           private :new
-          
+
         end
-        
-        
+
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # Details of axis. This is used to label axes, and indicate the
       # orientation.
-      
+
       class AxisInfo < Base
-        
+
         # :stopdoc:
         NAMES_BY_VALUE = ['OTHER', 'NORTH', 'SOUTH', 'EAST', 'WEST', 'UP', 'DOWN']
         # :startdoc:
-        
+
         def initialize(name_, orientation_)  # :nodoc:
           @name = name_
           case orientation_
@@ -289,28 +288,28 @@ module RGeo
             @orientation = orientation_.to_i
           end
         end
-        
-        
+
+
         # Human readable name for axis. Possible values are "X", "Y",
         # "Long", "Lat" or any other short string.
         attr_reader :name
-        
+
         # Gets enumerated value for orientation.
         attr_reader :orientation
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "AXIS"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           [NAMES_BY_VALUE[@orientation]]
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Creates an AxisInfo. you must pass the human readable name for
           # the axis (e.g. "X", "Y", "Long", "Lat", or other short string)
           # and either an integer orientation code or a string. Possible
@@ -318,73 +317,73 @@ module RGeo
           # "<tt>SOUTH</tt>", "<tt>EAST</tt>", "<tt>WEST</tt>",
           # "<tt>UP</tt>", and "<tt>DOWN</tt>", or the corresponding
           # integer values 0-5.
-          
+
           def create(name_, orientation_)
             new(name_, orientation_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # A named projection parameter value. The linear units of
       # parameters' values match the linear units of the containing
       # projected coordinate system. The angular units of parameter
       # values match the angular units of the geographic coordinate
       # system that the projected coordinate system is based on.
-      
+
       class ProjectionParameter < Base
-        
+
         def initialize(name_, value_)  # :nodoc:
           @name = name_
           @value = value_.to_f
         end
-        
-        
+
+
         # The parameter name.
         attr_reader :name
-        
+
         # The parameter value.
         attr_reader :value
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "PARAMETER"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           [@value]
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a parameter given the name and value.
-          
+
           def create(name_, value_)
             new(name_, value_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # Parameters for a geographic transformation into WGS84. The Bursa
       # Wolf parameters should be applied to geocentric coordinates, where
       # the X axis points towards the Greenwich Prime Meridian, the Y axis
       # points East, and the Z axis points North.
-      
+
       class WGS84ConversionInfo < Base
-        
+
         def initialize(dx_, dy_, dz_, ex_, ey_, ez_, ppm_)  # :nodoc:
           @dx = dx_.to_f
           @dy = dy_.to_f
@@ -394,65 +393,65 @@ module RGeo
           @ez = ez_.to_f
           @ppm = ppm_.to_f
         end
-        
-        
+
+
         # Bursa Wolf shift in meters.
         attr_reader :dx
-        
+
         # Bursa Wolf shift in meters.
         attr_reader :dy
-        
+
         # Bursa Wolf shift in meters.
         attr_reader :dz
-        
+
         # Bursa Wolf rotation in arc seconds.
         attr_reader :ex
-        
+
         # Bursa Wolf rotation in arc seconds.
         attr_reader :ey
-        
+
         # Bursa Wolf rotation in arc seconds.
         attr_reader :ez
-        
+
         # Bursa Wolf scaling in in parts per million.
         attr_reader :ppm
-        
-        
+
+
         def _to_wkt(open_, close_)  # :nodoc:
           "TOWGS84#{open_}#{@dx},#{@dy},#{@dz},#{@ex},#{@ey},#{@ez},#{@ppm}#{close_}"
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create the horizontal datum shift transformation into WGS84,
           # given the seven Bursa Wolf parameters.
           # The Bursa Wolf shift should be in meters, the rotation in arc
           # seconds, and the scaling in parts per million.
-          
+
           def create(dx_, dy_, dz_, ex_, ey_, ez_, ppm_)
             new(dx_, dy_, dz_, ex_, ey_, ez_, ppm_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # A base interface for metadata applicable to coordinate system
       # objects.
-      # 
+      #
       # The metadata items "Abbreviation"’, "Alias", "Authority",
       # "AuthorityCode", "Name" and "Remarks" were specified in the Simple
       # Features interfaces, so they have been kept here.
-      # 
+      #
       # This specification does not dictate what the contents of these
       # items should be. However, the following guidelines are suggested:
-      # 
+      #
       # When CS_CoordinateSystemAuthorityFactory is used to create an
       # object, the "Authority" and "AuthorityCode" values should be set
       # to the authority name of the factory object, and the authority
@@ -460,27 +459,27 @@ module RGeo
       # may not be set. (If the authority is EPSG, the implementer may
       # consider using the corresponding metadata values in the EPSG
       # tables.)
-      # 
+      #
       # When CS_CoordinateSystemFactory creates an object, the "Name"
       # should be set to the value supplied by the client. All of the
       # other metadata items should be left empty.
-      # 
+      #
       # == Notes
-      # 
+      #
       # This is a non-instantiable abstract class.
-      # 
+      #
       # Most subclasses will have a set of optional parameters in their
       # "create" method to set the metadata fields. These parameters are,
       # in order:
-      # 
+      #
       # * <b>authority</b>: authority name
       # * <b>authority_code</b>: authority-specific identification code
       # * <b>abbreviation</b>: an abbreviation
       # * <b>alias</b>: an alias
       # * <b>remarks</b>: provider-supplied remarks.
-      
+
       class Info < Base
-        
+
         def initialize(name_, authority_=nil, authority_code_=nil, abbreviation_=nil, alias_=nil, remarks_=nil)  # :nodoc:
           @name = name_
           @authority = authority_ ? authority_.to_s : nil
@@ -489,14 +488,14 @@ module RGeo
           @alias = alias_ ? alias_.to_s : nil
           @remarks = remarks_ ? remarks_.to_s : nil
         end
-        
-        
+
+
         # Gets the abbreviation.
         attr_reader :abbreviation
-        
+
         # Gets the alias.
         attr_reader :alias
-        
+
         # Gets the authority name.
         # An Authority is an organization that maintains definitions of
         # Authority Codes. For example the European Petroleum Survey Group
@@ -505,7 +504,7 @@ module RGeo
         # ID. For example, the EPSG code for a WGS84 Lat/Lon coordinate
         # system is "4326".
         attr_reader :authority
-        
+
         # Gets the authority-specific identification code.
         # The AuthorityCode is a compact string defined by an Authority to
         # reference a particular spatial reference object. For example,
@@ -514,187 +513,187 @@ module RGeo
         # consist of a few digits. The EPSG code for WGS84 Lat/Lon is
         # "4326".
         attr_reader :authority_code
-        
+
         # Gets the name.
         attr_reader :name
-        
+
         # Gets the provider-supplied remarks.
         attr_reader :remarks
-        
-        
+
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # Base interface for defining units.
-      # 
+      #
       # == Notes
-      # 
+      #
       # Normally, you will instantiate one of the subclasses LinearUnit or
       # AngularUnit. However, it is possible to instantiate Unit if it is
       # not clear whether the data refers to a LinearUnit or AngularUnit.
-      
+
       class Unit < Info
-        
+
         def initialize(name_, conversion_factor_, *optional_)  # :nodoc:
           super(name_, *optional_)
           @conversion_factor = conversion_factor_.to_f
         end
-        
-        
+
+
         # This field is not part of the OGC CT spec, but is part of the
         # SFS. It is an alias of the appropriate field in the subclass,
         # i.e. LinearUnit#meters_per_unit or AngularUnit#radians_per_unit.
         attr_reader :conversion_factor
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "UNIT"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           [@conversion_factor]
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a bare Unit that does not specify whether it is a
           # LinearUnit or an AngularUnit, given a unit name and a
           # conversion factor. You may also provide the optional
           # parameters specified by the Info interface.
-          
+
           def create(name_, conversion_factor_, *optional_)
             new(name_, conversion_factor_, *optional_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # Definition of linear units.
-      
+
       class LinearUnit < Unit
-        
-        
+
+
         # Returns the number of meters per LinearUnit.
         # Also available as Unit#conversion_factor.
-        
+
         def meters_per_unit
           @conversion_factor
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a LinearUnit given a unit name and a conversion factor
           # in meters per unit. You may also provide the optional
           # parameters specified by the Info interface.
-          
+
           def create(name_, meters_per_unit_, *optional_)
             new(name_, meters_per_unit_, *optional_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # Definition of angular units.
-      
+
       class AngularUnit < Unit
-        
-        
+
+
         # Returns the number of radians per AngularUnit.
         # Also available as Unit#conversion_factor.
-        
+
         def radians_per_unit
           @conversion_factor
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create an AngularUnit given a unit name and a conversion
           # factor in radians per unit. You may also provide the optional
           # parameters specified by the Info interface.
-          
+
           def create(name_, radians_per_unit_, *optional_)
             new(name_, radians_per_unit_, *optional_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # A meridian used to take longitude measurements from.
-      
+
       class PrimeMeridian < Info
-        
+
         def initialize(name_, angular_unit_, longitude_, *optional_)  # :nodoc:
           super(name_, *optional_)
           @angular_unit = angular_unit_
           @longitude = longitude_.to_f
         end
-        
-        
+
+
         # Returns the AngularUnits.
         attr_reader :angular_unit
-        
+
         # Returns the longitude value relative to the Greenwich Meridian.
         # The longitude is expressed in this objects angular units.
         attr_reader :longitude
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "PRIMEM"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           [@longitude]
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a PrimeMeridian given a name, AngularUnits, and the
           # longitude relative to the Greenwich Meridian, expressed in
           # the AngularUnits. You may also provide the optional parameters
           # specified by the Info interface.
-          
+
           def create(name_, angular_unit_, longitude_, *optional_)
             new(name_, angular_unit_, longitude_, *optional_)
           end
-          
-          
+
+
         end
-        
-        
+
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # An approximation of the Earth's surface as a squashed sphere.
-      
+
       class Ellipsoid < Info
-        
+
         def initialize(name_, semi_major_axis_, semi_minor_axis_, inverse_flattening_, ivf_definitive_, linear_unit_, *optional_)  # :nodoc:
           super(name_, *optional_)
           @semi_major_axis = semi_major_axis_.to_f
@@ -703,64 +702,64 @@ module RGeo
           @ivf_definitive = ivf_definitive_ ? true : false
           @linear_unit = linear_unit_
         end
-        
-        
+
+
         # Gets the equatorial radius. The returned length is expressed in
         # this object's axis units.
         attr_reader :semi_major_axis
-        
+
         # Gets the polar radius. The returned length is expressed in this
         # object's axis units.
         attr_reader :semi_minor_axis
-        
+
         # Returns the value of the inverse of the flattening constant. The
         # inverse flattening is related to the equatorial/polar radius by
         # the formula ivf=re/(re-rp). For perfect spheres, this formula
         # breaks down, and a special IVF value of zero is used.
         attr_reader :inverse_flattening
-        
+
         # Is the Inverse Flattening definitive for this ellipsoid? Some
         # ellipsoids use the IVF as the defining value, and calculate the
         # polar radius whenever asked. Other ellipsoids use the polar
         # radius to calculate the IVF whenever asked. This distinction can
         # be important to avoid floating-point rounding errors.
         attr_reader :ivf_definitive
-        
+
         # Returns the LinearUnit. The units of the semi-major and
         # semi-minor axis values.
         attr_reader :axis_unit
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "SPHEROID"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           [@semi_major_axis, @inverse_flattening]
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create an Ellipsoid given a name, semi-major and semi-minor
           # axes, the inverse flattening, a boolean indicating whether
           # the inverse flattening is definitive, and the LinearUnit
           # indicating the axis units. The LinearUnit is optional and
           # may be set to nil. You may also provide the optional parameters
           # specified by the Info interface.
-          
+
           def create(name_, semi_major_axis_, semi_minor_axis_, inverse_flattening_, ivf_definitive_, linear_unit_, *optional_)
             new(name_, semi_major_axis_, semi_minor_axis_, inverse_flattening_, ivf_definitive_, linear_unit_, *optional_)
           end
-          
-          
+
+
           # Create an Ellipsoid given a name, semi-major and semi-minor
           # axes, and the LinearUnit indicating the axis units. In the
           # resulting ellipsoid, the inverse flattening is not definitive.
           # The LinearUnit is optional and may be set to nil. You may also
           # provide the optional parameters specified by the Info interface.
-          
+
           def create_ellipsoid(name_, semi_major_axis_, semi_minor_axis_, linear_unit_, *optional_)
             semi_major_axis_ = semi_major_axis_.to_f
             semi_minor_axis_ = semi_minor_axis_.to_f
@@ -768,14 +767,14 @@ module RGeo
             inverse_flattening_ = 0.0 if inverse_flattening_.infinite?
             new(name_, semi_major_axis_, semi_minor_axis_, inverse_flattening_, false, linear_unit_, *optional_)
           end
-          
-          
+
+
           # Create an Ellipsoid given a name, semi-major axis, inverse
           # flattening, and the LinearUnit indicating the axis units. In
           # the resulting ellipsoid, the inverse flattening is definitive.
           # The LinearUnit is optional and may be set to nil. You may also
           # provide the optional parameters specified by the Info interface.
-          
+
           def create_flattened_sphere(name_, semi_major_axis_, inverse_flattening_, linear_unit_, *optional_)
             semi_major_axis_ = semi_major_axis_.to_f
             inverse_flattening_ = inverse_flattening_.to_f
@@ -783,15 +782,15 @@ module RGeo
             semi_minor_axis_ = semi_major_axis_ if semi_minor_axis_.infinite?
             new(name_, semi_major_axis_, semi_minor_axis_, inverse_flattening_, true, linear_unit_, *optional_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # A set of quantities from which other quantities are calculated.
       # For the OGC abstract model, it can be defined as a set of real
       # points on the earth that have coordinates. EG. A datum can be
@@ -803,344 +802,344 @@ module RGeo
       # as axis of spin). The definition of the datum may also include
       # the temporal behavior (such as the rate of change of the
       # orientation of the coordinate axes).
-      # 
+      #
       # == Notes
-      # 
+      #
       # This is a non-instantiable abstract class. You must instantiate
       # one of the subclasses HorizontalDatum, VerticalDatum, or
       # LocalDatum.
-      
+
       class Datum < Info
-        
+
         def initialize(name_, datum_type_, *optional_)  # :nodoc:
           super(name_, *optional_)
           @datum_type = datum_type_.to_i
         end
-        
-        
+
+
         # Gets the type of the datum as an enumerated code.
         attr_reader :datum_type
-        
-        
+
+
         def _wkt_content(open_, close_)  # :nodoc:
           []
         end
-        
-        
+
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # Procedure used to measure vertical distances.
-      
+
       class VerticalDatum < Datum
-        
+
         def _wkt_typename  # :nodoc:
           "VERT_DATUM"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           [@datum_type]
         end
-        
+
         class << self
-          
-          
+
+
           # Create a VerticalDatum given a name and a datum type code.
           # You may also provide the optional parameters specified by the
           # Info interface.
-          
+
           def create(name_, datum_type_, *optional_)
             new(name_, datum_type_, *optional_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # Local datum. If two local datum objects have the same datum type
       # and name, then they can be considered equal. This means that
       # coordinates can be transformed between two different local
       # coordinate systems, as long as they are based on the same local
       # datum.
-      
+
       class LocalDatum < Datum
-        
+
         def _wkt_typename  # :nodoc:
           "LOCAL_DATUM"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           [@datum_type]
         end
-        
+
         class << self
-          
-          
+
+
           # Create a LocalDatum given a name and a datum type code. You
           # may also provide the optional parameters specified by the
           # Info interface.
-          
+
           def create(name_, datum_type_, *optional_)
             new(name_, datum_type_, *optional_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # Procedure used to measure positions on the surface of the Earth.
-      
+
       class HorizontalDatum < Datum
-        
+
         def initialize(name_, datum_type_, ellipsoid_, wgs84_parameters_, *optional_)  # :nodoc:
           super(name_, datum_type_, *optional_)
           @ellipsoid = ellipsoid_
           @wgs84_parameters = wgs84_parameters_
         end
-        
-        
+
+
         # Returns the Ellipsoid.
         attr_reader :ellipsoid
-        
+
         # Gets preferred parameters for a Bursa Wolf transformation into
         # WGS84. The 7 returned values correspond to (dx,dy,dz) in meters,
         # (ex,ey,ez) in arc-seconds, and scaling in parts-per-million.
         attr_reader :wgs84_parameters
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "DATUM"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           array_ = [@ellipsoid._to_wkt(open_, close_)]
           array_ << @wgs84_parameters._to_wkt(open_, close_) if @wgs84_parameters
           array_
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a HorizontalDatum given a name, datum type code,
           # Ellipsoid, and WGS84ConversionInfo. The WGS84ConversionInfo
           # is optional and may be set to nil. You may also provide the
           # optional parameters specified by the Info interface.
-          
+
           def create(name_, datum_type_, ellipsoid_, wgs84_parameters_, *optional_)
             new(name_, datum_type_, ellipsoid_, wgs84_parameters_, *optional_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # A projection from geographic coordinates to projected coordinates.
-      
+
       class Projection < Info
-        
+
         def initialize(name_, class_name_, parameters_, *optional_)  # :nodoc:
           super(name_, *optional_)
           @class_name = class_name_.to_s
           @parameters = parameters_ ? parameters_.dup : []
         end
-        
-        
+
+
         # Gets the projection classification name
         # (e.g. "Transverse_Mercator").
         attr_reader :class_name
-        
-        
+
+
         # Gets number of parameters of the projection.
-        
+
         def num_parameters
           @parameters.size
         end
-        
-        
+
+
         # Gets an inexed parameter of the projection.
-        
+
         def get_parameter(index_)
           @parameters[index_]
         end
-        
-        
+
+
         # Iterates over the parameters of the projection.
-        
+
         def each_parameter(&block_)
           @parameters.each(&block_)
         end
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "PROJECTION"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           []
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a Projection given a name, a projection class, and an
           # array of ProjectionParameter. You may also provide the
           # optional parameters specified by the Info interface.
-          
+
           def create(name_, class_name_, parameters_, *optional_)
             new(name_, class_name_, parameters_, *optional_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # Base interface for all coordinate systems.
-      # 
+      #
       # A coordinate system is a mathematical space, where the elements
       # of the space are called positions. Each position is described by
       # a list of numbers. The length of the list corresponds to the
       # dimension of the coordinate system. So in a 2D coordinate system
       # each position is described by a list containing 2 numbers.
-      # 
+      #
       # However, in a coordinate system, not all lists of numbers
       # correspond to a position -- some lists may be outside the domain
       # of the coordinate system. For example, in a 2D Lat/Lon coordinate
       # system, the list (91,91) does not correspond to a position.
-      # 
+      #
       # Some coordinate systems also have a mapping from the mathematical
       # space into locations in the real world. So in a Lat/Lon coordinate
       # system, the mathematical position (lat, long) corresponds to a
       # location on the surface of the Earth. This mapping from the
       # mathematical space into real-world locations is called a Datum.
-      # 
+      #
       # == Notes
-      # 
+      #
       # This is a non-instantiable abstract class. You must instantiate
       # one of the subclasses GeocentricCoordinateSystem,
       # GeographicCoordinateSystem, ProjectedCoordinateSystem,
       # VerticalCoordinateSystem, LocalCoordinateSystem, or
       # CompoundCoordinateSystem.
-      
+
       class CoordinateSystem < Info
-        
+
         def initialize(name_, dimension_, *optional_)  # :nodoc:
           super(name_, *optional_)
           @dimension = dimension_.to_i
         end
-        
-        
+
+
         # Dimension of the coordinate system
         attr_reader :dimension
-        
-        
+
+
         # Gets axis details for dimension within coordinate system. Each
         # dimension in the coordinate system has a corresponding axis.
-        
+
         def get_axis(dimension_)
           nil
         end
-        
-        
+
+
         # Gets units for dimension within coordinate system. Each
         # dimension in the coordinate system has corresponding units.
-        
+
         def get_units(dimension_)
           nil
         end
-        
-        
+
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # An aggregate of two coordinate systems (CRS). One of these is
       # usually a CRS based on a two dimensional coordinate system such
       # as a geographic or a projected coordinate system with a horizontal
       # datum. The other is a vertical CRS which is a one-dimensional
       # coordinate system with a vertical datum.
-      
+
       class CompoundCoordinateSystem < CoordinateSystem
-        
+
         def initialize(name_, head_, tail_, *optional_)  # :nodoc:
           super(name_, head_.dimension + tail_.dimension, *optional_)
           @head = head_
           @tail = tail_
         end
-        
-        
+
+
         # Gets first sub-coordinate system.
         attr_reader :head
-        
+
         # Gets second sub-coordinate system.
         attr_reader :tail
-        
-        
+
+
         # Implements CoordinateSystem#get_axis
-        
+
         def get_axis(index_)
           hd_ = @head.dimension
           index_ < hd_ ? @head.get_axis(index_) : @tail.get_axis(index_ - hd_)
         end
-        
-        
+
+
         # Implements CoordinateSystem#get_units
-        
+
         def get_units(index_)
           hd_ = @head.dimension
           index_ < hd_ ? @head.get_units(index_) : @tail.get_units(index_ - hd_)
         end
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "COMPD_CS"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           [@head._to_wkt(open_, close_), @tail._to_wkt(open_, close_)]
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a CompoundCoordinateSystem given two sub-coordinate
           # systems. You may also provide the optional parameters
           # specified by the Info interface.
-          
+
           def create(name_, head_, tail_, *optional_)
             new(name_, head_, tail_, *optional_)
           end
-          
-          
+
+
         end
-        
-        
+
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # A local coordinate system, with uncertain relationship to the
       # world. In general, a local coordinate system cannot be related to
       # other coordinate systems. However, if two objects supporting this
@@ -1152,78 +1151,78 @@ module RGeo
       # have a mechanism for correlating local datums. (E.g. from a
       # database of transformations, which is created and maintained from
       # real-world measurements.)
-      # 
+      #
       # == Notes
-      # 
+      #
       # RGeo's implementation does not provide the Coordinate
       # Transformation (CT) package.
-      
+
       class LocalCoordinateSystem < CoordinateSystem
-        
+
         def initialize(name_, local_datum_, unit_, axes_, *optional_)  # :nodoc:
           super(name_, axes_.size, *optional_)
           @local_datum = local_datum_
           @unit = unit_
           @axes = axes_.dup
         end
-        
-        
+
+
         # Gets the local datum.
         attr_reader :local_datum
-        
-        
+
+
         # Implements CoordinateSystem#get_axis
-        
+
         def get_axis(index_)
           @axes[index_]
         end
-        
-        
+
+
         # Implements CoordinateSystem#get_units
-        
+
         def get_units(index_)
           @unit
         end
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "LOCAL_CS"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           [@local_datum._to_wkt(open_, close_), @unit._to_wkt(open_, close_)] + @axes.map{ |ax_| ax_._to_wkt(open_, close_) }
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a LocalCoordinateSystem given a name, a LocalDatum, a
           # Unit, and an array of at least one AxisInfo. You may also
           # provide the optional parameters specified by the Info
           # interface.
-          
+
           def create(name_, local_datum_, unit_, axes_, *optional_)
             new(name_, local_datum_, unit_, axes_, *optional_)
           end
-          
-          
+
+
         end
-        
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # A 3D coordinate system, with its origin at the centre of the
       # Earth. The X axis points towards the prime meridian. The Y axis
       # points East or West. The Z axis points North or South. By default
       # the Z axis will point North, and the Y axis will point East (e.g.
       # a right handed system), but you should check the axes for
       # non-default values.
-      
+
       class GeocentricCoordinateSystem < CoordinateSystem
-        
+
         def initialize(name_, horizontal_datum_, prime_meridian_, linear_unit_, axis0_, axis1_, axis2_, *optional_)  # :nodoc:
           super(name_, 3, *optional_)
           @horizontal_datum = horizontal_datum_
@@ -1233,39 +1232,39 @@ module RGeo
           @axis1 = axis1_
           @axis2 = axis2_
         end
-        
-        
+
+
         # Returns the HorizontalDatum. The horizontal datum is used to
         # determine where the centre of the Earth is considered to be.
         # All coordinate points will be measured from the centre of the
         # Earth, and not the surface.
         attr_reader :horizontal_datum
-        
+
         # Returns the PrimeMeridian.
         attr_reader :prime_meridian
-        
+
         # Gets the units used along all the axes.
         attr_reader :linear_unit
-        
-        
+
+
         # Implements CoordinateSystem#get_units
-        
+
         def get_units(index_)
           @linear_unit
         end
-        
-        
+
+
         # Implements CoordinateSystem#get_axis
-        
+
         def get_axis(index_)
           [@axis0, @axis1, @axis2][index_]
         end
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "GEOCCS"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           arr_ = [@horizontal_datum._to_wkt(open_, close_), @prime_meridian._to_wkt(open_, close_), @linear_unit._to_wkt(open_, close_)]
           arr_ << @axis0._to_wkt(open_, close_) if @axis0
@@ -1273,130 +1272,130 @@ module RGeo
           arr_ << @axis2._to_wkt(open_, close_) if @axis2
           arr_
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a GeocentricCoordinateSystem given a name, a
           # HorizontalDatum, a PrimeMeridian, a LinearUnit, and three
           # AxisInfo objects. The AxisInfo are optional and may be nil.
           # You may also provide the optional parameters specified by the
           # Info interface.
-          
+
           def create(name_, horizontal_datum_, prime_meridian_, linear_unit_, axis0_, axis1_, axis2_, *optional_)
             new(name_, horizontal_datum_, prime_meridian_, linear_unit_, axis0_, axis1_, axis2_, *optional_)
           end
-          
-          
+
+
         end
-        
-        
+
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # A one-dimensional coordinate system suitable for vertical
       # measurements.
-      
+
       class VerticalCoordinateSystem < CoordinateSystem
-        
+
         def initialize(name_, vertical_datum_, vertical_unit_, axis_, *optional_)  # :nodoc:
           super(name_, 1, *optional_)
           @vertical_datum = vertical_datum_
           @vertical_unit = vertical_unit_
           @axis = axis_
         end
-        
-        
+
+
         # Gets the vertical datum, which indicates the measurement method.
         attr_reader :vertical_datum
-        
+
         # Gets the units used along the vertical axis. The vertical units
         # must be the same as the CS_CoordinateSystem units.
         attr_reader :vertical_unit
-        
-        
+
+
         # Implements CoordinateSystem#get_units
-        
+
         def get_units(index_)
           @vertical_unit
         end
-        
-        
+
+
         # Implements CoordinateSystem#get_axis
-        
+
         def get_axis(index_)
           @axis
         end
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "VERT_CS"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           arr_ = [@vertical_datum._to_wkt(open_, close_), @vertical_unit._to_wkt(open_, close_)]
           arr_ << @axis._to_wkt(open_, close_) if @axis
           arr_
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a VerticalCoordinateSystem given a name, a
           # VerticalDatum, a LinearUnit, and an AxisInfo. The AxisInfo is
           # optional and may be nil. You may also provide the optional
           # parameters specified by the Info interface.
-          
+
           def create(name_, vertical_datum_, vertical_unit_, axis_, *optional_)
             new(name_, vertical_datum_, vertical_unit_, axis_, *optional_)
           end
-          
-          
+
+
         end
-        
-        
+
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # A 2D coordinate system suitable for positions on the Earth's surface.
-      # 
+      #
       # == Notes
-      # 
+      #
       # This is a non-instantiable abstract class. You must instantiate
       # one of the subclasses GeographicCoordinateSystem or
       # ProjectedCoordinateSystem.
-      
+
       class HorizontalCoordinateSystem < CoordinateSystem
-        
+
         def initialize(name_, horizontal_datum_, *optional_)  # :nodoc:
           super(name_, 2, *optional_)
           @horizontal_datum = horizontal_datum_
         end
-        
-        
+
+
         # Returns the HorizontalDatum.
         attr_reader :horizontal_datum
-        
-        
+
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # A coordinate system based on latitude and longitude. Some
       # geographic coordinate systems are Lat/Lon, and some are Lon/Lat.
       # You can find out which this is by examining the axes. You should
       # also check the angular units, since not all geographic coordinate
       # systems use degrees.
-      
+
       class GeographicCoordinateSystem < HorizontalCoordinateSystem
-        
+
         def initialize(name_, angular_unit_, horizontal_datum_, prime_meridian_, axis0_, axis1_, *optional_)  # :nodoc:
           super(name_, horizontal_datum_, *optional_)
           @prime_meridian = prime_meridian_
@@ -1404,86 +1403,86 @@ module RGeo
           @axis0 = axis0_
           @axis1 = axis1_
         end
-        
-        
+
+
         # Returns the PrimeMeridian.
         attr_reader :prime_meridian
-        
+
         # Returns the AngularUnit. The angular unit must be the same as
         # the CS_CoordinateSystem units.
         attr_reader :angular_unit
-        
-        
+
+
         # Implements CoordinateSystem#get_units
-        
+
         def get_units(index_)
           @angular_unit
         end
-        
+
 
         # Implements CoordinateSystem#get_axis
-        
+
         def get_axis(index_)
           index_ == 1 ? @axis1 : @axis0
         end
-        
-        
+
+
         # Gets the number of available conversions to WGS84 coordinates.
-        
+
         def num_conversion_to_wgs84
           @horizontal_datum.wgs84_parameters ? 1 : 0
         end
-        
-        
+
+
         # Gets details on a conversion to WGS84. Some geographic
         # coordinate systems provide several transformations into WGS84,
         # which are designed to provide good accuracy in different areas
         # of interest. The first conversion (with index=0) should provide
         # acceptable accuracy over the largest possible area of interest.
-        
+
         def get_wgs84_conversion_info(index_)
           @horizontal_datum.wgs84_parameters
         end
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "GEOGCS"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           arr_ = [@horizontal_datum._to_wkt(open_, close_), @prime_meridian._to_wkt(open_, close_), @angular_unit._to_wkt(open_, close_)]
           arr_ << @axis0._to_wkt(open_, close_) if @axis0
           arr_ << @axis1._to_wkt(open_, close_) if @axis1
           arr_
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a GeographicCoordinateSystem, given a name, an
           # AngularUnit, a HorizontalDatum, a PrimeMeridian, and two
           # AxisInfo objects. The AxisInfo objects are optional and may
           # be set to nil. You may also provide the optional parameters
           # specified by the Info interface.
-          
+
           def create(name_, angular_unit_, horizontal_datum_, prime_meridian_, axis0_, axis1_, *optional_)
             new(name_, angular_unit_, horizontal_datum_, prime_meridian_, axis0_, axis1_, *optional_)
           end
-          
-          
+
+
         end
-        
-        
+
+
       end
-      
-      
+
+
       # == OGC spec description
-      # 
+      #
       # A 2D cartographic coordinate system.
-      
+
       class ProjectedCoordinateSystem < HorizontalCoordinateSystem
-        
+
         def initialize(name_, geographic_coordinate_system_, projection_, linear_unit_, axis0_, axis1_, *optional_)  # :nodoc:
           super(name_, geographic_coordinate_system_.horizontal_datum, *optional_)
           @geographic_coordinate_system = geographic_coordinate_system_
@@ -1492,37 +1491,37 @@ module RGeo
           @axis0 = axis0_
           @axis1 = axis1_
         end
-        
-        
+
+
         # Returns the GeographicCoordinateSystem.
         attr_reader :geographic_coordinate_system
-        
+
         # Gets the projection.
         attr_reader :projection
-        
+
         # Returns the LinearUnits. The linear unit must be the same as
         # the CS_CoordinateSystem units.
         attr_reader :linear_unit
-        
-        
+
+
         # Implements CoordinateSystem#get_units
-        
+
         def get_units(index_)
           @linear_unit
         end
-        
-        
+
+
         # Implements CoordinateSystem#get_axis
-        
+
         def get_axis(index_)
           index_ == 1 ? @axis1 : @axis0
         end
-        
-        
+
+
         def _wkt_typename  # :nodoc:
           "PROJCS"
         end
-        
+
         def _wkt_content(open_, close_)  # :nodoc:
           arr_ = [@geographic_coordinate_system._to_wkt(open_, close_), @projection._to_wkt(open_, close_)]
           @projection.each_parameter{ |param_| arr_ << param_._to_wkt(open_, close_) }
@@ -1531,31 +1530,31 @@ module RGeo
           arr_ << @axis1._to_wkt(open_, close_) if @axis1
           arr_
         end
-        
-        
+
+
         class << self
-          
-          
+
+
           # Create a ProjectedCoordinateSystem given a name, a
           # GeographicCoordinateSystem, and Projection, a LinearUnit, and
           # two AxisInfo objects. The AxisInfo objects are optional and
           # may be set to nil. You may also provide the optional
           # parameters specified by the Info interface.
-          
+
           def create(name_, geographic_coordinate_system_, projection_, linear_unit_, axis0_, axis1_, *optional_)
             new(name_, geographic_coordinate_system_, projection_, linear_unit_, axis0_, axis1_, *optional_)
           end
-          
-          
+
+
         end
-        
-        
+
+
       end
-      
-      
+
+
     end
-    
-    
+
+
   end
-  
+
 end