adiwg-mdtranslator 2.3.5 → 2.4.0

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

@@ -0,0 +1,128 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc planar coordinate information
+
+# History:
+#  Stan Smith 2017-10-05 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module PlanarCoordinateInformation
+
+               def self.unpack(xPlanarCI, hResourceInfo, hResponseObj)
+
+                  # instance classes needed in script
+                  intMetadataClass = InternalMetadata.new
+
+                  # map projection 4.1.2.4.1 (plance) - planar coordinate encoding method
+                  # -> resourceInfo.spatialRepresentationTypes
+                  encoding = xPlanarCI.xpath('./plance').text
+                  unless encoding.empty?
+                     hResourceInfo[:spatialRepresentationTypes] << encoding
+                  end
+
+                  # map projection 4.1.2.4.2 (coordrep) - coordinate representation
+                  xCoordRep = xPlanarCI.xpath('./coordrep')
+                  unless xCoordRep.empty?
+
+                     hCoordResolution = intMetadataClass.newCoordinateResolution
+
+                     # map projection 4.1.2.4.2.1 (absres) - abscissa resolution
+                     # -> resourceInfo.spatialResolutions.spatialResolution.coordinateResolution.abscissaResolutionX
+                     abscissa = xCoordRep.xpath('./absres').text
+                     unless abscissa.empty?
+                        hCoordResolution[:abscissaResolutionX] = abscissa.to_f
+                     end
+
+                     # map projection 4.1.2.4.2.2 (ordres) - ordinate resolution
+                     # -> resourceInfo.spatialResolutions.spatialResolution.coordinateResolution.ordinateResolutionY
+                     ordinate = xCoordRep.xpath('./ordres').text
+                     unless ordinate.empty?
+                        hCoordResolution[:ordinateResolutionY] = ordinate.to_f
+                     end
+
+                     # map projection 4.1.2.4.4 (plandu) - planar distance units
+                     # -> resourceInfo.spatialResolutions.spatialResolution.coordinateResolution.unitOfMeasure
+                     distUnits = xPlanarCI.xpath('./plandu').text
+                     unless distUnits.empty?
+                        hCoordResolution[:unitOfMeasure] = distUnits
+                     end
+
+                     hResolution = intMetadataClass.newSpatialResolution
+                     hResolution[:coordinateResolution] = hCoordResolution
+                     hResourceInfo[:spatialResolutions] << hResolution
+
+                  end
+
+                  # map projection 4.1.2.4.3 (distbrep) - distance and bearing representation
+                  xBDRep = xPlanarCI.xpath('./distbrep')
+                  unless xBDRep.empty?
+
+                     hBDResolution = intMetadataClass.newBearingDistanceResolution
+
+                     # map projection 4.1.2.4.3.1 (distres) - distance resolution
+                     # -> resourceInfo.spatialResolutions.spatialResolution.bearingDistanceResolution.distanceResolution
+                     distRes = xBDRep.xpath('./distres').text
+                     unless distRes.empty?
+                        hBDResolution[:distanceResolution] = distRes.to_f
+                     end
+
+                     # map projection 4.1.2.4.4 (plandu) - planar distance units
+                     # -> resourceInfo.spatialResolutions.spatialResolution.bearingDistanceResolution.distanceUnitOfMeasure
+                     distUnits = xPlanarCI.xpath('./plandu').text
+                     unless distUnits.empty?
+                        hBDResolution[:distanceUnitOfMeasure] = distUnits
+                     end
+                     # map projection 4.1.2.4.3.2 (bearres) - bearing resolution
+                     # -> resourceInfo.spatialResolutions.spatialResolution.bearingDistanceResolution.bearingResolution
+                     bearingRes = xBDRep.xpath('./bearres').text
+                     unless bearingRes.empty?
+                        hBDResolution[:bearingResolution] = bearingRes.to_f
+                     end
+
+                     # map projection 4.1.2.4.3.3 (bearunit) - bearing units
+                     # -> resourceInfo.spatialResolutions.spatialResolution.bearingDistanceResolution.bearingUnitsOfMeasure
+                     bearingUnits = xBDRep.xpath('./bearunit').text
+                     unless bearingUnits.empty?
+                        hBDResolution[:bearingUnitsOfMeasure] = bearingUnits
+                     end
+
+                     # map projection 4.1.2.4.3.4 (bearrefd) - bearing reference direction
+                     # -> resourceInfo.spatialResolutions.spatialResolution.bearingDistanceResolution.bearingReferenceDirection
+                     bearingDirection = xBDRep.xpath('./bearrefd').text
+                     unless bearingDirection.empty?
+                        hBDResolution[:bearingReferenceDirection] = bearingDirection
+                     end
+
+                     # map projection 4.1.2.4.3.5 (bearrefm) - bearing reference meridian
+                     # -> resourceInfo.spatialResolutions.spatialResolution.bearingDistanceResolution.bearingReferenceMeridian
+                     bearingMeridian = xBDRep.xpath('./bearrefm').text
+                     unless bearingMeridian.empty?
+                        hBDResolution[:bearingReferenceMeridian] = bearingMeridian
+                     end
+
+                     hResolution = intMetadataClass.newSpatialResolution
+                     hResolution[:bearingDistanceResolution] = hBDResolution
+                     hResourceInfo[:spatialResolutions] << hResolution
+
+                  end
+
+                  # map projection 4.1.2.4.4 (plandu) - planar distance units
+                  # -> resourceInfo.spatialResolutions.spatialResolution.bearingDistanceResolution.distanceUnitOfMeasure
+                  # -> resourceInfo.spatialResolutions.spatialResolution.coordinateResolution.unitOfMeasure
+
+                  return nil
+
+               end
+
+            end
+
+         end
+      end
+   end
+end

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

