proj4 2.20.2 → 2.20.4

package/dist/proj4-src.js

@@ -259,7 +259,7 @@
 
   class PROJJSONBuilderBase {
     static getId(node) {
-      const idNode = node.find((child) => Array.isArray(child) && child[0] === "ID");
+      const idNode = node.find((child) => Array.isArray(child) && child[0] === 'ID');
       if (idNode && idNode.length >= 3) {
         return {
           authority: idNode[1],
@@ -269,20 +269,20 @@
       return null;
     }
 
-    static convertUnit(node, type = "unit") {
+    static convertUnit(node, type = 'unit') {
       if (!node || node.length < 3) {
-        return { type, name: "unknown", conversion_factor: null };
+        return { type, name: 'unknown', conversion_factor: null };
       }
 
       const name = node[1];
       const conversionFactor = parseFloat(node[2]) || null;
 
-      const idNode = node.find((child) => Array.isArray(child) && child[0] === "ID");
+      const idNode = node.find((child) => Array.isArray(child) && child[0] === 'ID');
       const id = idNode
         ? {
-            authority: idNode[1],
-            code: parseInt(idNode[2], 10),
-          }
+          authority: idNode[1],
+          code: parseInt(idNode[2], 10),
+        }
         : null;
 
       return {
@@ -294,7 +294,7 @@
     }
 
     static convertAxis(node) {
-      const name = node[1] || "Unknown";
+      const name = node[1] || 'Unknown';
 
       // Determine the direction
       let direction;
@@ -308,16 +308,16 @@
         else throw new Error(`Unknown axis abbreviation: ${abbreviation}`);
       } else {
         // Use the explicit direction provided in the AXIS node
-        direction = node[2]?.toLowerCase() || "unknown";
+        direction = node[2] ? node[2].toLowerCase() : 'unknown';
       }
 
-      const orderNode = node.find((child) => Array.isArray(child) && child[0] === "ORDER");
+      const orderNode = node.find((child) => Array.isArray(child) && child[0] === 'ORDER');
       const order = orderNode ? parseInt(orderNode[1], 10) : null;
 
       const unitNode = node.find(
         (child) =>
           Array.isArray(child) &&
-          (child[0] === "LENGTHUNIT" || child[0] === "ANGLEUNIT" || child[0] === "SCALEUNIT")
+          (child[0] === 'LENGTHUNIT' || child[0] === 'ANGLEUNIT' || child[0] === 'SCALEUNIT')
       );
       const unit = this.convertUnit(unitNode);
 
@@ -331,7 +331,7 @@
 
     static extractAxes(node) {
       return node
-        .filter((child) => Array.isArray(child) && child[0] === "AXIS")
+        .filter((child) => Array.isArray(child) && child[0] === 'AXIS')
         .map((axis) => this.convertAxis(axis))
         .sort((a, b) => (a.order || 0) - (b.order || 0)); // Sort by the "order" property
     }
@@ -339,17 +339,17 @@
     static convert(node, result = {}) {
 
       switch (node[0]) {
-        case "PROJCRS":
-          result.type = "ProjectedCRS";
+        case 'PROJCRS':
+          result.type = 'ProjectedCRS';
           result.name = node[1];
-          result.base_crs = node.find((child) => Array.isArray(child) && child[0] === "BASEGEOGCRS")
-            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "BASEGEOGCRS"))
+          result.base_crs = node.find((child) => Array.isArray(child) && child[0] === 'BASEGEOGCRS')
+            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === 'BASEGEOGCRS'))
             : null;
-          result.conversion = node.find((child) => Array.isArray(child) && child[0] === "CONVERSION")
-            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "CONVERSION"))
+          result.conversion = node.find((child) => Array.isArray(child) && child[0] === 'CONVERSION')
+            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === 'CONVERSION'))
             : null;
 
-          const csNode = node.find((child) => Array.isArray(child) && child[0] === "CS");
+          const csNode = node.find((child) => Array.isArray(child) && child[0] === 'CS');
           if (csNode) {
             result.coordinate_system = {
               type: csNode[1],
@@ -357,7 +357,7 @@
             };
           }
 
