Comment effectuer une géotransformation et une projection à 3 paramètres dans Sql Server 2012?

J'ai un tableau avec des colonnes de longitude latitude (NAD27). Je calcule deux autres colonnes, X et Y, représentant l'emplacement de Web Mercator (WGS84).

Actuellement, j'utilise un Arcmap pour ce faire, en appliquant la géotransformation recommandée pour la zone d'étude - la géotransformation à 3 paramètres (géocentrique) - pour passer de NAD27 à WGS84.

Je voudrais le faire entièrement dans Sql Server 2012. D'après ce que je peux dire, Sql Server ne prend pas en charge les transformations de données prêtes à l'emploi. Quelqu'un connaît-il une bibliothèque Sql qui prend en charge cette géotransformation? Je voudrais simplement utiliser les mêmes coefficients en SQL que j'utilise actuellement en Arcmap.

J'ai également besoin de projeter de WGS84 lat / long dans web mercator. Je vois cette formule implémentée en javascript , mais si quelqu'un a une procédure stockée SQL qui fait cela, ce serait génial.

En ce qui concerne le javascript à SQL, voici probablement comment vous géreriez cela:

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

Je pense que ce qui suit répondra à votre première question. Cela nécessitera un peu de vérification des erreurs. Pour vous aider, vous pouvez trouver l'équation originale ici: 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

Edit: quelques variables qui auraient dû être des noms de colonnes, et une virgule et des parenthèses manquantes.

Edit: une parenthèse de plus.

J'ai testé cette formule et elle fonctionne en utilisant des points aléatoires contre la transformation ArcGISs. N'oubliez pas que vos unités peuvent être en pieds / degrés. N'oubliez pas que ces résultats sont des deltas, vous devrez donc les ajouter à vos valeurs pour obtenir vos résultats finaux.

Ceci est un lien vers une question similaire:

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

Je pense que la première réponse peut être utile pour comprendre ce que vous pouvez et ce que vous ne pouvez pas faire dans SQL Server et pour connaître certaines méthodes pour résoudre votre problème.