val philly: Polygon = GeoJson.fromFile[Polygon]("src/test/resources/philly.json")

  val images: Seq[LandsatImage] =
    Landsat8Query()
      .withStartDate(ZonedDateTime.of(2012, 1, 12, 0, 0, 0, 0, ZoneOffset.UTC))
      .withEndDate(ZonedDateTime.of(2015, 11, 5, 0, 0, 0, 0, ZoneOffset.UTC))
      .withMaxCloudCoverage(80.0)
      .intersects(philly)
      .collect() match {
        case Success(r) => r
        case Failure(e) => throw e
      }

  val filtered: Seq[LandsatImage] = images.filter(_.imageExistsS3())

  filtered
    .foreach { image =>
      println(s"${image.sceneId} - ${image.cloudPercentage}% clouds")
      println(s"\tAquisition Date: ${image.aquisitionDate}")
      println(s"\tUSGS Thumbnail URL: ${image.thumbnailUrl}")
      println(s"\tAWS large thumbnail URL: ${image.largeThumbnail}")
      println(s"\tGoogle URL: ${image.googleUrl}")
      println(s"\tFootprint GeoJSON: ${image.footprint.toGeoJson}")
    }

  println(s"Results: ${filtered.size} images.")

  implicit val d = scala.concurrent.duration.Duration(30, scala.concurrent.duration.SECONDS)

  val result: Try[QueryResult] =
    Landsat8Query()
      .intersects(-75.65185546874999, 39.69701710019832, -74.718017578125, 40.24009510908543)
      .execute()

  result match {
    case Success(r) =>
      for (image <- r.images.take(1)) {
        println(image.thumbnailUrl)
        println(image.largeThumbnail)
        println(image.smallThumbnail)
      }
      println(s"RESULT COUNT: ${r.metadata.found}")
    case Failure(_) =>
      println("No results found!")
  }