# coding=utf-8
# Copyright 2023 The Chirp Authors.
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
#     http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.

"""Soundscape datasets."""

import dataclasses
import tempfile
from typing import Any, Callable, Sequence
import warnings

from absl import logging
from chirp import audio_utils
from chirp.data import tfds_features
from chirp.data.bird_taxonomy import bird_taxonomy
from chirp.data.soundscapes import dataset_fns
from chirp.data.soundscapes import soundscapes_lib
from etils import epath
import numpy as np
import pandas as pd
import tensorflow as tf
import tensorflow_datasets as tfds

_DESCRIPTION = """
Soundscape datasets.
"""

_CITATION = """
@inproceedings{kahl2019overview,
  title={Overview of BirdCLEF 2019: large-scale bird recognition in soundscapes},
  author={Kahl, Stefan and St{\"o}ter, Fabian-Robert and Go{\"e}au, Herv{\'e} and Glotin, Herv{\'e} and Planque, Robert and Vellinga, Willem-Pier and Joly, Alexis},
  booktitle={Working Notes of CLEF 2019-Conference and Labs of the Evaluation Forum},
  number={2380},
  pages={1--9},
  year={2019},
  organization={CEUR}
}
"""


@dataclasses.dataclass
class SoundscapesConfig(bird_taxonomy.BirdTaxonomyConfig):
  """Dataset configuration for Soundscape datasets.

  Attributes:
    audio_glob: Pattern to match to find audio files.
    class_list_name: The name of the ClassList to use for labels. This is
      typically a list of either all regionally feasible species in the area
      (for fully-annotated datasets) or the list of all species annotated (if
      only a subset has been labeled).
    metadata_load_fn: Because the metadata don't always appear under the same
      format, we specify for each config the way to load metadata. This function
      outputs a dataframe, where each row contains the metadata for some audio
      file. The column names of the dataframe should match the keys in
      metadata_fields.
    metadata_fields: Maps the fields of the metadata DataFrame to tfds.features
      datatypes.
    annotation_filename: Filename for the annotations file. Defaults to None, in
      which case the filename is derived from the config name.
    annotation_load_fn: Because the annotations don't always appear in the same
      format, we specify a function to load the annotations.
    keep_unknown_annotation: An "unknown" annotations appears in some datasets.
      This boolean decides whether it should keep this annotation (and
      therefore) add a species named "unknown" in the label space, or just scrub
      all "unknown" annotations.
    full_length_unknown_guard: If True, add an "unknown" annotation from the
      beginning of the recording to the beginning of the first annotation and
      another "unknown" annotation from the end of the last annotation to the
      end of the recording.
    supervised: Whether this is a supervised dataset. If so, any segment which
      overlaps an 'unknown' label will be dropped (to avoid downward bias on
      eval stats).
    audio_dir: Base directory for soundscapes data.
  """

  audio_glob: str = ''
  class_list_name: str = ''
  metadata_load_fn: soundscapes_lib.MetadataLoaderType | None = None
  metadata_fields: dict[str, soundscapes_lib.MetadataFeature] | None = None
  annotation_filename: str | None = None
  annotation_load_fn: Callable[[epath.Path], pd.DataFrame] | None = None
  keep_unknown_annotation: bool = False
  full_length_unknown_guard: bool = False
  supervised: bool = True
  audio_dir = epath.Path('gs://chirp-public-bucket/soundscapes')


