poly_partition

defmodule PolyPartition do
  alias PolyPartition.Helpers
  alias PolyPartition.Output
  alias PolyPartition.Geometry
  @moduledoc """
  Documentation for PolyBisector.
  """


  def split_list(list, bound, prior) do
    cond do
      list == prior -> list
      true ->
        result = []
        Enum.flat_map(list, fn(x) ->
          cond do
            Geometry.area(x) > bound -> result ++ Helpers.split(x, 0)
            true -> result ++ [x]
          end
        end)
        |> split_list(bound, list)
    end
  end

  def split_polys(list, bound) do
    list
    |> Enum.filter(fn(x) -> Geometry.area(x) > 0 end)
    |> split_list(bound, [])
  end

  @doc """
    Takes a MultiPolygon and returns a list of polygons (_not_ linear rings!)
  """
  def getPolys(input) do
    for n <- input.coordinates do
      cond do
        List.first(hd(n)) == List.last(hd(n)) -> tl(hd(n))
        true -> hd(n) #NO HOLES IN GONS
      end
    end
  end

  @doc """
    Takes a list of multipolygons, gives a list of polygons
  """
  def getAllPolys(input_list) do
    for n <- input_list do
      getPolys(n)
    end
    |> List.foldr([], fn(x, acc) -> acc ++ x end)
  end

  def get_split(input, bound) do
    getAllPolys(input)
    |> split_polys(bound)
    |> Enum.map(fn(x) -> Output.mPProducer(x) end)
  end

  def get_split_geojson(input, bound) do
    getAllPolys(input)
    |> split_polys(bound)
    |> Enum.map(fn(x) -> Output.mPProducer(x) |> Output.wrap_geoJSON_feature end)
    |> Output.wrap_geoJSON
    |> Poison.encode
  end
end