working test of percentage shared by walls

Solving a bug from ep workflow

See merge request Guille/hub!62

hub/catalog_factories/construction/nrel_catalog.py

@@ -40,7 +40,7 @@ class NrelCatalog(Catalog):
     _catalog_windows = []
     windows = self._constructions['library']['windows']['window']
     for window in windows:
-      frame_ratio = window['frame_ratio']['#text']
+      frame_ratio = float(window['frame_ratio']['#text'])
       g_value = window['shgc']
       overall_u_value = float(window['conductivity']['#text']) / float(window['thickness']['#text'])
       name = window['@name']
@@ -54,9 +54,9 @@ class NrelCatalog(Catalog):
     for material in materials:
       material_id = material['@id']
       name = material['@name']
-      solar_absorptance = material['solar_absorptance']['#text']
-      thermal_absorptance = material['thermal_absorptance']['#text']
-      visible_absorptance = material['visible_absorptance']['#text']
+      solar_absorptance = float(material['solar_absorptance']['#text'])
+      thermal_absorptance = float(material['thermal_absorptance']['#text'])
+      visible_absorptance = float(material['visible_absorptance']['#text'])
       no_mass = False
       thermal_resistance = None,
       conductivity = None,
@@ -64,11 +64,11 @@ class NrelCatalog(Catalog):
       specific_heat = None
       if 'no_mass' in material and material['no_mass'] == 'true':
         no_mass = True
-        thermal_resistance = material['thermal_resistance']['#text']
+        thermal_resistance = float(material['thermal_resistance']['#text'])
       else:
