geo_concerns 0.0.1

data/app/models/concerns/geo_concerns/vector_file_behavior.rb

@@ -0,0 +1,14 @@
+module GeoConcerns
+  # Attributes and methods for vector files
+  module VectorFileBehavior
+    extend ActiveSupport::Concern
+    include ::GeoConcerns::GeoreferencedBehavior
+    # Retrieve the Vector Work of which this Object is a member
+    # @return [GeoConcerns::VectorWork]
+    def vector_work
+      generic_works.select do |parent|
+        parent.class.included_modules.include?(::GeoConcerns::VectorWorkBehavior)
+      end.to_a
+    end
+  end
+end

data/app/models/concerns/geo_concerns/vector_work_behavior.rb

@@ -0,0 +1,78 @@
+module GeoConcerns
+  # Attributes and methods for vector works
+  module VectorWorkBehavior
+    extend ActiveSupport::Concern
+    include ::GeoConcerns::MetadataExtractionHelper
+
+    included do
+      type [Hydra::PCDM::Vocab::PCDMTerms.Object,
+            Hydra::Works::Vocab::WorksTerms.GenericWork,
+            ::GeoConcerns::GeoTerms.VectorWork]
+    end
+
+    def vector_files
+      members.select(&:vector_file?)
+    end
+
+    def metadata_files
+      members.select(&:external_metadata_file?)
+    end
+
+    # Defines type by what it is and isn't
+    # @return [Boolean]
+    def image_work?
+      false
+    end
+
+    def image_file?
+      false
+    end
+
+    def raster_work?
+      false
+    end
+
+    def raster_file?
+      false
+    end
+
+    def vector_work?
+      true
+    end
+
+    def vector_file?
+      false
+    end
+
+    def external_metadata_file?
+      false
+    end
+
+    # Retrieve all Raster Works for which this Vector Work can be extracted
+    # @return [Array]
+    def raster_works
+      ordered_by.select do |parent|
+        parent.class.included_modules.include?(::GeoConcerns::RasterWorkBehavior)
+      end
+    end
+
+    # Retrieve the only Raster Work for which feature extraction generates this Vector Work
+    # @return [GeoConcerns::RasterWork]
+    def raster_work
+      raster_works.first
+    end
+
+    def to_solr(solr_doc = {})
+      super.tap do |doc|
+        doc[solr_name("ordered_by", :symbol)] ||= []
+        doc[solr_name("ordered_by", :symbol)] += send(:ordered_by_ids)
+      end
+    end
+
+    private
+
+      def solr_name(*args)
+        ActiveFedora.index_field_mapper.solr_name(*args)
+      end
+  end
+end

data/app/presenters/geo_concerns/geo_concerns_show_presenter.rb

@@ -0,0 +1,30 @@
+module GeoConcerns
+  class GeoConcernsShowPresenter < CurationConcerns::WorkShowPresenter
+    delegate :has?, :first, to: :solr_document
+    delegate :spatial, :temporal, :issued, :coverage, :provenance, to: :solr_document
+
+    def members(presenter)
+      # TODO: member ids appear twice in member_ids_ssim.
+      # Figure out why instead of removing duplicates.
+      ids = solr_document.fetch('member_ids_ssim', [])
+      CurationConcerns::PresenterFactory.build_presenters(ids.uniq,
+                                                          presenter,
+                                                          current_ability)
+    end
+
+    def external_metadata_file_formats_presenters
+      # filter for external metadata files
+      members(::FileSetPresenter).select do |member|
+        MetadataFormatService.include? member.solr_document[:mime_type_ssi]
+      end
+    end
+
+    def attribute_to_html(field, options = {})
+      if field == :coverage
+        ::CoverageRenderer.new(field, send(field), options).render
+      else
+        super field, options
+      end
+    end
+  end
+end

data/app/presenters/geo_concerns/image_work_show_presenter.rb

@@ -0,0 +1,18 @@
+module GeoConcerns
+  class ImageWorkShowPresenter < GeoConcernsShowPresenter
+    def raster_work_presenters
+      # filter for raster works
+      members(::GeoConcerns::RasterWorkShowPresenter).select do |member|
+        format = member.solr_document[:has_model_ssim]
+        format ? format.first == 'RasterWork' : false
+      end
+    end
+
+    def image_file_presenters
+      # filter for image files
+      members(::FileSetPresenter).select do |member|
+        ImageFormatService.include? member.solr_document[:mime_type_ssi]
+      end
+    end
+  end
+end

data/app/presenters/geo_concerns/raster_work_show_presenter.rb

@@ -0,0 +1,18 @@
+module GeoConcerns
+  class RasterWorkShowPresenter < GeoConcernsShowPresenter
+    def vector_work_presenters
+      # filter for vector works
+      members(::GeoConcerns::VectorWorkShowPresenter).select do |member|
+        format = member.solr_document[:has_model_ssim]
+        format ? format.first == 'VectorWork' : false
+      end
+    end
+
+    def raster_file_presenters
+      # filter for raster files
+      members(::FileSetPresenter).select do |member|
+        RasterFormatService.include? member.solr_document[:mime_type_ssi]
+      end
+    end
+  end
+end

