@dataforsyningen/saul 3.0.1

package/.github/workflows/release-package.yml

@@ -0,0 +1,31 @@
+name: Node.js Package
+
+on:
+  release:
+    types: [created]
+
+jobs:
+  build:
+    runs-on: ubuntu-latest
+    steps:
+      - uses: actions/checkout@v3
+      - uses: actions/setup-node@v3
+        with:
+          node-version: 16
+      - run: npm ci
+      # - run: npm test
+
+  # Important note:
+  # This only works in github actions when you've added an "automation token" by the name of "NPM_TOKEN" to "repository secrets"
+  publish-npm:
+    needs: build
+    runs-on: ubuntu-latest
+    steps:
+      - uses: actions/checkout@v3
+      - uses: actions/setup-node@v3
+        with:
+          node-version: 16
+          registry-url: https://registry.npmjs.org
+      - run: npm publish --access public
+        env:
+          NODE_AUTH_TOKEN: ${{secrets.NPM_TOKEN}}

package/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2022 SDFI
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

package/README.md

@@ -0,0 +1,45 @@
+# SAUL
+
+This is a collection of stand-alone methods for photogrammetry calculations and API interactions that are primarily used by https://github.com/SDFIdk/skraafoto_frontend
+
+You can use these as building blocks for your own web application that wants to interface with STAC API for Danske Skråfotos.
+
+## Installation
+
+1. Install iva NPM:
+```
+npm i --save @dataforsyningen/saul
+```
+
+2. Create a configuration file or Javascript object with proper values and make sure to load it into your application before using any of the SAUL utilities. This package ships with an example configuration file at `/example/config.js.example`
+
+## How to use
+
+This is an example of how to import and use those utilities in Javascript:
+```
+import auth from './config.js' // This is the path to your configuration file
+import {getSTAC} from '@dataforsyningen/saul'
+
+let response = await getSTAC('/some-stac-endpoint', auth)
+console.log(response)
+```
+
+## API docs
+
+To read the docs, first build the documentation with JSDoc:
+```
+npm run docs
+```
+A `docs` directory will be created with a bunch of HTML files containing the documentation.
+
+## How to test
+
+From within `node_modules/@dataforsyningen/saul` you can run the `test` command to check that saul core features are still working.
+```
+npm run test
+```
+
+## Acknowledgements
+
+SAUL is made available under the MIT license by
+The Skråfoto Team @ [SDFI](https://sdfi.dk/)

package/config.js.example

@@ -0,0 +1,11 @@
+export default {
+  
+  // STAC TOKEN can be aquired from https://dataforsyningen.dk/
+  API_STAC_TOKEN: "",
+  API_STAC_BASEURL: "https://api.dataforsyningen.dk/skraafotoapi",
+
+  // DHM URL and auth can be aquired from https://datafordeler.dk/dataoversigt/danmarks-hoejdemodel-dhm/koter/
+  API_DHM_BASEURL: "https://services.datafordeler.dk/DHMTerraen/DHMKoter/1.0.0/GEOREST/HentKoter",
+  API_DHM_USERNAME: "",
+  API_DHM_PASSWORD: ""
+}

package/example.html

@@ -0,0 +1,58 @@
+<!DOCTYPE html>
+<html lang="en">
+  <head>
+    <meta charset="utf-8">
+    <meta name="viewport" content="width=device-width">
+    <title>Saul modules example page</title>
+  </head>
+  <body>
+
+    <h1>Saul modules example</h1>
+
+    <p>Just opening this example.html file might not work becuase of CORS restrictions when trying to access STAC API.</p>
+    <p>You should run this HTML file from a simple server like Python's <code>http.server</code></p>
+
+    <h2>This is the result</h2>
+    <div id="result"></div>
+
+    <script type="module">
+
+      // You should add your own config.js with proper tokens, passwords, etc.
+      // Check config.js.example for info on how to set it up 
+      import auth from './config.js'
+      
+      import {image2world, world2image, getZ, iterate} from './modules/saul-core.js'
+      import {getSTAC} from './modules/api.js'
+
+      let xcoor = 579782.75
+      let ycoor = 6131096.03
+
+      const result_display_element = document.getElementById('result')
+      let apidata = await getSTAC('/collections/skraafotos2019/items/2019_83_37_2_0046_00001113', auth)
+
+      // Let's call an imported method
+      let zcoor = await getZ(xcoor, ycoor, auth)
+      
+      let image_coor = world2image(apidata, xcoor, ycoor, zcoor)
+      
+      let col = image_coor[0]
+      let row = image_coor[1]
+      
+      let world_coor = image2world(apidata, col, row, zcoor)
+      
+      col = 7375
+      row = 5382
+      let ite = await iterate(apidata, col, row, auth, 0.05)
+
+      // Write out result for the user to see
+      if (image_coor) {
+        result.innerText = xcoor + ' ' + ycoor + ' ' + zcoor + '\n' + image_coor[0] + '  ' + image_coor[1] + '\n' +ite
+      } else {
+        result.innerText = 'There was no result'
+      }
+      
+
+    </script>
+
+  </body>
+</html>

package/index.js

@@ -0,0 +1,21 @@
+// Import/export modules with their public methods
+
+export {
+  image2world,
+  world2image,
+  getZ,
+  iterate
+} from './modules/saul-core.js'
+
+export {
+  postSTAC,
+  getSTAC,
+  getDHM,
+  get,
+  post
+} from './modules/api.js'
+
+export {
+  epsg25832proj,
+  createTranslator
+} from './modules/saul-projection.js'

package/modules/api.js

@@ -0,0 +1,186 @@
+/**
+ * SAUL API utilities 
+ */
+
+let error_msg
+let load_stack = []
+
+// These custom events can only be created in a browser environment
+const loadstart = isBrowser() ? new CustomEvent('loadstart') : null
+const loadend = isBrowser() ? new CustomEvent('loadend') : null
+const loaderror = isBrowser() ? new CustomEvent('loaderror', {
+  detail: {
+    name: getErrorMsg()
+  }
+}) : null
+
+
+/** Check if code is run in a Browser or Node enviroment */
+function isBrowser() {
+  // Check if the environment is Node.js
+  if (typeof process === "object" &&
+    typeof require === "function") {
+    return false
+  }
+  // Check if the environment is
+  // a Service worker
+  if (typeof importScripts === "function") {
+    return false
+  }
+  // Check if the environment is a Browser
+  if (typeof window === "object") {
+    return true
+  }
+}
+
+/** Sends 'loadstart' event if nothing else is currently loading */
+function startLoading() {
+  if (isBrowser() && load_stack.length < 1) {
+    document.dispatchEvent(loadstart)
+  }
+  load_stack.push(1)
+}
+
+/** Sends 'loadend' event if current loads have finished */
+function endLoading() {
+    load_stack.pop()
+    if (isBrowser() && load_stack.length < 1) {
+      document.dispatchEvent(loadend)
+    }
+}
+
+/** Sends 'loaderror' event and resets load status */
+function interruptLoading() {
+  load_stack = []
+  if (isBrowser()) {
+    document.dispatchEvent(loaderror)
+  }
+}
+
+/** Getter for error messages */
+function getErrorMsg() {
+  return error_msg
+}
+
+/** Returns JSON response data */
+function HttpResponseHandler(response) {
+  if (!response.ok) {
+    error_msg = response.status
+    interruptLoading()
+    throw new Error(`HTTP error! Status: ${ response.status }`)
+  }
+  // We assume the returned data is JSON
+  return response.json()
+}
+
+/** 
+ * GET HTTP responsee from API
+ * @param {string} url - API service URL, including endpoint paths and query parameters.
+ * @param {object} [config] - Custom request configs. See https://developer.mozilla.org/en-US/docs/Web/API/Fetch_API/Using_Fetch#supplying_request_options
+ * @returns {object} response object
+ */
+function get(url, config = {}) {
+  if (!url) {
+    console.error('Could not fetch data. Missing API URL')
+  } else {
+    startLoading()
+    return fetch( url, {
+      ...config,
+      method: 'GET'
+    })
+    .then((response) => {
+      return HttpResponseHandler(response)
+    })
+    .then((response) => {
+      // Finally, return the parsed JSON response
+      endLoading()
+      return response
+    })
+    .catch((error) => {
+      // ... unless something goes wrong
+      console.error(`Fetch error: ${error}`)
+      endLoading()
+      return error
+    })
+  }
+}
+
+/** 
+ * POST HTTP request to API
+ * @param {string} url - API service URL, including endpoint paths and query parameters.
+ * @param {object} requestbody - Request data
+ * @param {string} token - Authentication token from Dataforsyningen
+ * @returns {object} response object
+ */
+function post(url, requestbody, token) {
+  if (!url || !token || !requestbody) {
+    console.error('Could not fetch data. Missing API token, request body, or URL')
+  } else {
+    startLoading()
+    return fetch( url, {
+      method: 'POST',
+      headers: {
+        'token': token,
+        'Content-Type': 'application/json'
+      },
+      body: JSON.stringify(requestbody)
+    })
+    .then((response) => {
+      return HttpResponseHandler(response)
+    })
+    .then((response) => {
+      // Finally, return the parsed JSON response
+      endLoading()
+      return response
+    })
+    .catch((error) => {
+      // ... unless something goes wrong
+      console.error(`Fetch error: ${error}`)
+      endLoading()
+      return error
+    })
+  }
+}
+
+/** 
+ * API method to fetch data from DHM
+ * @param {string} query - DHM API query. Find details at https://datafordeler.dk/dataoversigt/danmarks-hoejdemodel-dhm/koter/
+ * @param {{API_DHM_BASEURL: string, API_DHM_USERNAME: string, API_DHM_PASSWORD: string}} auth - API autentication data. See ../config.js.example for reference.
+ * @returns {object} Response data
+ */
+function getDHM(query, auth) {
+  const auth_params = `&username=${auth.API_DHM_USERNAME}&password=${auth.API_DHM_PASSWORD}`
+  return get(encodeURI(auth.API_DHM_BASEURL + query + auth_params))
+  .then((data) => data)
+}
+
+/** 
+ * API method to GET data from STAC API
+ * @param {string} query - STAC API query.
+ * @param {{API_STAC_BASEURL: string, API_STAC_TOKEN: string}} auth - API autentication data. See ../config.js.example for reference.
+ * @returns {object} Response data
+ */
+function getSTAC(query, auth) {
+  return get(auth.API_STAC_BASEURL + query, {headers: {'token': auth.API_STAC_TOKEN}})
+  .then((data) => data)
+}
+
+/** 
+ * API method to POST data to STAC API
+ * @param {string} endpoint - STAC API endpoint. Ex. `/collections`
+ * @param {object} data - request body
+ * @param {{API_STAC_BASEURL: string, API_STAC_TOKEN: string}} auth - API autentication data. See ../config.js.example for reference.
+ * @returns {object} Reponse data
+ */
+function postSTAC(endpoint, data, auth) {
+  return post(auth.API_STAC_BASEURL + endpoint, data, auth.API_STAC_TOKEN)
+  .then((data) => data)
+}
+
+export {
+  get,
+  post,
+  getSTAC,
+  postSTAC,
+  getDHM
+}

package/modules/saul-core.js

@@ -0,0 +1,181 @@
+/**
+ * SAUL photogrammetry utilities 
+ */
+ 
+import { getDHM } from './api.js'
+
+/** 
+ * Converts x,y coordinates from an image to real world lat,lon coordinates
+ * @param {object} image_data - skraafoto-stac-api image data
+ * @param {number} col - Image x coordinate (from left to right). Should be a coordinate for the entire image, not just the part displayed in the viewport.
+ * @param {number} row - Image y coordinate (from top to bottom). Should be a coordinate for the entire image, not just the part displayed in the viewport.
+ * @param {number} [Z] - Elevation (geoide)
+ * @returns {array} [longitude, latitude, elevation] 
+ */
+function image2world(image_data, col, row, Z = 0) {
+
+  // constants pulled from image_data
+  const xx0  = image_data.properties['pers:interior_orientation'].principal_point_offset[0]
+  const yy0  = image_data.properties['pers:interior_orientation'].principal_point_offset[1]
+  const ci   = image_data.properties['pers:interior_orientation'].focal_length
+  const pix  = image_data.properties['pers:interior_orientation'].pixel_spacing[0]
+  const dimXi = image_data.properties['pers:interior_orientation'].sensor_array_dimensions[0]
+  const dimYi = image_data.properties['pers:interior_orientation'].sensor_array_dimensions[1]
+  const X0   = image_data.properties['pers:perspective_center'][0]
+  const Y0   = image_data.properties['pers:perspective_center'][1]
+  const Z0   = image_data.properties['pers:perspective_center'][2]
+  const Ome  = image_data.properties['pers:omega']
+  const Phi  = image_data.properties['pers:phi']
+  const Kap  = image_data.properties['pers:kappa']
+
+  // recalc values
+  var c = ci*(-1)
+  var dimX = dimXi*pix/2*(-1)
+  var dimY = dimYi*pix/2*(-1)
+  
+  // ... Do calculations ...
+  var o = radians(Ome)
+  var p = radians(Phi)
+  var k = radians(Kap)
+  var D11 =   Math.cos(p) * Math.cos(k)
+  var D12 = - Math.cos(p) * Math.sin(k)
+  var D13 =   Math.sin(p)
+  var D21 =   Math.cos(o) * Math.sin(k) + Math.sin(o) * Math.sin(p) * Math.cos(k)
+  var D22 =   Math.cos(o) * Math.cos(k) - Math.sin(o) * Math.sin(p) * Math.sin(k)
+  var D23 = - Math.sin(o) * Math.cos(p)
+  var D31 =   Math.sin(o) * Math.sin(k) - Math.cos(o) * Math.sin(p) * Math.cos(k)
+  var D32 =   Math.sin(o) * Math.cos(k) + Math.cos(o) * Math.sin(p) * Math.sin(k)
+  var D33 =   Math.cos(o) * Math.cos(p)
+
+  var x_dot = ((col*pix)-dimX*-1)-xx0
+  var y_dot = ((row*pix)-dimY*-1)-yy0
+
+  var kx = (D11*x_dot + D12*y_dot + D13*c)/(D31*x_dot + D32*y_dot + D33*c)
+  var ky = (D21*x_dot + D22*y_dot + D23*c)/(D31*x_dot + D32*y_dot + D33*c)
+
+  var X = (Z-Z0)*kx + X0
+  var Y = (Z-Z0)*ky + Y0
+
+  return[X,Y,Z]
+}
+
+/** 
+ * Converts lat,lon coordinates to x,y coordinates within a specific image
+ * @param {Object} image_data - skraafoto-stac-api image data
+ * @param {Number} Y - northing
+ * @param {Number} X - easting
+ * @param {Number} [Z] - elevation (geoide)
+ * @returns {array} [x,y] Column/row image coordinate 
+ */
+function world2image(image_data, X, Y, Z = 0) {
+
+  // constants pulled from image_data
+  const xx0  = image_data.properties['pers:interior_orientation'].principal_point_offset[0]
+  const yy0  = image_data.properties['pers:interior_orientation'].principal_point_offset[1]
+  const ci   = image_data.properties['pers:interior_orientation'].focal_length
+  const pix  = image_data.properties['pers:interior_orientation'].pixel_spacing[0]
+  const dimXi = image_data.properties['pers:interior_orientation'].sensor_array_dimensions[0]
+  const dimYi = image_data.properties['pers:interior_orientation'].sensor_array_dimensions[1]
+  const X0   = image_data.properties['pers:perspective_center'][0]
+  const Y0   = image_data.properties['pers:perspective_center'][1]
+  const Z0   = image_data.properties['pers:perspective_center'][2]
+  const Ome  = image_data.properties['pers:omega']
+  const Phi  = image_data.properties['pers:phi']
+  const Kap  = image_data.properties['pers:kappa']
+
+  // recalc values
+  var c = ci*(-1)
+  var dimX = dimXi*pix/2*(-1)
+  var dimY = dimYi*pix/2*(-1)
+  
+  // ... Do calculations ...
+  var o = radians(Ome)
+  var p = radians(Phi)
+  var k = radians(Kap)
+  var D11 =   Math.cos(p) * Math.cos(k)
+  var D12 = - Math.cos(p) * Math.sin(k)
+  var D13 =   Math.sin(p)
+  var D21 =   Math.cos(o) * Math.sin(k) + Math.sin(o) * Math.sin(p) * Math.cos(k)
+  var D22 =   Math.cos(o) * Math.cos(k) - Math.sin(o) * Math.sin(p) * Math.sin(k)
+  var D23 = - Math.sin(o) * Math.cos(p)
+  var D31 =   Math.sin(o) * Math.sin(k) - Math.cos(o) * Math.sin(p) * Math.cos(k)
+  var D32 =   Math.sin(o) * Math.cos(k) + Math.cos(o) * Math.sin(p) * Math.sin(k)
+  var D33 =   Math.cos(o) * Math.cos(p)
+
+  var x_dot = (-1)*c*((D11*(X-X0)+D21*(Y-Y0)+D31*(Z-Z0))/(D13*(X-X0)+D23*(Y-Y0)+D33*(Z-Z0)))
+  var y_dot = (-1)*c*((D12*(X-X0)+D22*(Y-Y0)+D32*(Z-Z0))/(D13*(X-X0)+D23*(Y-Y0)+D33*(Z-Z0)))
+
+  var col = ((x_dot-xx0)+(dimX))*(-1)/pix
+  var row = ((y_dot-yy0)+(dimY))*(-1)/pix
+
+  return [col, row]
+}
+
+/** Converts degress to radians */
+function radians(degrees) {
+  return degrees * (Math.PI / 180)
+}
+
+/** 
+ * Fetches elevation based on X,Y coordinates using DHM/Koter endpoint
+ * @param {number} xcoor - EPSG:25832 X coordinate
+ * @param {number} ycoor - EPSG:25832 Y coordinate
+ * @param {{API_DHM_BASEURL: string, API_DHM_USERNAME: string, API_DHM_PASSWORD: string}} auth - API autentication data. See ../config.js.example for reference.
+ * @returns {number} Elevation in meters 
+ */
+async function getZ(xcoor, ycoor, auth) {
+  let zcoor_data = await getDHM(`?geop=POINT(${xcoor} ${ycoor})&elevationmodel=dtm`, auth)
+  let z = zcoor_data.HentKoterRespons.data[0].kote
+  return z
+}
+
+/** Iterates guessing at a world coordinate using image coordinates and elevation info */
+function iterateRecursive(image_data, col, row, z, count, limit, auth, i) {
+
+  // Get coordinates based on current z value
+  const worldcoor = image2world(image_data,col,row,z)
+
+  i--
+  if (i <= 0) {
+    // max iterations reached
+    return [worldcoor, 0, count]
+  }
+
+  // Get new z value from coordinates
+  return getZ(worldcoor[0],worldcoor[1], auth).then((new_z) => {
+    
+    // How big is the different between z and new z?
+    const delta = Math.abs(new_z - z)
+    
+    if (delta >= limit) {
+      // If the difference is too big, try building coordinates with the new z
+      count = count + 1
+      return iterateRecursive(image_data, col, row, new_z, count, limit, auth, i)
+    } else {
+      // If the difference is small, return coordinates
+      return [worldcoor, delta, count]
+    }
+  })
+}
+
+/** 
+ * Tries to guess world coordinate for a pixel position within an image using STAC API image data and requests to DHM elevation data.
+ * @param {object} image_data - image item data from STAC API
+ * @param {number} col - image x coordinate (from left to right)
+ * @param {number} row - image y coordinate (from top to bottom)
+ * @param {{API_DHM_BASEURL: string, API_DHM_USERNAME: string, API_DHM_PASSWORD: string}} auth - API autentication data. See ../config.js.example for reference.
+ * @param {number} [limit] - result may be inaccurate within this limit. Default is 1.
+ * @returns {array} [world coordindates (array), elevation discrepancy, calculation iterations]
+ */
+function iterate(image_data, col, row, auth, limit = 1) {
+  let i = 20 // Maximum number of iterations allowed
+  return iterateRecursive(image_data, col, row, 0.5, 0, limit, auth, i)
+}
+
+export { 
+  image2world,
+  world2image,
+  getZ,
+  iterate
+}
+ 

package/modules/saul-projection.js

@@ -0,0 +1,33 @@
+/**
+ * SAUL map projection and transformation utilities 
+ */
+
+import proj4 from 'proj4'
+
+/** Defines an EPSG:25832 projection using proj4.
+ * Learn about EPSG:25832 at https://epsg.io/25832 and proj4 at http://proj4js.org/
+ * EPSG:25832 definition can be requested at Definition can be requested from https://epsg.io/?format=json&q=25832
+ * @param {object} proj4object - a proj4 instance
+ * @returns {object} Returns same proj4 instance with an added EPSG:25832 projection definition
+ */
+function epsg25832proj(proj4object) {
+  proj4object.defs('EPSG:25832', "+proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs")
+  return proj4object
+}
+
+/** Creates an object with methods for translating coordinates from one projection to the other
+ * @param {string} [projection1] - from Projection
+ * @param {string} [projection2] - to projection
+ * @returns {object} An object with two methods is returned, its methods are `forward` which projects from the first projection to the second and `inverse` which projects from the second to the first.
+ */
+function createTranslator(projection1 = 'WGS84', projection2 = 'EPSG:25832') {
+  if (!proj4.defs['EPSG:25832']) {
+    epsg25832proj(proj4) // Define EPSG:25832 projection if missing
+  }
+  return proj4(projection1, projection2)
+}
+
+export {
+  epsg25832proj,
+  createTranslator
+}

package/package.json

@@ -0,0 +1,32 @@
+{
+  "name": "@dataforsyningen/saul",
+  "version": "3.0.1",
+  "description": "This is a collection of stand-alone methods for photogrammetry calculations and API interactions for 'Danske Skråfotos'.",
+  "main": "index.js",
+  "type": "module",
+  "scripts": {
+    "test": "echo 'TESTS WILL ONLY RUN UNDER NODE v18' && node ./test/saul-projection.test.js && node ./test/saul-core.test.js",
+    "docs": "jsdoc ./modules/*.js --destination docs"
+  },
+  "repository": {
+    "type": "git",
+    "url": "git+https://github.com/SDFIdk/saul.git"
+  },
+  "keywords": [
+    "skråfoto",
+    "SDFI",
+    "saul"
+  ],
+  "author": "Styrelsen for Dataforsyning og Infrastruktur (https://sdfi.dk/)",
+  "license": "MIT",
+  "bugs": {
+    "url": "https://github.com/SDFIdk/saul/issues"
+  },
+  "homepage": "https://github.com/SDFIdk/saul#readme",
+  "devDependencies": {
+    "jsdoc": "^3.6.10"
+  },
+  "dependencies": {
+    "proj4": "^2.8.0"
+  }
+}

package/test/saul-core.test.js

@@ -0,0 +1,119 @@
+// Testing saul-core methods
+
+// You should add your own config.js with proper tokens, passwords, etc.
+// Check config.js.example for info on how to set it up
+import auth from '../config.js'
+
+import assert from 'assert'
+import {getSTAC} from '../modules/api.js'
+import {world2image, image2world, getZ, iterate} from '../modules/saul-core.js'
+
+/*
+const item = await getSTAC('/collections/skraafotos2019/items/2019_83_37_2_0046_00001113', auth)
+const world_x = 580341
+const world_y = 6130925
+const image_x = 7373
+const image_y = 3205
+const world_elevation = 32.7199
+*/
+
+const item = await getSTAC('/collections/skraafotos2021/items/2021_83_36_4_0013_00003824', auth)
+const world_x = 587739
+const world_y = 6130951
+const image_x = 41839.698939
+const image_y = 3204.326930
+const world_elevation = 38.874336
+
+function is_equalIsh(num1, num2) {
+  const deviation = 0.05
+  if (Math.abs(num1 - num2) > deviation) {
+    return false
+  } else {
+    return true
+  }
+} 
+
+// Test getZ
+try {
+
+  let elevation = await getZ(world_x, world_y, auth)
+
+  assert(is_equalIsh(elevation, world_elevation), "getZ fail: This is not the Z value we were hoping for.")
+  
+  console.log("Test getZ OK")
+
+} catch(error) {
+  console.error(error)
+}
+
+// Test iterate
+try {
+
+  let ite = await iterate(item, image_x, image_y, auth, 0.05)
+
+  assert(is_equalIsh(ite[0][0], world_x), "iterate fail: The coordinate values are not equal") //test if the values are equal
+  assert(is_equalIsh(ite[0][1], world_y), "iterate fail: The coordinate values are not equal") //test if the values are equal
+
+  console.log("Test iterate OK")
+
+} catch(error) {
+  console.error(error)
+}
+
+// Test world2image
+try {
+
+  let xy = world2image(item, world_x, world_y, world_elevation)
+
+  assert(is_equalIsh(xy[0], image_x), "world2image fail: The x coordinate values are not equal")
+  assert(is_equalIsh(xy[1], image_y), "world2image fail: The y coordinate values are not equal")
+  
+  console.log("Test world2image OK")
+
+} catch(error) {
+  console.error(error)
+}
+
+// Test image2world
+try {
+
+  let coords = image2world(item, image_x, image_y, world_elevation)
+
+  assert(is_equalIsh(coords[0], world_x), "image2world fail: The longitude values are not equal")
+  assert(is_equalIsh(coords[1], world_y), "image2world fail: The latitude values are not equal")
+
+  console.log("Test image2world OK")
+
+} catch(error) {
+  console.error(error)
+}
+
+// Test world2image > image2world 
+try {
+
+  let xy1 = world2image(item, world_x, world_y, world_elevation)
+  let coords1 = image2world(item, xy1[0], xy1[1], world_elevation)
+
+  assert(is_equalIsh(coords1[0], world_x), "world2image > image2world fail: The longitude values are not equal")
+  assert(is_equalIsh(coords1[1], world_y), "world2image > image2world fail: The latitude values are not equal")
+
+  console.log("Test world2image > image2world OK")
+
+} catch(error) {
+  console.error(error)
+}
+
+// Test world2image > iterate 
+try {
+
+  let xy2 = world2image(item, world_x, world_y, world_elevation)
+  let coords2 = iterate(item, xy2[0], xy2[1], auth)
+
+  assert(is_equalIsh(coords2[0], world_x), "world2image > iterate fail: The longitude values are not equal")
+  assert(is_equalIsh(coords2[1], world_y), "world2image > iterate fail: The latitude values are not equal")
+
+  console.log("Test world2image > iterate OK")
+
+} catch(error) {
+  console.error(error)
+}

package/test/saul-projection.test.js

@@ -0,0 +1,34 @@
+// Testing saul-projection methods
+
+import assert from 'assert'
+import proj4 from 'proj4'
+import {epsg25832proj, createTranslator} from '../modules/saul-projection.js'
+
+// Test epsg25832proj
+try {
+
+  const projection = epsg25832proj(proj4)
+
+  assert.ok(projection.defs['EPSG:25832'], "Projection not defined.")
+  
+  console.log("Test epsg25832proj Ok")
+
+} catch(error) {
+  console.error(error)
+}
+
+// Test createTranslator
+try {
+
+  const translator = createTranslator()
+  const forwarded_coors = translator.forward([10.31,55.37])
+  const inverse_coors = translator.inverse(forwarded_coors)
+
+  assert.strictEqual(forwarded_coors[0], 583023.3267001238, "Failed to forward coordinates")
+  assert.strictEqual(inverse_coors[0], 10.31, "Failed to inverse coordinates")
+  
+  console.log("Test createTranslator Ok")
+
+} catch(error) {
+  console.error(error)
+}