¿Cómo realizo una geotransformación y proyección de 3 parámetros en Sql Server 2012?

Tengo una tabla con columnas de latitud longitud (NAD27). Calculo otras dos columnas, X e Y, que representan la ubicación de Web Mercator (WGS84).

Actualmente estoy usando un Arcmap para hacer esto, aplicando la geotransformación recomendada para el área de estudio, la geotransformación de 3 parámetros (geocéntrica) , para pasar de NAD27 a WGS84.

Me gustaría hacer esto completamente dentro de Sql Server 2012. Por lo que puedo decir, Sql Server no admite transformaciones de datos de fábrica. ¿Alguien sabe de una biblioteca SQL que admite esta geotransformación? Me gustaría simplemente usar los mismos coeficientes en SQL que estoy usando actualmente en Arcmap.

También necesito proyectar desde WGS84 lat / long en web mercator. Veo esta fórmula implementada en javascript , pero si alguien tiene un procedimiento almacenado SQL que hace esto, sería genial.

Con respecto al javascript para SQL, esta es probablemente la forma en que manejarías eso:

SELECT  FromX, 
        FromY, 
        CASE WHEN FromX > 180 THEN NULL ELSE FromX * 0.017453292519943295 * 6378137.0 END AS mercatorX_lon2,
        CASE WHEN FromY > 90 THEN NULL ELSE 3189068.5 * LOG((1.0 + SIN(FromY * 0.017453292519943295)) / (1.0 - SIN(FromY * 0.017453292519943295))) END AS mercatorY_lat2
FROM TABLENAME

Creo que lo siguiente responderá a su primera pregunta. Requerirá bastante verificación de errores. Para ayudar, puede encontrar la ecuación original aquí: http://www.colorado.edu/geography/gcraft/notes/datum/gif/molodens.gif

--fromTheta :column --radians
--fromLamda :column --radians
--fromH     :column --meters

DECLARE @fromA float = 6378206.4        --radius of earth, meters
DECLARE @fromF float =1.0/294.9786982   --Flattening
DECLARE @toA float =6378137.0           --radius of earth, meters
DECLARE @toF float = 1.0/298.257223563  --Flattening
DECLARE @dA float = @toA - @fromA       --change in equatorial radius
DECLARE @dX float = -8.0                --change in X, meters
DECLARE @dY float = 160.0               --change in Y, meters
DECLARE @dZ float = 176.0               --change in Z, meters
DECLARE @dF float = @toF-@fromF         --change in flattening
DECLARE @fromES float = 2.0*@fromF - @fromF*@fromF --first eccentricity squared
DECLARE @bda float = 1.0-@fromF         --polar radius divided by equatorial radius

--RM = (@fromA*(1-@fromES)/POWER(1-@fromES*sin(fromTheta)*sin(fromTheta), 1.5))

--RN = (@fromA/SQRT(1.00-@fromES*sin(fromTheta)*sin(fromTheta)))

SELECT 

((((-@dX*sin(fromTheta)*cos(fromLamda)-@dY*sin(fromTheta)*sin(fromLamda))+@dZ*cos(fromTheta))+@dA*(@fromA/SQRT(1.00-@fromES*sin(fromTheta)*sin(fromTheta)))*@fromES*sin(fromTheta)*cos(fromTheta)/@fromA)+@df*((@fromA*(1-@fromES)/POWER(1-@fromES*sin(fromTheta)*sin(fromTheta), 1.5))/@bda+(@fromA/SQRT(1.00-@fromES*sin(fromTheta)*sin(fromTheta)))*@bda)*sin(fromTheta)*cos(fromTheta))/((@fromA*(1-@fromES)/POWER(1-@fromES*sin(fromTheta)*sin(fromTheta), 1.5)) + fromH) AS deltaTheta,
(-@dX*sin(fromLamda)+@dY*cos(fromLamda))/((((@fromA/SQRT(1.00-@fromES*sin(fromTheta)*sin(fromTheta))) +fromH) * cos(fromTheta)) AS deltaLamda,
@dX*cos(fromTheta)*cos(fromLamda)+@dY*cos(fromTheta)*sin(fromLamda)+@dZ*sin(fromTheta)-@da*@fromA/(@fromA/SQRT(1.00-@fromES*sin(fromTheta)*sin(fromTheta)))+@dF*@bda*(@fromA/SQRT(1.00-@fromES*sin(fromTheta)*sin(fromTheta)))*sin(fromTheta)*sin(fromTheta) AS deltaH

FROM TABLENAME

Editar: un par de variables que deberían haber sido nombres de columna, y una coma y paréntesis faltantes.

Editar: un paréntesis más.

He probado esta fórmula y funciona utilizando puntos aleatorios contra la transformación de ArcGIS. Recuerde que sus unidades pueden estar en pies / grados. Recuerde también que estos resultados son deltas, por lo que deberá agregarlos a sus valores para obtener sus resultados finales.

Este es un enlace a una pregunta similar:

http://sqlspatialtools.codeplex.com/discussions/286893

Creo que la primera respuesta puede ser útil para comprender lo que puede y lo que no puede hacer en SQL Server y conocer algunos métodos para resolver su problema.