Wie führe ich eine 3-Parameter-Geotransformation und -Projektion in SQL Server 2012 durch?

Ich habe eine Tabelle mit Spalten zum Längen- und Breitengrad (NAD27). Ich berechne zwei weitere Spalten, X und Y, die den Standort von Web Mercator (WGS84) darstellen.

Derzeit verwende ich dazu eine Arcmap, indem ich die empfohlene Geotransformation für das Untersuchungsgebiet - die 3-Parameter-Geotransformation (geozentrisch) - anwende, um von NAD27 zu WGS84 zu gelangen.

Ich möchte dies vollständig in SQL Server 2012 tun. Soweit ich weiß, unterstützt SQL Server keine sofort einsatzbereiten Datumstransformationen. Kennt jemand eine SQL-Bibliothek, die diese Geotransformation unterstützt? Ich möchte einfach die gleichen Koeffizienten in SQL verwenden, die ich derzeit in Arcmap verwende.

Ich muss auch von WGS84 lat / long in Web Mercator projizieren. Ich sehe diese Formel in Javascript implementiert , aber wenn jemand eine gespeicherte SQL-Prozedur hat, die dies tut, wäre es großartig.

In Bezug auf das Javascript zu SQL würden Sie wahrscheinlich so damit umgehen:

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

Ich denke, das Folgende wird Ihre erste Frage beantworten. Es wird einige Fehlerprüfungen erfordern. Zur Unterstützung finden Sie die ursprüngliche Gleichung hier: 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

Bearbeiten: ein paar Variablen, die Spaltennamen sein sollten, und ein fehlendes Komma und eine fehlende Klammer.

Bearbeiten: noch eine Klammer.

Ich habe diese Formel getestet und sie funktioniert mit zufälligen Punkten gegen die ArcGIS-Transformation. Denken Sie daran, dass Ihre Einheiten möglicherweise in Fuß / Grad angegeben sind. Denken Sie auch daran, dass es sich bei diesen Ergebnissen um Deltas handelt. Sie müssen sie daher zu Ihren Werten addieren, um die endgültigen Ergebnisse zu erhalten.

Dies ist ein Link zu einer ähnlichen Frage:

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

Ich denke, dass die erste Antwort hilfreich sein kann, um zu verstehen, was Sie in SQL Server können und was nicht, und um einige Methoden zur Lösung Ihres Problems zu kennen.