-          const lengthUnitNode = node.find((child) => Array.isArray(child) && child[0] === "LENGTHUNIT");
+          const lengthUnitNode = node.find((child) => Array.isArray(child) && child[0] === 'LENGTHUNIT');
           if (lengthUnitNode) {
             const unit = this.convertUnit(lengthUnitNode);
             result.coordinate_system.unit = unit; // Add unit to coordinate_system
@@ -366,23 +366,23 @@
           result.id = this.getId(node);
           break;
 
-        case "BASEGEOGCRS":
-        case "GEOGCRS":
-          result.type = "GeographicCRS";
+        case 'BASEGEOGCRS':
+        case 'GEOGCRS':
+          result.type = 'GeographicCRS';
           result.name = node[1];
         
           // Handle DATUM or ENSEMBLE
           const datumOrEnsembleNode = node.find(
-            (child) => Array.isArray(child) && (child[0] === "DATUM" || child[0] === "ENSEMBLE")
+            (child) => Array.isArray(child) && (child[0] === 'DATUM' || child[0] === 'ENSEMBLE')
           );
           if (datumOrEnsembleNode) {
             const datumOrEnsemble = this.convert(datumOrEnsembleNode);
-            if (datumOrEnsembleNode[0] === "ENSEMBLE") {
+            if (datumOrEnsembleNode[0] === 'ENSEMBLE') {
               result.datum_ensemble = datumOrEnsemble;
             } else {
               result.datum = datumOrEnsemble;
             }
-            const primem = node.find((child) => Array.isArray(child) && child[0] === "PRIMEM");
+            const primem = node.find((child) => Array.isArray(child) && child[0] === 'PRIMEM');
             if (primem && primem[1] !== 'Greenwich') {
               datumOrEnsemble.prime_meridian = {
                 name: primem[1],
@@ -392,42 +392,42 @@
           }
         
           result.coordinate_system = {
-            type: "ellipsoidal",
+            type: 'ellipsoidal',
             axis: this.extractAxes(node),
           };
         
           result.id = this.getId(node);
           break;
 
-        case "DATUM":
-          result.type = "GeodeticReferenceFrame";
+        case 'DATUM':
+          result.type = 'GeodeticReferenceFrame';
           result.name = node[1];
-          result.ellipsoid = node.find((child) => Array.isArray(child) && child[0] === "ELLIPSOID")
-            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "ELLIPSOID"))
+          result.ellipsoid = node.find((child) => Array.isArray(child) && child[0] === 'ELLIPSOID')
+            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === 'ELLIPSOID'))
             : null;
           break;
         
-        case "ENSEMBLE":
-          result.type = "DatumEnsemble";
+        case 'ENSEMBLE':
+          result.type = 'DatumEnsemble';
           result.name = node[1];
         
           // Extract ensemble members
           result.members = node
-            .filter((child) => Array.isArray(child) && child[0] === "MEMBER")
+            .filter((child) => Array.isArray(child) && child[0] === 'MEMBER')
             .map((member) => ({
-              type: "DatumEnsembleMember",
+              type: 'DatumEnsembleMember',
               name: member[1],
               id: this.getId(member), // Extract ID as { authority, code }
             }));
         
           // Extract accuracy
-          const accuracyNode = node.find((child) => Array.isArray(child) && child[0] === "ENSEMBLEACCURACY");
+          const accuracyNode = node.find((child) => Array.isArray(child) && child[0] === 'ENSEMBLEACCURACY');
           if (accuracyNode) {
             result.accuracy = parseFloat(accuracyNode[1]);
           }
         
           // Extract ellipsoid
-          const ellipsoidNode = node.find((child) => Array.isArray(child) && child[0] === "ELLIPSOID");
+          const ellipsoidNode = node.find((child) => Array.isArray(child) && child[0] === 'ELLIPSOID');
           if (ellipsoidNode) {
             result.ellipsoid = this.convert(ellipsoidNode); // Convert the ellipsoid node
           }
@@ -436,66 +436,66 @@
           result.id = this.getId(node);
           break;
 
-        case "ELLIPSOID":
-          result.type = "Ellipsoid";
+        case 'ELLIPSOID':
+          result.type = 'Ellipsoid';
           result.name = node[1];
           result.semi_major_axis = parseFloat(node[2]);
           result.inverse_flattening = parseFloat(node[3]);
