adiwg-mdtranslator 2.3.5 → 2.4.0

data/lib/adiwg/mdtranslator/readers/fgdc/modules/mapProjections/projection_spaceOblique.rb

@@ -0,0 +1,53 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc map projection - space oblique
+
+# History:
+#  Stan Smith 2017-10-18 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+require_relative 'projection_common'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module SpaceObliqueProjection
+
+               def self.unpack(xParams, hProjection, hResponseObj)
+
+                  # map projection 4.1.2.1.19 (spaceobq) - Space Oblique Mercator (Landsat)
+                  unless xParams.empty?
+                     paramCount = 0
+                     hProjection[:projectionName] = 'space oblique'
+
+                     # -> ReferenceSystemParameters.projection.landsatNumber
+                     paramCount += ProjectionCommon.unpackLandSat(xParams, hProjection, hResponseObj)
+
+                     # -> ReferenceSystemParameters.projection.landsatPath
+                     paramCount += ProjectionCommon.unpackLandSatPath(xParams, hProjection, hResponseObj)
+
+                     # -> ReferenceSystemParameters.projection.falseEasting
+                     # -> ReferenceSystemParameters.projection.falseNorthing
+                     paramCount += ProjectionCommon.unpackFalseNE(xParams, hProjection, hResponseObj)
+
+                     # verify parameter count
+                     if paramCount == 4
+                        return hProjection
+                     else
+                        hResponseObj[:readerExecutionMessages] << 'space oblique projection is missing one or more parameters'
+                        return nil
+                     end
+                  end
+
+                  return nil
+
+               end
+
+            end
+
+         end
+      end
+   end
+end

data/lib/adiwg/mdtranslator/readers/fgdc/modules/mapProjections/projection_stereographic.rb

@@ -0,0 +1,53 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc map projection - stereographic
+
+# History:
+#  Stan Smith 2017-10-18 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+require_relative 'projection_common'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module StereographicProjection
+
+               def self.unpack(xParams, hProjection, hResponseObj)
+
+                  # map projection 4.1.2.1.20 (stereo) - Stereographic
+                  unless xParams.empty?
+                     paramCount = 0
+                     hProjection[:projectionName] = 'stereographic'
+
+                     # -> ReferenceSystemParameters.projection.longitudeOfProjectionCenter
+                     paramCount += ProjectionCommon.unpackLongPC(xParams, hProjection, hResponseObj)
+
+                     # -> ReferenceSystemParameters.projection.latitudeOfProjectionCenter
+                     paramCount += ProjectionCommon.unpackLatPC(xParams, hProjection, hResponseObj)
+
+                     # -> ReferenceSystemParameters.projection.falseEasting
+                     # -> ReferenceSystemParameters.projection.falseNorthing
+                     paramCount += ProjectionCommon.unpackFalseNE(xParams, hProjection, hResponseObj)
+
+                     # verify parameter count
+                     if paramCount == 4
+                        return hProjection
+                     else
+                        hResponseObj[:readerExecutionMessages] << 'stereographic projection is missing one or more parameters'
+                        return nil
+                     end
+                  end
+
+                  return nil
+
+               end
+
+            end
+
+         end
+      end
+   end
+end

data/lib/adiwg/mdtranslator/readers/fgdc/modules/mapProjections/projection_transverseMercator.rb

@@ -0,0 +1,54 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc map projection - transverse mercator
+
+# History:
+#  Stan Smith 2017-10-04 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module TransverseMercatorProjection
+
+               def self.unpack(xParams, hProjection, hResponseObj)
+
+                  # map projection 4.1.2.1.21 (transmer) - Transverse Mercator
+                  unless xParams.empty?
+                     paramCount = 0
+                     hProjection[:projectionName] = 'transverse mercator'
+
+                     # -> ReferenceSystemParameters.projection.scaleFactorAtCentralMeridian
+                     paramCount += ProjectionCommon.unpackSFCM(xParams, hProjection, hResponseObj)
+
+                     # -> ReferenceSystemParameters.projection.longitudeOfCentralMeridian
+                     paramCount += ProjectionCommon.unpackLongCM(xParams, hProjection, hResponseObj)
+
+                     # -> ReferenceSystemParameters.projection.latitudeOfProjectionOrigin
+                     paramCount += ProjectionCommon.unpackLatPO(xParams, hProjection, hResponseObj)
+
+                     # -> ReferenceSystemParameters.projection.falseEasting
+                     # -> ReferenceSystemParameters.projection.falseNorthing
+                     paramCount += ProjectionCommon.unpackFalseNE(xParams, hProjection, hResponseObj)
+
+                     if paramCount == 5
+                        return hProjection
+                     else
+                        hResponseObj[:readerExecutionMessages] << 'transverse mercator projection is missing one or more parameters'
+                        return nil
+                     end
+                  end
+
+                  return nil
+
+               end
+
+            end
+
+         end
+      end
+   end
+end