@@ -0,0 +1,159 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc map grid coordinate system
+
+# History:
+#  Stan Smith 2017-10-04 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+require_relative 'mapProjections/projection_transverseMercator'
+require_relative 'mapProjections/projection_polarStereo'
+require_relative 'mapProjections/projection_lambertConic'
+require_relative 'mapProjections/projection_obliqueMercator'
+require_relative 'mapProjections/projection_polyconic'
+require_relative 'mapProjections/projection_equirectangular'
+require_relative 'mapProjections/projection_azimuthEquidistant'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module MapGridSystem
+
+               def self.unpack(xMapGrid, hResponseObj)
+
+                  # instance classes needed in script
+                  intMetadataClass = InternalMetadata.new
+                  hProjection = intMetadataClass.newProjection
+
+                  # grid system 4.1.2.2.1 (gridsysn) - grid coordinate system name
+                  # -> ReferenceSystemParameters.projection.projectionIdentifier.identifier
+                  name = xMapGrid.xpath('./gridsysn').text
+                  unless name.empty?
+                     hIdentifier = intMetadataClass.newIdentifier
+                     hIdentifier[:identifier] = name
+                     hProjection[:projectionIdentifier] = hIdentifier
+                  end
+
+                  # grid system 4.1.2.2.2 (utm) - universal transverse mercator
+                  xUTM = xMapGrid.xpath('./utm')
+                  unless xUTM.empty?
+
+                     # grid system 4.1.2.2.2.1 (utmzone) - utm zone number {-60..-1, 1..60}
+                     # -> ReferenceSystemParameters.projection.zone
+                     zone = xUTM.xpath('./utmzone').text
+                     unless zone.empty?
+                        hProjection[:zone] = zone
+                     end
+
+                     # + transverse mercator
+                     xTransMer = xUTM.xpath('./transmer')
+                     unless xTransMer.empty?
+                        return TransverseMercatorProjection.unpack(xTransMer, hProjection, hResponseObj)
+                     end
+
+                  end
+
+                  # grid system 4.1.2.2.3 (ups) - universal polar stereographic
+                  xUSP = xMapGrid.xpath('./ups')
+                  unless xUSP.empty?
+
+                     # grid system 4.1.2.2.3.1 (upszone) - utm zone number {-60..-1, 1..60}
+                     # -> ReferenceSystemParameters.projection.zone
+                     zone = xUSP.xpath('./upszone').text
+                     unless zone.empty?
+                        hProjection[:zone] = zone
+                     end
+
+                     # + polar stereographic
+                     xPolarS = xUSP.xpath('./polarst')
+                     unless xPolarS.empty?
+                        return PolarStereoProjection.unpack(xPolarS, hProjection, hResponseObj)
+                     end
+
+                  end
+
+                  # grid system 4.1.2.2.4 (spcs) - state plane coordinate system
+                  xStateP = xMapGrid.xpath('./spcs')
+                  unless xStateP.empty?
+
+                     # grid system 4.1.2.2.4.1 (spcszone) - state plane zone number {nnnn}
+                     # -> ReferenceSystemParameters.projection.zone
+                     zone = xStateP.xpath('./spcszone').text
+                     unless zone.empty?
+                        hProjection[:zone] = zone
+                     end
+
+                     # + [ lambert conformal conic | transverse mercator | oblique mercator | polyconic ]
+                     # + lambert conformal conic
+                     xLambert = xStateP.xpath('./lambertc')
+                     unless xLambert.empty?
+                        return LambertConicProjection.unpack(xLambert, hProjection, hResponseObj)
+                     end
+
+                     # + transverse mercator
+                     xTransMer = xStateP.xpath('./transmer')
+                     unless xTransMer.empty?
+                        return TransverseMercatorProjection.unpack(xTransMer, hProjection, hResponseObj)
+                     end
+
+                     # + oblique mercator
+                     xObliqueM = xStateP.xpath('./obqmerc')
+                     unless xObliqueM.empty?
+                        return ObliqueMercatorProjection.unpack(xObliqueM, hProjection, hResponseObj)
+                     end
+
+                     # + polyconic
+                     xPolyCon = xStateP.xpath('./polycon')
+                     unless xPolyCon.empty?
+                        return PolyconicProjection.unpack(xPolyCon, hProjection, hResponseObj)
+                     end
+
+                  end
+
+                  # grid system 4.1.2.2.5 (arcsys) - equal arc-second coordinate system
+                  xArc = xMapGrid.xpath('./arcsys')
+                  unless xArc.empty?
+
+                     # grid system 4.1.2.2.5.1 (arcszone) - state plane zone number {1..18}
+                     # -> ReferenceSystemParameters.projection.zone
+                     zone = xArc.xpath('./arczone').text
+                     unless zone.empty?
+                        hProjection[:zone] = zone
+                     end
+
+                     # + [ equirectangular | azimuthal equidistant ]
+                     # + equirectangular
+                     xEquiR = xArc.xpath('./equirect')
+                     unless xEquiR.empty?
+                        return EquirectangularProjection.unpack(xEquiR, hProjection, hResponseObj)
+                     end
+
+                     # + azimuthal equidistant
+                     xAzimuthE = xArc.xpath('./azimequi')
+                     unless xAzimuthE.empty?
+                        return AzimuthEquidistantProjection.unpack(xAzimuthE, hProjection, hResponseObj)
+                     end
+
+                  end
+
+                  # grid system 4.1.2.2.6 (othergrd) - other coordinate system {text}
+                  # -> ReferenceSystemParameters.projection.otherGridDescription
+                  otherG = xMapGrid.xpath('./othergrd').text
+                  unless otherG.empty?
+                     hProjection[:otherGridDescription] = otherG
+                     hProjection[:projectionName] = 'other grid system'
+                     return hProjection
+                  end
+
+                  return nil
+
+               end
+
+            end
+
+         end
+      end
+   end
+end

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

