武汉大学定量遥感作业
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
■绘制图1.4 本地观测天顶角和传感器扫描角的数值对比
1.数据处理
① 已知R 0=6378km,选取SPOT 卫星轨道高度h=832km,则
sin θ=
R0+ℎR0
sin α可以化为： sin θ=1.13 sin α
② 在Excel 中用函数=ROUND((RAND()*num ),0)&"°"生成随机α值30个，
利用=SIN(RADIANS())计算出sin α； 利用sinθ=1.13 sinα换算出sinθ； 利用=DEGREES(ASIN())计算出θ值。
2.绘图
观测天顶角和传感器扫描角的数值表
No. α sinα sinθ θ 1 11 0.19 0.22 12.45 2 2 0.03 0.04 2.26 3 6 0.10 0.12 6.78 4 60 0.87 0.98 78.13 5 56 0.83 0.94 69.52 6 1 0.02 0.02 1.13 7 10 0.17 0.20 11.32 8 49 0.75 0.85 58.52 9 26 0.44 0.50 29.69 10 26 0.44 0.50 29.69 11 51 0.78 0.88 61.42 12 59 0.86 0.97 75.60 13 59 0.86 0.97 75.60 14 7 0.12 0.14 7.92 15 50 0.77 0.87 59.95 16 28 0.47 0.53 32.04 17 56 0.83 0.94 69.52 18 59 0.86 0.97 75.60 19 56 0.83 0.94 69.52 20 56 0.83 0.94 69.52 21 36 0.59 0.66 41.62 22 5 0.09 0.10 5.65 23 34 0.56 0.63 39.19 24 17 0.29 0.33 19.29
25 47 0.73 0.83 55.73 26 29 0.48 0.55 33.22 27 36 0.59 0.66 41.62 28 22 0.37 0.42 25.04 29 36 0.59 0.66 41.62 30 22 0.37 0.42 25.04
R 绘图代码如下：
>data<-read.table("clipboard", header = T) > data
>library(ggplot2)
>p<-ggplot(data=data, mapping=aes(x=α, y=θ))
>p+stat_smooth()
>p+stat_smooth()+labs(y='观测天顶角/(°)',x='传感器扫描角度/(°)'，title ='图1.4 本地观测天顶角和传感器扫描角的数值对比')
>p+stat_smooth()+geom_abline(linetype=1，colour=2)
图像说明：从本地观测天顶角和传感器扫描角的数值对比图可以看出，很明显，地球表面一点的观测天顶角总是比传感器扫描角大，并且传感器扫描角越大，它们之间的差就越大。