股票翻转坐标有什么用？

股票翻转坐标有什么用

接下来具体说说

绕任意转轴的旋转矩阵与指数坐标解释

绕任意转轴的旋转矩阵可以使用罗德里格斯公式（Rodrigues' formula）来表示。该公式将旋转向量和旋转角度**为旋转矩阵。

假设旋转向量为v，旋转角度为θ，单位向量k表示旋转轴方向，则旋转矩阵R可以通过以下公式计算：

R = I + sin(θ)K + (1 - cos(θ))K^2

其中，I是单位矩阵，K是反对称矩阵，定义为：

| 0 -kz ky |

K = | kz 0 -kx || -ky kx 0 |

指数坐标（Exponential Coordinates）是描述刚体在空间中的位置和姿态的一种方法。它使用旋转矩阵和平移向量来表示刚体的变换。

对于旋转矩阵R和平移向量t，指数坐标可以表示为：

[ R t ][ 0 1 ]

其中，R是旋转矩阵，t是平移向量。

指数坐标可以通过指数映射（Exponential Map）和对数映射（Logarithmic Map）与旋转矩阵之间进行转换。指数映射将旋转矩阵映射到对应的指数坐标，而对数映射则将指数坐标映射回旋转矩阵。

这些概念和公式在刚体运动学和机器人学中经常使用，用于描述和计算刚体的姿态和位姿

在Python中，你可以使用NumPy库来进行矩阵运算和向量操作。下面是一个示例代码，演示如何计算绕任意转轴的旋转矩阵以及指数坐标：

import numpy as np def rotation_matrix (axis, theta) : """ 计算绕任意轴的旋转矩阵 :param axis: 旋转轴的单位向量 :param theta: 旋转角度（弧度） :return: 旋转矩阵 """ axis = axis / np.linalg.norm(axis) # 将轴向量单位化 a = np.cos(theta / 2.0 ) b, c, d = -axis * np.sin(theta / 2.0 ) rotation_mat = np.array([[a*a + b*b - c*c - d*d, 2 * (b*c - a*d), 2 * (b*d + a*c)], [ 2 * (b*c + a*d), a*a + c*c - b*b - d*d, 2 * (c*d - a*b)], [ 2 * (b*d - a*c), 2 * (c*d + a*b), a*a + d*d - b*b - c*c]]) return rotation_mat def exponential_coordinates (rotation_mat) : """ 计算旋转矩阵的指数坐标 :param rotation_mat: 旋转矩阵 :return: 指数坐标 """ theta = np.arccos((np.trace(rotation_mat) - 1 ) / 2.0 ) axis = np.array([rotation_mat[ 2 , 1 ] - rotation_mat[ 1 , 2 ], rotation_mat[ 0 , 2 ] - rotation_mat[ 2 , 0 ], rotation_mat[ 1 , 0 ] - rotation_mat[ 0 , 1 ]]) axis = axis / ( 2 * np.sin(theta)) return axis, theta # 示例使用 axis = np.array([ 1 , 1 , 1 ]) # 旋转轴 axis = axis / np.linalg.norm(axis) # 单位化 theta = np.pi / 2 # 旋转角度（弧度） # 计算绕任意轴的旋转矩阵 rotation_mat = rotation_matrix(axis, theta)print( "Rotation Matrix:" )print(rotation_mat) # 计算旋转矩阵的指数坐标 exponential_coords = exponential_coordinates(rotation_mat)print( "Exponential Coordinates:" )print( "Axis:" , exponential_coords[ 0 ])print( "Theta:" , exponential_coords[ 1 ])

这段代码定义了两个函数，rotation_matrix用于计算绕任意轴的旋转矩阵，exponential_coordinates用于计算旋转矩阵的指数坐标。你可以通过传入旋转轴和旋转角度来获取旋转矩阵和指数坐标。

