docs: def create_station(name, latlonalt, parent_frame=WGS84, orientation='N', mask=None):
      """Create a ground station instance
  
      Args:
          name (str): Name of the station
  
          latlonalt (tuple of float): coordinates of the station, as follow:
  
              * Latitude in degrees
              * Longitude in degrees
              * Altitude to sea level in meters
  
          parent_frame (Frame): Planetocentric rotating frame of reference of
              coordinates.
          orientation (str or float): Heading of the station
              Acceptable values are 'N', 'S', 'E', 'W' or any angle in radians
          mask: (2D array of float): First dimension is azimuth counterclockwise strictly increasing.
              Second dimension is elevation. Both in radians
  
      Return:
          TopocentricFrame
      """
  
      if isinstance(orientation, str):
          orient = {'N': np.pi, 'S': 0., 'E': np.pi / 2., 'W': 3 * np.pi / 2.}
          heading = orient[orientation]
      else:
          heading = orientation
  
      latlonalt = list(latlonalt)
      latlonalt[:2] = np.radians(latlonalt[:2])
      coordinates = TopocentricFrame._geodetic_to_cartesian(*latlonalt)
  
      mtd = '_to_%s' % parent_frame.__name__
      dct = {
          mtd: TopocentricFrame._to_parent_frame,
          'latlonalt': latlonalt,
          'coordinates': coordinates,
          'parent_frame': parent_frame,
          'heading': heading,
          'orientation': orientation,
          'mask': np.array(mask) if mask else None,
      }
      cls = _MetaFrame(name, (TopocentricFrame,), dct)
      cls + parent_frame
  
      return cls