def get_a_tree(self):
""" Retrieve a single tree and assign its species, standid, and plot; create a list of tuples for each year of its life. Gather the equations it needs to have its biomass computed.
"""
sql = self.tree_query.format(tid=self.tid)
self.cur.execute(sql)
for index,row in enumerate(self.cur):
if index == 0:
self.species = str(row[1]).strip().lower()
self.standid = str(row[2]).strip().lower()
self.plotid = str(row[3]).strip().lower()
else:
pass
# append to state to ( year, dbh, status, dbh_code )
try:
self.state.append( [int(str(row[6])), round(float(str(row[4])),3), str(row[5]), str(row[7])] )
except Exception:
self.state.append( [int(str(row[6])), None, str(row[5]), str(row[7])] )
# get the equation for that tree
sql_2 = self.eqn_query.format(species=self.species)
# note : f = lambda x : biomass_basis.which_fx('as_d2htcasc')(0.5, x, 5, 6, 7, 2, 3)
self.cur.execute(sql_2)
for row in self.cur:
# identify the form of the equation -- as long as its not comp bio we just extract it -- if it's comp bio we need to get each reference. in some cases the equation might cross the 'big' barrier so we need to index it by this
form = str(row[2]).strip().lower()
try:
woodden = round(float(str(row[11])),3)
except:
woodden = None
try:
h1 = round(float(str(row[3])),6)
except:
h1 = None
try:
h2 = round(float(str(row[4])),6)
except:
h2 = None
try:
h3 = round(float(str(row[5])),6)
except:
h3 = None
try:
b1 = round(float(str(row[6])),6)
except:
b1 = None
try:
b2 = round(float(str(row[7])),6)
except:
b2 = None
try:
b3 = round(float(str(row[8])),6)
except:
b3 = None
try:
j1 = round(float(str(row[9])),6)
except:
j1 = None
try:
j2 = round(float(str(row[10])),6)
except:
j2 = None
if form != 'as_compbio':
this_eqn = lambda x : biomass_basis.which_fx(form)(woodden, x, b1, b2, b3, j1, j2, h1, h2, h3)
self.eqns.update({str(row[1]):this_eqn})
elif form == 'as_compbio':
this_eqn = lambda x: biomass_basis.which_fx('as_biopak')(woodden, x, b1, b2, b3, j1, j2, h1, h2, h3)
self.eqns.update({str(row[12]):this_eqn})
def select_eqns(self):
""" Get only the equations you need based on the species on that plot
Queries the database for individual species that will be on this stand and makes an equation table
"""
list_species = []
self.cur.execute(self.species_list.format(standid = self.standid))
for row in self.cur:
list_species.append(str(row[0]).strip().lower())
for each_species in list_species:
sql2 = self.eqn_query.format(species=each_species)
self.cur.execute(sql2)
for row in self.cur:
form = str(row[2]).strip().lower()
try:
woodden = round(float(str(row[11])),3)
except:
woodden = None
self.woodden = woodden
try:
h1 = round(float(str(row[3])),6)
except:
h1 = None
try:
h2 = round(float(str(row[4])),6)
except:
h2 = None
try:
h3 = round(float(str(row[5])),6)
except:
h3 = None
try:
b1 = round(float(str(row[6])),6)
except:
b1 = None
try:
b2 = round(float(str(row[7])),6)
except:
b2 = None
try:
b3 = round(float(str(row[8])),6)
except:
b3 = None
try:
j1 = round(float(str(row[9])),6)
except:
j1 = None
try:
j2 = round(float(str(row[10])),6)
except:
j2 = None
if form != 'as_compbio':
this_eqn = lambda x : biomass_basis.which_fx(form)(woodden, x, b1, b2, b3, j1, j2, h1, h2, h3)
if each_species not in self.eqns:
self.eqns.update({each_species:{str(row[1]).strip():this_eqn}})
elif each_species in self.eqns:
# when there are 2 + sizes
self.eqns[each_species].update({str(row[1]).strip():this_eqn})
elif form == 'as_compbio':
this_eqn = lambda x: biomass_basis.which_fx('as_biopak')(woodden, x, b1, b2, b3, j1, j2, h1, h2, h3)
self.eqns.update({each_species: {str(row[12]):this_eqn}})