请注意，这里的旋转角度是以弧度为单位。另外，为了进行矩阵运算，我们使用了NumPy库来进行向量的归一化、矩阵的乘法等操作。如果你还没有安装NumPy，可以使用以下命令进行安装：

pip install numpy

希望对你有帮助！如果你有任何问题，请随时提问。

以下是使用C++编写的计算绕任意转轴的旋转矩阵和指数坐标的示例代码：

# include <iostream> # include <cmath> # include <Eigen/Dense> Eigen:: Matrix3d rotationMatrix ( const Eigen::Vector3d& axis, double theta) { Eigen::Matrix3d rotation; double cosTheta = std :: cos (theta); double sinTheta = std :: sin (theta); double x = axis.x(); double y = axis.y(); double z = axis.z(); rotation << cosTheta + ( 1 - cosTheta) * x * x, ( 1 - cosTheta) * x * y - sinTheta * z, ( 1 - cosTheta) * x * z + sinTheta * y, ( 1 - cosTheta) * x * y + sinTheta * z, cosTheta + ( 1 - cosTheta) * y * y, ( 1 - cosTheta) * y * z - sinTheta * x, ( 1 - cosTheta) * x * z - sinTheta * y, ( 1 - cosTheta) * y * z + sinTheta * x, cosTheta + ( 1 - cosTheta) * z * z; return rotation;}Eigen:: Matrix4d exponentialCoordinates ( const Eigen::Matrix3d& rotation, const Eigen::Vector3d& translation) { Eigen::Matrix4d exponential; exponential.block< 3 , 3 >( 0 , 0 ) = rotation; exponential.block< 3 , 1 >( 0 , 3 ) = translation; exponential.block< 1 , 4 >( 3 , 0 ) << 0 , 0 , 0 , 1 ; return exponential;} int main () { Eigen:: Vector3d axis ( 1.0 , 0.0 , 0.0 ) ; // 旋转轴的单位向量 double theta = M_PI / 2.0 ; // 旋转角度（弧度） Eigen::Matrix3d rotation = rotationMatrix(axis, theta); std :: cout << "Rotation Matrix:\n" << rotation << std :: endl ; Eigen:: Vector3d translation ( 1.0 , 2.0 , 3.0 ) ; // 平移向量 Eigen::Matrix4d exponential = exponentialCoordinates(rotation, translation); std :: cout << "Exponential Coordinates:\n" << exponential << std :: endl ; return 0 ;}

上述代码使用了Eigen库来进行矩阵运算和向量操作。它定义了两个函数，rotationMatrix用于计算绕任意轴的旋转矩阵，exponentialCoordinates用于计算指数坐标。在main函数中，我们指定旋转轴和角度，然后计算旋转矩阵和指数坐标，并将结果打印输出。

请确保在编译和运行代码之前，已经安装了Eigen库，并将其包含在代码中。

绕任意转轴的旋转矩阵和指数坐标在许多领域和场景中都有广泛的应用，特别是在计算机图形学、机器人学、物体姿态估计和模拟等领域。

以下是一些应用场景的示例：

1. 计算机图形学：在三维计算机图形学中，通过旋转矩阵可以实现物体的旋转变换，从而实现模型的旋转、动画和视角变换等效果。

2. 机器人学：在机器人学中，旋转矩阵和指数坐标常用于描述机器人的姿态和运动。通过旋转矩阵可以表示机器人末端执行器相对于基座的旋转姿态，而指数坐标可以表示机器人的运动轨迹。

3. 物体姿态估计：在计算机视觉和机器人感知中，通过旋转矩阵和指数坐标可以估计物体的姿态和旋转角度。这对于物体识别、跟踪和位姿估计等任务非常重要。

4. 模拟和仿真：在物理仿真和虚拟现实等领域，旋转矩阵和指数坐标可用于模拟物体的旋转和运动。通过对物体施加旋转矩阵，可以模拟物体的旋转行为，并将其用于物理仿真和虚拟环境中。

各类坐标转换不完全说明

01

定义

