code: def node_to_point_geometry(node):
      
      assert "pointGeometry" in node.tag
      for subnode in node.nodes:
          if "Point" in subnode.tag:
              # Position
              lon, lat = map(float, subnode.nodes[0].text.split())
              point = Point(lon, lat)
          elif "upperSeismoDepth" in subnode.tag:
              upper_depth = float_(subnode.text)
          elif "lowerSeismoDepth" in subnode.tag:
              lower_depth = float_(subnode.text)
          else:
              # Redundent
              pass
      assert lower_depth > upper_depth
      return point, upper_depth, lower_depth
docs: Reads the node and returns the point geometry, upper depth and lower depth