-        conductivity = material['conductivity']['#text']
-        density = material['density']['#text']
-        specific_heat = material['specific_heat']['#text']
+        conductivity = float(material['conductivity']['#text'])
+        density = float(material['density']['#text'])
+        specific_heat = float(material['specific_heat']['#text'])
       _material = Material(material_id,
                            name,
                            solar_absorptance,
@@ -96,7 +96,7 @@ class NrelCatalog(Catalog):
         material_id = layer['material'][0]
         thickness = 0
         if 'thickness' in layer:
-          thickness = layer['thickness']['#text']
+          thickness = float(layer['thickness']['#text'])
         for material in self._catalog_materials:
           if str(material_id) == str(material.id):
             layers.append(Layer(layer_id, layer_name, material, thickness))
@@ -114,18 +114,20 @@ class NrelCatalog(Catalog):
       climate_zone = archetype['@climate_zone']
       construction_period = \
         ConstructionHelper().reference_standard_to_construction_period[archetype['@reference_standard']]
-      average_storey_height = archetype['average_storey_height']['#text']
-      thermal_capacity = str(float(archetype['thermal_capacity']['#text']) * 1000)
-      extra_loses_due_to_thermal_bridges = archetype['extra_loses_due_to_thermal_bridges']['#text']
-      indirect_heated_ratio = archetype['indirect_heated_ratio']['#text']
-      infiltration_rate_for_ventilation_system_off = archetype['infiltration_rate_for_ventilation_system_off']['#text']
-      infiltration_rate_for_ventilation_system_on = archetype['infiltration_rate_for_ventilation_system_on']['#text']
+      average_storey_height = float(archetype['average_storey_height']['#text'])
+      thermal_capacity = float(archetype['thermal_capacity']['#text']) * 1000
+      extra_loses_due_to_thermal_bridges = float(archetype['extra_loses_due_to_thermal_bridges']['#text'])
+      indirect_heated_ratio = float(archetype['indirect_heated_ratio']['#text'])
+      infiltration_rate_for_ventilation_system_off = \
+        float(archetype['infiltration_rate_for_ventilation_system_off']['#text'])
+      infiltration_rate_for_ventilation_system_on = \
+        float(archetype['infiltration_rate_for_ventilation_system_on']['#text'])
 
       archetype_constructions = []
       for archetype_construction in archetype['constructions']['construction']:
         for construction in self._catalog_constructions:
           if construction.id == archetype_construction['@id']:
-            window_ratio = archetype_construction['window_ratio']['#text']
+            window_ratio = float(archetype_construction['window_ratio']['#text'])
             window_id = archetype_construction['window']
             _construction = None
             _window = None

hub/city_model_structure/attributes/plane.py

@@ -55,7 +55,7 @@ class Plane:
       self._equation = (a, b, c, d)
     return self._equation
 
-  def distance(self, point):
+  def distance_to_point(self, point):
     """
     Distance between the given point and the plane
     :return: float

hub/city_model_structure/building_demand/surface.py

@@ -32,7 +32,6 @@ class Surface:
     self._area = None
     self._lower_corner = None
     self._upper_corner = None
-    self._shared_surfaces = []
     self._global_irradiance = dict()
     self._perimeter_polygon = perimeter_polygon
     self._holes_polygons = holes_polygons
@@ -42,6 +41,7 @@ class Surface:
     self._inverse = None
     self._associated_thermal_boundaries = []
     self._vegetation = None
+    self._percentage_shared = None
 
   @property
   def name(self):
@@ -283,12 +283,6 @@ class Surface:
       self._inverse = Surface(new_solid_polygon, new_perimeter_polygon, new_holes_polygons, cte.VIRTUAL_INTERNAL)
     return self._inverse
 
-  def shared_surfaces(self):
-    """
-    Raises not implemented error
-    """
-    raise NotImplementedError
-
   def divide(self, z):
     """
     Divides a surface at Z plane
@@ -336,3 +330,19 @@ class Surface:
     :param value: None or [ThermalBoundary]
     """
     self._associated_thermal_boundaries = value
+
+  @property
+  def percentage_shared(self):
+    """
+    Get percentage of the wall shared with other walls
+    :return: float
+    """
+    return self._percentage_shared
+
+  @percentage_shared.setter
+  def percentage_shared(self, value):
+    """
+    Set percentage of the wall shared with other walls
+    :param value: float
+    """
+    self._percentage_shared = value

hub/city_model_structure/building_demand/thermal_zone.py

@@ -510,10 +510,14 @@ class ThermalZone:
           _schedule.values = values[:day]
           _schedules.append(_schedule)
 
+      _internal_gain.average_internal_gain = _average_internal_gain
+      _internal_gain.convective_fraction = 0
+      _internal_gain.radiative_fraction = 0
+      _internal_gain.latent_fraction = 0
+      if _average_internal_gain != 0:
         _internal_gain.convective_fraction = _convective_fraction / _average_internal_gain
         _internal_gain.radiative_fraction = _radiative_fraction / _average_internal_gain
         _internal_gain.latent_fraction = _latent_fraction / _average_internal_gain
-      _internal_gain.average_internal_gain = _average_internal_gain
       _internal_gain.type = 'mean_value'
       _internal_gain.schedules = _schedules
       self._internal_gains = [_internal_gain]

hub/city_model_structure/building_demand/usage.py

@@ -91,6 +91,10 @@ class Usage:
                               + self.occupancy.sensible_radiative_internal_gain
                               + self.occupancy.latent_internal_gain)
           _internal_gain.average_internal_gain = _total_heat_gain
+          _internal_gain.latent_fraction = 0
+          _internal_gain.radiative_fraction = 0
+          _internal_gain.convective_fraction = 0
+          if _total_heat_gain != 0:
             _internal_gain.latent_fraction = self.occupancy.latent_internal_gain / _total_heat_gain
             _internal_gain.radiative_fraction = self.occupancy.sensible_radiative_internal_gain / _total_heat_gain
             _internal_gain.convective_fraction = self.occupancy.sensible_convective_internal_gain / _total_heat_gain

hub/city_model_structure/building_demand/usage_zone.py

@@ -91,6 +91,10 @@ class UsageZone:
                               + self.occupancy.sensible_radiative_internal_gain
                               + self.occupancy.latent_internal_gain)
           _internal_gain.average_internal_gain = _total_heat_gain
+          _internal_gain.latent_fraction = 0
+          _internal_gain.radiative_fraction = 0
+          _internal_gain.convective_fraction = 0
+          if _total_heat_gain != 0:
             _internal_gain.latent_fraction = self.occupancy.latent_internal_gain / _total_heat_gain
             _internal_gain.radiative_fraction = self.occupancy.sensible_radiative_internal_gain / _total_heat_gain
             _internal_gain.convective_fraction = self.occupancy.sensible_convective_internal_gain / _total_heat_gain

hub/city_model_structure/city_object.py

@@ -24,11 +24,15 @@ class CityObject:
     self._surfaces = surfaces
     self._type = None
     self._city_object_lower_corner = None
+    self._city_object_upper_corner = None
     self._detailed_polyhedron = None
     self._simplified_polyhedron = None
     self._min_x = ConfigurationHelper().max_coordinate
     self._min_y = ConfigurationHelper().max_coordinate
     self._min_z = ConfigurationHelper().max_coordinate
+    self._max_x = ConfigurationHelper().min_coordinate
+    self._max_y = ConfigurationHelper().min_coordinate
+    self._max_z = ConfigurationHelper().min_coordinate
     self._centroid = None
     self._external_temperature = dict()
     self._global_horizontal = dict()
@@ -212,6 +216,16 @@ class CityObject:
       self._city_object_lower_corner = [self._min_x, self._min_y, self._min_z]
     return self._city_object_lower_corner
 
+  @property
+  def upper_corner(self):
+    """
+    Get city object upper corner coordinates [x, y, z]
+    :return: [x,y,z]
+    """
+    if self._city_object_upper_corner is None:
+      self._city_object_upper_corner = [self._max_x, self._max_y, self._max_z]
+    return self._city_object_upper_corner
+
   @property
   def sensors(self) -> List[Sensor]:
     """

hub/exports/building_energy/idf.py

@@ -130,10 +130,7 @@ class Idf:
       self._idf.newidfobject(self._MATERIAL_NOMASS,
                              Name=layer.material.name,
                              Roughness=self._ROUGHNESS,
-                             Thermal_Resistance=layer.material.thermal_resistance,
-                             Thermal_Absorptance=layer.material.thermal_absorptance,
-                             Solar_Absorptance=layer.material.solar_absorptance,
-                             Visible_Absorptance=layer.material.visible_absorptance
+                             Thermal_Resistance=layer.material.thermal_resistance
                              )
     else:
       self._idf.newidfobject(self._MATERIAL,
@@ -323,9 +320,12 @@ class Idf:
 
   def _add_occupancy(self, thermal_zone, zone_name):
     number_of_people = thermal_zone.occupancy.occupancy_density * thermal_zone.total_floor_area
-    fraction_radiant = thermal_zone.occupancy.sensible_radiative_internal_gain / \
-                       (thermal_zone.occupancy.sensible_radiative_internal_gain +
-                        thermal_zone.occupancy.sensible_convective_internal_gain)
+    fraction_radiant = 0
+    total_sensible = thermal_zone.occupancy.sensible_radiative_internal_gain + \
+                     thermal_zone.occupancy.sensible_convective_internal_gain
+    if total_sensible != 0:
+      fraction_radiant = thermal_zone.occupancy.sensible_radiative_internal_gain / total_sensible
+
     self._idf.newidfobject(self._PEOPLE,
                            Name=f'{zone_name}_occupancy',
                            Zone_or_ZoneList_Name=zone_name,
@@ -377,7 +377,6 @@ class Idf:
     self._rename_building(self._city.name)
     self._lod = self._city.level_of_detail.geometry
     for building in self._city.buildings:
-
       for internal_zone in building.internal_zones:
         for thermal_zone in internal_zone.thermal_zones:
           for thermal_boundary in thermal_zone.thermal_boundaries:

hub/exports/building_energy/insel/insel_monthly_energy_balance.py

@@ -117,6 +117,8 @@ class InselMonthlyEnergyBalance(Insel):
 
     for thermal_boundary in thermal_zone.thermal_boundaries:
       type_code = _CONSTRUCTION_CODE[thermal_boundary.type]
+      window_area = 0
+      if thermal_boundary.window_ratio < 1:
         window_area = thermal_boundary.opaque_area * thermal_boundary.window_ratio / (1 - thermal_boundary.window_ratio)
 
       parameters.append(type_code)

hub/helpers/geometry_helper.py

@@ -49,6 +49,7 @@ class GeometryHelper:
   """
   Geometry helper class
   """
+  # todo: complete dictionary
   srs_transformations = {
     'urn:adv:crs:ETRS89_UTM32*DE_DHHN92_NH': 'epsg:25832'
   }
@@ -63,15 +64,25 @@ class GeometryHelper:
 
   @staticmethod
   def city_mapping(city, building_names=None, plot=False):
+    """
+    Returns a shared_information dictionary like
+
+    {
+     "building_name" : [{line: 0 coordinate_1: [x,y,z], coordinate_2:[x, y, z], points: 0}]
+    }
+    """
+    shared_information = {}
     if building_names is None:
       building_names = [b.name for b in city.buildings]
     x = int((city.upper_corner[0] - city.lower_corner[0]) * 0.5) + 1
     y = int((city.upper_corner[1] - city.lower_corner[1]) * 0.5) + 1
-    city_map = [['' for _ in range(y + 1)] for _ in range(x + 1)]
+    city_map = [[{} for _ in range(y + 1)] for _ in range(x + 1)]
     img = Image.new('RGB', (x + 1, y + 1), "black")  # create a new black image
     city_image = img.load()  # create the pixel map
     for building_name in building_names:
       building = city.city_object(building_name)
+      shared_information[building_name]: []
+      line = 0
       for ground in building.grounds:
         length = len(ground.perimeter_polygon.coordinates) - 1
         for i, coordinate in enumerate(ground.perimeter_polygon.coordinates):
@@ -79,6 +90,7 @@ class GeometryHelper:
           if i == length:
             j = 0
           next_coordinate = ground.perimeter_polygon.coordinates[j]
+          line_dictionary = {"line": line, "coordinate_1": coordinate, "coordinate_2":next_coordinate, "points": 0}
           point = GeometryHelper.coordinate_to_map_point(coordinate, city)
           distance = GeometryHelper.distance_between_points(coordinate, next_coordinate)
           if distance == 0:
@@ -88,7 +100,7 @@ class GeometryHelper:
           for k in range(0, int(distance)):
             x = MapPoint(point.x + (delta_x * k), point.y + (delta_y * k)).x
             y = MapPoint(point.x + (delta_x * k), point.y + (delta_y * k)).y
-            if city_map[x][y] == '':
+            if city_map[x][y] == {}:
               city_map[x][y] = building.name
               city_image[x, y] = (100, 0, 0)
             elif city_map[x][y] != building.name:
@@ -101,6 +113,7 @@ class GeometryHelper:
                 neighbour.neighbours = [building]
               elif building not in neighbour.neighbours:
                 neighbour.neighbours.append(building)
+          line += 1
     if plot:
       img.show()
 

hub/imports/construction/nrcan_physics_parameters.py

@@ -51,7 +51,7 @@ class NrcanPhysicsParameters:
             for thermal_zone in internal_zone.thermal_zones:
               thermal_zone.total_floor_area = thermal_zone.footprint_area
         else:
-          number_of_storeys = int(float(building.eave_height) / float(building.average_storey_height))
+          number_of_storeys = int(building.eave_height / building.average_storey_height)
           thermal_zone = building.internal_zones[0].thermal_zones[0]
           thermal_zone.total_floor_area = thermal_zone.footprint_area * number_of_storeys
       else:
@@ -69,7 +69,7 @@ class NrcanPhysicsParameters:
     nrcan_archetypes = nrcan_catalog.entries('archetypes')
     for building_archetype in nrcan_archetypes:
       construction_period_limits = building_archetype.construction_period.split('_')
-      if int(construction_period_limits[0]) <= int(year_of_construction) < int(construction_period_limits[1]):
+      if int(construction_period_limits[0]) <= year_of_construction < int(construction_period_limits[1]):
         if (str(function) == str(building_archetype.function)) and \
            (climate_zone == str(building_archetype.climate_zone)):
           return building_archetype
@@ -135,12 +135,12 @@ class NrcanPhysicsParameters:
         # The agreement is that the layers are defined from outside to inside
         external_layer = construction_archetype.layers[0]
         external_surface = thermal_boundary.parent_surface
-        external_surface.short_wave_reflectance = 1 - float(external_layer.material.solar_absorptance)
-        external_surface.long_wave_emittance = 1 - float(external_layer.material.solar_absorptance)
+        external_surface.short_wave_reflectance = 1 - external_layer.material.solar_absorptance
+        external_surface.long_wave_emittance = 1 - external_layer.material.solar_absorptance
         internal_layer = construction_archetype.layers[len(construction_archetype.layers) - 1]
         internal_surface = thermal_boundary.internal_surface
-        internal_surface.short_wave_reflectance = 1 - float(internal_layer.material.solar_absorptance)
-        internal_surface.long_wave_emittance = 1 - float(internal_layer.material.solar_absorptance)
+        internal_surface.short_wave_reflectance = 1 - internal_layer.material.solar_absorptance
+        internal_surface.long_wave_emittance = 1 - internal_layer.material.solar_absorptance
 
         for thermal_opening in thermal_boundary.thermal_openings:
           if construction_archetype.window is not None:

hub/imports/construction/nrel_physics_parameters.py

@@ -58,7 +58,7 @@ class NrelPhysicsParameters:
             for thermal_zone in internal_zone.thermal_zones:
               thermal_zone.total_floor_area = thermal_zone.footprint_area
         else:
-          number_of_storeys = int(float(building.eave_height) / float(building.average_storey_height))
+          number_of_storeys = int(building.eave_height / building.average_storey_height)
           thermal_zone = building.internal_zones[0].thermal_zones[0]
           thermal_zone.total_floor_area = thermal_zone.footprint_area * number_of_storeys
       else:
@@ -78,7 +78,7 @@ class NrelPhysicsParameters:
       construction_period_limits = building_archetype.construction_period.split(' - ')
       if construction_period_limits[1] == 'PRESENT':
         construction_period_limits[1] = 3000
-      if int(construction_period_limits[0]) <= int(year_of_construction) < int(construction_period_limits[1]):
+      if int(construction_period_limits[0]) <= year_of_constructionF < int(construction_period_limits[1]):
         if (str(function) == str(building_archetype.function)) and \
            (climate_zone == str(building_archetype.climate_zone)):
           return building_archetype
@@ -130,12 +130,12 @@ class NrelPhysicsParameters:
         # The agreement is that the layers are defined from outside to inside
         external_layer = construction_archetype.layers[0]
         external_surface = thermal_boundary.parent_surface
-        external_surface.short_wave_reflectance = 1 - float(external_layer.material.solar_absorptance)
-        external_surface.long_wave_emittance = 1 - float(external_layer.material.solar_absorptance)
+        external_surface.short_wave_reflectance = 1 - external_layer.material.solar_absorptance
+        external_surface.long_wave_emittance = 1 - external_layer.material.solar_absorptance
         internal_layer = construction_archetype.layers[len(construction_archetype.layers) - 1]
         internal_surface = thermal_boundary.internal_surface
-        internal_surface.short_wave_reflectance = 1 - float(internal_layer.material.solar_absorptance)
-        internal_surface.long_wave_emittance = 1 - float(internal_layer.material.solar_absorptance)
+        internal_surface.short_wave_reflectance = 1 - internal_layer.material.solar_absorptance
+        internal_surface.long_wave_emittance = 1 - internal_layer.material.solar_absorptance
 
         for thermal_opening in thermal_boundary.thermal_openings:
           if construction_archetype.window is not None:

hub/imports/geometry/rhino.py

@@ -124,7 +124,7 @@ class Rhino:
           # todo: this is a hack for dompark project it should not be done this way windows should be correctly modeled
           # if the distance between the wall plane and the window is less than 2m
           # and the window Z coordinate it's between the wall Z, it's a window of that wall
-          if plane.distance(corner) <= 2:
+          if plane.distance_to_point(corner) <= 2:
             # check if the window is in the right high.
             if surface.upper_corner[2] >= corner[2] >= surface.lower_corner[2]:
               if surface.holes_polygons is None:

hub/imports/usage/comnet_usage_parameters.py

@@ -189,6 +189,11 @@ class ComnetUsageParameters:
         _schedule_values[v, day] += value * archetype.appliances.density
         _sum += value * archetype.appliances.density * _number_of_days_per_type[day]
 
+    _latent_fraction = 0
+    _radiative_fraction = 0
+    _convective_fraction = 0
+    _average_internal_gain = 0
+    if _total_heat_gain != 0:
       _latent_fraction = _latent_heat_gain / _total_heat_gain
       _radiative_fraction = _radiative_heat_gain / _total_heat_gain
       _convective_fraction = _convective_heat_gain / _total_heat_gain

hub/unittests/tests_data/neighbours.geojson

@@ -0,0 +1,594 @@
+{
+  "type":"FeatureCollection",
+  "features":[
+    {
+      "type":"Feature",
+      "id":1,
+      "geometry":{
+        "type":"Polygon",
+        "coordinates":[
+          [
+            [
+              -73.580414175680588,
+              45.497641136608358
+            ],
+            [
+              -73.581414175680588,
+              45.497641136608358
+            ],
+            [
+              -73.581414175680588,
+              45.498641136608358
+            ],
+            [
+              -73.580414175680588,
+              45.498641136608358
+            ],
+            [
+              -73.580414175680588,
+              45.497641136608358
+            ]
+          ]
+        ]
+      },
+      "properties":{
+        "OBJECTID_12":1,
+        "gml_id":"1340908",
+        "gml_parent":"fme-gen-5fa2a82b-c38e-4bf0-9e8f-10a47b9f64f7",
+        "citygml_ta":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_fe":"cityObjectMember",
+        "citygml__1":" ",
+        "citygml__2":" ",
+        "gml_descri":" ",
+        "gml_name":" ",
+        "citygml_cr":" ",
+        "citygml_te":" ",
+        "externalRe":" ",
+        "external_1":" ",
+        "external_2":" ",
+        "citygml_ge":" ",
+        "citygml_re":" ",
+        "citygml__3":" ",
+        "citygml_ap":" ",
+        "citygml_cl":" ",
+        "citygml__4":" ",
+        "citygml_fu":" ",
+        "citygml__5":" ",
+        "citygml_us":" ",
+        "citygml__6":" ",
+        "citygml_ye":" ",
+        "citygml__7":" ",
+        "citygml_ro":" ",
+        "citygml__8":" ",
+        "citygml_me":21.824000000000002,
+        "citygml__9":"#m",
+        "citygml_st":" ",
+        "citygml_10":" ",
+        "citygml_11":" ",
+        "citygml_12":" ",
+        "citygml_13":" ",
+        "citygml_14":" ",
+        "citygml_ou":" ",
+        "citygml_in":" ",
+        "citygml_bo":" ",
+        "citygml_le":" ",
+        "citygml_15":" ",
+        "citygml_co":" ",
+        "citygml_ad":" ",
+        "Volume":"2783.169",
+        "parcelle":" ",
+        "OBJECTID":778,
+        "gml_id_1":"ebc7f916-d094-4de0-8c35-fc18eb8622f2",
+        "gml_pare_1":"1340908",
+        "citygml_16":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_17":"boundedBy",
+        "citygml_18":" ",
+        "citygml_19":" ",
+        "gml_desc_1":" ",
+        "gml_name_1":" ",
+        "citygml_20":" ",
+        "citygml_21":" ",
+        "external_3":" ",
+        "external_4":" ",
+        "external_5":" ",
+        "citygml_22":" ",
+        "citygml_23":" ",
+        "citygml_24":" ",
+        "citygml_25":" ",
+        "citygml_26":" ",
+        "citygml_op":" ",
+        "Area":"229.287",
+        "FID_":0,
+        "Join_Count":2,
+        "TARGET_FID":779,
+        "gml_id_12":"1340908",
+        "gml_pare_2":"fme-gen-5fa2a82b-c38e-4bf0-9e8f-10a47b9f64f7",
+        "citygml_27":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_28":"cityObjectMember",
+        "citygml_29":" ",
+        "citygml_30":" ",
+        "gml_desc_2":" ",
+        "gml_name_2":" ",
+        "citygml_31":" ",
+        "citygml_32":" ",
+        "external_6":" ",
+        "external_7":" ",
+        "external_8":" ",
+        "citygml_33":" ",
+        "citygml_34":" ",
+        "citygml_35":" ",
+        "citygml_36":" ",
+        "citygml_37":" ",
+        "citygml_38":" ",
+        "citygml_39":" ",
+        "citygml_40":" ",
+        "citygml_41":" ",
+        "citygml_42":" ",
+        "citygml_43":" ",
+        "citygml_44":" ",
+        "citygml_45":" ",
+        "citygml_46":" ",
+        "citygml_47":21.824000000000002,
+        "citygml_48":"#m",
+        "citygml_49":" ",
+        "citygml_50":" ",
+        "citygml_51":" ",
+        "citygml_52":" ",
+        "citygml_53":" ",
+        "citygml_54":" ",
+        "citygml_55":" ",
+        "citygml_56":" ",
+        "citygml_57":" ",
+        "citygml_58":" ",
+        "citygml_59":" ",
+        "citygml_60":" ",
+        "citygml_61":" ",
+        "Volume_1":"2783.169",
+        "Field":0,
+        "Field1":0,
+        "OBJECTID_1":778,
+        "gml_id_12_":"ebc7f916-d094-4de0-8c35-fc18eb8622f2",
+        "gml_pare_3":"1340908",
+        "citygml_62":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_63":"boundedBy",
+        "citygml_64":" ",
+        "citygml_65":" ",
+        "gml_desc_3":" ",
+        "gml_name_3":" ",
+        "citygml_66":" ",
+        "citygml_67":" ",
+        "external_9":" ",
+        "externa_10":" ",
+        "externa_11":" ",
+        "citygml_68":" ",
+        "citygml_69":" ",
+        "citygml_70":" ",
+        "citygml_71":" ",
+        "citygml_72":" ",
+        "citygml_73":" ",
+        "Area_1":"229.287",
+        "cityGML_hi":0,
+        "Z_Min":49.0745,
+        "Z_Max":69.165000000000006,
+        "Shape_Leng":59.532834838799999,
+        "ID_UEV":"01002777",
+        "CIVIQUE_DE":" 1460",
+        "CIVIQUE_FI":" 1460",
+        "NOM_RUE":"rue Sherbrooke Ouest  (MTL+MTO+WMT)",
+        "MUNICIPALI":"50",
+        "ETAGE_HORS":3,
+        "NOMBRE_LOG":1,
+        "ANNEE_CONS":1885,
+        "CODE_UTILI":"5010",
+        "LIBELLE_UT":"Immeuble commercial",
+        "CATEGORIE_":"Régulier",
+        "MATRICULE8":"9839-57-1941-6-000-0000",
+        "SUPERFICIE":193,
+        "SUPERFIC_1":609,
+        "NO_ARROND_":"REM19",
+        "Shape_Le_1":0.00076452447366199996,
+        "Shape_Ar_1":2.2162879886799998e-08,
+        "Z_Min_1":null,
+        "Z_Max_1":null,
+        "Shape_Length":59.532834838827348,
+        "Shape_Area":161.83671944596372
+      }
+    },
+    {
+      "type":"Feature",
+      "id":2,
+      "geometry":{
+        "type":"Polygon",
+        "coordinates":[
+          [
+            [
+              -73.581414175680588,
+              45.497641136608358
+            ],
+            [
+              -73.582214175680588,
+              45.497641136608358
+            ],
+            [
+              -73.582214175680588,
+              45.498441136608358
+            ],
+            [
+              -73.581414175680588,
+              45.498441136608358
+            ],
+            [
+              -73.581414175680588,
+              45.497641136608358
+            ]
+          ]
+        ]
+  
+      },
+      "properties":{
+        "OBJECTID_12":2,
+        "gml_id":"1340974",
+        "gml_parent":"fme-gen-5fa2a82b-c38e-4bf0-9e8f-10a47b9f64f7",
+        "citygml_ta":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_fe":"cityObjectMember",
+        "citygml__1":" ",
+        "citygml__2":" ",
+        "gml_descri":" ",
+        "gml_name":" ",
+        "citygml_cr":" ",
+        "citygml_te":" ",
+        "externalRe":" ",
+        "external_1":" ",
+        "external_2":" ",
+        "citygml_ge":" ",
+        "citygml_re":" ",
+        "citygml__3":" ",
+        "citygml_ap":" ",
+        "citygml_cl":" ",
+        "citygml__4":" ",
+        "citygml_fu":" ",
+        "citygml__5":" ",
+        "citygml_us":" ",
+        "citygml__6":" ",
+        "citygml_ye":" ",
+        "citygml__7":" ",
+        "citygml_ro":" ",
+        "citygml__8":" ",
+        "citygml_me":21.643999999999998,
+        "citygml__9":"#m",
+        "citygml_st":" ",
+        "citygml_10":" ",
+        "citygml_11":" ",
+        "citygml_12":" ",
+        "citygml_13":" ",
+        "citygml_14":" ",
+        "citygml_ou":" ",
+        "citygml_in":" ",
+        "citygml_bo":" ",
+        "citygml_le":" ",
+        "citygml_15":" ",
+        "citygml_co":" ",
+        "citygml_ad":" ",
+        "Volume":"8410.522",
+        "parcelle":" ",
+        "OBJECTID":779,
+        "gml_id_1":"96e73b07-262d-43a8-84ce-608133b39f16",
+        "gml_pare_1":"1340974",
+        "citygml_16":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_17":"boundedBy",
+        "citygml_18":" ",
+        "citygml_19":" ",
+        "gml_desc_1":" ",
+        "gml_name_1":" ",
+        "citygml_20":" ",
+        "citygml_21":" ",
+        "external_3":" ",
+        "external_4":" ",
+        "external_5":" ",
+        "citygml_22":" ",
+        "citygml_23":" ",
+        "citygml_24":" ",
+        "citygml_25":" ",
+        "citygml_26":" ",
+        "citygml_op":" ",
+        "Area":"553.859",
+        "FID_":0,
+        "Join_Count":3,
+        "TARGET_FID":780,
+        "gml_id_12":"1340974",
+        "gml_pare_2":"fme-gen-5fa2a82b-c38e-4bf0-9e8f-10a47b9f64f7",
+        "citygml_27":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_28":"cityObjectMember",
+        "citygml_29":" ",
+        "citygml_30":" ",
+        "gml_desc_2":" ",
+        "gml_name_2":" ",
+        "citygml_31":" ",
+        "citygml_32":" ",
+        "external_6":" ",
+        "external_7":" ",
+        "external_8":" ",
+        "citygml_33":" ",
+        "citygml_34":" ",
+        "citygml_35":" ",
+        "citygml_36":" ",
+        "citygml_37":" ",
+        "citygml_38":" ",
+        "citygml_39":" ",
+        "citygml_40":" ",
+        "citygml_41":" ",
+        "citygml_42":" ",
+        "citygml_43":" ",
+        "citygml_44":" ",
+        "citygml_45":" ",
+        "citygml_46":" ",
+        "citygml_47":21.643999999999998,
+        "citygml_48":"#m",
+        "citygml_49":" ",
+        "citygml_50":" ",
+        "citygml_51":" ",
+        "citygml_52":" ",
+        "citygml_53":" ",
+        "citygml_54":" ",
+        "citygml_55":" ",
+        "citygml_56":" ",
+        "citygml_57":" ",
+        "citygml_58":" ",
+        "citygml_59":" ",
+        "citygml_60":" ",
+        "citygml_61":" ",
+        "Volume_1":"8410.522",
+        "Field":0,
+        "Field1":0,
+        "OBJECTID_1":779,
+        "gml_id_12_":"96e73b07-262d-43a8-84ce-608133b39f16",
+        "gml_pare_3":"1340974",
+        "citygml_62":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_63":"boundedBy",
+        "citygml_64":" ",
+        "citygml_65":" ",
+        "gml_desc_3":" ",
+        "gml_name_3":" ",
+        "citygml_66":" ",
+        "citygml_67":" ",
+        "external_9":" ",
+        "externa_10":" ",
+        "externa_11":" ",
+        "citygml_68":" ",
+        "citygml_69":" ",
+        "citygml_70":" ",
+        "citygml_71":" ",
+        "citygml_72":" ",
+        "citygml_73":" ",
+        "Area_1":"553.859",
+        "cityGML_hi":0,
+        "Z_Min":47.817900000000002,
+        "Z_Max":69.462000000000003,
+        "Shape_Leng":124.143194192,
+        "ID_UEV":"01002773",
+        "CIVIQUE_DE":" 1438",
+        "CIVIQUE_FI":" 1438",
+        "NOM_RUE":"rue Sherbrooke Ouest  (MTL+MTO+WMT)",
+        "MUNICIPALI":"50",
+        "ETAGE_HORS":3,
+        "NOMBRE_LOG":2,
+        "ANNEE_CONS":1885,
+        "CODE_UTILI":"1000",
+        "LIBELLE_UT":"Logement",
+        "CATEGORIE_":"Régulier",
+        "MATRICULE8":"9839-57-4570-0-000-0000",
+        "SUPERFICIE":249,
+        "SUPERFIC_1":506,
+        "NO_ARROND_":"REM19",
+        "Shape_Le_1":0.00099703639048799998,
+        "Shape_Ar_1":2.8543276304299999e-08,
+        "Z_Min_1":null,
+        "Z_Max_1":null,
+        "Shape_Length":124.143194192441,
+        "Shape_Area":464.30094602931189
+      }
+    },
+    {
+      "type":"Feature",
+      "id":3,
+      "geometry":{
+        "type":"Polygon",
+        "coordinates":[
+          [
+            [
+              -73.581914175680588,
+              45.498441136608358
+            ],
+            [
+              -73.581914175680588,
+              45.499641136608358
+            ],
+            [
+              -73.580914175680588,
+              45.499641136608358
+            ],
+            [
+              -73.580914175680588,
+              45.498641136608358
+            ],
+            [
+              -73.581414175680588,
+              45.498641136608358
+            ],
+            [
+              -73.581414175680588,
+              45.498441136608358
+            ],
+            [
+              -73.581914175680588,
+              45.498441136608358
+            ]
+
+          ]
+        ]
+      },
+      "properties":{
+        "OBJECTID_12":3,
+        "gml_id":"1340910",
+        "gml_parent":"fme-gen-5fa2a82b-c38e-4bf0-9e8f-10a47b9f64f7",
+        "citygml_ta":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_fe":"cityObjectMember",
+        "citygml__1":" ",
+        "citygml__2":" ",
+        "gml_descri":" ",
+        "gml_name":" ",
+        "citygml_cr":" ",
+        "citygml_te":" ",
+        "externalRe":" ",
+        "external_1":" ",
+        "external_2":" ",
+        "citygml_ge":" ",
+        "citygml_re":" ",
+        "citygml__3":" ",
+        "citygml_ap":" ",
+        "citygml_cl":" ",
+        "citygml__4":" ",
+        "citygml_fu":" ",
+        "citygml__5":" ",
+        "citygml_us":" ",
+        "citygml__6":" ",
+        "citygml_ye":" ",
+        "citygml__7":" ",
+        "citygml_ro":" ",
+        "citygml__8":" ",
+        "citygml_me":21.916,
+        "citygml__9":"#m",
+        "citygml_st":" ",
+        "citygml_10":" ",
+        "citygml_11":" ",
+        "citygml_12":" ",
+        "citygml_13":" ",
+        "citygml_14":" ",
+        "citygml_ou":" ",
+        "citygml_in":" ",
+        "citygml_bo":" ",
+        "citygml_le":" ",
+        "citygml_15":" ",
+        "citygml_co":" ",
+        "citygml_ad":" ",
+        "Volume":"2257.436",
+        "parcelle":" ",
+        "OBJECTID":780,
+        "gml_id_1":"8222a1c7-e161-421a-8478-22d2a116e0b4",
+        "gml_pare_1":"1340910",
+        "citygml_16":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_17":"boundedBy",
+        "citygml_18":" ",
+        "citygml_19":" ",
+        "gml_desc_1":" ",
+        "gml_name_1":" ",
+        "citygml_20":" ",
+        "citygml_21":" ",
+        "external_3":" ",
+        "external_4":" ",
+        "external_5":" ",
+        "citygml_22":" ",
+        "citygml_23":" ",
+        "citygml_24":" ",
+        "citygml_25":" ",
+        "citygml_26":" ",
+        "citygml_op":" ",
+        "Area":"144.697",
+        "FID_":0,
+        "Join_Count":2,
+        "TARGET_FID":781,
+        "gml_id_12":"1340910",
+        "gml_pare_2":"fme-gen-5fa2a82b-c38e-4bf0-9e8f-10a47b9f64f7",
+        "citygml_27":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_28":"cityObjectMember",
+        "citygml_29":" ",
+        "citygml_30":" ",
+        "gml_desc_2":" ",
+        "gml_name_2":" ",
+        "citygml_31":" ",
+        "citygml_32":" ",
+        "external_6":" ",
+        "external_7":" ",
+        "external_8":" ",
+        "citygml_33":" ",
+        "citygml_34":" ",
+        "citygml_35":" ",
+        "citygml_36":" ",
+        "citygml_37":" ",
+        "citygml_38":" ",
+        "citygml_39":" ",
+        "citygml_40":" ",
+        "citygml_41":" ",
+        "citygml_42":" ",
+        "citygml_43":" ",
+        "citygml_44":" ",
+        "citygml_45":" ",
+        "citygml_46":" ",
+        "citygml_47":21.916,
+        "citygml_48":"#m",
+        "citygml_49":" ",
+        "citygml_50":" ",
+        "citygml_51":" ",
+        "citygml_52":" ",
+        "citygml_53":" ",
+        "citygml_54":" ",
+        "citygml_55":" ",
+        "citygml_56":" ",
+        "citygml_57":" ",
+        "citygml_58":" ",
+        "citygml_59":" ",
+        "citygml_60":" ",
+        "citygml_61":" ",
+        "Volume_1":"2257.436",
+        "Field":0,
+        "Field1":0,
+        "OBJECTID_1":780,
+        "gml_id_12_":"8222a1c7-e161-421a-8478-22d2a116e0b4",
+        "gml_pare_3":"1340910",
+        "citygml_62":"http://www.opengis.net/citygml/building/2.0",
+        "citygml_63":"boundedBy",
+        "citygml_64":" ",
+        "citygml_65":" ",
+        "gml_desc_3":" ",
+        "gml_name_3":" ",
+        "citygml_66":" ",
+        "citygml_67":" ",
+        "external_9":" ",
+        "externa_10":" ",
+        "externa_11":" ",
+        "citygml_68":" ",
+        "citygml_69":" ",
+        "citygml_70":" ",
+        "citygml_71":" ",
+        "citygml_72":" ",
+        "citygml_73":" ",
+        "Area_1":"144.697",
+        "cityGML_hi":0,
+        "Z_Min":48.983400000000003,
+        "Z_Max":67.617000000000004,
+        "Shape_Leng":52.283656634099998,
+        "ID_UEV":"01002775",
+        "CIVIQUE_DE":" 1448",
+        "CIVIQUE_FI":" 1448",
+        "NOM_RUE":"rue Sherbrooke Ouest  (MTL+MTO+WMT)",
+        "MUNICIPALI":"50",
+        "ETAGE_HORS":3,
+        "NOMBRE_LOG":1,
+        "ANNEE_CONS":1885,
+        "CODE_UTILI":"5010",
+        "LIBELLE_UT":"Immeuble commercial",
+        "CATEGORIE_":"Régulier",
+        "MATRICULE8":"9839-57-3057-9-000-0000",
+        "SUPERFICIE":167,
+        "SUPERFIC_1":354,
+        "NO_ARROND_":"REM19",
+        "Shape_Le_1":0.00074417728924999998,
+        "Shape_Ar_1":1.92186900974e-08,
+        "Z_Min_1":null,
+        "Z_Max_1":null,
+        "Shape_Length":52.283656634094768,
+        "Shape_Area":123.24449716965384
+      }
+    }
+  ]
+}