首先在讲述本文前，需要明确地理坐标系和投影坐标系的区别。两者定义以及联系如下：

地理坐标系： 是使用三维球面来定义地球表面位置，是基于球体或旋转椭球体的坐标系。使用基于经纬度的坐标系统来描述地球上某个点所处的位置。

投影坐标系： 在二维平面中进行定义。与地理坐标系不同，使用基于X，Y值的坐标系统来描述地球上某个点所处的位置。投影坐标系始终基于地理坐标系，而后者则是基于球体或旋转椭球体的。

两者联系： 投影坐标系始终基于地理坐标系，它对应于某个唯一的地理坐标系，通常地理坐标系和投影坐标系是一对多的关系。而地理坐标系是基于球体或旋转椭球体的。

02

火星、百度、高德等坐标介绍

其次我们要搞懂什么是 火星坐标系 ，百度以及高德等互联网地图采用的坐标系跟火星坐标系的关系。

GCJ-02 坐标系，即火星坐标系。 这是由国测局02年发布的坐标系，所以又称为国测局坐标，它是对经纬度数据加入随机的偏差进行加密的算法。互联网地图都是基于GCJ-02进行一次或者多次加密的，不允许使用WGS-84坐标。

BD09坐标系： 百度地图采用的百度坐标系，在GCJ-02基础上二次加密而成。高德、腾讯都是使用GCJ-02 坐标系。

03

GCJ-02 坐标系、BD09坐标系和WGS84 坐标系互转

WGS84 坐标系： Google地图使用，也是大多数科研人员使用到的地理坐标系。

本文参考GitHub的代码，在此基础上改进一点点。实现GCJ-02 坐标系、BD09坐标系和WGS84 坐标系之间的互转

使用前记得改路径，代码如下：