data/app/presenters/geo_concerns/vector_work_show_presenter.rb

@@ -0,0 +1,10 @@
+module GeoConcerns
+  class VectorWorkShowPresenter < GeoConcernsShowPresenter
+    def vector_file_presenters
+      # filter for vector files
+      members(::FileSetPresenter).select do |member|
+        VectorFormatService.include? member.solr_document[:mime_type_ssi]
+      end
+    end
+  end
+end

data/app/processors/geo_concerns/processors/base_geo_processor.rb

@@ -0,0 +1,86 @@
+module GeoConcerns
+  module Processors
+    module BaseGeoProcessor
+      extend ActiveSupport::Concern
+
+      included do
+        # Returns a formatted gdal_translate command. Used to translate a raster
+        # format into a different format. Also used in generating thumbnails
+        # from vector data.
+        #
+        # @param in_path [String] file input path
+        # #param options [Hash] creation options
+        # @param out_path [String] processor output file path
+        # @return [String] command for transforming a raster dataset
+        def self.translate(in_path, options, out_path)
+          "gdal_translate -outsize #{options[:output_size]} -q -ot Byte "\
+            "-of #{options[:output_format]} \"#{in_path}\" #{out_path}"
+        end
+
+        # Returns a path to an intermediate temp file.
+        #
+        # @param path [String] input file path to base temp path on
+        # @return [String] tempfile path
+        def self.intermediate_file_path(path)
+          ext = File.extname(path)
+          "#{File.dirname(path)}/#{File.basename(path, ext)}_temp#{ext}"
+        end
+      end
+
+      def options_for(_format)
+        {
+          label: label,
+          output_format: output_format,
+          output_size: output_size,
+          output_srid: output_srid,
+          basename: basename
+        }
+      end
+
+      # Returns the label directive or an empty string.
+      #
+      # @return [Sting] output label
+      def label
+        directives.fetch(:label, '')
+      end
+
+      # Tranforms the format directive into a GDAL output format.
+      #
+      # @return [Sting] derivative output format
+      def output_format
+        format = directives.fetch(:format, '').upcase
+        case format
+        when 'JPG'
+          'JPEG'
+        when 'TIF'
+          'GTiff'
+        else
+          format
+        end
+      end
+
+      # Tranforms the size directive into a GDAL size parameter.
+      #
+      # @return [String] derivative size
+      def output_size
+        return unless directives[:size]
+        directives[:size].tr('x', ' ')
+      end
+
+      # Gets srid for reprojection derivative or returns WGS 84.
+      #
+      # @return [String] spatial reference code
+      def output_srid
+        directives.fetch(:srid, 'EPSG:4326')
+        # directives[:srid] ? directives[:srid] : 'EPSG:4326'
+      end
+
+      # Extracts the base file name (without extension) from the source file path.
+      #
+      # @return [String] base file name for source
+      def basename
+        File.basename(source_path, File.extname(source_path))
+      end
+    end
+  end
+end

data/app/processors/geo_concerns/processors/raster.rb

@@ -0,0 +1,11 @@
+module GeoConcerns
+  module Processors
+    module Raster
+      extend ActiveSupport::Autoload
+
+      eager_autoload do
+        autoload :Processor
+      end
+    end
+  end
+end

data/app/processors/geo_concerns/processors/raster/aig.rb

@@ -0,0 +1,44 @@
+module GeoConcerns
+  module Processors
+    module Raster
+      class Aig < GeoConcerns::Processors::Raster::Base
+        include GeoConcerns::Processors::Zip
+
+        def self.encode(path, options, output_file)
+          unzip(path, output_file) do |zip_path|
+            info = gdalinfo(zip_path)
+            options[:min_max] = get_raster_min_max(info)
+            case options[:label]
+            when :thumbnail
+              encode_raster(zip_path, options, output_file)
+            when :display_raster
+              reproject_raster(zip_path, options, output_file)
+            end
+          end
+        end
+
+        def self.reproject_raster(in_path, options, out_path)
+          options[:output_size] = '100% 100%'
+          intermediate_file = intermediate_file_path(out_path)
+          execute warp(in_path, options, intermediate_file)
+          execute translate(intermediate_file, options, out_path)
+          FileUtils.rm_rf(intermediate_file)
+        end
+
+        # Returns a formatted gdal_translate command to translate a raster
+        # format into a different format with a scaling options. This command
+        # scales the min and max values of the raster into the 0 to 255 range.
+        # Scale is inverted (255 to 0) to create a better visualization.
+        #
+        # @param in_path [String] file input path
+        # #param options [Hash] creation options
+        # @param out_path [String] processor output file path
+        # @return [String] command for tranforming a usgs dem dataset
+        def self.translate(in_path, options, out_path)
+          "gdal_translate -scale #{options[:min_max]} 255 0 -outsize #{options[:output_size]} "\
+             "-q -ot Byte -of #{options[:output_format]} \"#{in_path}\" #{out_path}"
+        end
+      end
+    end
+  end
+end

