- 1、本文档共3页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
如何绘制日食界限图
如何绘制日食界限图?
张金龙
Jinlongzhang01@
日食界限图是表示日食可见区域的一种地图。包括日食的南北界限,日出日落时刻线,日出日落时初亏线,日出日落时复圆线。如2011年1月4日将发生一次日偏食,在欧洲大部分,非洲撒哈拉及其以北地区,中亚,西亚以及我国的西部可以见到日偏食。可见区域由若干条曲线围成。各曲线代表的含义分别为,
I, 日出时复圆。
II, 日出时食甚
III, 日出时初亏
IV, 日落时复圆
V, 日落时食甚
VI, 日落时初亏
还有一条线段表示日食南界。在这个区域内部,可以见到日偏食的完整过程。区域之外则不可见。由于月亮自西向东运行,其速度远远超过太阳的向东运行的速度。因此,日食总是从西侧开始,在东侧结束。
公元前2000年到公元3000年每一次日食kmz文件, 可以在 http://xjubier.free.fr/en/site_pages/SolarEclipsesGoogleEarth.html 下载。kmz文件可以被google earth读取。
为了生成清晰度更高的日食界限图,笔者编写了相应的R函数,以绘制日食界限图。该函数基于maptools、mapproj、mapdata。使用前请先下载这三个程序包。
由于绘制地图必须采用一定的投影方式,这里默认的投影方式为斜轴等面积投影。当然,高手可以通过阅读R的帮助信息,改变相应的地图投影方式,线段的颜色等,以达到满意的效果。
图2 - 6
plotkml绘制的日食界限图
plotkml <-
function(file, projection = "azequalarea" , parameter = NULL, grid = TRUE, info = TRUE,... ){
require(maptools)
require(mapproj)
require(mapdata)
#match.arg(projection)
zz <- getKMLcoordinates(file, ignoreAltitude = TRUE)
x11(11,9)
m <- map('world',proj= projection, parameter = parameter, orient=c(zz[[1]][2], zz[[1]][1], 0), ...)
for(i in 1:length(zz)){ ## Add lines to graph
if(is.data.frame(zz[[i]])){
if(ncol(zz[[i]]) > 3){
i = i + 1
} ## Avoid ncol > 3
}
if(is.list(zz[[i]])){
i = i + 1 ## Avoid list > 3
}
if(is.vector(zz[[i]])){
points(mapproject(list(y = zz[[i]][2], x= zz[[i]][1])), pch = "*", cex = 2, col = 2)
}
if(is.matrix(zz[[i]])){
lines(mapproject(list(y = zz[[i]][,2], x= zz[[i]][,1])), col = "blue", ...)
}
## print(i) used in debugging
}
if(grid){ ## Add map grid
map.grid(m, lty = 1, col = "grey")
}
file1 <- readLines(file)
if(title){ ## Add title
#file <- "c:/ASE_1800_04_24__4cdcf4f5330fb.kml"
daterecord <- file1[grepl("when", file1)]
date <- substring(gsub(" ", "", daterecord), 7, 16)
type <- file1[grepl("Type:", file1)]
st
您可能关注的文档
- 四川省机动车维修业户.doc
- 团队人员周鑫煜张杰波.doc
- 地点香港大会堂低座音乐厅.doc
- 城市园林绿化企业资质公示-.doc
- 地址邮编天津市河西区珠江道25号300222招生咨询022.doc
- 天任TP300系列在丝网印刷机中的应用.doc
- 孟津县行政审批项目及服务性项目登记表(一).doc
- 天姿盘点婚纱摄影中那些超梦幻的主题.doc
- 学号0145102英班.doc
- 审批签字.doc
- 地理-内蒙古锡林郭勒盟2023~2024学年第一学期高三年级全盟统考试题和答案.docx
- 历史|湖北元月调考暨湖北省部分市州元月高三期末联考历史试卷及答案.docx
- 历史-内蒙古锡林郭勒盟2023~2024学年第一学期高三年级全盟统考试题和答案.docx
- 数学(理科)-内蒙古锡林郭勒盟2023~2024学年第一学期高三年级全盟统考试题和答案.docx
- 历史丨金太阳24~239C湖北省十堰市2024届高三上学期1月调研考试历史试卷及答案.docx
- 青岛某高层办公楼地下防水施工方案(sbs卷材).docx
- 土木工程实习报告表格.doc
- 工程流体力学教学课件ppt作者闻建龙工程流体力学习题答案部分.doc
- 围手术期管理制度(1)_0.doc
- 2023年公务员多省联考《申论》题(内蒙古旗县卷)及参考答案.doc
文档评论(0)