From d84be88a91f25c337234f5a60f336a1785569080 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Marc-Andr=C3=A9=20Barbeau?=
 <marc-andre.barbeau@msp.gouv.qc.ca>
Date: Thu, 31 Oct 2019 13:11:11 -0400
Subject: [PATCH] feat(projection): add mtm projection

---
 packages/geo/src/lib/map/shared/map.utils.ts | 32 +++++++++++++++-----
 1 file changed, 25 insertions(+), 7 deletions(-)

diff --git a/packages/geo/src/lib/map/shared/map.utils.ts b/packages/geo/src/lib/map/shared/map.utils.ts
index 3b568eacc3..1d276d13c4 100644
--- a/packages/geo/src/lib/map/shared/map.utils.ts
+++ b/packages/geo/src/lib/map/shared/map.utils.ts
@@ -35,7 +35,7 @@ export function stringToLonLat(str: string, mapProjection: string): {lonLat: [nu
   let directionLat: string;
   let decimalLat: string;
   let pattern: string;
-  let timeZone: string;
+  let zone: string;
   let radius: string;
   let conf: string;
   let lon: any;
@@ -57,6 +57,9 @@ export function stringToLonLat(str: string, mapProjection: string): {lonLat: [nu
   const patternUtm = '(UTM)\-?(\\d{1,2})[\\s,.]*(\\d+[\\s.,]?\\d+)[\\s,.]+(\\d+[\\s.,]?\\d+)';
   const utmRegex =  new RegExp(`^${patternUtm}`, 'gi');
 
+  const patternMtm = '(MTM)\-?(\\d{1,2})[\\s,.]*(\\d+[\\s.,]?\\d+)[\\s,.]+(\\d+[\\s.,]?\\d+)';
+  const mtmRegex =  new RegExp(`^${patternMtm}`, 'gi');
+
   const ddCoord = '([-+])?(\\d{1,3})[,.](\\d+)';
   const patternDd = `${ddCoord}[,.]?\\s*${ddCoord}`;
   const ddRegex =  new RegExp(`^${patternDd}`, 'g');
@@ -112,12 +115,27 @@ export function stringToLonLat(str: string, mapProjection: string): {lonLat: [nu
       lon = convertDMSToDD(parseFloat(degreesLon), parseFloat(minutesLon), parseFloat(secondsLon), directionLon);
       lat = convertDMSToDD(parseFloat(degreesLat), parseFloat(minutesLat), parseFloat(secondsLat), directionLat);
 
-    } else if (utmRegex.test(coordStr)) {
-      isXYCoords = true;
-      [, pattern, timeZone, lon, lat] = coordStr.match(patternUtm);
-      const utm = '+proj=' + pattern + ' +zone=' + timeZone;
-      const wgs84 = '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs';
-      [lon, lat] = proj4(utm.toLocaleLowerCase(), wgs84, [parseFloat(lon), parseFloat(lat)]);
+  } else if (utmRegex.test(coordStr)) {
+    isXYCoords = true;
+    [, pattern, zone, lon, lat] = coordStr.match(patternUtm);
+    const utm = '+proj=' + pattern + ' +zone=' + zone;
+    const wgs84 = '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs';
+    [lon, lat] = proj4(utm.toLocaleLowerCase(), wgs84, [parseFloat(lon), parseFloat(lat)]);
+
+  } else if (mtmRegex.test(coordStr)) {
+    isXYCoords = true;
+    [, pattern, zone, lon, lat] = coordStr.match(patternMtm);
+    let lon0;
+    if (Number(zone) <=2) {
+      lon0 = -50 - Number(zone) * 3;
+    } else if (Number(zone) >=12) {
+      lon0 = -81 - (Number(zone)-12) * 3;
+    } else {
+      lon0 = -49.5 - Number(zone) * 3;
+    }
+    const mtm = `+proj=tmerc +lat_0=0 +lon_0=${lon0} +k=0.9999 +x_0=304800 +y_0=0 +ellps=GRS80 +units=m +no_defs`;
+    const wgs84 = '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs';
+    [lon, lat] = proj4(mtm, wgs84, [parseFloat(lon), parseFloat(lat)]);
 
   } else if (dmdRegex.test(coordStr)) {
     [,