Class: URBANopt::RNM::OhUgRate

Inherits:
Object
  • Object
show all
Defined in:
lib/urbanopt/rnm/oh_ug_rate.rb

Instance Attribute Summary collapse

Instance Method Summary collapse

Constructor Details

#initialize ⇒ OhUgRate

Returns a new instance of OhUgRate.



13
14
15
16
17
 
# File 'lib/urbanopt/rnm/oh_ug_rate.rb', line 13

def initialize
  @av_height = av_height
  @type = type
  @n_buildings_street = n_buildings_street
end

Instance Attribute Details

#av_height ⇒ Object

Returns the value of attribute av_height.



11
12
13
 
# File 'lib/urbanopt/rnm/oh_ug_rate.rb', line 11

def av_height
  @av_height
end

#n_buildings_street ⇒ Object

Returns the value of attribute n_buildings_street.



11
12
13
 
# File 'lib/urbanopt/rnm/oh_ug_rate.rb', line 11

def n_buildings_street
  @n_buildings_street
end

#type ⇒ Object

Returns the value of attribute type.



11
12
13
 
# File 'lib/urbanopt/rnm/oh_ug_rate.rb', line 11

def type
  @type
end

Instance Method Details

#classify_street_type(street_type, ug_ratio) ⇒ Object

defining a method that decides to place each street either OH or UG, according to its average height and the threshold height, if the street has a height below the threshold height it has to modelled has OH, otherwise UG



84
85
86
87
88
89
90
91
 
# File 'lib/urbanopt/rnm/oh_ug_rate.rb', line 84

def classify_street_type(street_type, ug_ratio)
  h = threshold_height(street_type, ug_ratio)
  if @av_height < h
    @type = 'A' # 'OH'
  else
    @type = 'S' # 'UG'
  end
end

#height_building(coordinates_building, coordinates_street, floors) ⇒ Object

for each street calculate the average height given the buildings in the streets and it calculates the number of buildings in each street



21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
 
# File 'lib/urbanopt/rnm/oh_ug_rate.rb', line 21

def height_building(coordinates_building, coordinates_street, floors)
  height = []
  id_building = []
  height_sum = 0
  n_buildings_street = 0
  for i in 0..coordinates_street.length - 1
    dist_min = 5000
    for j in 0..floors.length - 1
      for k in 0..coordinates_building[j].length - 1
        y = coordinates_building[j][k][1] - coordinates_street[i][1]
        x = coordinates_building[j][k][0] - coordinates_street[i][0]
        distance = (x + y) / 2 # finding the distance of each building node with each street node
        if distance < dist_min
          distance = (x**2 + y**2)**0.5
          if distance < dist_min
            dist_min = distance
            height[i] = floors[j]
            id_building[i] = coordinates_building[j][k][3].split('_')[0] # for future implementations set as attributes :id, :coordinates
          end
        end
      end
    end
    ii = 0
    while ii != i && id_building[i] != id_building[ii] # id_building[ii] != id_bui
      ii += 1
      if ii == i
        n_buildings_street += 1
      end
    end
    height_sum += height[i]
  end
  @n_buildings_street = n_buildings_street
  @av_height = (height_sum / (i + 1)).to_f.round(2)
end

#threshold_height(street_type, ug_ratio) ⇒ Object

defining a method which defines the “threshold height”, in an iterative way and adding 0.1m of height until the threhold limit is reached when the % of streets in the district above the threshold is equal to the UG rate defined by the user



58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
 
# File 'lib/urbanopt/rnm/oh_ug_rate.rb', line 58

def threshold_height(street_type, ug_ratio)
  h_threshold = 0
  tot_build_in_streets = 0
  n_street_oh = []
  for i in 0..street_type.length - 1
    tot_build_in_streets += street_type[i].n_buildings_street
  end
  # puts tot_build_in_streets
  street_set_oh = (tot_build_in_streets * (1 - ug_ratio)).round
  ii = 0
  while ii == 0 || n_street_oh[ii - 1] < street_set_oh
    n_street_oh[ii] = 0
    for i in 0..street_type.length - 1
      if street_type[i].av_height < h_threshold
        n_street_oh[ii] += street_type[i].n_buildings_street
      end
    end
    ii += 1
    h_threshold += 0.10
  end
  h_threshold -= 0.10
  return h_threshold
end