adiwg-mdtranslator 2.3.5 → 2.4.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 711b22c64a7c6acc742500ca82ceadfabfc979c3
-  data.tar.gz: 359f7b98507e8f2420f4c4cce21bff8719a55f28
+  metadata.gz: 20d12fe84ddfccbe18ea8ee3a2d4c28aedc471fd
+  data.tar.gz: 47d8c57d4668555bcf4b15ca67099675a9a8e9f4
 SHA512:
-  metadata.gz: f573f0909635e78ad98b365cdcb35c083ed9733617cc62575d64ed1c3b73413164fd1761bbb15782de42d4563c3f96a3d406e27c9da9655ddf811157de1ce4e8
-  data.tar.gz: 824d2cfc455565074996d874cff4f9b82f65d9ccb5b743dbb4d7daecfe5618aa88dd8c83a9d2c1ecf5df7e189532c54b94408989b7b5449fd42ba5efc16bbc15
+  metadata.gz: a3b11672f637c332273913e9d88005b88427c98cec0dc5d96275ec6f6823f5ceaffd0a1940988a14829ccec53476d19d6936fe68f65fde25af8dbc8e2a9d127f
+  data.tar.gz: d00a5b73bd938e76f396e996eca3ddf0a555b3cdb806950f50d8d29e127ce292a40dcb646721643cb1b64019f65603d00724a337780575a5e666c16de6ea0fd4

data/lib/adiwg/mdtranslator/internal/internal_metadata_obj.rb

@@ -780,21 +780,22 @@ class InternalMetadata
       {
          systemType: nil,
          systemIdentifier: {},
-         systemParameters: {}
+         systemParameterSet: {}
       }
    end
 
-   def newReferenceSystemParameters
+   def newReferenceSystemParameterSet
       {
          projection: {},
          ellipsoid: {},
-         datumIdentifier: {}
+         verticalDatum: {}
       }
    end
 
    def newProjection
       {
          projectionIdentifier: {},
+         projectionName: nil,
          zone: nil,
          standardParallel1: nil,
          standardParallel2: nil,
@@ -808,11 +809,18 @@ class InternalMetadata
          longitudeOfProjectionCenter: nil,
          latitudeOfProjectionCenter: nil,
          scaleFactorAtCenterLine: nil,
+         scaleFactorAtCentralMeridian: nil,
          straightVerticalLongitudeFromPole: nil,
          scaleFactorAtProjectionOrigin: nil,
          azimuthAngle: nil,
          azimuthMeasurePointLongitude: nil,
-         obliqueLinePoint: []
+         obliqueLinePoints: [],
+         landsatNumber: nil,
+         landsatPath: nil,
+         localPlanarDescription: nil,
+         localPlanarGeoreference: nil,
+         otherProjectionDescription: nil,
+         otherGridDescription: nil
       }
    end
 
@@ -826,6 +834,7 @@ class InternalMetadata
    def newEllipsoid
       {
          ellipsoidIdentifier: {},
+         ellipsoidName: nil,
          semiMajorAxis: nil,
          axisUnits: nil,
          denominatorOfFlatteningRatio: nil
@@ -836,7 +845,9 @@ class InternalMetadata
       {
          scaleFactor: nil,
          measure: {},
-         geographicMeasure: {},
+         coordinateResolution: {},
+         bearingDistanceResolution: {},
+         geographicResolution: {},
          levelOfDetail: nil
       }
    end
@@ -849,14 +860,44 @@ class InternalMetadata
       }
    end
 
-   def newGeographicMeasure
+   def newCoordinateResolution
       {
-         latitudeMeasure: nil,
-         longitudeMeasure: nil,
+         abscissaResolutionX: nil,
+         ordinateResolutionY: nil,
          unitOfMeasure: nil
       }
    end
 