data/lib/adiwg/mdtranslator/readers/fgdc/modules/mapProjections/projection_vanDerGrinten.rb

@@ -0,0 +1,48 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc map projection - van der grinten
+
+# History:
+#  Stan Smith 2017-10-18 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module VanDerGrintenProjection
+
+               def self.unpack(xParams, hProjection, hResponseObj)
+
+                  # map projection 4.1.2.1.22 (vdgrin) - van der Grinten
+                  unless xParams.empty?
+                     paramCount = 0
+                     hProjection[:projectionName] = 'van der grinten'
+
+                     # -> ReferenceSystemParameters.projection.longitudeOfCentralMeridian
+                     paramCount += ProjectionCommon.unpackLongCM(xParams, hProjection, hResponseObj)
+
+                     # -> ReferenceSystemParameters.projection.falseEasting
+                     # -> ReferenceSystemParameters.projection.falseNorthing
+                     paramCount += ProjectionCommon.unpackFalseNE(xParams, hProjection, hResponseObj)
+
+                     if paramCount == 3
+                        return hProjection
+                     else
+                        hResponseObj[:readerExecutionMessages] << 'van der grinten projection is missing one or more parameters'
+                        return nil
+                     end
+                  end
+
+                  return nil
+
+               end
+
+            end
+
+         end
+      end
+   end
+end

data/lib/adiwg/mdtranslator/readers/fgdc/modules/module_fgdc.rb

@@ -69,10 +69,10 @@ module ADIWG
                   end
 
                   # metadata (spref 4) - spatial reference
-                  # xSpatialRef = xMetadata.xpath('./spref')
-                  # unless xSpatialRef.empty?
-                  #    SpatialReference.unpack(xSpatialRef, hResponseObj)
-                  # end
+                  xSpatialRef = xMetadata.xpath('./spref')
+                  unless xSpatialRef.empty?
+                     SpatialReference.unpack(xSpatialRef, hResourceInfo, hResponseObj)
+                  end
 
                   # metadata (eainfo 5) - entity and attribute
                   xEntity = xMetadata.xpath('./eainfo')

data/lib/adiwg/mdtranslator/readers/fgdc/modules/module_horizontalPlanar.rb

