python计算方位角_利用python计算太阳天顶角、方位角、高度角

在遥感计算中一般都会用到天顶角、方位角、高度角。之前都是直接在excel中输入公式，这种方式输入公式比较麻烦，而且容易出错。后来在网上看到吉林大学汪自军博士的计算程序。链接：http://blog.sciencenet.cn/home.php?mod=space&uid=43777&do=blog&id=238552

非常好用，但是我看了一下，汪博士主要用matlab、C++、FORTRAN三种语言编写的，而且输入文件是txt格式。txt文件的操作性我认为不如excel。

说一下主要改进的地方：

1、编程语言改为python

2、输入和输出文件格式均改为excel

3、以为太阳高度角与天顶角互余，这次也加进去了。

欢迎讨论，下面附上主要代码：

table = data.sheets()[0] # 打开第一张表

nrows = table.nrows # 获取表的行数

ncols = table.ncols

station=table.col_values(0)[1:]

year=table.col_values(1)[1:]

month=table.col_values(2)[1:]

day=table.col_values(3)[1:]

hour=table.col_values(4)[1:]

min=table.col_values(5)[1:]

sec=table.col_values(6)[1:]

lon=table.col_values(7)[1:]

lat =tab