# -*- coding: utf-8 -*-
"""
Editor Spyder

Este é um arquivo de script temporário.
"""

import geopandas as gp
import matplotlib.pyplot as plt
from shapely.geometry import Point, Polygon


def extratic_latlon(list):
    area=[(-48,-28), (-48,-22),(-41,-22),(-41,-28)]
    area = Polygon(shell=area)
    target_points=[]
    for coords in list:
        for coord in coords:
            point = Point(coord)
            if area.contains(point):
                print(point)
                target_points.append(point)
    return target_points


def extratic_polygon(multipolygon):
    return list(multipolygon)

#Armadilha muito importante ao lidar com shape file,
#todos precisam estar juntos na mesma pasta
brasil = gp.read_file('/home/ladsin/Área de Trabalho/GisPython/Brasil/UFEBRASIL.shp')
coords = []
lat_lon = []
for polygon in brasil.geometry:
    if polygon.type == 'MultiPolygon':
        for i in range(len(polygon)):
            coords.append(polygon[i])
    else:
        coords.append(polygon)
for coord in range(len(coords)):
    lat_lon.append(coords[coord].exterior.coords.xy)
coords = []
for i in lat_lon:
    coords.append(extratic_latlon(i))