@@ -0,0 +1,81 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc horizontal planar data reference
+
+# History:
+#  Stan Smith 2017-10-03 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+require_relative 'module_mapProjection'
+require_relative 'module_mapGridSystem'
+require_relative 'module_mapLocalSystem'
+require_relative 'module_mapCoordinateInfo'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module PlanarReference
+
+               def self.unpack(xPlanar, hResourceInfo, hResponseObj)
+
+                  # instance classes needed in script
+                  intMetadataClass = InternalMetadata.new
+                  hProjection = nil
+
+                  # get distance unit of measure for use in map projections and grid systems
+                  distanceUnits = nil
+                  xPlanCI = xPlanar.xpath('./planci')
+                  unless xPlanCI.empty?
+                     planDU = xPlanCI.xpath('./plandu').text
+                     unless planDU.empty?
+                        distanceUnits = planDU
+                     end
+                  end
+
+                  # planar reference 4.1.2.1 (mapproj) - map projection
+                  xMapProjection = xPlanar.xpath('./mapproj')
+                  unless xMapProjection.empty?
+                     hProjection = MapProjection.unpack(xMapProjection, hResponseObj)
+                  end
+
+                  # planar reference 4.1.2.2 (gridsys) - grid coordinate system
+                  xMapGrid = xPlanar.xpath('./gridsys')
+                  unless xMapGrid.empty?
+                     hProjection = MapGridSystem.unpack(xMapGrid, hResponseObj)
+                  end
+
+                  # planar reference 4.1.2.3 (localp) - local planar
+                  xMapLocal = xPlanar.xpath('./localp')
+                  unless xMapLocal.empty?
+                     hProjection = MapLocalSystem.unpack(xMapLocal, hResponseObj)
+                  end
+
+                  unless hProjection.nil?
+                     hReferenceSystem = intMetadataClass.newSpatialReferenceSystem
+                     hSystemParameters = intMetadataClass.newReferenceSystemParameterSet
+                     unless distanceUnits.nil?
+                        hProjection[:falseEastingNorthingUnits] = distanceUnits
+                     end
+                     hSystemParameters[:projection] = hProjection
+                     hReferenceSystem[:systemParameterSet] = hSystemParameters
+                     hResourceInfo[:spatialReferenceSystems] << hReferenceSystem
+                  end
+
+                  # planar reference 4.1.2.4 (planci) - planar coordinate information
+                  xPlanarCI = xPlanar.xpath('./planci')
+                  unless xPlanarCI.empty?
+                     PlanarCoordinateInformation.unpack(xPlanarCI, hResourceInfo, hResponseObj)
+                  end
+
+                  return hResourceInfo
+
+               end
+
+            end
+
+         end
+      end
+   end
+end

data/lib/adiwg/mdtranslator/readers/fgdc/modules/module_horizontalReference.rb