1. # - *- coding: utf-8 -* - 2. import math 3. import pandas as pd 4. 5. x_pi = 3.14159265358979324 * 3000.0 / 180.0 6. pi = 3.1415926535897932384626 # π 7. a = 6378245.0 # 长半轴 8. ee = 0.00669342162296594323 # 偏心率平方 9. 10. 11. def gcj02 _to_ bd09(lng, lat): 12. """ 13. 火星坐标系(GCJ-02)转百度坐标系(BD-09) 14. 谷歌、高德——>百度 15. :param lng:火星坐标经度 16. :param lat:火星坐标纬度 17. :return: 18. """ 19. z = math.sqrt(lng * lng + lat * lat) + 0.00002 * math.sin(lat * x_pi) 20. theta = math.atan2(lat, lng) + 0.000003 * math.cos(lng * x_pi) 21. bd_lng = z * math.cos(theta) + 0.0065 22. bd_lat = z * math.sin(theta) + 0.006 23. return [bd _lng, bd_ lat] 24. 25. 26. def bd09 _to_ gcj02(bd _lon, bd_ lat): 27. """ 28. 百度坐标系(BD-09)转火星坐标系(GCJ-02) 29. 百度——>谷歌、高德 30. :param bd_lat:百度坐标纬度 31. :param bd_lon:百度坐标经度 32. :return:转换后的坐标列表形式 33. """ 34. x = bd_lon - 0.0065 35. y = bd_lat - 0.006 36. z = math.sqrt(x * x + y * y) - 0.00002 * math.sin(y * x_pi) 37. theta = math.atan2(y, x) - 0.000003 * math.cos(x * x_pi) 38. gg_lng = z * math.cos(theta) 39. gg_lat = z * math.sin(theta) 40. return [gg _lng, gg_ lat] 41. 42. 43. def wgs84 _to_ gcj02(lng, lat): 44. """ 45. WGS84转GCJ02(火星坐标系) 46. :param lng:WGS84坐标系的经度 47. :param lat:WGS84坐标系的纬度 48. :return: 49. """ 50. if out _of_ china(lng, lat): # 判断是否在国内 51. return [lng, lat] 52. dlat = _transformlat(lng - 105.0, lat - 35.0) 53. dlng = _transformlng(lng - 105.0, lat - 35.0) 54. radlat = lat / 180.0 * pi 55. magic = math.sin(radlat) 56. magic = 1 - ee * magic * magic 57. sqrtmagic = math.sqrt(magic) 58. dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi) 59. dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi) 60. mglat = lat + dlat 61. mglng = lng + dlng 62. return [mglng, mglat] 63. 64. 65. def gcj02 _to_ wgs84(lng, lat): 66. """ 67. GCJ02(火星坐标系)转GPS84 68. :param lng:火星坐标系的经度 69. :param lat:火星坐标系纬度 70. :return: 71. """ 72. 73. dlat = _transformlat(lng - 105.0, lat - 35.0) 74. dlng = _transformlng(lng - 105.0, lat - 35.0) 75. radlat = lat / 180.0 * pi 76. magic = math.sin(radlat) 77. magic = 1 - ee * magic * magic 78. sqrtmagic = math.sqrt(magic) 79. dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi) 80. dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi) 81. mglat = lat + dlat 82. mglng = lng + dlng 83. return [lng * 2 - mglng, lat * 2 - mglat] 84. 85. 86. def bd09 _to_ wgs84(bd _lon, bd_ lat): 87. lon, lat = bd09 _to_ gcj02(bd _lon, bd_ lat) 88. return gcj02 _to_ wgs84(lon, lat) 89. 90. 91. def wgs84 _to_ bd09(lon, lat): 92. lon, lat = wgs84 _to_ gcj02(lon, lat) 93. return gcj02 _to_ bd09(lon, lat) 94. 95. 96. def _transformlat(lng, lat): 97. ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + \ 98. 0.1 * lng * lat + 0.2 * math.sqrt(math.fabs(lng)) 99. ret += (20.0 * math.sin(6.0 * lng * pi) + 20.0 * 100. math.sin(2.0 * lng * pi)) * 2.0 / 3.0 101. ret += (20.0 * math.sin(lat * pi) + 40.0 * 102. math.sin(lat / 3.0 * pi)) * 2.0 / 3.0 103. ret += (160.0 * math.sin(lat / 12.0 * pi) + 320 * 104. math.sin(lat * pi / 30.0)) * 2.0 / 3.0 105. return ret 106. 107. 108. def _transformlng(lng, lat): 109. ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + \ 110. 0.1 * lng * lat + 0.1 * math.sqrt(math.fabs(lng)) 111. ret += (20.0 * math.sin(6.0 * lng * pi) + 20.0 * 112. math.sin(2.0 * lng * pi)) * 2.0 / 3.0 113. ret += (20.0 * math.sin(lng * pi) + 40.0 * 114. math.sin(lng / 3.0 * pi)) * 2.0 / 3.0 115. ret += (150.0 * math.sin(lng / 12.0 * pi) + 300.0 * 116. math.sin(lng / 30.0 * pi)) * 2.0 / 3.0 117. return ret 118. 119. if __name__ == ' __main__ ': 120. # 读取csv格式数据 121. df = pd.read_csv('./ original.csv') 122. # 新建两列用于储存转换后的坐标 123. df['lng_new'] = 0.000000 124. df['lat_new'] = 0.000000 125. for i in range(0, len(df)): 126. # 根据需要选择转换方法，这里以百度坐标转国测局坐标 127. df[ 'lng_new' ][ i ] = bd09 _to_ gcj02(df[ 'lng' ][ i ], df[ 'lat' ][ i ])[0] 128. df[ 'lat_new' ][ i ] = bd09 _to_ gcj02(df[ 'lng' ][ i ], df[ 'lat' ][ i ])[1] 129. # 将转换后的数据另存为新文件 130. df.to_csv('./new.csv')

04

结果示意

原始数据：

以上就是股票翻转坐标有什么用？的详细内容，希望通过阅读小编的文章之后能够有所收获！