Liên quan đến javascript sang SQL, đây có lẽ là cách bạn sẽ xử lý việc đó:
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
Tôi nghĩ rằng sau đây sẽ trả lời câu hỏi đầu tiên của bạn. Nó sẽ yêu cầu khá nhiều kiểm tra lỗi. Để hỗ trợ, bạn có thể tìm phương trình ban đầu tại đây: 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
Chỉnh sửa: một vài biến nên là tên cột và dấu phẩy và dấu ngoặc đơn bị thiếu.
Chỉnh sửa: thêm một dấu ngoặc đơn.
Tôi đã thử nghiệm công thức này và nó hoạt động bằng cách sử dụng các điểm ngẫu nhiên chống lại biến đổi ArcGIS. Hãy nhớ rằng các đơn vị của bạn có thể được tính bằng feet / độ. Ngoài ra, hãy nhớ những kết quả này là đồng bằng, vì vậy bạn sẽ phải thêm chúng vào giá trị của bạn để có được kết quả cuối cùng.