@@ -0,0 +1,167 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc horizontal data reference
+
+# History:
+#  Stan Smith 2017-10-02 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+require_relative 'module_horizontalPlanar'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module HorizontalReference
+
+               def self.unpack(xHorizontalRef, hResourceInfo, hResponseObj)
+
+                  # instance classes needed in script
+                  intMetadataClass = InternalMetadata.new
+
+                  # horizontal reference 4.1.1 (geograph) - geographic coordinate system
+                  xGeographic = xHorizontalRef.xpath('./geograph')
+                  unless xGeographic.empty?
+
+                     hResolution = intMetadataClass.newSpatialResolution
+                     hGeoResolution = intMetadataClass.newGeographicResolution
+                     hasValue = false
+
+                     # geographic reference 4.1.1.1 (latres) - latitude resolution
+                     # -> spatialResolution.geographicResolution.latitudeResolution
+                     latResolution = xGeographic.xpath('./latres').text
+                     unless latResolution.empty?
+                        hGeoResolution[:latitudeResolution] = latResolution.to_f
+                        hasValue = true
+                     end
+
+                     # geographic reference 4.1.1.2 (longres) - longitude resolution
+                     # -> spatialResolution.geographicResolution.longitudeResolution
+                     longResolution = xGeographic.xpath('./longres').text
+                     unless longResolution.empty?
+                        hGeoResolution[:longitudeResolution] = longResolution.to_f
+                        hasValue = true
+                     end
+
+                     # geographic reference 4.1.1.3 (geogunit) - latitude/longitude units
+                     # -> spatialResolution.geographicResolution.unitOfMeasure
+                     unitMeasure = xGeographic.xpath('./geogunit').text
+                     unless unitMeasure.empty?
+                        hGeoResolution[:unitOfMeasure] = unitMeasure
+                        hasValue = true
+                     end
+
+                     if hasValue
+                        hResolution[:geographicResolution] = hGeoResolution
+                        hResourceInfo[:spatialResolutions] << hResolution
+                     end
+
+                  end
+
+                  # horizontal reference 4.1.2 (planar) - planar coordinate system []
+                  axPlanar = xHorizontalRef.xpath('./planar')
+                  unless axPlanar.empty?
+                     axPlanar.each do |xPlanar|
+                        PlanarReference.unpack(xPlanar, hResourceInfo, hResponseObj)
+                     end
+                  end
+
+                  # horizontal reference 4.1.3 (local) - local coordinate system
+                  xLocal = xHorizontalRef.xpath('./local')
+                  unless xLocal.empty?
+
+                     hProjection = intMetadataClass.newProjection
+                     hProjection[:projectionName] = 'local'
+
+                     # local planar 4.1.3.1 (localdes) - local description
+                     # -> referenceSystemParameters.projection.localPlanarDescription
+                     description = xLocal.xpath('./localdes').text
+                     unless description.empty?
+                        hProjection[:localPlanarDescription] = description
+                     end
+
+                     # local planar 4.1.3.2 (localgeo) - local georeference information
+                     # -> referenceSystemParameters.projection.localPlanarGeoreference
+                     georeference = xLocal.xpath('./localgeo').text
+                     unless georeference.empty?
+                        hProjection[:localPlanarGeoreference] = georeference
+                     end
+
+                     hReferenceSystem = intMetadataClass.newSpatialReferenceSystem
+                     hSystemParameters = intMetadataClass.newReferenceSystemParameterSet
+                     hSystemParameters[:projection] = hProjection
+                     hReferenceSystem[:systemParameterSet] = hSystemParameters
+                     hResourceInfo[:spatialReferenceSystems] << hReferenceSystem
+
+                  end
+
+                  # horizontal reference 4.1.4 (geodetic) - parameters for shape of earth
+                  xGeodetic = xHorizontalRef.xpath('./geodetic')
+                  unless xGeodetic.empty?
+
+                     hEllipsoid = intMetadataClass.newEllipsoid
+
+                     # geodetic model 4.1.4.1 (horizdn) - horizontal datum name
+                     # -> referenceSystemParameters.datumIdentifier.identifier.identifier
+                     datumName = xGeodetic.xpath('./horizdn').text
+                     unless datumName.empty?
+                        hIdentifier = intMetadataClass.newIdentifier
+                        hIdentifier[:identifier] = datumName
+                        hEllipsoid[:ellipsoidIdentifier] = hIdentifier
+                     end
+
+                     # geodetic model 4.1.4.2 (ellips) - ellipsoid name
+                     # -> referenceSystemParameters.ellipsoid.ellipsoidName
+                     ellipsoidName = xGeodetic.xpath('./ellips').text
+                     unless ellipsoidName.empty?
+                        hEllipsoid[:ellipsoidName] = ellipsoidName
+                     end
+
+                     # geodetic model 4.1.4.3 (semiaxis) - semi-major axis
+                     # -> referenceSystemParameters.ellipsoid.semiMajorAxis
+                     semiAxis = xGeodetic.xpath('./semiaxis').text
+                     unless semiAxis.empty?
+                        hEllipsoid[:semiMajorAxis] = semiAxis.to_f
+                     end
+
+                     # geodetic model 4.1.2.4.4 (plandu) - distance units
+                     # take value from the first 'planar' section with 'plandu' specified
+                     axPlanar = xHorizontalRef.xpath('./planar')
+                     unless axPlanar.empty?
+                        axPlanar.each do |xPlanar|
+                           xPlanCI = xPlanar.xpath('./planci')
+                           unless xPlanCI.empty?
+                              units = xPlanCI.xpath('./plandu').text
+                              unless units.empty?
+                                 hEllipsoid[:axisUnits] = units
+                              end
+                           end
+                        end
+                     end
+
+                     # geodetic model 4.1.4.4 (denflat) - denominator of flattening ratio
+                     # -> referenceSystemParameters.ellipsoid.denominatorOfFlatteningRatio
+                     flattening = xGeodetic.xpath('./denflat').text
+                     unless flattening.empty?
+                        hEllipsoid[:denominatorOfFlatteningRatio] = flattening.to_f
+                     end
+
+                     hReferenceSystem = intMetadataClass.newSpatialReferenceSystem
+                     hSystemParameters = intMetadataClass.newReferenceSystemParameterSet
+                     hSystemParameters[:ellipsoid] = hEllipsoid
+                     hReferenceSystem[:systemParameterSet] = hSystemParameters
+                     hResourceInfo[:spatialReferenceSystems] << hReferenceSystem
+
+                  end
+
+                  return hResourceInfo
+
+               end
+
+            end
+
+         end
+      end
+   end
+end