@@ -0,0 +1,59 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc map grid coordinate system
+
+# History:
+#  Stan Smith 2017-10-04 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+require_relative 'mapProjections/projection_transverseMercator'
+require_relative 'mapProjections/projection_polarStereo'
+require_relative 'mapProjections/projection_lambertConic'
+require_relative 'mapProjections/projection_obliqueMercator'
+require_relative 'mapProjections/projection_polyconic'
+require_relative 'mapProjections/projection_equirectangular'
+require_relative 'mapProjections/projection_azimuthEquidistant'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module MapLocalSystem
+
+               def self.unpack(xMapLocal, hResponseObj)
+
+                  # instance classes needed in script
+                  intMetadataClass = InternalMetadata.new
+                  hProjection = intMetadataClass.newProjection
+
+                  # -> ReferenceSystemParameters.projection.projectionIdentifier.identifier
+                  hIdentifier = intMetadataClass.newIdentifier
+                  hIdentifier[:identifier] = 'local planar coordinate system'
+                  hProjection[:projectionIdentifier] = hIdentifier
+                  hProjection[:projectionName] = 'local planar'
+
+                  # local planar 4.1.2.3.1 (localpd) - local planar description
+                  # -> ReferenceSystemParameters.projection.localPlanarDescription
+                  description = xMapLocal.xpath('./localpd').text
+                  unless description.empty?
+                     hProjection[:localPlanarDescription] = description
+                  end
+
+                  # local planar 4.1.2.3.2 (localpgi) - local planar georeference information
+                  # -> ReferenceSystemParameters.projection.localPlanarGeoreference
+                  georeference = xMapLocal.xpath('./localpgi').text
+                  unless georeference.empty?
+                     hProjection[:localPlanarGeoreference] = georeference
+                  end
+
+                  return hProjection
+
+               end
+
+            end
+
+         end
+      end
+   end
+end

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

