proj4 2.15.0 → 2.16.1

package/lib/projections/tmerc.js

@@ -6,7 +6,7 @@ import pj_mlfn from '../common/pj_mlfn';
 import pj_inv_mlfn from '../common/pj_inv_mlfn';
 import adjust_lon from '../common/adjust_lon';
 
-import {EPSLN, HALF_PI} from '../constants/values';
+import { EPSLN, HALF_PI } from '../constants/values';
 import sign from '../common/sign';
 
 export function init() {
@@ -40,8 +40,7 @@ export function forward(p) {
 
     if ((Math.abs(Math.abs(b) - 1)) < EPSLN) {
       return (93);
-    }
-    else {
+    } else {
       x = 0.5 * this.a * this.k0 * Math.log((1 + b) / (1 - b)) + this.x0;
       y = cos_phi * Math.cos(delta_lon) / Math.sqrt(1 - Math.pow(b, 2));
       b = Math.abs(y);
@@ -49,12 +48,10 @@ export function forward(p) {
       if (b >= 1) {
         if ((b - 1) > EPSLN) {
           return (93);
-        }
-        else {
+        } else {
           y = 0;
         }
-      }
-      else {
+      } else {
         y = Math.acos(y);
       }
 
@@ -64,8 +61,7 @@ export function forward(p) {
 
       y = this.a * this.k0 * (y - this.lat0) + this.y0;
     }
-  }
-  else {
+  } else {
     var al = cos_phi * delta_lon;
     var als = Math.pow(al, 2);
     var c = this.ep2 * Math.pow(cos_phi, 2);
@@ -77,18 +73,18 @@ export function forward(p) {
     al = al / Math.sqrt(con);
     var ml = pj_mlfn(lat, sin_phi, cos_phi, this.en);
 
-    x = this.a * (this.k0 * al * (1 +
-      als / 6 * (1 - t + c +
-      als / 20 * (5 - 18 * t + ts + 14 * c - 58 * t * c +
-      als / 42 * (61 + 179 * ts - ts * t - 479 * t))))) +
-      this.x0;
-
-    y = this.a * (this.k0 * (ml - this.ml0 +
-      sin_phi * delta_lon * al / 2 * (1 +
-      als / 12 * (5 - t + 9 * c + 4 * cs +
-      als / 30 * (61 + ts - 58 * t + 270 * c - 330 * t * c +
-      als / 56 * (1385 + 543 * ts - ts * t - 3111 * t)))))) +
-      this.y0;
+    x = this.a * (this.k0 * al * (1
+      + als / 6 * (1 - t + c
+        + als / 20 * (5 - 18 * t + ts + 14 * c - 58 * t * c
+          + als / 42 * (61 + 179 * ts - ts * t - 479 * t)))))
+        + this.x0;
+
+    y = this.a * (this.k0 * (ml - this.ml0
+      + sin_phi * delta_lon * al / 2 * (1
+        + als / 12 * (5 - t + 9 * c + 4 * cs
+          + als / 30 * (61 + ts - 58 * t + 270 * c - 330 * t * c
+            + als / 56 * (1385 + 543 * ts - ts * t - 3111 * t))))))
+          + this.y0;
   }
 
   p.x = x;
@@ -120,12 +116,10 @@ export function inverse(p) {
 
     if ((g === 0) && (h === 0)) {
       lon = 0;
-    }
-    else {
+    } else {
       lon = adjust_lon(Math.atan2(g, h) + this.long0);
     }
-  }
-  else { // ellipsoidal form
+  } else { // ellipsoidal form
     con = this.ml0 + y / this.k0;
     phi = pj_inv_mlfn(con, this.es, this.en);
 
@@ -142,17 +136,16 @@ export function inverse(p) {
       var ds = Math.pow(d, 2);
       con = con * tan_phi;
 
-      lat = phi - (con * ds / (1 - this.es)) * 0.5 * (1 -
-        ds / 12 * (5 + 3 * t - 9 * c * t + c - 4 * cs -
-        ds / 30 * (61 + 90 * t - 252 * c * t + 45 * ts + 46 * c -
-        ds / 56 * (1385 + 3633 * t + 4095 * ts + 1574 * ts * t))));
+      lat = phi - (con * ds / (1 - this.es)) * 0.5 * (1
+        - ds / 12 * (5 + 3 * t - 9 * c * t + c - 4 * cs
+          - ds / 30 * (61 + 90 * t - 252 * c * t + 45 * ts + 46 * c
+            - ds / 56 * (1385 + 3633 * t + 4095 * ts + 1574 * ts * t))));
 
-      lon = adjust_lon(this.long0 + (d * (1 -
-        ds / 6 * (1 + 2 * t + c -
-        ds / 20 * (5 + 28 * t + 24 * ts + 8 * c * t + 6 * c -
-        ds / 42 * (61 + 662 * t + 1320 * ts + 720 * ts * t)))) / cos_phi));
-    }
-    else {
+      lon = adjust_lon(this.long0 + (d * (1
+        - ds / 6 * (1 + 2 * t + c
+          - ds / 20 * (5 + 28 * t + 24 * ts + 8 * c * t + 6 * c
+            - ds / 42 * (61 + 662 * t + 1320 * ts + 720 * ts * t)))) / cos_phi));
+    } else {
       lat = HALF_PI * sign(y);
       lon = 0;
     }
@@ -164,7 +157,7 @@ export function inverse(p) {
   return p;
 }
 
-export var names = ["Fast_Transverse_Mercator", "Fast Transverse Mercator"];
+export var names = ['Fast_Transverse_Mercator', 'Fast Transverse Mercator'];
 export default {
   init: init,
   forward: forward,

package/lib/projections/tpers.js

@@ -1,4 +1,3 @@
-
 var mode = {
   N_POLE: 0,
   S_POLE: 1,
@@ -6,23 +5,23 @@ var mode = {
   OBLIQ: 3
 };
 
-import { D2R, HALF_PI, EPSLN } from "../constants/values";
-import hypot from "../common/hypot";
+import { D2R, HALF_PI, EPSLN } from '../constants/values';
+import hypot from '../common/hypot';
 
 var params = {
-  h:     { def: 100000, num: true },           // default is Karman line, no default in PROJ.7
-  azi:   { def: 0, num: true, degrees: true }, // default is North
-  tilt:  { def: 0, num: true, degrees: true }, // default is Nadir
-  long0: { def: 0, num: true },                // default is Greenwich, conversion to rad is automatic
-  lat0:  { def: 0, num: true }                 // default is Equator, conversion to rad is automatic
+  h: { def: 100000, num: true }, // default is Karman line, no default in PROJ.7
+  azi: { def: 0, num: true, degrees: true }, // default is North
+  tilt: { def: 0, num: true, degrees: true }, // default is Nadir
+  long0: { def: 0, num: true }, // default is Greenwich, conversion to rad is automatic
+  lat0: { def: 0, num: true } // default is Equator, conversion to rad is automatic
 };
 
 export function init() {
   Object.keys(params).forEach(function (p) {
-    if (typeof this[p] === "undefined") {
+    if (typeof this[p] === 'undefined') {
       this[p] = params[p].def;
     } else if (params[p].num && isNaN(this[p])) {
-      throw new Error("Invalid parameter value, must be numeric " + p + " = " + this[p]);
+      throw new Error('Invalid parameter value, must be numeric ' + p + ' = ' + this[p]);
     } else if (params[p].num) {
       this[p] = parseFloat(this[p]);
     }
@@ -41,12 +40,12 @@ export function init() {
     this.cosph0 = Math.cos(this.lat0);
   }
 
-  this.pn1 = this.h / this.a;  // Normalize relative to the Earth's radius
+  this.pn1 = this.h / this.a; // Normalize relative to the Earth's radius
 
   if (this.pn1 <= 0 || this.pn1 > 1e10) {
-    throw new Error("Invalid height");
+    throw new Error('Invalid height');
   }
-  
+
   this.p = 1 + this.pn1;
   this.rp = 1 / this.p;
   this.h1 = 1 / this.pn1;
@@ -99,7 +98,7 @@ export function forward(p) {
       break;
   }
 
-  // Tilt 
+  // Tilt
   var yt, ba;
   yt = y * this.cg + x * this.sg;
   ba = 1 / (yt * this.sw * this.h1 + this.cw);
@@ -160,7 +159,7 @@ export function inverse(p) {
   return p;
 }
 
-export var names = ["Tilted_Perspective", "tpers"];
+export var names = ['Tilted_Perspective', 'tpers'];
 export default {
   init: init,
   forward: forward,

package/lib/projections/utm.js

@@ -1,8 +1,7 @@
 import adjust_zone from '../common/adjust_zone';
 import etmerc from './etmerc';
 export var dependsOn = 'etmerc';
-import {D2R} from '../constants/values';
-
+import { D2R } from '../constants/values';
 
 export function init() {
   var zone = adjust_zone(this.zone, this.long0);
@@ -10,7 +9,7 @@ export function init() {
     throw new Error('unknown utm zone');
   }
   this.lat0 = 0;
-  this.long0 =  ((6 * Math.abs(zone)) - 183) * D2R;
+  this.long0 = ((6 * Math.abs(zone)) - 183) * D2R;
   this.x0 = 500000;
   this.y0 = this.utmSouth ? 10000000 : 0;
   this.k0 = 0.9996;
@@ -20,7 +19,7 @@ export function init() {
   this.inverse = etmerc.inverse;
 }
 
-export var names = ["Universal Transverse Mercator System", "utm"];
+export var names = ['Universal Transverse Mercator System', 'utm'];
 export default {
   init: init,
   names: names,

package/lib/projections/vandg.js

@@ -1,23 +1,22 @@
 import adjust_lon from '../common/adjust_lon';
 
-import {HALF_PI, EPSLN} from '../constants/values';
+import { HALF_PI, EPSLN } from '../constants/values';
 
 import asinz from '../common/asinz';
 
 /* Initialize the Van Der Grinten projection
-  ----------------------------------------*/
+  ---------------------------------------- */
 export function init() {
-  //this.R = 6370997; //Radius of earth
+  // this.R = 6370997; //Radius of earth
   this.R = this.a;
 }
 
 export function forward(p) {
-
   var lon = p.x;
   var lat = p.y;
 
   /* Forward equations
-    -----------------*/
+    ----------------- */
   var dlon = adjust_lon(lon - this.long0);
   var x, y;
 
@@ -30,8 +29,7 @@ export function forward(p) {
     x = this.x0;
     if (lat >= 0) {
       y = this.y0 + Math.PI * this.R * Math.tan(0.5 * theta);
-    }
-    else {
+    } else {
       y = this.y0 + Math.PI * this.R * -Math.tan(0.5 * theta);
     }
     //  return(OK);
@@ -50,15 +48,14 @@ export function forward(p) {
     con = -con;
   }
   x = this.x0 + con;
-  //con = Math.abs(con / (Math.PI * this.R));
+  // con = Math.abs(con / (Math.PI * this.R));
   var q = asq + g;
   con = Math.PI * this.R * (m * q - al * Math.sqrt((msq + asq) * (asq + 1) - q * q)) / (msq + asq);
   if (lat >= 0) {
-    //y = this.y0 + Math.PI * this.R * Math.sqrt(1 - con * con - 2 * al * con);
+    // y = this.y0 + Math.PI * this.R * Math.sqrt(1 - con * con - 2 * al * con);
     y = this.y0 + con;
-  }
-  else {
-    //y = this.y0 - Math.PI * this.R * Math.sqrt(1 - con * con - 2 * al * con);
+  } else {
+    // y = this.y0 - Math.PI * this.R * Math.sqrt(1 - con * con - 2 * al * con);
     y = this.y0 - con;
   }
   p.x = x;
@@ -67,7 +64,7 @@ export function forward(p) {
 }
 
 /* Van Der Grinten inverse equations--mapping x,y to lat/long
-  ---------------------------------------------------------*/
+  --------------------------------------------------------- */
 export function inverse(p) {
   var lon, lat;
   var xx, yy, xys, c1, c2, c3;
@@ -78,7 +75,7 @@ export function inverse(p) {
   var d;
 
   /* inverse equations
-    -----------------*/
+    ----------------- */
   p.x -= this.x0;
   p.y -= this.y0;
   con = Math.PI * this.R;
@@ -95,23 +92,20 @@ export function inverse(p) {
   if (Math.abs(con) > 1) {
     if (con >= 0) {
       con = 1;
-    }
-    else {
+    } else {
       con = -1;
     }
   }
   th1 = Math.acos(con) / 3;
   if (p.y >= 0) {
     lat = (-m1 * Math.cos(th1 + Math.PI / 3) - c2 / 3 / c3) * Math.PI;
-  }
-  else {
+  } else {
     lat = -(-m1 * Math.cos(th1 + Math.PI / 3) - c2 / 3 / c3) * Math.PI;
   }
 
   if (Math.abs(xx) < EPSLN) {
     lon = this.long0;
-  }
-  else {
+  } else {
     lon = adjust_lon(this.long0 + Math.PI * (xys - 1 + Math.sqrt(1 + 2 * (xx * xx - yy * yy) + xys * xys)) / 2 / xx);
   }
 
@@ -120,7 +114,7 @@ export function inverse(p) {
   return p;
 }
 
-export var names = ["Van_der_Grinten_I", "VanDerGrinten", "vandg"];
+export var names = ['Van_der_Grinten_I', 'VanDerGrinten', 'Van_der_Grinten', 'vandg'];
 export default {
   init: init,
   forward: forward,

package/lib/projections.js

@@ -1,5 +1,5 @@
-import merc from "./projections/merc";
-import longlat from "./projections/longlat";
+import merc from './projections/merc';
+import longlat from './projections/longlat';
 var projs = [merc, longlat];
 var names = {};
 var projStore = [];
@@ -11,13 +11,17 @@ function add(proj, i) {
     return true;
   }
   projStore[len] = proj;
-  proj.names.forEach(function(n) {
+  proj.names.forEach(function (n) {
     names[n.toLowerCase()] = len;
   });
   return this;
 }
 
-export {add};
+export { add };
+
+export function getNormalizedProjName(n) {
+  return n.replace(/[-\(\)\s]+/g, ' ').trim().replace(/ /g, '_');
+}
 
 export function get(name) {
   if (!name) {
@@ -27,6 +31,10 @@ export function get(name) {
   if (typeof names[n] !== 'undefined' && projStore[names[n]]) {
     return projStore[names[n]];
   }
+  n = getNormalizedProjName(n);
+  if (n in names && projStore[names[n]]) {
+    return projStore[names[n]];
+  }
 }
 
 export function start() {

package/lib/transform.js

@@ -1,4 +1,4 @@
-import {D2R, R2D, PJD_3PARAM, PJD_7PARAM, PJD_GRIDSHIFT} from './constants/values';
+import { D2R, R2D, PJD_3PARAM, PJD_7PARAM, PJD_GRIDSHIFT } from './constants/values';
 import datum_transform from './datum_transform';
 import adjust_axis from './adjust_axis';
 import proj from './Proj';
@@ -7,8 +7,8 @@ import checkSanity from './checkSanity';
 
 function checkNotWGS(source, dest) {
   return (
-    (source.datum.datum_type === PJD_3PARAM || source.datum.datum_type === PJD_7PARAM || source.datum.datum_type === PJD_GRIDSHIFT) && dest.datumCode !== 'WGS84') ||
-    ((dest.datum.datum_type === PJD_3PARAM || dest.datum.datum_type === PJD_7PARAM || dest.datum.datum_type === PJD_GRIDSHIFT) && source.datumCode !== 'WGS84');
+    (source.datum.datum_type === PJD_3PARAM || source.datum.datum_type === PJD_7PARAM || source.datum.datum_type === PJD_GRIDSHIFT) && dest.datumCode !== 'WGS84')
+  || ((dest.datum.datum_type === PJD_3PARAM || dest.datum.datum_type === PJD_7PARAM || dest.datum.datum_type === PJD_GRIDSHIFT) && source.datumCode !== 'WGS84');
 }
 
 export default function transform(source, dest, point, enforceAxis) {

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "proj4",
-  "version": "2.15.0",
+  "version": "2.16.1",
   "description": "Proj4js is a JavaScript library to transform point coordinates from one coordinate system to another, including datum transformations.",
   "homepage": "https://proj4js.github.io/proj4js/",
   "main": "dist/proj4-src.js",
@@ -10,10 +10,18 @@
     "doc": "docs"
   },
   "scripts": {
-    "prepare": "grunt",
-    "build": "grunt",
-    "build:tmerc": "grunt build:tmerc",
-    "test": "npm run build && istanbul test node_modules/mocha/bin/_mocha test/test.js"
+    "prepare": "npm run build && npm run test:all",
+    "build": "grunt && npm run rollup",
+    "build:tmerc": "grunt build:tmerc && npm run rollup",
+    "rollup": "rollup -c",
+    "pretest": "npm run lint",
+    "lint": "eslint *.js *.mjs scripts lib test",
+    "format": "npm run lint -- --fix",
+    "test": "npm run build && npm run test:all",
+    "test:coverage": "nyc npm run test:ci",
+    "test:all": "npm run test:ci && npm run test:browser",
+    "test:browser": "node test/puppeteer-tests.mjs",
+    "test:ci": "npx -y mocha test/test-ci.mjs --reporter dot"
   },
   "repository": {
     "type": "git",
@@ -22,25 +30,24 @@
   "author": "",
   "license": "MIT",
   "devDependencies": {
-    "chai": "~4.1.2",
-    "curl-amd": "^0.8.12",
-    "grunt": "^1.0.1",
-    "grunt-cli": "~1.2.0",
-    "grunt-contrib-connect": "~1.0.2",
-    "grunt-contrib-jshint": "~3.2.0",
-    "grunt-contrib-uglify": "~3.1.0",
-    "grunt-mocha-phantomjs": "~4.0.0",
-    "grunt-rollup": "^6.0.0",
-    "istanbul": "~0.4.5",
-    "mocha": "~4.0.0",
-    "rollup": "^0.50.0",
-    "rollup-plugin-json": "^2.3.0",
-    "rollup-plugin-node-resolve": "^3.0.0",
-    "rollup-plugin-replace": "^2.0.0",
-    "tin": "~0.5.0"
+    "@rollup/plugin-commonjs": "^28.0.3",
+    "@rollup/plugin-json": "^6.1.0",
+    "@rollup/plugin-node-resolve": "^15.2.3",
+    "@rollup/plugin-replace": "^5.0.5",
+    "@rollup/plugin-terser": "^0.4.4",
+    "@stylistic/eslint-plugin": "^4.2.0",
+    "chai": "^5.0.0",
+    "eslint": "^9.25.1",
+    "grunt": "^1.6.1",
+    "grunt-cli": "^1.4.3",
+    "mocha": "^10.2.0",
+    "nyc": "^17.1.0",
+    "puppeteer": "^24.2.0",
+    "rollup": "^4.9.6"
   },
   "dependencies": {
+    "geographiclib-geodesic": "^2.1.1",
     "mgrs": "1.0.0",
-    "wkt-parser": "^1.4.0"
+    "wkt-parser": "^1.5.0"
   }
 }

package/projs.js

@@ -29,7 +29,7 @@ import tpers from './lib/projections/tpers';
 import geos from './lib/projections/geos';
 import eqearth from './lib/projections/eqearth';
 import bonne from './lib/projections/bonne';
-export default function(proj4){
+export default function (proj4) {
   proj4.Proj.projections.add(tmerc);
   proj4.Proj.projections.add(etmerc);
   proj4.Proj.projections.add(utm);
@@ -61,4 +61,4 @@ export default function(proj4){
   proj4.Proj.projections.add(geos);
   proj4.Proj.projections.add(eqearth);
   proj4.Proj.projections.add(bonne);
-}
+}

package/publish.sh

@@ -5,13 +5,13 @@ VERSION=$(node -e "console.log(require('./package.json').version)")
 
 # Build
 git checkout -b build
-node_modules/.bin/grunt
+npm run build
 git add dist -f
 git commit -m "build $VERSION"
 
 # Tag and push
 git tag -f v$VERSION -m "$VERSION"
-git push --tags git@github.com:proj4js/proj4js.git $VERSION
+git push --tags git@github.com:proj4js/proj4js.git v$VERSION
 
 # Publish
 npm publish

package/rollup.config.mjs

@@ -0,0 +1,29 @@
+import json from '@rollup/plugin-json';
+import nodeResolve from '@rollup/plugin-node-resolve';
+import terser from '@rollup/plugin-terser';
+import replace from '@rollup/plugin-replace';
+import commonjs from '@rollup/plugin-commonjs';
+
+export default {
+  input: './lib/index.js',
+  output: [{
+    file: './dist/proj4-src.js',
+    format: 'umd',
+    name: 'proj4'
+  },
+  {
+    file: './dist/proj4.js',
+    format: 'umd',
+    name: 'proj4',
+    plugins: [terser()]
+  }],
+  plugins: [
+    replace({
+      preventAssignment: true,
+      __VERSION__: process.env.npm_package_version
+    }),
+    json(),
+    nodeResolve(),
+    commonjs()
+  ]
+};

package/scripts/readme.md

@@ -0,0 +1,10 @@
+The script in this directory can be used to update `lib/constants/Datum.js` with data from a [`proj.db`](https://proj.org/en/stable/resource_files.html#proj-db) file placed in the root of the repository.
+
+After running the script, the formatting needs to be fixed.
+
+A complete command sequence (issued from the root of the repository) would look like this:
+
+    npm install sqlite3
+    node scripts/updateDatums.mjs
+    npm uninstall sqlite3
+    npm run format

package/scripts/updateDatums.mjs

@@ -0,0 +1,119 @@
+import sqlite3 from 'sqlite3';
+import fs from 'fs';
+import path from 'path';
+import { fileURLToPath } from 'url';
+import datums from '../lib/constants/Datum.js'; // Import the datums object directly
+
+const DATUM_OVERRIDES = {
+  EPSG_4149: {
+    towgs84: '674.374,15.056,405.346'
+  },
+  EPSG_4267: {
+    towgs84: '-8.0,160.0,176.0'
+  },
+  EPSG_4818: {
+    towgs84: '589,76,480'
+  }
+};
+
+// Get the current file's directory
+const __filename = fileURLToPath(import.meta.url);
+const __dirname = path.dirname(__filename);
+
+// Open proj.db
+const dbPath = path.resolve(__dirname, '../proj.db');
+const db = new sqlite3.Database(dbPath);
+
+for (const key in datums) {
+  if (key === datums[key].datumName && datums[datums[key].datumName] === datums[key]) {
+    delete datums[key];
+  }
+}
+
+const databaseDatumNames = new Set();
+
+// Query proj.db for Helmert transforms with method_code = 9606
+db.all(
+  `SELECT ht.source_crs_auth_name, 
+          ht.source_crs_code, 
+          ht.tx, ht.ty, ht.tz, 
+          ht.rx, ht.ry, ht.rz, 
+          ht.scale_difference,
+          ht.accuracy, 
+          gd.name AS datum_name,
+          cv.name AS datum_code,
+          e.name AS ellipse_name,
+          ht.method_code
+   FROM helmert_transformation ht
+   JOIN crs_view cv 
+     ON ht.source_crs_auth_name = cv.auth_name 
+    AND ht.source_crs_code = cv.code
+   JOIN geodetic_crs gcrs
+     ON cv.table_name = 'geodetic_crs'
+    AND cv.auth_name = gcrs.auth_name
+    AND cv.code = gcrs.code
+   JOIN geodetic_datum gd
+     ON gcrs.datum_auth_name = gd.auth_name
+    AND gcrs.datum_code = gd.code
+   JOIN ellipsoid e
+     ON gd.ellipsoid_auth_name = e.auth_name
+    AND gd.ellipsoid_code = e.code
+   WHERE ht.deprecated = 0 
+     AND ht.method_code IN (9606, 9607, 9603)
+     AND (ht.target_crs_auth_name = 'EPSG' AND ht.target_crs_code IN ('4326', '7019', '4258')) -- WGS84, GRS80, ETRS89
+     AND cv.type != 'vertical' -- Exclude vertical datums
+   ORDER BY
+     ht.accuracy ASC,
+     CASE ht.method_code 
+       WHEN 9606 THEN 1
+       WHEN 9607 THEN 2
+       WHEN 9603 THEN 3
+       ELSE 4
+     END ASC`, // Sort by method_code (9606 first, then 9607), then by accuracy`, // Assuming lower accuracy values are better
+  (err, rows) => {
+    if (err) {
+      console.error('Error querying proj.db:', err);
+      return;
+    }
+
+    rows.forEach((row) => {
+      const normalizedDatumName = row.datum_name.toLowerCase();
+
+      // Skip if the datumName already exists (case-insensitive)
+      if (databaseDatumNames.has(normalizedDatumName)) {
+        return;
+      }
+
+      if (row.method_code === 9607 && row.rx) {
+        row.rx = -row.rx;
+        row.ry = -row.ry;
+        row.rz = -row.rz;
+      }
+
+      // Construct the towgs84 string from tx, ty, tz, rx, ry, rz, and scale_difference
+      const towgs84 = row.rx ? `${row.tx},${row.ty},${row.tz},${row.rx},${row.ry},${row.rz},${row.scale_difference}` : `${row.tx},${row.ty},${row.tz}`;
+
+      datums[`${row.source_crs_auth_name}_${row.source_crs_code}`] = {
+        towgs84
+      };
+
+      databaseDatumNames.add(normalizedDatumName);
+    });
+
+    Object.assign(datums, DATUM_OVERRIDES);
+
+    // Write updated datums back to Datum.js
+    const datumFilePath = path.resolve(__dirname, '../lib/constants/Datum.js');
+    const updatedContent = `var datums = ${JSON.stringify(datums, null, 2)};\n\n`
+      + `for (var key in datums) {\n`
+      + `  var datum = datums[key];\n`
+      + `  if (!datum.datumName) {\n`
+      + `    continue;\n`
+      + `  }\n`
+      + `  datums[datum.datumName] = datum;\n`
+      + `}\n\n`
+      + `export default datums;`;
+    fs.writeFileSync(datumFilePath, updatedContent, 'utf-8');
+    console.log('Datum.js updated successfully!');
+  }
+);