-          node.find((child) => Array.isArray(child) && child[0] === "LENGTHUNIT")
-            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "LENGTHUNIT"), result)
+          node.find((child) => Array.isArray(child) && child[0] === 'LENGTHUNIT')
+            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === 'LENGTHUNIT'), result)
             : null;
           break;
 
-        case "CONVERSION":
-          result.type = "Conversion";
+        case 'CONVERSION':
+          result.type = 'Conversion';
           result.name = node[1];
-          result.method = node.find((child) => Array.isArray(child) && child[0] === "METHOD")
-            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "METHOD"))
+          result.method = node.find((child) => Array.isArray(child) && child[0] === 'METHOD')
+            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === 'METHOD'))
             : null;
           result.parameters = node
-            .filter((child) => Array.isArray(child) && child[0] === "PARAMETER")
+            .filter((child) => Array.isArray(child) && child[0] === 'PARAMETER')
             .map((param) => this.convert(param));
           break;
 
-        case "METHOD":
-          result.type = "Method";
+        case 'METHOD':
+          result.type = 'Method';
           result.name = node[1];
           result.id = this.getId(node);
           break;
 
-        case "PARAMETER":
-          result.type = "Parameter";
+        case 'PARAMETER':
+          result.type = 'Parameter';
           result.name = node[1];
           result.value = parseFloat(node[2]);
           result.unit = this.convertUnit(
             node.find(
               (child) =>
                 Array.isArray(child) &&
-                (child[0] === "LENGTHUNIT" || child[0] === "ANGLEUNIT" || child[0] === "SCALEUNIT")
+                (child[0] === 'LENGTHUNIT' || child[0] === 'ANGLEUNIT' || child[0] === 'SCALEUNIT')
             )
           );
           result.id = this.getId(node);
           break;
 
-        case "BOUNDCRS":
-          result.type = "BoundCRS";
+        case 'BOUNDCRS':
+          result.type = 'BoundCRS';
 
           // Process SOURCECRS
-          const sourceCrsNode = node.find((child) => Array.isArray(child) && child[0] === "SOURCECRS");
+          const sourceCrsNode = node.find((child) => Array.isArray(child) && child[0] === 'SOURCECRS');
           if (sourceCrsNode) {
             const sourceCrsContent = sourceCrsNode.find((child) => Array.isArray(child));
             result.source_crs = sourceCrsContent ? this.convert(sourceCrsContent) : null;
           }
 
           // Process TARGETCRS
-          const targetCrsNode = node.find((child) => Array.isArray(child) && child[0] === "TARGETCRS");
+          const targetCrsNode = node.find((child) => Array.isArray(child) && child[0] === 'TARGETCRS');
           if (targetCrsNode) {
             const targetCrsContent = targetCrsNode.find((child) => Array.isArray(child));
             result.target_crs = targetCrsContent ? this.convert(targetCrsContent) : null;
           }
 
           // Process ABRIDGEDTRANSFORMATION
-          const transformationNode = node.find((child) => Array.isArray(child) && child[0] === "ABRIDGEDTRANSFORMATION");
+          const transformationNode = node.find((child) => Array.isArray(child) && child[0] === 'ABRIDGEDTRANSFORMATION');
           if (transformationNode) {
             result.transformation = this.convert(transformationNode);
           } else {
@@ -503,25 +503,25 @@
           }
           break;
 
-        case "ABRIDGEDTRANSFORMATION":
-          result.type = "Transformation";
+        case 'ABRIDGEDTRANSFORMATION':
+          result.type = 'Transformation';
           result.name = node[1];
-          result.method = node.find((child) => Array.isArray(child) && child[0] === "METHOD")
-            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "METHOD"))
+          result.method = node.find((child) => Array.isArray(child) && child[0] === 'METHOD')
+            ? this.convert(node.find((child) => Array.isArray(child) && child[0] === 'METHOD'))
             : null;
 
           result.parameters = node