+   def newBearingDistanceResolution
+      {
+         distanceResolution: nil,
+         distanceUnitOfMeasure: nil,
+         bearingResolution: nil,
+         bearingUnitOfMeasure: nil,
+         bearingReferenceDirection: nil,
+         bearingReferenceMeridian: nil
+      }
+   end
+
+   def newGeographicResolution
+      {
+         latitudeResolution: nil,
+         longitudeResolution: nil,
+         unitOfMeasure: nil
+      }
+   end
+
+   def newVerticalDatum
+      {
+         datumIdentifier: {},
+         isDepthSystem: false,
+         encodingMethod: nil,
+         verticalResolution: nil,
+         unitOfMeasure: nil
+      }
+
+   end
+
    def newSpatialRepresentation
       {
          gridRepresentation: {},

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

@@ -0,0 +1,57 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc map projection - albers
+
+# History:
+#  Stan Smith 2017-10-04 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+require_relative 'projection_common'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module AlbersProjection
+
+               def self.unpack(xParams, hProjection, hResponseObj)
+
+                  # map projection 4.1.2.1.2 (albers) - Albers Conical Equal Area
+                  unless xParams.empty?
+                     paramCount = 0
+                     hProjection[:projectionName] = 'albers'
+
+                     # -> ReferenceSystemParameters.projection.standardParallel1
+                     # -> ReferenceSystemParameters.projection.standardParallel2
+                     paramCount += ProjectionCommon.unpackStandParallel(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)
+
+                     # verify parameter count
+                     if paramCount == 6
+                        return hProjection
+                     else
+                        hResponseObj[:readerExecutionMessages] << 'albers 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_azimuthEquidistant.rb

@@ -0,0 +1,53 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc map projection - azimuth equidistant
+
+# History:
+#  Stan Smith 2017-10-16 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+require_relative 'projection_common'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module AzimuthEquidistantProjection
+
+               def self.unpack(xParams, hProjection, hResponseObj)
+
+                  # map projection 4.1.2.1.3 (azimequi) - Azimuthal Equidistant
+                  unless xParams.empty?
+                     paramCount = 0
+                     hProjection[:projectionName] = 'azimuthal equidistant'
+
+                     # -> 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)
+
+                     # verify parameter count
+                     if paramCount == 4
+                        return hProjection
+                     else
+                        hResponseObj[:readerExecutionMessages] << 'azimuth equidistant 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_common.rb

@@ -0,0 +1,252 @@
+# Reader - fgdc to internal data structure
+# unpack fgdc map projection - map projection parent
+
+# History:
+#  Stan Smith 2017-10-16 original script
+
+require 'nokogiri'
+require 'adiwg/mdtranslator/internal/internal_metadata_obj'
+
+module ADIWG
+   module Mdtranslator
+      module Readers
+         module Fgdc
+
+            module ProjectionCommon
+
+               # instance classes needed in script
+               @intMetadataClass = InternalMetadata.new
+
+               # standard parallels
+               def self.unpackStandParallel(xParams, hProjection, hResponseObj)
+                  paramCount = 0
+                  aStandP = xParams.xpath('./stdparll')
+                  unless aStandP.empty?
+                     aStandP.each_with_index do |xStandP, index|
+                        standP = xStandP.text
+                        unless standP.empty?
+                           if index == 0
+                              hProjection[:standardParallel1] = standP.to_f
+                              paramCount += 1
+                           else
+                              hProjection[:standardParallel2] = standP.to_f
+                              paramCount += 1
+                           end
+                        end
+                     end
+                  end
+                  return paramCount
+               end
+
+               # longitude of central meridian
+               def self.unpackLongCM(xParams, hProjection, hResponseObj)
+                  longCM = xParams.xpath('./longcm').text
+                  unless longCM.empty?
+                     hProjection[:longitudeOfCentralMeridian] = longCM.to_f
+                     return 1
+                  end
+                  return 0
+               end
+
+               # latitude of projection origin
+               def self.unpackLatPO(xParams, hProjection, hResponseObj)
+                  latPO = xParams.xpath('./latprjo').text
+                  unless latPO.empty?
+                     hProjection[:latitudeOfProjectionOrigin] = latPO.to_f
+                     return 1
+                  end
+                  return 0
+               end
+
+               # false northing and easting(xParams, hProjection, hResponseObj)
+               def self.unpackFalseNE(xParams, hProjection, hResponseObj)
+                  paramCount = 0
+                  falseE = xParams.xpath('./feast').text
+                  unless falseE.empty?
+                     hProjection[:falseEasting] = falseE.to_f
+                     paramCount += 1
+                  end
+                  falseN = xParams.xpath('./fnorth').text
+                  unless falseN.empty?
+                     hProjection[:falseNorthing] = falseN.to_f
+                     paramCount += 1
+                  end
+                  return paramCount
+               end
+
+               # height of perspective point above surface
+               def self.unpackHeightAS(xParams, hProjection, hResponseObj)
+                  heightAS = xParams.xpath('./heightpt').text
+                  unless heightAS.empty?
+                     hProjection[:heightOfProspectivePointAboveSurface] = heightAS.to_f
+                     return 1
+                  end
+                  return 0
+               end
+
+               # longitude of projection center
+               def self.unpackLongPC(xParams, hProjection, hResponseObj)
+                  longPC = xParams.xpath('./longpc').text
+                  unless longPC.empty?
+                     hProjection[:longitudeOfProjectionCenter] = longPC.to_f
+                     return 1
+                  end
+                  return 0
+               end
+
+               # latitude of projection center
+               def self.unpackLatPC(xParams, hProjection, hResponseObj)
+                  latPC = xParams.xpath('./latprjc').text
+                  unless latPC.empty?
+                     hProjection[:latitudeOfProjectionCenter] = latPC.to_f
+                     return 1
+                  end
+                  return 0
+               end
+
+               # scale factor at equator
+               def self.unpackSFEquator(xParams, hProjection, hResponseObj)
+                  sFEquator = xParams.xpath('./sfequat').text
+                  unless sFEquator.empty?
+                     hProjection[:scaleFactorAtEquator] = sFEquator.to_f
+                     return 1
+                  end
+                  return 0
+               end
+
+               # scale factor at center line
+               def self.unpackSFCenter(xParams, hProjection, hResponseObj)
+                  sFCenter = xParams.xpath('./sfctrlin').text
+                  unless sFCenter.empty?
+                     hProjection[:scaleFactorAtCenterLine] = sFCenter.to_f
+                     return 1
+                  end
+                  return 0
+               end
+
+               # scale factor at central meridian
+               def self.unpackSFCM(xParams, hProjection, hResponseObj)
+                  sFCM = xParams.xpath('./sfctrmer').text
+                  unless sFCM.empty?
+                     hProjection[:scaleFactorAtCentralMeridian] = sFCM.to_f
+                     return 1
+                  end
+                  return 0
+               end
+
+               # scale factor at projection origin
+               def self.unpackSFPO(xParams, hProjection, hResponseObj)
+                  sFProjectO = xParams.xpath('./sfprjorg').text
+                  unless sFProjectO.empty?
+                     hProjection[:scaleFactorAtProjectionOrigin] = sFProjectO.to_f
+                     return 1
+                  end
+                  return 0
+               end
+
+               # oblique line azimuth
+               def self.unpackObliqueLA(xParams, hProjection, hResponseObj)
+                  xObliqueLA = xParams.xpath('./obqlazim')
+                  unless xObliqueLA.empty?
+                     paramCount = 0
+                     lineAzimuth = xObliqueLA.xpath('./asimangl').text
+                     unless lineAzimuth.empty?
+                        hProjection[:azimuthAngle] = lineAzimuth.to_f
+                        paramCount += 1
+                     end
+                     lineLong = xObliqueLA.xpath('./asimptl').text
+                     unless lineLong.empty?
+                        hProjection[:azimuthMeasurePointLongitude] = lineLong.to_f
+                        paramCount += 1
+                     end
+                     return paramCount
+                  end
+                  return 0
+               end
+
+               # oblique line point
+               def self.unpackObliqueLP(xParams, hProjection, hResponseObj)
+                  xObliqueLP = xParams.xpath('./obqlpt')
+                  unless xObliqueLP.empty?
+                     paramCount = 0
+                     hPoint1 = @intMetadataClass.newObliqueLinePoint
+                     hPoint2 = @intMetadataClass.newObliqueLinePoint
+
+                     lat1 = xObliqueLP.xpath('./obqllat[1]').text
+                     unless lat1.empty?
+                        hPoint1[:azimuthLineLatitude] = lat1.to_f
+                        paramCount += 1
+                     end
+                     long1 = xObliqueLP.xpath('./obqllong[1]').text
+                     unless long1.empty?
+                        hPoint1[:azimuthLineLongitude] = long1.to_f
+                        paramCount += 1
+                     end
+
+                     lat2 = xObliqueLP.xpath('./obqllat[2]').text
+                     unless lat2.empty?
+                        hPoint2[:azimuthLineLatitude] = lat2.to_f
+                        paramCount += 1
+                     end
+                     long2 = xObliqueLP.xpath('./obqllong[2]').text
+                     unless long2.empty?
+                        hPoint2[:azimuthLineLongitude] = long2.to_f
+                        paramCount += 1
+                     end
+
+                     if paramCount == 4
+                        hProjection[:obliqueLinePoints] << hPoint1
+                        hProjection[:obliqueLinePoints] << hPoint2
+                        return hProjection[:obliqueLinePoints].length
+                     end
+
+                  end
+                  return 0
+               end
+
+               # straight vertical longitude from pole
+               def self.unpackVSLong(xParams, hProjection, hResponseObj)
+                  sVLongP = xParams.xpath('./svlong').text
+                  unless sVLongP.empty?
+                     hProjection[:straightVerticalLongitudeFromPole] = sVLongP.to_f
+                     return 1
+                  end
+                  return 0
+               end
+
+               # landsat number
+               def self.unpackLandSat(xParams, hProjection, hResponseObj)
+                  landsatN = xParams.xpath('./landsat').text
+                  unless landsatN.empty?
+                     hProjection[:landsatNumber] = landsatN.to_i
+                     return 1
+                  end
+                  return 0
+               end
+
+               # landsat path number
+               def self.unpackLandSatPath(xParams, hProjection, hResponseObj)
+                  landsatP = xParams.xpath('./pathnum').text
+                  unless landsatP.empty?
+                     hProjection[:landsatPath] = landsatP.to_i
+                     return 1
+                  end
+                  return 0
+               end
+
+               # other projection
+               def self.unpackOtherProjection(xParams, hProjection, hResponseObj)
+                  other = xParams.xpath('./otherprj').text
+                  unless other.empty?
+                     hProjection[:otherProjectionDescription] = other
+                     return 1
+                  end
+                  return 0
+               end
+
+            end
+
+         end
+      end
+   end
+end