def load_horizon_array(array, loca = None):
if loca: place = th.location(loca)
else: place = th.location('default')
#place = s.location('akb')
gamma = np.deg2rad(np.array([v[0] for v in array]))
alpha = np.deg2rad(np.array([v[1] for v in array]))
return horizon(gamma,alpha,place)
def load_horizon_array(array, loca=None):
if loca: place = th.location(loca)
else: place = th.location('default')
#place = s.location('akb')
gamma = np.deg2rad(np.array([v[0] for v in array]))
alpha = np.deg2rad(np.array([v[1] for v in array]))
return horizon(gamma, alpha, place)
def load_horizon(fname,loca):
'''Open SolarShading app's csv output, and returns horizon data arrays'''
fin = open(fname,'r')
h = np.array(fin.readline().rstrip(',\n').split(','),dtype=float) # header line
gamma = np.deg2rad(np.array(fin.readline().rstrip(',\n').split(','),dtype=float)) # azimuth angle [deg]
alpha = np.deg2rad(np.array(fin.readline().rstrip(',\n').split(','),dtype=float)) # altitude angle [deg]
pname = re.sub('\..+','',fname)
if loca: place = th.location(loca)
elif len(h)>3: # for android .csv files.
place = s.location('default',h[0],-h[1],h[2])
#gamma = gamma - pi
else: place = th.location('default')
fin.close()
return horizon(gamma,alpha,place) #place, hor
def load_horizon(fname, loca):
"""Open SolarShading app's csv output, and returns horizon data arrays"""
fin = open(fname, "r")
h = np.array(fin.readline().rstrip(",\n").split(","), dtype=float) # header line
gamma = np.deg2rad(np.array(fin.readline().rstrip(",\n").split(","), dtype=float)) # azimuth angle [deg]
alpha = np.deg2rad(np.array(fin.readline().rstrip(",\n").split(","), dtype=float)) # altitude angle [deg]
pname = re.sub("\..+", "", fname)
if loca:
place = th.location(loca)
elif len(h) > 3: # for android .csv files.
place = s.location("default", h[0], -h[1], h[2])
# gamma = gamma - pi
else:
place = th.location("default")
fin.close()
return horizon(gamma, alpha, place) # place, hor
def load_horizon(fname, loca):
'''Open SolarShading app's csv output, and returns horizon data arrays'''
fin = open(fname, 'r')
h = np.array(fin.readline().rstrip(',\n').split(','),
dtype=float) # header line
gamma = np.deg2rad(
np.array(fin.readline().rstrip(',\n').split(','),
dtype=float)) # azimuth angle [deg]
alpha = np.deg2rad(
np.array(fin.readline().rstrip(',\n').split(','),
dtype=float)) # altitude angle [deg]
pname = re.sub('\..+', '', fname)
if loca: place = th.location(loca)
elif len(h) > 3: # for android .csv files.
place = s.location('default', h[0], -h[1], h[2])
#gamma = gamma - pi
else:
place = th.location('default')
fin.close()
return horizon(gamma, alpha, place) #place, hor
C = th.collectorArray(th.chuanghui_H50(), 1, slope=20, azimuth=0, DFR=True, rho_g=0.2) # Change
C.beta = 0 # Calculations for horizontal radiation
gamma = gamma - pi # S = 0 deg for solar calculations & plotting
r = 90 - np.rad2deg(alpha) # r in degrees
fig = plt.figure(1)
ax = fig.add_axes([0.1, 0.1, 0.8, 0.8], polar=True)
plot_sunpath(ax, L, C)
plot_horizon(ax, gamma, r)
fig.savefig("test.png")
return fig
if __name__ == "__main__":
args = get_args()
L = th.location() # Default location for starters
C = th.collectorArray(th.chuanghui_H50(), 1, slope=20, azimuth=0, DFR=True, rho_g=0.2) # Change
# C = s.collector('akb1') # Collector object from solarfun.py
C.beta = 0 # Calculations for horizontal radiation
plt.ion()
fig = plt.figure(1)
if args.image:
axi = fig.add_axes([0.1, 0.1, 0.8, 0.8], frame_on=False)
axi = make_invisible(axi)
pic = plt.imread(args.image)
axi.imshow(pic)
ax = fig.add_axes([0.1, 0.1, 0.8, 0.8], polar=True)
ax.patch.set_alpha(0)
if args.horizon:
H = load_horizon(args.horizon, args.location) # Horizon data object from .csv file
rho_g=0.2) # Change
C.beta = 0 # Calculations for horizontal radiation
gamma = gamma - pi # S = 0 deg for solar calculations & plotting
r = 90 - np.rad2deg(alpha) # r in degrees
fig = plt.figure(1)
ax = fig.add_axes([0.1, 0.1, 0.8, 0.8], polar=True)
plot_sunpath(ax, L, C)
plot_horizon(ax, gamma, r)
fig.savefig('test.png')
return fig
if __name__ == "__main__":
args = get_args()
L = th.location() # Default location for starters
C = th.collectorArray(th.chuanghui_H50(),
1,
slope=20,
azimuth=0,
DFR=True,
rho_g=0.2) # Change
#C = s.collector('akb1') # Collector object from solarfun.py
C.beta = 0 # Calculations for horizontal radiation
plt.ion()
fig = plt.figure(1)
if args.image:
axi = fig.add_axes([0.1, 0.1, 0.8, 0.8], frame_on=False)
axi = make_invisible(axi)
pic = plt.imread(args.image)
