Skip to content

convertTrajectory bugs

@fhoueto.amz CC - @gehrmann @debasisc

  • "unitZ": "osdu:reference-data--UnitOfMeasure:ft[US]" does not work and neither does "unitZ": "osdu:reference-data--UnitOfMeasure:ft%5BUS%5D" (there is reference data with record id: "id": "osdu:reference-data--UnitOfMeasure:ft%5BUS%5D" in the system.
  • unitXY in api spec is said to be optional but in API call it is mandatory. I believe it should be optional because the unit of the CRS should be retrieved from the CRS PR. There appears to be a bug that unitXY is assumed to be the unit of the projCRS. The bug should be fixed or as workaround it must be said the unitXY is set to the unit of the projected CRS.
  • unitXY in api spec is said to only control dx,dy input but it appears to somehow relate to projected coordinates. The unit of the projected CRS is not respected. The main issue I think is the unit of projected CRS is ignored (i.e., it should not be needed to know the unit of the CRS to compute N and E).

The following are simple test cases where from a starting point we "drill" due north horizontally. So we expect the Y coordinates to increase by 100m if the MD is 100m.

I tried it with record-id and with PR for UOM. Both have the bug.

Case 1: projCRS is in ft[US]; unitZ in m; unitXY set to m

  • Output Y at station 2 is expected to be around 9,000,000 ftUS + 100m = 9,000,328.08 ftUS.
  • The actual output is 9,000,100. (FAIL) The local min. curvature output of 100m unitZ was converted to XY unit as 100m.
Request 1 
{
  "azimuthReference": "GN",
  "interpolate": false,
  "referencePoint": {
    "x": 1000000.0,
    "y": 9000000.0,
    "z": 0
  },
  "unitZ": "osdu:reference-data--UnitOfMeasure:m",
  "unitXY": "osdu:reference-data--UnitOfMeasure:m",
  "inputStations": [
    {
      "md": 0,
      "azimuth": 0,
      "inclination": 90
    },
    {
      "md": 100,
      "azimuth": 0,
      "inclination": 90
    }
  ],
  "trajectoryCRS": "{{NAMESPACE}}:reference-data--CoordinateReferenceSystem:BoundProjected:EPSG::32066_EPSG::15851",
  "inputKind": "MD_Incl_Azim",
    "method": "AzimuthalEquidistant"
}
Response 1 
{
    "trajectoryCRS": "{\"authCode\":{\"auth\":\"OSDU\",\"code\":\"32066079\"},\"lateBoundCRS\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"32066\"},\"name\":\"NAD_1927_BLM_Zone_16N\",\"type\":\"LBC\",\"ver\":\"PE_10_9_1\",\"wkt\":\"PROJCS[\\\"NAD_1927_BLM_Zone_16N\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],PROJECTION[\\\"Transverse_Mercator\\\"],PARAMETER[\\\"False_Easting\\\",1640416.666666667],PARAMETER[\\\"False_Northing\\\",0.0],PARAMETER[\\\"Central_Meridian\\\",-87.0],PARAMETER[\\\"Scale_Factor\\\",0.9996],PARAMETER[\\\"Latitude_Of_Origin\\\",0.0],UNIT[\\\"Foot_US\\\",0.3048006096012192],AUTHORITY[\\\"EPSG\\\",32066]]\"},\"name\":\"NAD27 * OGP-Usa Conus / BLM zone 16N (US survey feet) [32066,15851]\",\"singleCT\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"15851\"},\"name\":\"NAD_1927_To_WGS_1984_79_CONUS\",\"type\":\"ST\",\"ver\":\"PE_10_9_1\",\"wkt\":\"GEOGTRAN[\\\"NAD_1927_To_WGS_1984_79_CONUS\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],GEOGCS[\\\"GCS_WGS_1984\\\",DATUM[\\\"D_WGS_1984\\\",SPHEROID[\\\"WGS_1984\\\",6378137.0,298.257223563]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],METHOD[\\\"NADCON\\\"],PARAMETER[\\\"Dataset_conus\\\",0.0],OPERATIONACCURACY[5.0],AUTHORITY[\\\"EPSG\\\",15851]]\"},\"type\":\"EBC\",\"ver\":\"PE_10_9_1\"}",
    "unitXY": "{\"scaleOffset\":{\"scale\":1.0,\"offset\":0.0},\"symbol\":\"m\",\"baseMeasurement\":{\"ancestry\":\"Length\",\"type\":\"UM\"},\"type\":\"USO\"}",
    "unitZ": "{\"scaleOffset\":{\"scale\":1.0,\"offset\":0.0},\"symbol\":\"m\",\"baseMeasurement\":{\"ancestry\":\"Length\",\"type\":\"UM\"},\"type\":\"USO\"}",
    "unitDls": "{\"scaleOffset\":{\"scale\":5.81776417331443E-4,\"offset\":0.0},\"symbol\":\"deg/30m\",\"baseMeasurement\":{\"ancestry\":\"Rotation_Per_Length\",\"type\":\"UM\"},\"type\":\"USO\"}",
    "stations": [
        {
            "md": 0.0,
            "inclination": 90.0,
            "azimuthTN": 359.19005152801196,
            "azimuthGN": 0.0,
            "dxTN": 0.0,
            "dyTN": 0.0,
            "point": {
                "x": 1000000.000000025,
                "y": 8999999.999999689,
                "z": 0.0
            },
            "wgs84Longitude": -88.93078989815679,
            "wgs84Latitude": 24.79293652254376,
            "dls": 0.0,
            "original": true,
            "dz": 0.0
        },
        {
            "md": 100.0,
            "inclination": 90.0,
            "azimuthTN": 359.19005152801196,
            "azimuthGN": 0.0,
            "dxTN": -1.4135796794430462,
            "dyTN": 99.9900084632953,
            "point": {
                "x": 1000000.0000073737,
                "y": 9000100.007059755,
                "z": -6.123233995736766E-15
            },
            "wgs84Longitude": -88.93079417227369,
            "wgs84Latitude": 24.793211666962428,
            "dls": 0.0,
            "original": true,
            "dz": 6.123233995736766E-15
        }
    ],
    "localCRS": "{\"lateBoundCRS\":{\"name\":\"Azimuthal Equidistant\",\"type\":\"LBC\",\"ver\":\"PE_10_9_1\",\"wkt\":\"PROJCS[\\\"Azimuthal Equidistant Lng=-88.93089504;Lat=24.79259766\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],PROJECTION[\\\"Modified Azimuthal_Equidistant\\\"],PARAMETER[\\\"False_Easting\\\",0.0],PARAMETER[\\\"False_Northing\\\",0.0],PARAMETER[\\\"Central_Meridian\\\",-88.93089504402704],PARAMETER[\\\"Latitude_Of_Origin\\\",24.792597655454724],UNIT[\\\"Foot_US\\\",0.3048006096012192]]\"},\"name\":\"Azimuthal Equidistant - NAD_1927_To_WGS_1984_79_CONUS\",\"singleCT\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"15851\"},\"name\":\"NAD_1927_To_WGS_1984_79_CONUS\",\"type\":\"ST\",\"ver\":\"PE_10_9_1\",\"wkt\":\"GEOGTRAN[\\\"NAD_1927_To_WGS_1984_79_CONUS\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],GEOGCS[\\\"GCS_WGS_1984\\\",DATUM[\\\"D_WGS_1984\\\",SPHEROID[\\\"WGS_1984\\\",6378137.0,298.257223563]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],METHOD[\\\"NADCON\\\"],PARAMETER[\\\"Dataset_conus\\\",0.0],OPERATIONACCURACY[5.0],AUTHORITY[\\\"EPSG\\\",15851]]\"},\"type\":\"EBC\",\"ver\":\"PE_10_9_1\"}",
    "method": "AzimuthalEquidistant",
    "operationsApplied": [
        "derived TN from GN azimuth by grid convergence 359.190052",
        "computed deflections via minimum curvature method",
        "computation method: AzimuthalEquidistant",
        "conversion from 'Azimuthal Equidistant' to 'GCS_North_American_1927'",
        "conversion from 'GCS_North_American_1927' to 'NAD_1927_BLM_Zone_16N'",
        "to WGS 84: conversion from NAD_1927_BLM_Zone_16N to GCS_North_American_1927; 2 points converted",
        "to WGS 84: transformation GCS_North_American_1927 to GCS_WGS_1984 using NAD_1927_To_WGS_1984_79_CONUS; 2 points successfully transformed"
    ],
    "inputKind": "MD_Incl_Azim"
}

Case 2: projCRS is in ft[US]; unitZ in ft; unitXY set to ft

  • Output Y at station 2 is expected to be around 9,000,000 ftUS + 100 ft = 9,000,100.0 ftUS.
  • The actual output is 9,000,100. (PASS) The local min. curvature output of 100ft unitZ was converted to XY unit as 100 ft.
Request 2 
{
  "azimuthReference": "GN",
  "interpolate": false,
  "referencePoint": {
    "x": 1000000.0,
    "y": 9000000.0,
    "z": 0
  },
  "unitZ": "osdu:reference-data--UnitOfMeasure:ft",
  "unitXY": "osdu:reference-data--UnitOfMeasure:ft",
  "inputStations": [
    {
      "md": 0,
      "azimuth": 0,
      "inclination": 90
    },
    {
      "md": 100,
      "azimuth": 0,
      "inclination": 90
    }
  ],
  "trajectoryCRS": "{{NAMESPACE}}:reference-data--CoordinateReferenceSystem:BoundProjected:EPSG::32066_EPSG::15851",
  "inputKind": "MD_Incl_Azim",
  "method": "AzimuthalEquidistant"
}
Response 2 
{
    "trajectoryCRS": "{\"authCode\":{\"auth\":\"OSDU\",\"code\":\"32066079\"},\"lateBoundCRS\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"32066\"},\"name\":\"NAD_1927_BLM_Zone_16N\",\"type\":\"LBC\",\"ver\":\"PE_10_9_1\",\"wkt\":\"PROJCS[\\\"NAD_1927_BLM_Zone_16N\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],PROJECTION[\\\"Transverse_Mercator\\\"],PARAMETER[\\\"False_Easting\\\",1640416.666666667],PARAMETER[\\\"False_Northing\\\",0.0],PARAMETER[\\\"Central_Meridian\\\",-87.0],PARAMETER[\\\"Scale_Factor\\\",0.9996],PARAMETER[\\\"Latitude_Of_Origin\\\",0.0],UNIT[\\\"Foot_US\\\",0.3048006096012192],AUTHORITY[\\\"EPSG\\\",32066]]\"},\"name\":\"NAD27 * OGP-Usa Conus / BLM zone 16N (US survey feet) [32066,15851]\",\"singleCT\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"15851\"},\"name\":\"NAD_1927_To_WGS_1984_79_CONUS\",\"type\":\"ST\",\"ver\":\"PE_10_9_1\",\"wkt\":\"GEOGTRAN[\\\"NAD_1927_To_WGS_1984_79_CONUS\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],GEOGCS[\\\"GCS_WGS_1984\\\",DATUM[\\\"D_WGS_1984\\\",SPHEROID[\\\"WGS_1984\\\",6378137.0,298.257223563]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],METHOD[\\\"NADCON\\\"],PARAMETER[\\\"Dataset_conus\\\",0.0],OPERATIONACCURACY[5.0],AUTHORITY[\\\"EPSG\\\",15851]]\"},\"type\":\"EBC\",\"ver\":\"PE_10_9_1\"}",
    "unitXY": "{\"abcd\":{\"a\":0.0,\"b\":0.3048,\"c\":1.0,\"d\":0.0},\"symbol\":\"ft\",\"baseMeasurement\":{\"ancestry\":\"L\",\"type\":\"UM\"},\"type\":\"UAD\"}",
    "unitZ": "{\"abcd\":{\"a\":0.0,\"b\":0.3048,\"c\":1.0,\"d\":0.0},\"symbol\":\"ft\",\"baseMeasurement\":{\"ancestry\":\"L\",\"type\":\"UM\"},\"type\":\"UAD\"}",
    "unitDls": "{\"scaleOffset\":{\"scale\":5.72614583987641E-4,\"offset\":0.0},\"symbol\":\"deg/100ft\",\"baseMeasurement\":{\"ancestry\":\"Rotation_Per_Length\",\"type\":\"UM\"},\"type\":\"USO\"}",
    "stations": [
        {
            "md": 0.0,
            "inclination": 90.0,
            "azimuthTN": 359.19005152801196,
            "azimuthGN": 0.0,
            "dxTN": 0.0,
            "dyTN": 0.0,
            "point": {
                "x": 1000000.000000025,
                "y": 8999999.999999689,
                "z": 0.0
            },
            "wgs84Longitude": -88.93078989815679,
            "wgs84Latitude": 24.79293652254376,
            "dls": 0.0,
            "original": true,
            "dz": 0.0
        },
        {
            "md": 99.99999999999999,
            "inclination": 90.0,
            "azimuthTN": 359.19005152801196,
            "azimuthGN": 0.0,
            "dxTN": -1.413579679443046,
            "dyTN": 99.99000846329528,
            "point": {
                "x": 1000000.0000073737,
                "y": 9000100.007059755,
                "z": -6.123233995736765E-15
            },
            "wgs84Longitude": -88.93079417227369,
            "wgs84Latitude": 24.793211666962428,
            "dls": 0.0,
            "original": true,
            "dz": 6.123233995736765E-15
        }
    ],
    "localCRS": "{\"lateBoundCRS\":{\"name\":\"Azimuthal Equidistant\",\"type\":\"LBC\",\"ver\":\"PE_10_9_1\",\"wkt\":\"PROJCS[\\\"Azimuthal Equidistant Lng=-88.93089504;Lat=24.79259766\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],PROJECTION[\\\"Modified Azimuthal_Equidistant\\\"],PARAMETER[\\\"False_Easting\\\",0.0],PARAMETER[\\\"False_Northing\\\",0.0],PARAMETER[\\\"Central_Meridian\\\",-88.93089504402704],PARAMETER[\\\"Latitude_Of_Origin\\\",24.792597655454724],UNIT[\\\"Foot_US\\\",0.3048006096012192]]\"},\"name\":\"Azimuthal Equidistant - NAD_1927_To_WGS_1984_79_CONUS\",\"singleCT\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"15851\"},\"name\":\"NAD_1927_To_WGS_1984_79_CONUS\",\"type\":\"ST\",\"ver\":\"PE_10_9_1\",\"wkt\":\"GEOGTRAN[\\\"NAD_1927_To_WGS_1984_79_CONUS\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],GEOGCS[\\\"GCS_WGS_1984\\\",DATUM[\\\"D_WGS_1984\\\",SPHEROID[\\\"WGS_1984\\\",6378137.0,298.257223563]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],METHOD[\\\"NADCON\\\"],PARAMETER[\\\"Dataset_conus\\\",0.0],OPERATIONACCURACY[5.0],AUTHORITY[\\\"EPSG\\\",15851]]\"},\"type\":\"EBC\",\"ver\":\"PE_10_9_1\"}",
    "method": "AzimuthalEquidistant",
    "operationsApplied": [
        "derived TN from GN azimuth by grid convergence 359.190052",
        "computed deflections via minimum curvature method",
        "computation method: AzimuthalEquidistant",
        "conversion from 'Azimuthal Equidistant' to 'GCS_North_American_1927'",
        "conversion from 'GCS_North_American_1927' to 'NAD_1927_BLM_Zone_16N'",
        "to WGS 84: conversion from NAD_1927_BLM_Zone_16N to GCS_North_American_1927; 2 points converted",
        "to WGS 84: transformation GCS_North_American_1927 to GCS_WGS_1984 using NAD_1927_To_WGS_1984_79_CONUS; 2 points successfully transformed"
    ],
    "inputKind": "MD_Incl_Azim"
}

Case 3: projCRS is in ft[US]; unitZ in m; unitXY set to ft

  • Output Y at station 2 is expected to be around 9,000,000 ftUS + 100 m = 9,000,328.08 ftUS.
  • The actual output is 9,000,328.11. (PASS) The local min. curvature output of 100m unitZ was converted to XY unit as 328 ft.
Request 3 
{
  "azimuthReference": "GN",
  "interpolate": false,
  "referencePoint": {
    "x": 1000000.0,
    "y": 9000000.0,
    "z": 0
  },
  "unitZ": "osdu:reference-data--UnitOfMeasure:m",
  "unitXY": "osdu:reference-data--UnitOfMeasure:ft",
  "inputStations": [
    {
      "md": 0,
      "azimuth": 0,
      "inclination": 90
    },
    {
      "md": 100,
      "azimuth": 0,
      "inclination": 90
    }
  ],
  "trajectoryCRS": "{{NAMESPACE}}:reference-data--CoordinateReferenceSystem:BoundProjected:EPSG::32066_EPSG::15851",
  "inputKind": "MD_Incl_Azim",
  "method": "AzimuthalEquidistant"
}
Response 3 
{
    "trajectoryCRS": "{\"authCode\":{\"auth\":\"OSDU\",\"code\":\"32066079\"},\"lateBoundCRS\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"32066\"},\"name\":\"NAD_1927_BLM_Zone_16N\",\"type\":\"LBC\",\"ver\":\"PE_10_9_1\",\"wkt\":\"PROJCS[\\\"NAD_1927_BLM_Zone_16N\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],PROJECTION[\\\"Transverse_Mercator\\\"],PARAMETER[\\\"False_Easting\\\",1640416.666666667],PARAMETER[\\\"False_Northing\\\",0.0],PARAMETER[\\\"Central_Meridian\\\",-87.0],PARAMETER[\\\"Scale_Factor\\\",0.9996],PARAMETER[\\\"Latitude_Of_Origin\\\",0.0],UNIT[\\\"Foot_US\\\",0.3048006096012192],AUTHORITY[\\\"EPSG\\\",32066]]\"},\"name\":\"NAD27 * OGP-Usa Conus / BLM zone 16N (US survey feet) [32066,15851]\",\"singleCT\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"15851\"},\"name\":\"NAD_1927_To_WGS_1984_79_CONUS\",\"type\":\"ST\",\"ver\":\"PE_10_9_1\",\"wkt\":\"GEOGTRAN[\\\"NAD_1927_To_WGS_1984_79_CONUS\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],GEOGCS[\\\"GCS_WGS_1984\\\",DATUM[\\\"D_WGS_1984\\\",SPHEROID[\\\"WGS_1984\\\",6378137.0,298.257223563]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],METHOD[\\\"NADCON\\\"],PARAMETER[\\\"Dataset_conus\\\",0.0],OPERATIONACCURACY[5.0],AUTHORITY[\\\"EPSG\\\",15851]]\"},\"type\":\"EBC\",\"ver\":\"PE_10_9_1\"}",
    "unitXY": "{\"abcd\":{\"a\":0.0,\"b\":0.3048,\"c\":1.0,\"d\":0.0},\"symbol\":\"ft\",\"baseMeasurement\":{\"ancestry\":\"L\",\"type\":\"UM\"},\"type\":\"UAD\"}",
    "unitZ": "{\"scaleOffset\":{\"scale\":1.0,\"offset\":0.0},\"symbol\":\"m\",\"baseMeasurement\":{\"ancestry\":\"Length\",\"type\":\"UM\"},\"type\":\"USO\"}",
    "unitDls": "{\"scaleOffset\":{\"scale\":5.81776417331443E-4,\"offset\":0.0},\"symbol\":\"deg/30m\",\"baseMeasurement\":{\"ancestry\":\"Rotation_Per_Length\",\"type\":\"UM\"},\"type\":\"USO\"}",
    "stations": [
        {
            "md": 0.0,
            "inclination": 90.0,
            "azimuthTN": 359.19005152801196,
            "azimuthGN": 0.0,
            "dxTN": 0.0,
            "dyTN": 0.0,
            "point": {
                "x": 1000000.000000025,
                "y": 8999999.999999689,
                "z": 0.0
            },
            "wgs84Longitude": -88.93078989815679,
            "wgs84Latitude": 24.79293652254376,
            "dls": 0.0,
            "original": true,
            "dz": 0.0
        },
        {
            "md": 100.0,
            "inclination": 90.0,
            "azimuthTN": 359.19005152801196,
            "azimuthGN": 0.0,
            "dxTN": -4.637728607096608,
            "dyTN": 328.051208869079,
            "point": {
                "x": 1000000.0000791344,
                "y": 9000328.107152136,
                "z": -6.123233995736766E-15
            },
            "wgs84Longitude": -88.9308039209204,
            "wgs84Latitude": 24.79383922729075,
            "dls": 0.0,
            "original": true,
            "dz": 6.123233995736766E-15
        }
    ],
    "localCRS": "{\"lateBoundCRS\":{\"name\":\"Azimuthal Equidistant\",\"type\":\"LBC\",\"ver\":\"PE_10_9_1\",\"wkt\":\"PROJCS[\\\"Azimuthal Equidistant Lng=-88.93089504;Lat=24.79259766\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],PROJECTION[\\\"Modified Azimuthal_Equidistant\\\"],PARAMETER[\\\"False_Easting\\\",0.0],PARAMETER[\\\"False_Northing\\\",0.0],PARAMETER[\\\"Central_Meridian\\\",-88.93089504402704],PARAMETER[\\\"Latitude_Of_Origin\\\",24.792597655454724],UNIT[\\\"Foot_US\\\",0.3048006096012192]]\"},\"name\":\"Azimuthal Equidistant - NAD_1927_To_WGS_1984_79_CONUS\",\"singleCT\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"15851\"},\"name\":\"NAD_1927_To_WGS_1984_79_CONUS\",\"type\":\"ST\",\"ver\":\"PE_10_9_1\",\"wkt\":\"GEOGTRAN[\\\"NAD_1927_To_WGS_1984_79_CONUS\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],GEOGCS[\\\"GCS_WGS_1984\\\",DATUM[\\\"D_WGS_1984\\\",SPHEROID[\\\"WGS_1984\\\",6378137.0,298.257223563]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],METHOD[\\\"NADCON\\\"],PARAMETER[\\\"Dataset_conus\\\",0.0],OPERATIONACCURACY[5.0],AUTHORITY[\\\"EPSG\\\",15851]]\"},\"type\":\"EBC\",\"ver\":\"PE_10_9_1\"}",
    "method": "AzimuthalEquidistant",
    "operationsApplied": [
        "derived TN from GN azimuth by grid convergence 359.190052",
        "computed deflections via minimum curvature method",
        "computation method: AzimuthalEquidistant",
        "conversion from 'Azimuthal Equidistant' to 'GCS_North_American_1927'",
        "conversion from 'GCS_North_American_1927' to 'NAD_1927_BLM_Zone_16N'",
        "to WGS 84: conversion from NAD_1927_BLM_Zone_16N to GCS_North_American_1927; 2 points converted",
        "to WGS 84: transformation GCS_North_American_1927 to GCS_WGS_1984 using NAD_1927_To_WGS_1984_79_CONUS; 2 points successfully transformed"
    ],
    "inputKind": "MD_Incl_Azim"
}

Case 4: projCRS is in ft[US]; unitZ in ft; unitXY set to m

  • Output Y at station 2 is expected to be around 9,000,000 ftUS + 100 ft = 9,000,100.0 ftUS.
  • The actual output is 9,000,030.48. (FAIL) The local coordinates ~100ft were converted to the XY unit for some reason.
Request 4 
{
  "azimuthReference": "GN",
  "interpolate": false,
  "referencePoint": {
    "x": 1000000.0,
    "y": 9000000.0,
    "z": 0
  },
  "unitZ": "osdu:reference-data--UnitOfMeasure:ft",
  "unitXY": "osdu:reference-data--UnitOfMeasure:m",
  "inputStations": [
    {
      "md": 0,
      "azimuth": 0,
      "inclination": 90
    },
    {
      "md": 100,
      "azimuth": 0,
      "inclination": 90
    }
  ],
  "trajectoryCRS": "{{NAMESPACE}}:reference-data--CoordinateReferenceSystem:BoundProjected:EPSG::32066_EPSG::15851",
  "inputKind": "MD_Incl_Azim",
  "method": "AzimuthalEquidistant"
}
Response 4 
{
    "trajectoryCRS": "{\"authCode\":{\"auth\":\"OSDU\",\"code\":\"32066079\"},\"lateBoundCRS\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"32066\"},\"name\":\"NAD_1927_BLM_Zone_16N\",\"type\":\"LBC\",\"ver\":\"PE_10_9_1\",\"wkt\":\"PROJCS[\\\"NAD_1927_BLM_Zone_16N\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],PROJECTION[\\\"Transverse_Mercator\\\"],PARAMETER[\\\"False_Easting\\\",1640416.666666667],PARAMETER[\\\"False_Northing\\\",0.0],PARAMETER[\\\"Central_Meridian\\\",-87.0],PARAMETER[\\\"Scale_Factor\\\",0.9996],PARAMETER[\\\"Latitude_Of_Origin\\\",0.0],UNIT[\\\"Foot_US\\\",0.3048006096012192],AUTHORITY[\\\"EPSG\\\",32066]]\"},\"name\":\"NAD27 * OGP-Usa Conus / BLM zone 16N (US survey feet) [32066,15851]\",\"singleCT\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"15851\"},\"name\":\"NAD_1927_To_WGS_1984_79_CONUS\",\"type\":\"ST\",\"ver\":\"PE_10_9_1\",\"wkt\":\"GEOGTRAN[\\\"NAD_1927_To_WGS_1984_79_CONUS\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],GEOGCS[\\\"GCS_WGS_1984\\\",DATUM[\\\"D_WGS_1984\\\",SPHEROID[\\\"WGS_1984\\\",6378137.0,298.257223563]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],METHOD[\\\"NADCON\\\"],PARAMETER[\\\"Dataset_conus\\\",0.0],OPERATIONACCURACY[5.0],AUTHORITY[\\\"EPSG\\\",15851]]\"},\"type\":\"EBC\",\"ver\":\"PE_10_9_1\"}",
    "unitXY": "{\"scaleOffset\":{\"scale\":1.0,\"offset\":0.0},\"symbol\":\"m\",\"baseMeasurement\":{\"ancestry\":\"Length\",\"type\":\"UM\"},\"type\":\"USO\"}",
    "unitZ": "{\"abcd\":{\"a\":0.0,\"b\":0.3048,\"c\":1.0,\"d\":0.0},\"symbol\":\"ft\",\"baseMeasurement\":{\"ancestry\":\"L\",\"type\":\"UM\"},\"type\":\"UAD\"}",
    "unitDls": "{\"scaleOffset\":{\"scale\":5.72614583987641E-4,\"offset\":0.0},\"symbol\":\"deg/100ft\",\"baseMeasurement\":{\"ancestry\":\"Rotation_Per_Length\",\"type\":\"UM\"},\"type\":\"USO\"}",
    "stations": [
        {
            "md": 0.0,
            "inclination": 90.0,
            "azimuthTN": 359.19005152801196,
            "azimuthGN": 0.0,
            "dxTN": 0.0,
            "dyTN": 0.0,
            "point": {
                "x": 1000000.000000025,
                "y": 8999999.999999689,
                "z": 0.0
            },
            "wgs84Longitude": -88.93078989815679,
            "wgs84Latitude": 24.79293652254376,
            "dls": 0.0,
            "original": true,
            "dz": 0.0
        },
        {
            "md": 99.99999999999999,
            "inclination": 90.0,
            "azimuthTN": 359.19005152801196,
            "azimuthGN": 0.0,
            "dxTN": -0.43085908629424047,
            "dyTN": 30.476954579612407,
            "point": {
                "x": 1000000.0000007066,
                "y": 9000030.4821516,
                "z": -6.123233995736765E-15
            },
            "wgs84Longitude": -88.93079120090563,
            "wgs84Latitude": 24.793020386563672,
            "dls": 0.0,
            "original": true,
            "dz": 6.123233995736765E-15
        }
    ],
    "localCRS": "{\"lateBoundCRS\":{\"name\":\"Azimuthal Equidistant\",\"type\":\"LBC\",\"ver\":\"PE_10_9_1\",\"wkt\":\"PROJCS[\\\"Azimuthal Equidistant Lng=-88.93089504;Lat=24.79259766\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],PROJECTION[\\\"Modified Azimuthal_Equidistant\\\"],PARAMETER[\\\"False_Easting\\\",0.0],PARAMETER[\\\"False_Northing\\\",0.0],PARAMETER[\\\"Central_Meridian\\\",-88.93089504402704],PARAMETER[\\\"Latitude_Of_Origin\\\",24.792597655454724],UNIT[\\\"Foot_US\\\",0.3048006096012192]]\"},\"name\":\"Azimuthal Equidistant - NAD_1927_To_WGS_1984_79_CONUS\",\"singleCT\":{\"authCode\":{\"auth\":\"EPSG\",\"code\":\"15851\"},\"name\":\"NAD_1927_To_WGS_1984_79_CONUS\",\"type\":\"ST\",\"ver\":\"PE_10_9_1\",\"wkt\":\"GEOGTRAN[\\\"NAD_1927_To_WGS_1984_79_CONUS\\\",GEOGCS[\\\"GCS_North_American_1927\\\",DATUM[\\\"D_North_American_1927\\\",SPHEROID[\\\"Clarke_1866\\\",6378206.4,294.9786982]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],GEOGCS[\\\"GCS_WGS_1984\\\",DATUM[\\\"D_WGS_1984\\\",SPHEROID[\\\"WGS_1984\\\",6378137.0,298.257223563]],PRIMEM[\\\"Greenwich\\\",0.0],UNIT[\\\"Degree\\\",0.0174532925199433]],METHOD[\\\"NADCON\\\"],PARAMETER[\\\"Dataset_conus\\\",0.0],OPERATIONACCURACY[5.0],AUTHORITY[\\\"EPSG\\\",15851]]\"},\"type\":\"EBC\",\"ver\":\"PE_10_9_1\"}",
    "method": "AzimuthalEquidistant",
    "operationsApplied": [
        "derived TN from GN azimuth by grid convergence 359.190052",
        "computed deflections via minimum curvature method",
        "computation method: AzimuthalEquidistant",
        "conversion from 'Azimuthal Equidistant' to 'GCS_North_American_1927'",
        "conversion from 'GCS_North_American_1927' to 'NAD_1927_BLM_Zone_16N'",
        "to WGS 84: conversion from NAD_1927_BLM_Zone_16N to GCS_North_American_1927; 2 points converted",
        "to WGS 84: transformation GCS_North_American_1927 to GCS_WGS_1984 using NAD_1927_To_WGS_1984_79_CONUS; 2 points successfully transformed"
    ],
    "inputKind": "MD_Incl_Azim"
}
Edited by Puneet Bhardwaj