Merge pull request #1218 from cbworden/sta2feature

Sta2feature

Author: cbworden

code.json

@@ -45,6 +45,49 @@
         "date": {
           "metadataLastUpdated": "2018-08-23"
         }
+      },
+      {
+        "name": "shakemap",
+        "organization": "U.S. Geological Survey",
+        "description": "Near-real-time maps of ground motion and shaking intensity",
+        "version": "4.1.0",
+        "status": "Production",
+
+        "permissions": {
+          "usageType": "openSource",
+          "licenses": [
+            {
+              "name": "Public Domain, CC0-1.0",
+              "URL": "https://github.com/usgs/shakemap/blob/master/LICENSE.md"
+            }
+          ]
+        },
+
+        "homepageURL": "https://usgs.github.io/shakemap/sm4_index.html",
+        "downloadURL": "https://github.com/usgs/shakemap/archive/master.zip",
+        "disclaimerURL": "https://github.com/usgs/shakemap/blob/master/LICENSE.md",
+        "repositoryURL": "https://github.com/usgs/shakemap.git",
+        "vcs": "git",
+
+        "laborHours": null,
+
+        "tags": [
+          "earthquake",
+          "shakemap"
+        ],
+
+        "languages": [
+          "python"
+        ],
+
+        "contact": {
+          "name": "C. Bruce Worden",
+          "email": "[email protected] "
+        },
+
+        "date": {
+          "metadataLastUpdated": "2022-04-28"
+        }
       }
     ]
 }

shakemap/mapping/mapmaker.py

@@ -1418,7 +1418,7 @@ def arclen(path):
             props = contour_object["properties"]
             multi_lines = sShape(contour_object["geometry"])
             pmulti_lines = proj.project_geometry(multi_lines, src_crs=geoproj)
-            for multi_line in pmulti_lines:
+            for multi_line in pmulti_lines.geoms:
                 pmulti_line = mapping(multi_line)["coordinates"]
                 x, y = zip(*pmulti_line)
                 contour_lens[props["value"]].append(arclen(pmulti_line))
@@ -1444,7 +1444,7 @@ def arclen(path):
             # isovalue)
             min_len = np.array(contour_lens[props["value"]]).mean()
 
-            for multi_line in pmulti_lines:
+            for multi_line in pmulti_lines.geoms:
                 pmulti_line = mapping(multi_line)["coordinates"]
                 x, y = zip(*pmulti_line)
                 # color = imt_cmap.getDataColor(props['value'])