制作天文软件，总是要根坐标旋转打交道，这两天，转得有点晕。界线图计算及输出，要在赤道坐标、贝塞尔坐标、视图正投影坐标、最大食为中心的坐标，以及一些中间坐标（如直角坐标、各种交线起算的坐标）之间转来转去，搞得很麻烦。总希望能一步就转到位，可总是不顺意。

谢谢许老师能提供这么好的程序并且能公开代码，也谢谢浪-淘-沙兄的指点。

我主要还是想重现星像，其他的应用因需而定，解决的思路还是跟浪兄的相似。

目前只是业余兴趣，顺便通过这个课题学习一些知识，编程、天文学、数学等方面的知识，其实从许老师的贴子中也摸索到了正确的方向和方法，而且许老师还提供了代码，在下已经非常感谢了，不敢再打扰老师。未来如果有这方面的代码提供，在下密切关注。

这种描点法得到的地图，确实会出现非线性缩放的问题。

还有，我用的地图像素如果与你用的不是相同大小的，我描点得到的数据，可能无法直接使用。需要一个变换比例尺。

 

所以我在想，最好使用的经纬度坐标，不是像素坐标。

有了经纬度坐标，不管是什么地图，都可以按比例投影到画面上。

这样，不管将来的地图如何放大，边界上每个点的位置都是固定的。

 

地图有各种投影展开法。采用象素法，非线性问题有点难解决。

我用的是时区图，经度是线性的。纬度是非线性的。我上面讲的像素描点法，正是用像素坐标表示经纬度，二者是线性的比例关系。

 

如果你描点的点是经纬度数据，基本可以直接使用，只是一个起算点的问题，只需在vml.js程序中改动下面这个语句：

 

     J*=2*Math.PI/2009, W=Math.PI/2-W*Math.PI/970;

 

这句正是比例变换语句

 

 

我知道许兄近日来为了日食南界和地图的算法忙得非常辛苦，常常忙到半夜。我以下这几点建议您可以在以后考虑一下是否增加。

本人对《寿星天文历》的“天象”部分有以下几点建议：

对于南极地图，因为展开的面积比较小，且以南极点为中心对称的。

经线是直的，纬线正好是同心圆。

采用极坐标来描点，很容易得到边界的经纬度。

所以要得到经纬度数据还是比较容易的。

 

其它地方的经纬度数据，用描点法都会产生非线性变形的。（经线可能是弯的，纬线也可能不是同心圆）。

已发布v4.72版

1、增加食界图按指定的步长进行步进切换，以方便比较不同日食的区别。

2、修正食界连线的Bug，因切入时序的判别条件有错误，造成个别食界图连线失常

3、中心线、本影南北界线、半影南北界线使用椭圆求解，以提高精度，将来有更精确地图，就很有必要提高中心线及本影精度。（不过这里仍采用中精度星历）

4、南北界计算中加入地球自转速度(地表速度)对观测点的切点计算的影响

5、改正地表速度的一个Bug，原选对轻微超出地面的地方，地表速度置为0，后发现这样会造成速度场不连续，影响迭代计算，导至半影南北界在地球边缘时有轻微的不连续

 

到今天为止，日食计算基本完整。当然，还有一些计算过程的中间数据报表没有输出（比如中心线、南北界的各时刻位置数据表、日食根数表、最初切入点、影轴切入点数据等等），这些数据基时在计算过程中都已在内存中产生，但考虑到一般用户看不懂这些数据，暂不提供。

 

凌日等天象，确实很重要。不过得等有空的时候来编写。这几天一直在计算日月食。

 

日月食、凌日，一般需要图形来支持，一想到作图，就害怕，麻烦得不得了。

我利用AutoCAD描南极洲地图，东经为负数，西经为正数，南纬为负数。

利用ATAN2函数得到经度，再判断比较计算出正（西）、负（东）数的。

 

请许兄试着绘制到日食界限图中。

如果所需要格式不合要求，请提出改正方案。

 

南极有几个小岛，数据是单独的。

 

我不知道m l是如何加到数据中的，怕加错了，引起麻烦。请许兄自行添加。

按许兄前面的示范，可能要删除结尾的一个逗号。

也请许兄辛苦一下吧。

数据见打包文件。

 

其中一个小岛的数据：

161.341,-79.951,
161.681,-79.964,
162.48,-79.924,
163.639,-79.782,
164.451,-79.565,
164.724,-79.229,
164.359,-79,
163.63,-78.784,
162.835,-78.786,
162.338,-79.003,
162.145,-79.339,
161.443,-79.58,
160.835,-79.735,
160.8,-79.876,
161.341,-79.951,

下载信息  [文件大小：   下载次数： ]

点击浏览该文件:南极洲地图数据表6.txt

你的精度比我原来的高出十倍以上，我原来的差不多只有1度(最多只有0.5度)的精度，以致于我都没有动力制作精密的日食界线图。看来有必要做一个地图的放缩功能，以显示出更精密的日食情况。

 

我今天下乡检查工作刚回来，要连续下乡一星期，所以这星期不会很有空，今天回来看到这个数据表，一个惊喜。我以前，搞这个地图，花了好多时间，却搞不出一个象样的。浪版主是专业级的，谢谢。

 

等一下，我做一个地图数据接口附助程序，方便你加入新的数据，以便调试。