-            .filter((child) => Array.isArray(child) && (child[0] === "PARAMETER" || child[0] === "PARAMETERFILE"))
+            .filter((child) => Array.isArray(child) && (child[0] === 'PARAMETER' || child[0] === 'PARAMETERFILE'))
             .map((param) => {
-              if (param[0] === "PARAMETER") {
+              if (param[0] === 'PARAMETER') {
                 return this.convert(param);
-              } else if (param[0] === "PARAMETERFILE") {
+              } else if (param[0] === 'PARAMETERFILE') {
                 return {
                   name: param[1],
                   value: param[2],
                   id: {
-                    "authority": "EPSG",
-                    "code": 8656
+                    'authority': 'EPSG',
+                    'code': 8656
                   }
                 };
               }
@@ -530,7 +530,7 @@
           // Adjust the Scale difference parameter if present
           if (result.parameters.length === 7) {
             const scaleDifference = result.parameters[6];
-            if (scaleDifference.name === "Scale difference") {
+            if (scaleDifference.name === 'Scale difference') {
               scaleDifference.value = Math.round((scaleDifference.value - 1) * 1e12) / 1e6;
             }
           }
@@ -538,14 +538,14 @@
           result.id = this.getId(node);
           break;
         
-        case "AXIS":
+        case 'AXIS':
           if (!result.coordinate_system) {
-            result.coordinate_system = { type: "unspecified", axis: [] };
+            result.coordinate_system = { type: 'unspecified', axis: [] };
           }
           result.coordinate_system.axis.push(this.convertAxis(node));
           break;
         
-        case "LENGTHUNIT":
+        case 'LENGTHUNIT':
           const unit = this.convertUnit(node, 'LinearUnit');
           if (result.coordinate_system && result.coordinate_system.axis) {
             result.coordinate_system.axis.forEach((axis) => {
@@ -578,7 +578,7 @@
       super.convert(node, result);
 
       // Skip `CS` and `USAGE` nodes for WKT2-2015
-      if (result.coordinate_system?.subtype === "Cartesian") {
+      if (result.coordinate_system && result.coordinate_system.subtype === 'Cartesian') {
         delete result.coordinate_system;
       }
       if (result.usage) {
@@ -594,7 +594,7 @@
       super.convert(node, result);
 
       // Handle `CS` node for WKT2-2019
-      const csNode = node.find((child) => Array.isArray(child) && child[0] === "CS");
+      const csNode = node.find((child) => Array.isArray(child) && child[0] === 'CS');
       if (csNode) {
         result.coordinate_system = {
           subtype: csNode[1],
@@ -603,13 +603,21 @@
       }
 
       // Handle `USAGE` node for WKT2-2019
-      const usageNode = node.find((child) => Array.isArray(child) && child[0] === "USAGE");
+      const usageNode = node.find((child) => Array.isArray(child) && child[0] === 'USAGE');
       if (usageNode) {
-        result.usage = {
-          scope: usageNode.find((child) => Array.isArray(child) && child[0] === "SCOPE")?.[1],
-          area: usageNode.find((child) => Array.isArray(child) && child[0] === "AREA")?.[1],
-          bbox: usageNode.find((child) => Array.isArray(child) && child[0] === "BBOX")?.slice(1),
-        };
+        const scope = usageNode.find((child) => Array.isArray(child) && child[0] === 'SCOPE');
+        const area = usageNode.find((child) => Array.isArray(child) && child[0] === 'AREA');
+        const bbox = usageNode.find((child) => Array.isArray(child) && child[0] === 'BBOX');
+        result.usage = {};
+        if (scope) {
+          result.usage.scope = scope[1];
+        }
+        if (area) {
+          result.usage.area = area[1];
+        }
+        if (bbox) {
+          result.usage.bbox = bbox.slice(1);
+        }
       }
 
       return result;
@@ -623,21 +631,21 @@
    */
   function detectWKT2Version(root) {
     // Check for WKT2-2019-specific nodes
-    if (root.find((child) => Array.isArray(child) && child[0] === "USAGE")) {
-      return "2019"; // `USAGE` is specific to WKT2-2019
+    if (root.find((child) => Array.isArray(child) && child[0] === 'USAGE')) {
+      return '2019'; // `USAGE` is specific to WKT2-2019
     }
 
     // Check for WKT2-2015-specific nodes
-    if (root.find((child) => Array.isArray(child) && child[0] === "CS")) {
-      return "2015"; // `CS` is valid in both, but default to 2015 unless `USAGE` is present
+    if (root.find((child) => Array.isArray(child) && child[0] === 'CS')) {
+      return '2015'; // `CS` is valid in both, but default to 2015 unless `USAGE` is present
     }
 
-    if (root[0] === "BOUNDCRS" || root[0] === "PROJCRS" || root[0] === "GEOGCRS") {
-      return "2015"; // These are valid in both, but default to 2015
+    if (root[0] === 'BOUNDCRS' || root[0] === 'PROJCRS' || root[0] === 'GEOGCRS') {
+      return '2015'; // These are valid in both, but default to 2015
     }
 
     // Default to WKT2-2015 if no specific indicators are found
-    return "2015";
+    return '2015';
   }
 
   /**
@@ -647,7 +655,7 @@
    */
   function buildPROJJSON(root) {
     const version = detectWKT2Version(root);
-    const builder = version === "2019" ? PROJJSONBuilder2019 : PROJJSONBuilder2015;
+    const builder = version === '2019' ? PROJJSONBuilder2019 : PROJJSONBuilder2015;
     return builder.convert(root);
   }
 
@@ -662,30 +670,30 @@
 
     // Check for WKT2-specific keywords
     if (
-      normalizedWKT.includes("PROJCRS") ||
-      normalizedWKT.includes("GEOGCRS") ||
-      normalizedWKT.includes("BOUNDCRS") ||
-      normalizedWKT.includes("VERTCRS") ||
-      normalizedWKT.includes("LENGTHUNIT") ||
-      normalizedWKT.includes("ANGLEUNIT") ||
-      normalizedWKT.includes("SCALEUNIT")
+      normalizedWKT.includes('PROJCRS') ||
+      normalizedWKT.includes('GEOGCRS') ||
+      normalizedWKT.includes('BOUNDCRS') ||
+      normalizedWKT.includes('VERTCRS') ||
+      normalizedWKT.includes('LENGTHUNIT') ||
+      normalizedWKT.includes('ANGLEUNIT') ||
+      normalizedWKT.includes('SCALEUNIT')
     ) {
-      return "WKT2";
+      return 'WKT2';
     }
 
     // Check for WKT1-specific keywords
     if (
-      normalizedWKT.includes("PROJCS") ||
-      normalizedWKT.includes("GEOGCS") ||
-      normalizedWKT.includes("LOCAL_CS") ||
-      normalizedWKT.includes("VERT_CS") ||
-      normalizedWKT.includes("UNIT")
+      normalizedWKT.includes('PROJCS') ||
+      normalizedWKT.includes('GEOGCS') ||
+      normalizedWKT.includes('LOCAL_CS') ||
+      normalizedWKT.includes('VERT_CS') ||
+      normalizedWKT.includes('UNIT')
     ) {
-      return "WKT1";
+      return 'WKT1';
     }
 
     // Default to WKT1 if no specific indicators are found
-    return "WKT1";
+    return 'WKT1';
   }
 
   var NEUTRAL = 1;
@@ -1010,7 +1018,7 @@
       if (result.units === 'meter') {
         result.to_meter = 1; // Only set to_meter if units are 'meter'
       }
-    } else if (unit?.name) {
+    } else if (unit && unit.name) {
       result.units = unit.name.toLowerCase();
       if (result.units === 'metre') {
         result.units = 'meter'; // Normalize 'metre' to 'meter'
@@ -1051,7 +1059,7 @@
       transformPROJJSON(projjson.source_crs, result);
 
       if (projjson.transformation) {
-        if (projjson.transformation.method?.name === 'NTv2') {
+        if (projjson.transformation.method && projjson.transformation.method.name === 'NTv2') {
           // Set nadgrids to the filename from the parameterfile
           result.nadgrids = projjson.transformation.parameters[0].value;
         } else {
@@ -1081,8 +1089,8 @@
         case 'type':
           if (value === 'GeographicCRS') {
             result.projName = 'longlat';
-          } else if (value === 'ProjectedCRS') {
-            result.projName = projjson.conversion?.method?.name; // Retain original capitalization
+          } else if (value === 'ProjectedCRS' && projjson.conversion && projjson.conversion.method) {
+            result.projName = projjson.conversion.method.name; // Retain original capitalization
           }
           break;
 
@@ -1124,7 +1132,7 @@
               const { units, to_meter } = processUnit(value.unit);
               result.units = units;
               result.to_meter = to_meter;
-            } else if (value.axis[0]?.unit) {
+            } else if (value.axis[0] && value.axis[0].unit) {
               const { units, to_meter } = processUnit(value.axis[0].unit);
               result.units = units;
               result.to_meter = to_meter;
@@ -1590,6 +1598,9 @@
   }
   var codes = ['3857', '900913', '3785', '102113'];
   function checkMercator(item) {
+    if (item.title) {
+      return item.title.toLowerCase().indexOf('epsg:') === 0 && codes.indexOf(item.title.substr(5)) > -1;
+    }
     var auth = match(item, 'authority');
     if (!auth) {
       return;
@@ -1612,31 +1623,27 @@
    * @returns {import('./defs').ProjectionDefinition}
    */
   function parse(code) {
+    let out;
     if (testObj(code)) {
       // check to see if this is a WKT string
       if (testDef(code)) {
-        return defs[code];
-      }
-      if (testWKT(code)) {
-        var out = wkt(code);
-        // test of spetial case, due to this being a very common and often malformed
-        if (checkMercator(out)) {
-          return defs['EPSG:3857'];
-        }
+        out = defs[code];
+      } else if (testWKT(code)) {
+        out = wkt(code);
         var maybeProjStr = checkProjStr(out);
         if (maybeProjStr) {
-          return projStr(maybeProjStr);
+          out = projStr(maybeProjStr);
         }
-        return out;
-      }
-      if (testProj(code)) {
-        return projStr(code);
+      } else if (testProj(code)) {
+        out = projStr(code);
       }
     } else if (!('projName' in code)) {
-      return wkt(code);
+      out = wkt(code);
     } else {
-      return code;
+      out = code;
     }
+    // test for special Web Mercator case, due to this being a very common and often malformed
+    return out && checkMercator(out) ? defs['EPSG:3857'] : out;
   }
 
   function extend (destination, source) {
@@ -2260,7 +2267,7 @@
       towgs84: '-178.3,-316.7,-131.5,5.278,6.077,10.979,19.166'
     },
     EPSG_4231: {
-      towgs84: '-83.11,-97.38,-117.22,0.0276,-0.2167,0.2147,0.1218'
+      towgs84: '-83.11,-97.38,-117.22,0.005693,-0.044698,0.044285,0.1218'
     },
     EPSG_4274: {
       towgs84: '-230.994,102.591,25.199,0.633,-0.239,0.9,1.95'
@@ -2290,7 +2297,7 @@
       towgs84: '8.846,-4.394,-1.122,-0.00237,-0.146528,0.130428,0.783926'
     },
     EPSG_4289: {
-      towgs84: '565.7381,50.4018,465.2904,-1.91514,1.60363,-9.09546,4.07244'
+      towgs84: '565.7381,50.4018,465.2904,-0.395026,0.330772,-1.876073,4.07244'
     },
     EPSG_4230: {
       towgs84: '-68.863,-134.888,-111.49,-0.53,-0.14,0.57,-3.4'
@@ -2320,7 +2327,7 @@
       towgs84: '-45,417,-3.5'
     },
     EPSG_4611: {
-      towgs84: '-162.619,-276.959,-161.764,0.067753,-2.243649,-1.158827,-1.094246'
+      towgs84: '-162.619,-276.959,-161.764,0.067753,-2.243648,-1.158828,-1.094246'
     },
     EPSG_4633: {
       towgs84: '137.092,131.66,91.475,-1.9436,-11.5993,-4.3321,-7.4824'
@@ -2338,7 +2345,7 @@
       towgs84: '482.5,-130.6,564.6,-1.042,-0.214,-0.631,8.15'
     },
     EPSG_4660: {
-      towgs84: '982.6087,552.753,-540.873,32.39344,-153.25684,-96.2266,16.805'
+      towgs84: '982.6087,552.753,-540.873,6.681627,-31.611492,-19.848161,16.805'
     },
     EPSG_4662: {
       towgs84: '97.295,-263.247,310.882,-1.5999,0.8386,3.1409,13.3259'
@@ -2368,7 +2375,7 @@
       towgs84: '217.109,86.452,23.711,0.0183,-0.0003,0.007,-0.0093'
     },
     EPSG_9333: {
-      towgs84: '0,0,0,-8.393,0.749,-10.276,0'
+      towgs84: '0,0,0,-0.008393,0.000749,-0.010276,0'
     },
     EPSG_9059: {
       towgs84: '0,0,0'
@@ -2524,7 +2531,7 @@
       towgs84: '674.374,15.056,405.346'
     },
     EPSG_4617: {
-      towgs84: '-0.991,1.9072,0.5129,1.25033e-7,4.6785e-8,5.6529e-8,0'
+      towgs84: '-0.991,1.9072,0.5129,0.02579,0.00965,0.01166,0'
     },
     EPSG_4663: {
       towgs84: '-210.502,-66.902,-48.476,2.094,-15.067,-5.817,0.485'
@@ -2587,10 +2594,10 @@
       towgs84: '275.57,676.78,229.6'
     },
     EPSG_4283: {
-      towgs84: '61.55,-10.87,-40.19,39.4924,32.7221,32.8979,-9.994'
+      towgs84: '0.06155,-0.01087,-0.04019,0.039492,0.032722,0.032898,-0.009994'
     },
     EPSG_4317: {
-      towgs84: '2.3287,-147.0425,-92.0802,-0.3092483,0.32482185,0.49729934,5.68906266'
+      towgs84: '2.3287,-147.0425,-92.0802,-0.309248,0.324822,0.497299,5.689063'
     },
     EPSG_4272: {
       towgs84: '59.47,-5.04,187.44,0.47,-0.1,1.024,-4.5993'
@@ -2713,7 +2720,7 @@
       towgs84: '-208.4058,-109.8777,-2.5764'
     },
     ESRI_104101: {
-      towgs84: '374,150,588'
+      towgs84: '372.87,149.23,585.29'
     },
     EPSG_4693: {
       towgs84: '0,-0.15,0.68'
@@ -3398,6 +3405,87 @@
     },
     EPSG_4818: {
       towgs84: '589,76,480'
+    },
+    EPSG_10328: {
+      towgs84: '0,0,0'
+    },
+    EPSG_9782: {
+      towgs84: '0,0,0'
+    },
+    EPSG_9777: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10690: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10639: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10739: {
+      towgs84: '0,0,0'
+    },
+    EPSG_7686: {
+      towgs84: '0,0,0'
+    },
+    EPSG_8900: {
+      towgs84: '0,0,0'
+    },
+    EPSG_5886: {
+      towgs84: '0,0,0'
+    },
+    EPSG_7683: {
+      towgs84: '0,0,0'
+    },
+    EPSG_6668: {
+      towgs84: '0,0,0'
+    },
+    EPSG_20046: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10299: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10310: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10475: {
+      towgs84: '0,0,0'
+    },
+    EPSG_4742: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10671: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10762: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10725: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10791: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10800: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10305: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10941: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10968: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10875: {
+      towgs84: '0,0,0'
+    },
+    EPSG_6318: {
+      towgs84: '0,0,0'
+    },
+    EPSG_10910: {
+      towgs84: '0,0,0'
     }
   };
 
@@ -6916,7 +7004,7 @@
     return (p);
   }
 
-  var names$n = ['Krovak', 'krovak'];
+  var names$n = ['Krovak', 'Krovak Modified', 'Krovak (North Orientated)', 'Krovak Modified (North Orientated)', 'krovak'];
   var krovak = {
     init: init$n,
     forward: forward$m,
@@ -8910,8 +8998,8 @@
 
     /* Place parameters in static storage for common use
         ------------------------------------------------- */
-    this.sin_p14 = Math.sin(this.lat0);
-    this.cos_p14 = Math.cos(this.lat0);
+    this.sin_p14 = Math.sin(this.lat0 || 0);
+    this.cos_p14 = Math.cos(this.lat0 || 0);
   }
 
   /* Orthographic forward equations--mapping lat,long to x,y
@@ -8926,7 +9014,7 @@
     var lat = p.y;
     /* Forward equations
         ----------------- */
-    dlon = adjust_lon(lon - this.long0, this.over);
+    dlon = adjust_lon(lon - (this.long0 || 0), this.over);
 
     sinphi = Math.sin(lat);
     cosphi = Math.cos(lat);
@@ -8936,7 +9024,7 @@
     ksp = 1;
     if ((g > 0) || (Math.abs(g) <= EPSLN)) {
       x = this.a * ksp * cosphi * Math.sin(dlon);
-      y = this.y0 + this.a * ksp * (this.cos_p14 * sinphi - this.sin_p14 * cosphi * coslon);
+      y = (this.y0 || 0) + this.a * ksp * (this.cos_p14 * sinphi - this.sin_p14 * cosphi * coslon);
     }
     p.x = x;
     p.y = y;
@@ -8949,36 +9037,39 @@
     var sinz, cosz; /* sin of z and cos of z      */
     var con;
     var lon, lat;
+    var long0, lat0;
     /* Inverse equations
         ----------------- */
-    p.x -= this.x0;
-    p.y -= this.y0;
+    p.x -= this.x0 || 0;
+    p.y -= this.y0 || 0;
     rh = Math.sqrt(p.x * p.x + p.y * p.y);
     z = asinz(rh / this.a);
 
     sinz = Math.sin(z);
     cosz = Math.cos(z);
 
-    lon = this.long0;
+    long0 = this.long0 || 0;
+    lat0 = this.lat0 || 0;
+    lon = long0;
     if (Math.abs(rh) <= EPSLN) {
-      lat = this.lat0;
+      lat = lat0;
       p.x = lon;
       p.y = lat;
       return p;
     }
     lat = asinz(cosz * this.sin_p14 + (p.y * sinz * this.cos_p14) / rh);
-    con = Math.abs(this.lat0) - HALF_PI;
+    con = Math.abs(lat0) - HALF_PI;
     if (Math.abs(con) <= EPSLN) {
-      if (this.lat0 >= 0) {
-        lon = adjust_lon(this.long0 + Math.atan2(p.x, -p.y), this.over);
+      if (lat0 >= 0) {
+        lon = adjust_lon(long0 + Math.atan2(p.x, -p.y), this.over);
       } else {
-        lon = adjust_lon(this.long0 - Math.atan2(-p.x, p.y), this.over);
+        lon = adjust_lon(long0 - Math.atan2(-p.x, p.y), this.over);
       }
       p.x = lon;
       p.y = lat;
       return p;
     }
-    lon = adjust_lon(this.long0 + Math.atan2((p.x * sinz), rh * this.cos_p14 * cosz - p.y * this.sin_p14 * sinz), this.over);
+    lon = adjust_lon(long0 + Math.atan2((p.x * sinz), rh * this.cos_p14 * cosz - p.y * this.sin_p14 * sinz), this.over);
     p.x = lon;
     p.y = lat;
     return p;
@@ -9119,7 +9210,7 @@
       if (this.face === FACE_ENUM.RIGHT) {
         lon = qsc_shift_lon_origin(lon, +HALF_PI);
       } else if (this.face === FACE_ENUM.BACK) {
-        lon = qsc_shift_lon_origin(lon, 3.14159265359);
+        lon = qsc_shift_lon_origin(lon, +SPI);
       } else if (this.face === FACE_ENUM.LEFT) {
         lon = qsc_shift_lon_origin(lon, -HALF_PI);
       }
@@ -9299,7 +9390,7 @@
       if (this.face === FACE_ENUM.RIGHT) {
         lp.lam = qsc_shift_lon_origin(lp.lam, -HALF_PI);
       } else if (this.face === FACE_ENUM.BACK) {
-        lp.lam = qsc_shift_lon_origin(lp.lam, -3.14159265359);
+        lp.lam = qsc_shift_lon_origin(lp.lam, -SPI);
       } else if (this.face === FACE_ENUM.LEFT) {
         lp.lam = qsc_shift_lon_origin(lp.lam, +HALF_PI);
       }
@@ -9353,9 +9444,9 @@
   /* Helper function: shift the longitude. */
   function qsc_shift_lon_origin(lon, offset) {
     var slon = lon + offset;
-    if (slon < -3.14159265359) {
+    if (slon < -SPI) {
       slon += TWO_PI;
-    } else if (slon > 3.14159265359) {
+    } else if (slon > +SPI) {
       slon -= TWO_PI;
     }
     return slon;
@@ -10552,7 +10643,7 @@
     nadgrid,
     transform,
     mgrs,
-    version: '2.20.2'
+    version: '2.20.4'
   });
   includedProjections(proj4);