class Soundscapes(bird_taxonomy.BirdTaxonomy):
  """DatasetBuilder for soundscapes data."""

  VERSION = tfds.core.Version('1.3.0')
  RELEASE_NOTES = {
      '1.0.0': (
          'Initial release. The label set corresponds to the full '
          'set of ~11 000 Xeno-Canto species.'
      ),
      '1.0.1': (
          'The label set is now restricted to the species present in each'
          'dataset.'
      ),
      '1.0.2': (
          'Streamlines data handling, and adds handling for a new '
          'Sapsucker Woods dataset.'
      ),
      '1.0.3': (
          'Adds handling for the new Cornell Sierra Nevadas dataset and '
          'the KitzesLab Powdermill dataset.'
      ),
      '1.0.4': 'Adds a unique recording ID and a segment ID to all samples.',
      '1.0.5': (
          'Adds Peru dataset and moves to new version of SSW annotations. '
          'Supervised segments with the "unknown" label are now dropped.'
      ),
      '1.0.6': 'Updates the dataset following Int16AsFloatTensor refactoring.',
      '1.0.7': 'Fix some dropped annotations in the Hawaii dataset.',
      '1.1.0': 'Adds full-length variants.',
      '1.2.0': 'Updated ebird2021 taxonomy.',
      '1.3.0': (
          'Switch to Zenodo versions of soundscape datasets from Cornell. '
          'Add the Colombia+Costa Rica dataset, remove the Colombia-only '
          'BirdCLEF dataset.'
      ),
  }
  BUILDER_CONFIGS = [
      # pylint: disable=unexpected-keyword-arg
      SoundscapesConfig(
          name='caples',  # TODO(mboudiaf) Try to interface caples metadata.
          class_list_name='caples',
          audio_glob='caples/audio/*',
          interval_length_s=5.0,
          localization_fn=audio_utils.slice_peaked_audio,
          annotation_load_fn=dataset_fns.load_caples_annotations,
          description='Annotated Caples recordings from 2018/2019.',
      ),
      SoundscapesConfig(
          name='caples_full_length',
          class_list_name='caples',
          audio_glob='caples/audio/*',
          annotation_filename='caples.csv',
          annotation_load_fn=dataset_fns.load_caples_annotations,
          keep_unknown_annotation=True,
          # Some recordings in Caples are only partially-annotated, so to avoid
          # scoring legitimate model predictions as false positives we pad with
          # "unknown" annotations before the first annotation and after the last
          # annotation.
          full_length_unknown_guard=True,
          description='Full-length annotated Caples recordings from 2018/2019.',
      ),
      SoundscapesConfig(
          name='hawaii',
          audio_glob='hawaii/audio/*.flac',
          interval_length_s=5.0,
          localization_fn=audio_utils.slice_peaked_audio,
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'Fully annotated Hawaii recordings. '
              'https://zenodo.org/record/7078499'
          ),
          class_list_name='hawaii',
      ),
      SoundscapesConfig(
          name='hawaii_full_length',
          audio_glob='hawaii/audio/*.flac',
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'Full-length, fully annotated Hawaii recordings. '
              'https://zenodo.org/record/7078499'
          ),
          class_list_name='hawaii',
      ),
      SoundscapesConfig(
          name='ssw',
          audio_glob='ssw/audio/*.flac',
          interval_length_s=5.0,
          localization_fn=audio_utils.slice_peaked_audio,
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          description=(
              'Annotated Sapsucker Woods recordings. '
              'https://zenodo.org/record/7018484'
          ),
          class_list_name='ssw',
      ),
      SoundscapesConfig(
          name='ssw_full_length',
          audio_glob='ssw/audio/*.flac',
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'Full-length, annotated Sapsucker Woods recordings. '
              'https://zenodo.org/record/7018484'
          ),
          class_list_name='ssw',
      ),
      SoundscapesConfig(
          name='coffee_farms',
          audio_glob='coffee_farms/audio/*.flac',
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          interval_length_s=5.0,
          localization_fn=audio_utils.slice_peaked_audio,
          keep_unknown_annotation=True,
          description=(
              'Colombian and Costa Rican coffee farm recordings. '
              'https://zenodo.org/record/7525349'
          ),
          class_list_name='coffee_farms',
      ),
      SoundscapesConfig(
          name='coffee_farms_full_length',
          audio_glob='coffee_farms/audio/*.flac',
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'Full-length Colombian and Costa Rican coffee farm recordings. '
              'https://zenodo.org/record/7525349'
          ),
          class_list_name='coffee_farms',
      ),
      SoundscapesConfig(
          name='high_sierras',
          audio_glob='high_sierras/audio/*.flac',
          interval_length_s=5.0,
          localization_fn=audio_utils.slice_peaked_audio,
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'High Sierras recordings. https://zenodo.org/record/7525805'
          ),
          class_list_name='high_sierras',
      ),
      SoundscapesConfig(
          name='high_sierras_full_length',
          audio_glob='high_sierras/audio/*.flac',
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'Full-length High Sierras recordings. '
              'https://zenodo.org/record/7525805'
          ),
          class_list_name='high_sierras',
      ),
      SoundscapesConfig(
          name='sierras_kahl',
          audio_glob='sierras_kahl/audio/*.flac',
          interval_length_s=5.0,
          localization_fn=audio_utils.slice_peaked_audio,
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'Sierra Nevada recordings. https://zenodo.org/record/7050014'
          ),
          class_list_name='sierras_kahl',
      ),
      SoundscapesConfig(
          name='sierras_kahl_full_length',
          audio_glob='sierras_kahl/audio/*.flac',
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'Full-length Sierra Nevada recordings. '
              'https://zenodo.org/record/7050014'
          ),
          class_list_name='sierras_kahl',
      ),
      SoundscapesConfig(
          name='powdermill',
          audio_glob='powdermill/*/*.wav',
          interval_length_s=5.0,
          localization_fn=audio_utils.slice_peaked_audio,
          annotation_load_fn=dataset_fns.load_powdermill_annotations,
          description=(
              'New England recordings from Powdermill Nature Reserve, '
              'Rector, PA. https://doi.org/10.1002/ecy.3329'
          ),
          class_list_name='powdermill',
      ),
      SoundscapesConfig(
          name='powdermill_full_length',
          audio_glob='powdermill/*/*.wav',
          annotation_filename='powdermill.csv',
          annotation_load_fn=dataset_fns.load_powdermill_annotations,
          keep_unknown_annotation=True,
          description=(
              'Full-length New England recordings from Powdermill '
              'Nature Reserve, Rector, PA. '
              'https://doi.org/10.1002/ecy.3329'
          ),
          class_list_name='powdermill',
      ),
      SoundscapesConfig(
          name='peru',
          audio_glob='peru/audio/*.flac',
          interval_length_s=5.0,
          localization_fn=audio_utils.slice_peaked_audio,
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'Soundscapes from the SW Amazon basin. '
              'https://zenodo.org/record/7079124#.YypL8-xufhM'
          ),
          class_list_name='peru',
      ),
      SoundscapesConfig(
          name='peru_full_length',
          audio_glob='peru/audio/*.flac',
          annotation_load_fn=dataset_fns.load_cornell_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'Full-length soundscapes from the SW Amazon basin. '
              'https://zenodo.org/record/7079124#.YypL8-xufhM'
          ),
          class_list_name='peru',
      ),
      SoundscapesConfig(
          name='weldy_calltype_full_length',
          audio_glob='weldy_calltype/annotated_recordings/*.wav',
          annotation_load_fn=dataset_fns.load_weldy_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'Full-length annotated bird call types from the PNW. '
              'https://zenodo.org/record/8047850'
          ),
          class_list_name='weldy_calltype',
      ),
      SoundscapesConfig(
          name='anuraset_full_length',
          audio_glob='anuraset/raw_data/*/*.wav',
          annotation_load_fn=dataset_fns.load_anuraset_annotations,
          annotation_filename='annotations.csv',
          keep_unknown_annotation=True,
          description=(
              'Full-length annotated frog vocalizations. '
              'https://zenodo.org/record/8056090'
          ),
          class_list_name='anuraset',
      ),
  ]

  def _info(self) -> tfds.core.DatasetInfo:
    dataset_class_list = soundscapes_lib.load_class_list(
        self.builder_config.class_list_name,
        self.builder_config.keep_unknown_annotation,
    )
    logging.info(
        'Currently considering a total of %s species.',
        len(dataset_class_list.classes),
    )
    full_length = self.builder_config.localization_fn is None
    if full_length:
      audio_feature_shape = [None]
    else:
      audio_feature_shape = [
          int(
              self.builder_config.sample_rate_hz
              * self.builder_config.interval_length_s
          )
      ]
    common_features = {
        'audio': tfds_features.Int16AsFloatTensor(
            shape=audio_feature_shape,
            sample_rate=self.builder_config.sample_rate_hz,
            encoding=tfds.features.Encoding.ZLIB,
        ),
        'label': tfds.features.Sequence(
            tfds.features.ClassLabel(names=dataset_class_list.classes)
        ),
        'filename': tfds.features.Text(),
        'recording_id': tfds.features.Scalar(dtype=tf.uint64),
        'segment_id': tfds.features.Scalar(dtype=tf.int64),
        'segment_start': tfds.features.Scalar(dtype=tf.uint64),
        'segment_end': tfds.features.Scalar(dtype=tf.uint64),
    }
    if full_length:
      common_features.update({
          'annotation_start': tfds.features.Sequence(
              tfds.features.Scalar(dtype=tf.uint64)
          ),
          'annotation_end': tfds.features.Sequence(
              tfds.features.Scalar(dtype=tf.uint64)
          ),
      })
    if self.builder_config.metadata_load_fn is not None:
      if self.builder_config.metadata_fields is None:
        raise ValueError(
            "If a 'metadata_load_fn' is specified, then the"
            "'metadata_fields' mapping must also be specied."
        )
      additional_features = {
          k.target_key: k.feature_type
          for k in self.builder_config.metadata_fields.values()
      }
      common_features.update(additional_features)
    return tfds.core.DatasetInfo(
        builder=self,
        features=tfds.features.FeaturesDict(common_features),
        supervised_keys=('audio', 'label'),
        homepage='https://github.com/google-research/chirp',
        citation=_CITATION,
    )

  def _split_generators(self, dl_manager: tfds.download.DownloadManager):
    # Defined as part of the tfds API for dividing the dataset into splits.
    # https://www.tensorflow.org/datasets/add_dataset#specifying_dataset_splits
    dl_manager._force_checksums_validation = (  # pylint: disable=protected-access
        False
    )

    # Get the state from the dl_manager which we'll use to create segments.
    all_audio_filenames = self.builder_config.audio_dir.glob(
        self.builder_config.audio_glob
    )
    if self.builder_config.supervised:
      # For supervised data, we first grab the annotated segments.
      filename = (
          self.builder_config.annotation_filename
          or f'{self.builder_config.class_list_name}.csv'
      )
      annotations_path = dl_manager.download_and_extract({
          'segments': (
              self.builder_config.audio_dir
              / self.builder_config.class_list_name
              / filename
          ).as_posix(),
      })['segments']
      annotations_df = self.builder_config.annotation_load_fn(annotations_path)
    else:
      annotations_df = None

    segments = soundscapes_lib.create_segments_df(
        all_audio_filenames,
        annotations_df,
        self.builder_config.supervised,
        self.builder_config.audio_dir,
        self.builder_config.metadata_fields,
        self.builder_config.metadata_load_fn,
    )
    return {
        'train': self._generate_examples(segments=segments),
    }

  def _generate_examples(self, segments: pd.DataFrame):
    """Generate examples from the dataframe of segments.

    Args:
      segments: Dataframe of segments. Each row (=segment) must minimally
        contain the following fields: ['filename', 'url', 'label',
        'start_time_s', 'end_time_s'].

    Returns:
      List of valid segments.
    """
    beam = tfds.core.lazy_imports.apache_beam
    librosa = tfds.core.lazy_imports.librosa
    info = self._info()
    full_length = self.builder_config.localization_fn is None
    # Drop any extraneous columns.
    for k in segments.columns.values:
      if k not in info.features and k not in [
          'url',
          'start_time_s',
          'end_time_s',
      ]:
        segments = segments.drop(k, axis=1)

    def _process_group(
        group: tuple[int, tuple[str, pd.DataFrame]]
    ) -> Sequence[tuple[str, dict[str, Any]]]:
      # Each filename gets a unique ID
      recording_id, (filename, segment_group) = group

      # Each segment in segment_group will generate a tf.Example. A lot of
      # fields, especially metadata ones will be shared between segments.
      # Therefore, we create a template.
      recording_template = segment_group.iloc[0].copy()
      recording_template['recording_id'] = recording_id

      url = recording_template['url']

      # Remove all the fields we don't need from the recording_template. We
      # set errors='ignore' as some fields to be dropped may already not
      # exist.
      recording_template = recording_template.drop(
          ['url', 'start_time_s', 'end_time_s'], errors='ignore'
      ).to_dict()

      # Load the audio associated with this group of segments
      with tempfile.NamedTemporaryFile(mode='w+b', suffix=url.suffix) as f:
        f.write(url.read_bytes())
        # librosa outputs lots of warnings which we can safely ignore when
        # processing all Xeno-Canto files and PySoundFile is unavailable.
        with warnings.catch_warnings():
          warnings.simplefilter('ignore')
          sr = self.builder_config.sample_rate_hz
          try:
            audio, _ = librosa.load(
                f.name, sr=sr, res_type=self.builder_config.resampling_method
            )
          except Exception as inst:  # pylint: disable=broad-except
            # We have no idea what can go wrong in librosa, so we catch a broad
            # exception here.
            logging.warning(
                (
                    'The audio at %s could not be loaded. Following'
                    'exception occured: %s'
                ),
                url,
                inst,
            )
            return []
          if not full_length:
            # We remove all short audios. These short audios are only observed
            # among caples_2020 unlabelled recordings.
            target_length = int(sr * self.builder_config.interval_length_s)
            if len(audio) < target_length:
              logging.warning('Skipping audio at %s because too short.', url)
              return []

          # Resampling can introduce artifacts that push the signal outside the
          # [-1, 1) interval.
          audio = np.clip(audio, -1.0, 1.0 - (1.0 / float(1 << 15)))

      class_list = soundscapes_lib.load_class_list(
          self.builder_config.class_list_name,
          self.builder_config.keep_unknown_annotation,
      )

      if full_length:
        annotations = soundscapes_lib.get_full_length_annotations(
            audio,
            segment_group,
            class_list,
            self.builder_config.sample_rate_hz,
            unknown_guard=self.builder_config.full_length_unknown_guard,
        )
        if annotations.empty:
          return []
        beam.metrics.Metrics.counter('soundscapes', 'examples').inc()
        return [(
            filename,
            {
                **recording_template,
                'label': annotations['label'].tolist(),
                'audio': audio,
                'segment_start': 0,
                'segment_end': len(audio),
                'segment_id': 0,
                'annotation_start': annotations['annotation_start'].tolist(),
                'annotation_end': annotations['annotation_end'].tolist(),
            },
        )]
      else:
        labeled_intervals = soundscapes_lib.get_labeled_intervals(
            audio,
            segment_group,
            class_list,
            self.builder_config.sample_rate_hz,
            self.builder_config.interval_length_s,
            self.builder_config.localization_fn,
            self.builder_config.supervised,
        )

        # Create a tf.Example for every segment.
        valid_segments = []
        for index, ((start, end), segment_labels) in enumerate(
            labeled_intervals.items()
        ):
          key = f'{filename}_{index}'
          valid_segments.append((
              key,
              {
                  **recording_template,
                  'label': list(segment_labels),
                  'audio': audio[start:end],
                  'segment_start': start,
                  'segment_end': end,
                  'segment_id': index,
              },
          ))
          beam.metrics.Metrics.counter('soundscapes', 'examples').inc()
        return valid_segments

    for group in enumerate(segments.groupby('filename')):
      for key, example in _process_group(group):
        yield key, example
