将球面投影到平面上,有多种投影方式。R中的mapproj程序包包含了大部分的投影,适用于地图的投影变换。 而星图中最常用的投影方式为球极平面投影(Stereographic projection), 这种投影方式有时也用于地图中。
这里介绍在R中如何实现球极平面投影。
定义函数stereo_proj
- lambda: 其中lambda为赤经
- phi: 为赤纬
- lambda0 =0: lambda0为地图(星图)的投影中心点 赤经
- phi0 =0: phio为地图(星图)投影中心点的 赤纬
- R =1: 为比例半径, 一般取值为1
函数定义如下:
1 | stereo_proj <-function(lambda,phi, lambda0 =0, phi0 =0, R =1){ |
1 | plot(stereo_proj(lambda = aaa$x, phi = aaa$y), |