本文整理汇总了Python中astral.Astral.solar_azimuth方法的典型用法代码示例。如果您正苦于以下问题：Python Astral.solar_azimuth方法的具体用法？Python Astral.solar_azimuth怎么用？Python Astral.solar_azimuth使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类astral.Astral的用法示例。
在下文中一共展示了Astral.solar_azimuth方法的2个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: testAzimuth
# 需要导入模块: from astral import Astral [as 别名]
# 或者: from astral.Astral import solar_azimuth [as 别名]
def testAzimuth():
city_name = "Jubail"
dd = Astral()
city = dd[city_name]
print("Latitude: %f, Longitude: %f" % (city.latitude, city.longitude))
dt = datetime.datetime.now(tz=city.tz)
print("Date & time: %s" % dt)
print("Date & time (UTC): %s" % dt.astimezone(pytz.utc))
print("Azimuth: %.02f" % dd.solar_azimuth(dt, city.latitude, city.longitude))示例2: calcFiniteSlopes
# 需要导入模块: from astral import Astral [as 别名]
# 或者: from astral.Astral import solar_azimuth [as 别名]
# Convert angular measurements to radians
azRad, elevRad = (360. - az + 90.)*np.pi/180., (90.-elev)*np.pi/180.
Sx, Sy = calcFiniteSlopes(elevGrid, dx) # Calculate slope in X and Y directions
AspectRad = np.arctan2(Sy, Sx) # Angle of aspect
SmagRad = np.arctan(np.sqrt(Sx**2. + Sy**2.)) # magnitude of slope in radians
return (((np.cos(elevRad) * np.cos(SmagRad)) + (np.sin(elevRad)* np.sin(SmagRad) * np.cos(azRad - AspectRad))))*255
a = Astral()
Latitude = 47.074531
longitude = 12.846210
b = datetime.datetime(2015,10,3,7,30)
c= a.solar_azimuth(b,Latitude,longitude)
d= a.solar_elevation(b,Latitude,longitude)
print d, a.solar_zenith(b, Latitude, longitude)
print c
headinfo,dtm = ascii.read_ascii("C:\Master\settings/vernagtferner14-16/dgm_vernagtferner.txt")
test = calcHillshade(dtm,headinfo[-2],c,d)
print np.min(test), np.max(test)
plt.imshow(test, cmap="Greys")
plt.show()
ascii.write_ascii("test.asc",headinfo,test,format="%i")