@@ -0,0 +1,195 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc map projection
+
+# History:
+#  Stan Smith 2017-10-03 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+require_relative 'mapProjections/projection_albers'
+require_relative 'mapProjections/projection_azimuthEquidistant'
+require_relative 'mapProjections/projection_equidistantConic'
+require_relative 'mapProjections/projection_equirectangular'
+require_relative 'mapProjections/projection_generalVertical'
+require_relative 'mapProjections/projection_gnomonic'
+require_relative 'mapProjections/projection_lambertEqualArea'
+require_relative 'mapProjections/projection_lambertConic'
+require_relative 'mapProjections/projection_mercator'
+require_relative 'mapProjections/projection_modifiedAlaska'
+require_relative 'mapProjections/projection_millerCylinder'
+require_relative 'mapProjections/projection_obliqueMercator'
+require_relative 'mapProjections/projection_orthographic'
+require_relative 'mapProjections/projection_polarStereo'
+require_relative 'mapProjections/projection_polyconic'
+require_relative 'mapProjections/projection_robinson'
+require_relative 'mapProjections/projection_sinusoidal'
+require_relative 'mapProjections/projection_spaceOblique'
+require_relative 'mapProjections/projection_stereographic'
+require_relative 'mapProjections/projection_transverseMercator'
+require_relative 'mapProjections/projection_vanDerGrinten'
+require_relative 'mapProjections/projection_parameters'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module MapProjection
+
+               def self.unpack(xMapProjection, hResponseObj)
+
+                  # instance classes needed in script
+                  intMetadataClass = InternalMetadata.new
+                  hProjection = intMetadataClass.newProjection
+
+                  # map projection 4.1.2.1.1 (mapprojn) - map projection name
+                  # -> ReferenceSystemParameters.projection.projectionIdentifier.identifier
+                  name = xMapProjection.xpath('./mapprojn').text
+                  unless name.empty?
+                     hIdentifier = intMetadataClass.newIdentifier
+                     hIdentifier[:identifier] = name
+                     hProjection[:projectionIdentifier] = hIdentifier
+                  end
+
+                  # map projection 4.1.2.1.2 (albers) - Albers Conical Equal Area
+                  xAlbers = xMapProjection.xpath('./albers')
+                  unless xAlbers.empty?
+                     return AlbersProjection.unpack(xAlbers, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.3 (azimequi) - Azimuthal Equidistant
+                  xAzimequi = xMapProjection.xpath('./azimequi')
+                  unless xAzimequi.empty?
+                     return AzimuthEquidistantProjection.unpack(xAzimequi, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.4 (equicon) - Equidistant Conic
+                  xEquiCon = xMapProjection.xpath('./equicon')
+                  unless xEquiCon.empty?
+                     return EquidistantConicProjection.unpack(xEquiCon, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.5 (equirect) - Equirectangular
+                  xEquirec = xMapProjection.xpath('./equirect')
+                  unless xEquirec.empty?
+                     return EquirectangularProjection.unpack(xEquirec, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.6 (gvnsp) - General Vertical Near-sided Perspective
+                  xGenVert = xMapProjection.xpath('./gvnsp')
+                  unless xGenVert.empty?
+                     return GeneralVerticalProjection.unpack(xGenVert, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.7 (gnomonic) - Gnomonic
+                  xGeomonic = xMapProjection.xpath('./gnomonic')
+                  unless xGeomonic.empty?
+                     return GnomonicProjection.unpack(xGeomonic, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.8 (lamberta) - Lambert Azimuthal Equal Area
+                  xLambertA = xMapProjection.xpath('./lamberta')
+                  unless xLambertA.empty?
+                     return LambertEqualAreaProjection.unpack(xLambertA, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.9 (lambertc) - Lambert Conformal Conic
+                  xLambertC = xMapProjection.xpath('./lambertc')
+                  unless xLambertC.empty?
+                     return LambertConicProjection.unpack(xLambertC, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.10 (mercator) - Mercator
+                  xMercator = xMapProjection.xpath('./mercator')
+                  unless xMercator.empty?
+                     return MercatorProjection.unpack(xMercator, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.11 (modsak) - Modified Stereographic for Alaska
+                  xAlaska = xMapProjection.xpath('./modsak')
+                  unless xAlaska.empty?
+                     return ModifiedAlaskaProjection.unpack(xAlaska, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.12 (miller) - Miller Cylindrical
+                  xMiller = xMapProjection.xpath('./miller')
+                  unless xMiller.empty?
+                     return MillerCylinderProjection.unpack(xMiller, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.13 (obqmerc) - Oblique Mercator
+                  xObliqueM = xMapProjection.xpath('./obqmerc')
+                  unless xObliqueM.empty?
+                     return ObliqueMercatorProjection.unpack(xObliqueM, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.14 (orthogr) - Orthographic
+                  xOrtho = xMapProjection.xpath('./orthogr')
+                  unless xOrtho.empty?
+                     return OrthographicProjection.unpack(xOrtho, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.15 (polarst) - Polar Stereographic
+                  xPolarStereo = xMapProjection.xpath('./polarst')
+                  unless xPolarStereo.empty?
+                     return PolarStereoProjection.unpack(xPolarStereo, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.16 (polycon) - Polyconic
+                  xPolycon = xMapProjection.xpath('./polycon')
+                  unless xPolycon.empty?
+                     return PolyconicProjection.unpack(xPolycon, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.17 (robinson) - Robinson
+                  xRobin = xMapProjection.xpath('./robinson')
+                  unless xRobin.empty?
+                     return RobinsonProjection.unpack(xRobin, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.18 (sinusoid) - Sinusoidal
+                  xSinu = xMapProjection.xpath('./sinusoid')
+                  unless xSinu.empty?
+                     return SinusoidalProjection.unpack(xSinu, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.19 (spaceobq) - Space Oblique Mercator (Landsat)
+                  xSpaceO = xMapProjection.xpath('./spaceobq')
+                  unless xSpaceO.empty?
+                     return SpaceObliqueProjection.unpack(xSpaceO, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.20 (stereo) - Stereographic
+                  xStereo = xMapProjection.xpath('./stereo')
+                  unless xStereo.empty?
+                     return StereographicProjection.unpack(xStereo, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.21 (transmer) - Transverse Mercator
+                  xTransMer = xMapProjection.xpath('./transmer')
+                  unless xTransMer.empty?
+                     return TransverseMercatorProjection.unpack(xTransMer, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.22 (vdgrin) - van der Grinten
+                  xVanDerG = xMapProjection.xpath('./vdgrin')
+                  unless xVanDerG.empty?
+                     return VanDerGrintenProjection.unpack(xVanDerG, hProjection, hResponseObj)
+                  end
+
+                  # map projection 4.1.2.1.23 (mapprojp) - projection parameter set
+                  xParamSet = xMapProjection.xpath('./mapprojp')
+                  unless xParamSet.empty?
+                     return ProjectionParameters.unpack(xParamSet, hProjection, hResponseObj)
+                  end
+
+                  return nil
+
+               end
+
+            end
+
+         end
+      end
+   end
+end