data/app/processors/geo_concerns/processors/raster/base.rb

@@ -0,0 +1,74 @@
+module GeoConcerns
+  module Processors
+    module Raster
+      class Base < Hydra::Derivatives::Processors::Processor
+        include Hydra::Derivatives::Processors::ShellBasedProcessor
+        include GeoConcerns::Processors::BaseGeoProcessor
+
+        def self.encode(path, options, output_file)
+          case options[:label]
+          when :thumbnail
+            encode_raster(path, options, output_file)
+          when :display_raster
+            reproject_raster(path, options, output_file)
+          end
+        end
+
+        def self.encode_raster(in_path, options, out_path)
+          execute translate(in_path, options, out_path)
+          File.unlink("#{out_path}.aux.xml")
+        end
+
+        def self.reproject_raster(in_path, options, out_path)
+          intermediate_file = intermediate_file_path(out_path)
+          execute warp(in_path, options, intermediate_file)
+          execute compress(intermediate_file, options, out_path)
+          FileUtils.rm_rf(intermediate_file)
+        end
+
+        # Returns a formatted gdalwarp. Used to transform a raster
+        # from one projection into another.
+        #
+        # @param in_path [String] file input path
+        # #param options [Hash] creation options
+        # @param out_path [String] processor output file path
+        # @return [String] command for reprojecting a raster
+        def self.warp(in_path, options, out_path)
+          "gdalwarp -q -r bilinear -t_srs #{options[:output_srid]} "\
+            " #{in_path} #{out_path} -co 'COMPRESS=NONE'"
+        end
+
+        # Returns a formatted gdal_translate command. Used compress
+        # an previously uncompressed raster.
+        #
+        # @param in_path [String] file input path
+        # #param options [Hash] creation options
+        # @param out_path [String] processor output file path
+        # @return [String] command for compressing a raster dataset
+        def self.compress(in_path, options, out_path)
+          "gdal_translate -q -ot Byte -a_srs #{options[:output_srid]} "\
+            "\"#{in_path}\" #{out_path} -co 'COMPRESS=LZW'"
+        end
+
+        # Given an output string from the gdalinfo command, returns
+        # a formatted string for the computed min and max values.
+        #
+        # @param info_string [String] ouput from gdalinfo
+        # @return [String] computed min and max values
+        def self.get_raster_min_max(info_string)
+          match = %r{(?<=Computed Min/Max=).*?(?=\s)}.match(info_string)
+          match ? match[0].tr(',', ' ') : ''
+        end
+
+        # Runs the gdalinfo command and returns the result as a string.
+        #
+        # @param path [String] path to raster file
+        # @return [String] output of gdalinfo
+        def self.gdalinfo(path)
+          stdout, _stderr, _status = Open3.capture3("gdalinfo -mm #{path}")
+          stdout
+        end
+      end
+    end
+  end
+end

data/app/processors/geo_concerns/processors/raster/dem.rb

@@ -0,0 +1,46 @@
+module GeoConcerns
+  module Processors
+    module Raster
+      class Dem < GeoConcerns::Processors::Raster::Base
+        def self.encode_raster(in_path, options, out_path)
+          intermediate_file = intermediate_file_path(out_path)
+          execute translate(in_path, options, intermediate_file)
+          execute hillshade(intermediate_file, options, out_path)
+          FileUtils.rm_rf(intermediate_file)
+          File.unlink("#{out_path}.aux.xml")
+        end
+
+        def self.reproject_raster(in_path, options, out_path)
+          intermediate_file = intermediate_file_path(out_path)
+          execute hillshade(in_path, options, intermediate_file)
+          execute warp(intermediate_file, options, out_path)
+          FileUtils.rm_rf(intermediate_file)
+        end
+
+        # Returns a formatted gdal_translate command to translate a vector
+        # in USGS DEM format into a different format.
+        #
+        # @param in_path [String] file input path
+        # #param options [Hash] creation options
+        # @param out_path [String] processor output file path
+        # @return [String] command for tranforming a usgs dem dataset
+        def self.translate(in_path, options, out_path)
+          "gdal_translate -outsize #{options[:output_size]} -q -ot Byte "\
+            "-of USGSDEM #{in_path} #{out_path}"
+        end
+
+        # Returns a formatted gdal hillshade command. Calculates hillshade
+        # on a raster that contains elevation data.
+        #
+        # @param in_path [String] file input path
+        # #param options [Hash] creation options
+        # @param output_file [String] processor output file path
+        # @return [String] command for generating a hillshade
+        def self.hillshade(in_path, options, out_path)
+          "gdaldem hillshade -q "\
+            "-of #{options[:output_format]} \"#{in_path}\" #{out_path}"
+        end
+      end
+    end
+  end
+end