Better display multiple sample locations

Author: lmrodriguezr

app/models/genome.rb

@@ -252,15 +252,18 @@ def source_attribute_groups
     @source_attribute_groups
   end
 
+  ##
+  # Finds the locations of all source samples associated to this genome, and
+  # returns them as an Array of 2-element Arrays ([lat, lon]) or +nil+
   def source_sample_locations
     coord = /([-+] *)?(\d+(?:[\.\,]\d+)?|\d+°(?:\d+['"])*)( *[NSEW])?/
     keys = {
       lat: %i[lat geographic_location_latitude latitude_start latitude_end],
       lon: %i[lon geographic_location_longitude longitude_start longitude_end]
     }
 
+    coords = { lat: nil, lon: nil }
     @_source_sample_locations ||=
-      coords = { lat: nil, lon: nil }
       source_cannonical_samples.map do |sample|
         # Try joint keys
         if sample[:lat_lon]
@@ -306,6 +309,29 @@ def source_sample_locations
       end
   end
 
+  ##
+  # Finds the rectangular bounds of all sample locations, with a minimum range
+  # of latitudes of +minlat+ and longitudes of +minlon+, and returns it as an
+  # array in the [south, west, north, east] order
+  def source_sample_area(minlat = 0.1, minlon = 0.1)
+    loc = source_sample_locations.compact
+    return unless loc.present?
+
+    min = { lat: minlat, lon: minlon }
+    rng = { lat: loc.map { |i| i[0] }.minmax, lon: loc.map { |i| i[1] }.minmax }
+
+    rng.each do |k, v|
+      width = v.inject(:-).abs
+      if width < min[k]
+        pad = (min[k] - width) / 2
+        rng[k][0] -= pad
+        rng[k][1] += pad
+      end
+    end
+
+    [rng[:lat][0], rng[:lon][0], rng[:lat][1], rng[:lon][1]]
+  end
+
   ##
   # TODO
   # Use source_cannonical_samples instead!

app/views/genomes/sample_map.html.erb

@@ -19,15 +19,15 @@
   <% end %>
 
   <script>
-    function load_map(lat, lon) {
+    function load_map(lat1, lon1, lat2, lon2) {
       return(
         new maplibregl.Map({
           container: 'map',
           style:
             'https://api.maptiler.com/maps/hybrid/style.json?' +
-              'key=<%= ENV["RAILS_SEQCODE_MAPTILER_KEY"] %>',
-          center: [lon, lat], // starting position [lng, lat]
-          zoom: 9 // starting zoom
+              'key=<%= ENV["RAILS_SEQCODE_MAPTILER_KEY"] ||
+                         "get_your_own_OpIi9ZULNHzrESv6T2vL" %>',
+          bounds: [lon1, lat1, lon2, lat2]
         })
       );
     }
@@ -38,7 +38,8 @@
           .addTo(map)
       );
     }
-    var map = load_map(<%= locs.first[0] %>, <%= locs.first[1] %>);
+    <% bounds = @genome.source_sample_area %>
+    var map = load_map(<%= bounds.join(', ') %>);
     <% locs.each do |lat_lon| %>
       var marker = add_marker(map, <%= lat_lon[0] %>, <%= lat_lon[1] %>);
     <% end %>