博文
如何绘制日食路线图(界限图):附R程序？
||||
如何绘制日食界限图-路线图(附R程序)？
张金龙
日食界限图是表示日食可见区域的一种地图。包括日食的南北界限，日出日落时刻线，日出日落时初亏线，日出日落时复圆线。如2011年1月4日将发生一次日偏食，在欧洲大部分，非洲撒哈拉及其以北地区，中亚，西亚以及我国的西部可以见到日偏食。可见区域由若干条曲线围成。
图1 日食界限图
各曲线代表的含义分别为，
I，日出时复圆。
II，日出时食甚
III， 日出时初亏
IV，日落时复圆
V，日落时食甚
VI， 日落时初亏
还有一条线段表示日食南界。在这个区域内部，可以见到日偏食的完整过程。区域之外则不可见。
由于月亮自西向东运行，其速度远远超过太阳的向东运行的速度。因此，日食总是从西侧开始，在东侧结束。
公元前2000年到公元3000年每一次日食kmz文件， 可以在 http://xjubier.free.fr/en/site_pages/SolarEclipsesGoogleEarth.html下载。kmz文件可以被google earth读取。
为了生成清晰度更高的日食界限图，笔者编写了相应的R函数plotkml，以绘制日食界限图。
点击下载 如何绘制日食路线图R程序
该函数基于maptools、mapproj、mapdata。使用前请先下载这三个程序包。由于绘制地图必须采用一定的投影方式，这里默认的投影方式为斜轴等面积投影。当然，高手可以通过阅读R的帮助信息，改变相应的地图投影方式，线段的颜色等，以达到满意的效果。
图2 - 6
plotkml绘制的日食界限图
图2 2009年7月22日日全食路线图
图3 2020年6月21日日环食
图4 2021年6月10日日环食
图5 2024年10月2日日环食
图6 2026年2月17日日环食
图7 2026年8月12日日全食
图8 2023年4月20日混合食（全环食）
图9 2035年9月2日日全食
图10 2024年4月20日日全食
https://m.sciencenet.cn/blog-255662-384353.html
上一篇：如何计算日食的时间和类型（附R源代码）？
下一篇：Phylomatic生成进化树及物种功能性状聚类
