Processamento de Dados.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Os códigos apresentados abaixo fazem parte da série de reportagens sobre grilagem de terras no litoral do Ceará, disponível [neste link](https://mais.opovo.com.br/reportagens-especiais/itarema-especulacao-imobiliaria/):\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import geopandas as gpd\n",
    "from shapely.geometry import Point, Polygon, LineString\n",
    "from pyproj import Geod\n",
    "from shapely import wkt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "pd.options.display.float_format = '{:.2f}'.format"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/alexandre/anaconda3/envs/credos/lib/python3.7/site-packages/geopandas/array.py:93: ShapelyDeprecationWarning: __len__ for multi-part geometries is deprecated and will be removed in Shapely 2.0. Check the length of the `geoms` property instead to get the  number of parts of a multi-part geometry.\n",
      "  aout[:] = out\n"
     ]
    }
   ],
   "source": [
    "assentamentos_federais_ce = gpd.read_file('assentamento_federal_ce.geojson', crs='EPSG:4326')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>cd_sipra</th>\n",
       "      <th>uf</th>\n",
       "      <th>nome_proje</th>\n",
       "      <th>municipio</th>\n",
       "      <th>area_hecta</th>\n",
       "      <th>capacidade</th>\n",
       "      <th>num_famili</th>\n",
       "      <th>fase</th>\n",
       "      <th>data_de_cr</th>\n",
       "      <th>forma_obte</th>\n",
       "      <th>data_obten</th>\n",
       "      <th>area_calc_</th>\n",
       "      <th>sr</th>\n",
       "      <th>descricao_</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>185</th>\n",
       "      <td>CE0092000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA PACHICU</td>\n",
       "      <td>ITAREMA</td>\n",
       "      <td>619.6443</td>\n",
       "      <td>34</td>\n",
       "      <td>28</td>\n",
       "      <td>5</td>\n",
       "      <td>23/11/1995</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>24/03/1995</td>\n",
       "      <td>613.49</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.72711 -2.99860, -39.72305 ...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>348</th>\n",
       "      <td>CE0097000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA PATOS</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>1598.6423</td>\n",
       "      <td>30</td>\n",
       "      <td>30</td>\n",
       "      <td>5</td>\n",
       "      <td>12/12/1995</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>26/12/1994</td>\n",
       "      <td>1599.67</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.68994 -3.03521, -39.69042 ...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      cd_sipra  uf  nome_proje municipio area_hecta  capacidade  num_famili  \\\n",
       "185  CE0092000  CE  PA PACHICU   ITAREMA   619.6443          34          28   \n",
       "348  CE0097000  CE    PA PATOS  AMONTADA  1598.6423          30          30   \n",
       "\n",
       "     fase  data_de_cr      forma_obte  data_obten  area_calc_    sr  \\\n",
       "185     5  23/11/1995  Desapropriação  24/03/1995      613.49  None   \n",
       "348     5  12/12/1995  Desapropriação  26/12/1994     1599.67  None   \n",
       "\n",
       "    descricao_                                           geometry  \n",
       "185       None  MULTIPOLYGON (((-39.72711 -2.99860, -39.72305 ...  \n",
       "348       None  MULTIPOLYGON (((-39.68994 -3.03521, -39.69042 ...  "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "assentamentos = assentamentos_federais_ce.loc[(assentamentos_federais_ce.nome_proje == 'PA PATOS') | (assentamentos_federais_ce.nome_proje == 'PA PACHICU')]\n",
    "assentamentos"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Marco</th>\n",
       "      <th>X_Este</th>\n",
       "      <th>Y_Norte</th>\n",
       "      <th>Lat</th>\n",
       "      <th>Lng</th>\n",
       "      <th>Confrontando</th>\n",
       "      <th>Azimute</th>\n",
       "      <th>Distância</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>418.610,0029</td>\n",
       "      <td>9.669.427,0674</td>\n",
       "      <td>-2,9905242189790</td>\n",
       "      <td>-39,73240215</td>\n",
       "      <td>ASSOCIAÇÃO DOS PRODUTORES RURAIS DA FAZENDA PATOS</td>\n",
       "      <td>213°32'49\"</td>\n",
       "      <td>368,429</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>418.406,4010</td>\n",
       "      <td>9669120,00660</td>\n",
       "      <td>-2,9933007588199</td>\n",
       "      <td>-39,73423604</td>\n",
       "      <td>ASSOCIAÇÃO DOS PRODUTORES RURAIS DA FAZENDA PATOS</td>\n",
       "      <td>120°43'53\"</td>\n",
       "      <td>634,036</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>418.951,4011</td>\n",
       "      <td>9668796,00660</td>\n",
       "      <td>-2,9962350574347</td>\n",
       "      <td>-39,72933396</td>\n",
       "      <td>ASSOCIAÇÃO DOS PRODUTORES RURAIS DA FAZENDA PATOS</td>\n",
       "      <td>191°57'24\"</td>\n",
       "      <td>521,31</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4</td>\n",
       "      <td>418.843,4010</td>\n",
       "      <td>9668286,00660</td>\n",
       "      <td>-3,0008480327824</td>\n",
       "      <td>-39,73030883</td>\n",
       "      <td>FRANCISCO DAS CHARGAS SALES</td>\n",
       "      <td>114°07'56\"</td>\n",
       "      <td>1.446,96</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>5</td>\n",
       "      <td>420.163,9070</td>\n",
       "      <td>9667694,42710</td>\n",
       "      <td>-3,0062075737095</td>\n",
       "      <td>-39,71843005</td>\n",
       "      <td>MANGUEZAL</td>\n",
       "      <td>13°39'00\"</td>\n",
       "      <td>474,989</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>6</td>\n",
       "      <td>420.276,0000</td>\n",
       "      <td>9668156,00000</td>\n",
       "      <td>-3,0020326899996</td>\n",
       "      <td>-39,717418669</td>\n",
       "      <td>MANGUEZAL</td>\n",
       "      <td>100°24'28\"</td>\n",
       "      <td>49,82</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>7</td>\n",
       "      <td>420.325,0000</td>\n",
       "      <td>9668147,00000</td>\n",
       "      <td>-3,0021143977676</td>\n",
       "      <td>-39,716977805</td>\n",
       "      <td>MANGUEZAL</td>\n",
       "      <td>140°05'57\"</td>\n",
       "      <td>230,723</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>8</td>\n",
       "      <td>420.473,0000</td>\n",
       "      <td>9667970,00000</td>\n",
       "      <td>-3,0037164790995</td>\n",
       "      <td>-39,715647098</td>\n",
       "      <td>ÁREA PRESUMIVEL DE MARINHA</td>\n",
       "      <td>92°49'51\"</td>\n",
       "      <td>567,61</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>9</td>\n",
       "      <td>421.039,9172</td>\n",
       "      <td>9667941,96860</td>\n",
       "      <td>-3,0039734062162</td>\n",
       "      <td>-39,710545955</td>\n",
       "      <td>ÁREA PRESUMIVEL DE MARINHA</td>\n",
       "      <td>302°43'23\"</td>\n",
       "      <td>965,825</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>10</td>\n",
       "      <td>420.227,3738</td>\n",
       "      <td>9668464,07150</td>\n",
       "      <td>-2,9992454788453</td>\n",
       "      <td>-39,717854404</td>\n",
       "      <td>ÁREA PRESUMIVEL DE MARINHA</td>\n",
       "      <td>300°46'12\"</td>\n",
       "      <td>1.882,35</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Marco        X_Este         Y_Norte               Lat            Lng  \\\n",
       "0      1  418.610,0029  9.669.427,0674  -2,9905242189790   -39,73240215   \n",
       "1      2  418.406,4010   9669120,00660  -2,9933007588199   -39,73423604   \n",
       "2      3  418.951,4011   9668796,00660  -2,9962350574347   -39,72933396   \n",
       "3      4  418.843,4010   9668286,00660  -3,0008480327824   -39,73030883   \n",
       "4      5  420.163,9070   9667694,42710  -3,0062075737095   -39,71843005   \n",
       "5      6  420.276,0000   9668156,00000  -3,0020326899996  -39,717418669   \n",
       "6      7  420.325,0000   9668147,00000  -3,0021143977676  -39,716977805   \n",
       "7      8  420.473,0000   9667970,00000  -3,0037164790995  -39,715647098   \n",
       "8      9  421.039,9172   9667941,96860  -3,0039734062162  -39,710545955   \n",
       "9     10  420.227,3738   9668464,07150  -2,9992454788453  -39,717854404   \n",
       "\n",
       "                                        Confrontando       Azimute Distância  \n",
       "0  ASSOCIAÇÃO DOS PRODUTORES RURAIS DA FAZENDA PATOS    213°32'49\"   368,429  \n",
       "1  ASSOCIAÇÃO DOS PRODUTORES RURAIS DA FAZENDA PATOS   120°43'53\"    634,036  \n",
       "2  ASSOCIAÇÃO DOS PRODUTORES RURAIS DA FAZENDA PATOS    191°57'24\"    521,31  \n",
       "3                        FRANCISCO DAS CHARGAS SALES   114°07'56\"   1.446,96  \n",
       "4                                          MANGUEZAL    13°39'00\"    474,989  \n",
       "5                                          MANGUEZAL    100°24'28\"     49,82  \n",
       "6                                          MANGUEZAL    140°05'57\"   230,723  \n",
       "7                         ÁREA PRESUMIVEL DE MARINHA     92°49'51\"    567,61  \n",
       "8                         ÁREA PRESUMIVEL DE MARINHA    302°43'23\"   965,825  \n",
       "9                         ÁREA PRESUMIVEL DE MARINHA    300°46'12\"  1.882,35  "
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "imovel_146 = pd.read_csv('imovel_irregular_146.csv')\n",
    "imovel_146"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/alexandre/anaconda3/envs/credos/lib/python3.7/site-packages/pandas/core/dtypes/cast.py:112: ShapelyDeprecationWarning: The array interface is deprecated and will no longer work in Shapely 2.0. Convert the '.coords' to a numpy array instead.\n",
      "  values = construct_1d_object_array_from_listlike(values)\n"
     ]
    }
   ],
   "source": [
    "imovel_146['Lat'] = imovel_146['Lat'].apply(lambda x: float(x.replace(',','.')))\n",
    "imovel_146['Lng'] = imovel_146['Lng'].apply(lambda x: float(x.replace(',','.')))\n",
    "imovel_146['geometry'] = [Point(xy) for xy in zip(imovel_146.Lng, imovel_146.Lat)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"100.0\" height=\"100.0\" viewBox=\"-39.7351836434 -3.0071551771095 0.0255852917999988 0.017578561530499748\" preserveAspectRatio=\"xMinYMin meet\"><g transform=\"matrix(1,0,0,-1,0,-5.9967317926885)\"><path fill-rule=\"evenodd\" fill=\"#66cc99\" stroke=\"#555555\" stroke-width=\"0.0005117058359999759\" opacity=\"0.6\" d=\"M -39.73240215,-2.990524218979 L -39.73423604,-2.9933007588199 L -39.72933396,-2.9962350574347 L -39.73030883,-3.0008480327824 L -39.71843005,-3.0062075737095 L -39.717418669,-3.0020326899996 L -39.716977805,-3.0021143977676 L -39.715647098,-3.0037164790995 L -39.710545955,-3.0039734062162 L -39.717854404,-2.9992454788453 L -39.73240215,-2.990524218979 z\" /></g></svg>"
      ],
      "text/plain": [
       "<shapely.geometry.polygon.Polygon at 0x7f07c0e07c10>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "lat_list = imovel_146['Lat'].to_list()\n",
    "lng_list = imovel_146['Lng'].to_list()\n",
    "polygon_geom = Polygon(zip(lng_list, lat_list))\n",
    "polygon_geom"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"100.0\" height=\"100.0\" viewBox=\"-39.74440086484413 -3.043104256772639 0.030166692380944937 0.04621374779590104\" preserveAspectRatio=\"xMinYMin meet\"><g transform=\"matrix(1,0,0,-1,0,-6.039994765749377)\"><g><path fill-rule=\"evenodd\" fill=\"#66cc99\" stroke=\"#555555\" stroke-width=\"0.0009242749559180207\" opacity=\"0.6\" d=\"M -39.72711044931032,-2.998602129265475 L -39.723053934257884,-3.000098109155095 L -39.718383198832484,-3.002310320053523 L -39.715945792751924,-3.003916931295058 L -39.71787406274631,-3.007388283613164 L -39.72556430957628,-3.021232571996849 L -39.72816493671804,-3.027147665971369 L -39.73013100221043,-3.03316829484196 L -39.73257055897797,-3.041392636483902 L -39.73539269607216,-3.039261115697143 L -39.73685690509973,-3.037079500029984 L -39.738704357252175,-3.03503150611387 L -39.74147321897293,-3.032804936007722 L -39.74268924455539,-3.032096496782009 L -39.736940521109254,-3.018734428581933 L -39.73638183567188,-3.017435758237033 L -39.73659728523881,-3.017343610662421 L -39.73652863564981,-3.017088370395955 L -39.736636363877004,-3.01684694146638 L -39.73689138135424,-3.016717857991663 L -39.737186978516924,-3.016851542222681 L -39.73737450881464,-3.017124341768274 L -39.73770719975062,-3.017554990947939 L -39.73828162177483,-3.017502402155769 L -39.738529975931954,-3.017633494712282 L -39.73888179836505,-3.017611633959561 L -39.73895511579634,-3.017449564527703 L -39.73921448586685,-3.017503484382652 L -39.7393580586843,-3.017442052372302 L -39.73949810598334,-3.017489902150166 L -39.73976576180618,-3.017421419670393 L -39.73999380146504,-3.017584730584499 L -39.74022446568323,-3.017504151227759 L -39.74050647929881,-3.017510200315433 L -39.738640927914616,-3.014969190406454 L -39.7308847718396,-3.003889534156097 L -39.72949401207949,-3.002288091512774 L -39.72841597139564,-3.001267754195892 L -39.72828351767418,-3.000997357475639 L -39.72711044931032,-2.998602129265475 z\" /></g></g></svg>"
      ],
      "text/plain": [
       "<shapely.geometry.multipolygon.MultiPolygon at 0x7f07c0943610>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "assentamentos.iloc[0].geometry"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "True"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "polygon_geom.intersects(assentamentos.iloc[0].geometry)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.73240 -2.99052, -39.73424 -2.993...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry\n",
       "0  POLYGON ((-39.73240 -2.99052, -39.73424 -2.993..."
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_imovel_146 = gpd.GeoDataFrame(index=[0],crs='EPSG:4326', geometry=[polygon_geom])\n",
    "gdf_imovel_146"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "resultado =  gpd.sjoin(assentamentos, gdf_imovel_146)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.001301363320369968"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "resultado.iloc[1].geometry.area"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>cd_sipra</th>\n",
       "      <th>uf</th>\n",
       "      <th>nome_proje</th>\n",
       "      <th>municipio</th>\n",
       "      <th>area_hecta</th>\n",
       "      <th>capacidade</th>\n",
       "      <th>num_famili</th>\n",
       "      <th>fase</th>\n",
       "      <th>data_de_cr</th>\n",
       "      <th>forma_obte</th>\n",
       "      <th>data_obten</th>\n",
       "      <th>area_calc_</th>\n",
       "      <th>sr</th>\n",
       "      <th>descricao_</th>\n",
       "      <th>geometry</th>\n",
       "      <th>index_right</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>185</th>\n",
       "      <td>CE0092000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA PACHICU</td>\n",
       "      <td>ITAREMA</td>\n",
       "      <td>619.6443</td>\n",
       "      <td>34</td>\n",
       "      <td>28</td>\n",
       "      <td>5</td>\n",
       "      <td>23/11/1995</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>24/03/1995</td>\n",
       "      <td>613.49</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.72711 -2.99860, -39.72305 ...</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>348</th>\n",
       "      <td>CE0097000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA PATOS</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>1598.6423</td>\n",
       "      <td>30</td>\n",
       "      <td>30</td>\n",
       "      <td>5</td>\n",
       "      <td>12/12/1995</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>26/12/1994</td>\n",
       "      <td>1599.67</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.68994 -3.03521, -39.69042 ...</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      cd_sipra  uf  nome_proje municipio area_hecta  capacidade  num_famili  \\\n",
       "185  CE0092000  CE  PA PACHICU   ITAREMA   619.6443          34          28   \n",
       "348  CE0097000  CE    PA PATOS  AMONTADA  1598.6423          30          30   \n",
       "\n",
       "     fase  data_de_cr      forma_obte  data_obten  area_calc_    sr  \\\n",
       "185     5  23/11/1995  Desapropriação  24/03/1995      613.49  None   \n",
       "348     5  12/12/1995  Desapropriação  26/12/1994     1599.67  None   \n",
       "\n",
       "    descricao_                                           geometry  index_right  \n",
       "185       None  MULTIPOLYGON (((-39.72711 -2.99860, -39.72305 ...            0  \n",
       "348       None  MULTIPOLYGON (((-39.68994 -3.03521, -39.69042 ...            0  "
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "resultado"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.73240 -2.99052, -39.73424 -2.993...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry\n",
       "0  POLYGON ((-39.73240 -2.99052, -39.73424 -2.993..."
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_imovel_146"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0    POLYGON ((-39.71843 -3.00621, -39.72873 -3.001...\n",
       "dtype: geometry"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_imovel_146.geometry.intersection(assentamentos.iloc[0].geometry)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0    POLYGON ((-39.73240 -2.99052, -39.73424 -2.993...\n",
       "Name: geometry, dtype: geometry"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_imovel_146.geometry"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"100.0\" height=\"100.0\" viewBox=\"-39.74440086484413 -3.043104256772639 0.030166692380944937 0.04621374779590104\" preserveAspectRatio=\"xMinYMin meet\"><g transform=\"matrix(1,0,0,-1,0,-6.039994765749377)\"><g><path fill-rule=\"evenodd\" fill=\"#66cc99\" stroke=\"#555555\" stroke-width=\"0.0009242749559180207\" opacity=\"0.6\" d=\"M -39.72711044931032,-2.998602129265475 L -39.723053934257884,-3.000098109155095 L -39.718383198832484,-3.002310320053523 L -39.715945792751924,-3.003916931295058 L -39.71787406274631,-3.007388283613164 L -39.72556430957628,-3.021232571996849 L -39.72816493671804,-3.027147665971369 L -39.73013100221043,-3.03316829484196 L -39.73257055897797,-3.041392636483902 L -39.73539269607216,-3.039261115697143 L -39.73685690509973,-3.037079500029984 L -39.738704357252175,-3.03503150611387 L -39.74147321897293,-3.032804936007722 L -39.74268924455539,-3.032096496782009 L -39.736940521109254,-3.018734428581933 L -39.73638183567188,-3.017435758237033 L -39.73659728523881,-3.017343610662421 L -39.73652863564981,-3.017088370395955 L -39.736636363877004,-3.01684694146638 L -39.73689138135424,-3.016717857991663 L -39.737186978516924,-3.016851542222681 L -39.73737450881464,-3.017124341768274 L -39.73770719975062,-3.017554990947939 L -39.73828162177483,-3.017502402155769 L -39.738529975931954,-3.017633494712282 L -39.73888179836505,-3.017611633959561 L -39.73895511579634,-3.017449564527703 L -39.73921448586685,-3.017503484382652 L -39.7393580586843,-3.017442052372302 L -39.73949810598334,-3.017489902150166 L -39.73976576180618,-3.017421419670393 L -39.73999380146504,-3.017584730584499 L -39.74022446568323,-3.017504151227759 L -39.74050647929881,-3.017510200315433 L -39.738640927914616,-3.014969190406454 L -39.7308847718396,-3.003889534156097 L -39.72949401207949,-3.002288091512774 L -39.72841597139564,-3.001267754195892 L -39.72828351767418,-3.000997357475639 L -39.72711044931032,-2.998602129265475 z\" /></g></g></svg>"
      ],
      "text/plain": [
       "<shapely.geometry.multipolygon.MultiPolygon at 0x7f07c0943610>"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "assentamentos.iloc[0].geometry"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "      <th>intersecao_cd_sipra</th>\n",
       "      <th>intersecao_nome</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...</td>\n",
       "      <td>CE0092000</td>\n",
       "      <td>PA PACHICU</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry intersecao_cd_sipra  \\\n",
       "0  POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...           CE0092000   \n",
       "\n",
       "  intersecao_nome  \n",
       "0      PA PACHICU  "
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cruzamento1 = assentamentos.iloc[0].geometry.intersection(gdf_imovel_146.iloc[0].geometry)\n",
    "cruzamento1 = gpd.GeoDataFrame(index=[0],crs='EPSG:4326', geometry=[cruzamento1])\n",
    "cruzamento1['intersecao_cd_sipra'] = assentamentos.iloc[0].cd_sipra\n",
    "cruzamento1['intersecao_nome'] = assentamentos.iloc[0].nome_proje\n",
    "cruzamento1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "      <th>intersecao_cd_sipra</th>\n",
       "      <th>intersecao_nome</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...</td>\n",
       "      <td>CE0097000</td>\n",
       "      <td>PA PATOS</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry intersecao_cd_sipra  \\\n",
       "0  POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...           CE0097000   \n",
       "\n",
       "  intersecao_nome  \n",
       "0        PA PATOS  "
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cruzamento2 = assentamentos.iloc[1].geometry.intersection(gdf_imovel_146.iloc[0].geometry)\n",
    "cruzamento2 = gpd.GeoDataFrame(index=[0],crs='EPSG:4326', geometry=[cruzamento2])\n",
    "cruzamento2['intersecao_cd_sipra'] = assentamentos.iloc[1].cd_sipra\n",
    "cruzamento2['intersecao_nome'] = assentamentos.iloc[1].nome_proje\n",
    "cruzamento2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "      <th>intersecao_cd_sipra</th>\n",
       "      <th>intersecao_nome</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...</td>\n",
       "      <td>CE0092000</td>\n",
       "      <td>PA PACHICU</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...</td>\n",
       "      <td>CE0097000</td>\n",
       "      <td>PA PATOS</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry intersecao_cd_sipra  \\\n",
       "0  POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...           CE0092000   \n",
       "0  POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...           CE0097000   \n",
       "\n",
       "  intersecao_nome  \n",
       "0      PA PACHICU  \n",
       "0        PA PATOS  "
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "intersecoes = pd.concat([cruzamento1,cruzamento2])\n",
    "intersecoes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "geod = Geod(ellps=\"WGS84\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "def getArea(pol):\n",
    "    return abs(geod.geometry_area_perimeter(pol)[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0   47.17\n",
       "0    0.32\n",
       "Name: geometry, dtype: float64"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "intersecoes.geometry.apply(lambda x: getArea(x)/10000)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "intersecoes['intersecao_area_km2'] = intersecoes.geometry.apply(lambda x: getArea(x)/10000)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "      <th>intersecao_cd_sipra</th>\n",
       "      <th>intersecao_nome</th>\n",
       "      <th>intersecao_area_km2</th>\n",
       "      <th>intersecao_imovel</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...</td>\n",
       "      <td>CE0092000</td>\n",
       "      <td>PA PACHICU</td>\n",
       "      <td>47.17</td>\n",
       "      <td>146</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...</td>\n",
       "      <td>CE0097000</td>\n",
       "      <td>PA PATOS</td>\n",
       "      <td>0.32</td>\n",
       "      <td>146</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry intersecao_cd_sipra  \\\n",
       "0  POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...           CE0092000   \n",
       "0  POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...           CE0097000   \n",
       "\n",
       "  intersecao_nome  intersecao_area_km2 intersecao_imovel  \n",
       "0      PA PACHICU                47.17               146  \n",
       "0        PA PATOS                 0.32               146  "
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "intersecoes['intersecao_imovel'] = '146'\n",
    "intersecoes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>ponto</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>421296</td>\n",
       "      <td>9667851</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>421399</td>\n",
       "      <td>9667416</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>421292</td>\n",
       "      <td>9667416</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4</td>\n",
       "      <td>421134</td>\n",
       "      <td>9666927</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>5</td>\n",
       "      <td>420869</td>\n",
       "      <td>9667016</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>6</td>\n",
       "      <td>419569</td>\n",
       "      <td>9667596</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>7</td>\n",
       "      <td>418781</td>\n",
       "      <td>9667950</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>8</td>\n",
       "      <td>418889</td>\n",
       "      <td>9668460</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>9</td>\n",
       "      <td>418344</td>\n",
       "      <td>9668784</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>10</td>\n",
       "      <td>418651</td>\n",
       "      <td>9669247</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>11</td>\n",
       "      <td>418999</td>\n",
       "      <td>9668974</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>12</td>\n",
       "      <td>419604</td>\n",
       "      <td>9668624</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    ponto       x        y\n",
       "0       1  421296  9667851\n",
       "1       2  421399  9667416\n",
       "2       3  421292  9667416\n",
       "3       4  421134  9666927\n",
       "4       5  420869  9667016\n",
       "5       6  419569  9667596\n",
       "6       7  418781  9667950\n",
       "7       8  418889  9668460\n",
       "8       9  418344  9668784\n",
       "9      10  418651  9669247\n",
       "10     11  418999  9668974\n",
       "11     12  419604  9668624"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "imovel_297 = pd.read_csv('imovel_irregular_297.csv')\n",
    "imovel_297"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "lat_list = imovel_297['y'].to_list()\n",
    "lng_list = imovel_297['x'].to_list()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"300\" height=\"300\" viewBox=\"418221.8 9666804.8 3299.4000000000233 2564.39999999851\" preserveAspectRatio=\"xMinYMin meet\"><g transform=\"matrix(1,0,0,-1,0,19336174.0)\"><path fill-rule=\"evenodd\" fill=\"#66cc99\" stroke=\"#555555\" stroke-width=\"21.996000000000155\" opacity=\"0.6\" d=\"M 421296.0,9667851.0 L 421399.0,9667416.0 L 421292.0,9667416.0 L 421134.0,9666927.0 L 420869.0,9667016.0 L 419569.0,9667596.0 L 418781.0,9667950.0 L 418889.0,9668460.0 L 418344.0,9668784.0 L 418651.0,9669247.0 L 418999.0,9668974.0 L 419604.0,9668624.0 L 421296.0,9667851.0 z\" /></g></svg>"
      ],
      "text/plain": [
       "<shapely.geometry.polygon.Polygon at 0x7f07c09ada50>"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "polygon_297 = Polygon(zip(lng_list, lat_list))\n",
    "polygon_297"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [],
   "source": [
    "gdf_imovel_297 = gpd.GeoDataFrame(index=[0],crs='EPSG:31984', geometry=[polygon_297]).to_crs('epsg:4326')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.70824 -3.00480, -39.70732 -3.008...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry\n",
       "0  POLYGON ((-39.70824 -3.00480, -39.70732 -3.008..."
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_imovel_297"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "      <th>intersecao_cd_sipra</th>\n",
       "      <th>intersecao_nome</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...</td>\n",
       "      <td>CE0092000</td>\n",
       "      <td>PA PACHICU</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry intersecao_cd_sipra  \\\n",
       "0  POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...           CE0092000   \n",
       "\n",
       "  intersecao_nome  \n",
       "0      PA PACHICU  "
      ]
     },
     "execution_count": 31,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cruzamento3 = assentamentos.iloc[0].geometry.intersection(gdf_imovel_297.iloc[0].geometry)\n",
    "cruzamento3 = gpd.GeoDataFrame(index=[0],crs='EPSG:4326', geometry=[cruzamento3])\n",
    "cruzamento3['intersecao_cd_sipra'] = assentamentos.iloc[0].cd_sipra\n",
    "cruzamento3['intersecao_nome'] = assentamentos.iloc[0].nome_proje\n",
    "cruzamento3"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "      <th>intersecao_cd_sipra</th>\n",
       "      <th>intersecao_nome</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.71787 -3.00739, -39.71595 -3.003...</td>\n",
       "      <td>CE0097000</td>\n",
       "      <td>PA PATOS</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry intersecao_cd_sipra  \\\n",
       "0  POLYGON ((-39.71787 -3.00739, -39.71595 -3.003...           CE0097000   \n",
       "\n",
       "  intersecao_nome  \n",
       "0        PA PATOS  "
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cruzamento4 = assentamentos.iloc[1].geometry.intersection(gdf_imovel_297.iloc[0].geometry)\n",
    "cruzamento4 = gpd.GeoDataFrame(index=[0],crs='EPSG:4326', geometry=[cruzamento4])\n",
    "cruzamento4['intersecao_cd_sipra'] = assentamentos.iloc[1].cd_sipra\n",
    "cruzamento4['intersecao_nome'] = assentamentos.iloc[1].nome_proje\n",
    "cruzamento4"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "      <th>intersecao_cd_sipra</th>\n",
       "      <th>intersecao_nome</th>\n",
       "      <th>intersecao_area_km2</th>\n",
       "      <th>intersecao_imovel</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...</td>\n",
       "      <td>CE0092000</td>\n",
       "      <td>PA PACHICU</td>\n",
       "      <td>100.43</td>\n",
       "      <td>297</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.71787 -3.00739, -39.71595 -3.003...</td>\n",
       "      <td>CE0097000</td>\n",
       "      <td>PA PATOS</td>\n",
       "      <td>87.65</td>\n",
       "      <td>297</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry intersecao_cd_sipra  \\\n",
       "0  POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...           CE0092000   \n",
       "0  POLYGON ((-39.71787 -3.00739, -39.71595 -3.003...           CE0097000   \n",
       "\n",
       "  intersecao_nome  intersecao_area_km2 intersecao_imovel  \n",
       "0      PA PACHICU               100.43               297  \n",
       "0        PA PATOS                87.65               297  "
      ]
     },
     "execution_count": 33,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "resultado_intersecoes_297 = pd.concat([cruzamento3,cruzamento4])\n",
    "resultado_intersecoes_297['intersecao_area_km2'] = resultado_intersecoes_297.geometry.apply(lambda x: getArea(x)/10000)\n",
    "resultado_intersecoes_297['intersecao_imovel'] = '297'\n",
    "resultado_intersecoes_297"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "      <th>intersecao_cd_sipra</th>\n",
       "      <th>intersecao_nome</th>\n",
       "      <th>intersecao_area_km2</th>\n",
       "      <th>intersecao_imovel</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...</td>\n",
       "      <td>CE0092000</td>\n",
       "      <td>PA PACHICU</td>\n",
       "      <td>47.17</td>\n",
       "      <td>146</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...</td>\n",
       "      <td>CE0097000</td>\n",
       "      <td>PA PATOS</td>\n",
       "      <td>0.32</td>\n",
       "      <td>146</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...</td>\n",
       "      <td>CE0092000</td>\n",
       "      <td>PA PACHICU</td>\n",
       "      <td>100.43</td>\n",
       "      <td>297</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.71787 -3.00739, -39.71595 -3.003...</td>\n",
       "      <td>CE0097000</td>\n",
       "      <td>PA PATOS</td>\n",
       "      <td>87.65</td>\n",
       "      <td>297</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry intersecao_cd_sipra  \\\n",
       "0  POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...           CE0092000   \n",
       "0  POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...           CE0097000   \n",
       "0  POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...           CE0092000   \n",
       "0  POLYGON ((-39.71787 -3.00739, -39.71595 -3.003...           CE0097000   \n",
       "\n",
       "  intersecao_nome  intersecao_area_km2 intersecao_imovel  \n",
       "0      PA PACHICU                47.17               146  \n",
       "0        PA PATOS                 0.32               146  \n",
       "0      PA PACHICU               100.43               297  \n",
       "0        PA PATOS                87.65               297  "
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "intersecoes = pd.concat([intersecoes,resultado_intersecoes_297])\n",
    "intersecoes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [],
   "source": [
    "intersecao_imv146_297 = gdf_imovel_146.iloc[0].geometry.intersection(gdf_imovel_297.iloc[0].geometry)\n",
    "cruzamento5 = gpd.GeoDataFrame(index=[0],crs='EPSG:4326', geometry=[intersecao_imv146_297])\n",
    "cruzamento5['categoria'] = 'Área de interseção'\n",
    "cruzamento5['descricao'] = f'Área de interseção ({round(getArea(intersecao_imv146_297)/10000,2)}) entre os imóveis 146 e 297'\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>descricao</th>\n",
       "      <th>geometry</th>\n",
       "      <th>categoria</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Interseção do imóvel 146 com o assentamento PA...</td>\n",
       "      <td>POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...</td>\n",
       "      <td>Área de interseção</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Interseção do imóvel 146 com o assentamento PA...</td>\n",
       "      <td>POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...</td>\n",
       "      <td>Área de interseção</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Interseção do imóvel 297 com o assentamento PA...</td>\n",
       "      <td>POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...</td>\n",
       "      <td>Área de interseção</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Interseção do imóvel 297 com o assentamento PA...</td>\n",
       "      <td>POLYGON ((-39.71787 -3.00739, -39.71595 -3.003...</td>\n",
       "      <td>Área de interseção</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Área de interseção (113.25) entre os imóveis 1...</td>\n",
       "      <td>POLYGON ((-39.72933 -2.99624, -39.73320 -2.993...</td>\n",
       "      <td>Área de interseção</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                           descricao  \\\n",
       "0  Interseção do imóvel 146 com o assentamento PA...   \n",
       "0  Interseção do imóvel 146 com o assentamento PA...   \n",
       "0  Interseção do imóvel 297 com o assentamento PA...   \n",
       "0  Interseção do imóvel 297 com o assentamento PA...   \n",
       "0  Área de interseção (113.25) entre os imóveis 1...   \n",
       "\n",
       "                                            geometry           categoria  \n",
       "0  POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...  Área de interseção  \n",
       "0  POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...  Área de interseção  \n",
       "0  POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...  Área de interseção  \n",
       "0  POLYGON ((-39.71787 -3.00739, -39.71595 -3.003...  Área de interseção  \n",
       "0  POLYGON ((-39.72933 -2.99624, -39.73320 -2.993...  Área de interseção  "
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "area_int = intersecoes.copy()\n",
    "area_int['descricao'] = area_int.apply(lambda x: f'Interseção do imóvel {x[\"intersecao_imovel\"]} com o assentamento {x[\"intersecao_nome\"]}: {x[\"intersecao_area_km2\"]} hectares', axis=1)\n",
    "area_int = area_int[['descricao','geometry']]\n",
    "area_int['categoria'] = 'Área de interseção'\n",
    "\n",
    "area_int = pd.concat([area_int,cruzamento5])\n",
    "\n",
    "area_int.to_file('resultados/interseções.geojson', driver ='GeoJSON')\n",
    "\n",
    "area_int"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>descricao</th>\n",
       "      <th>geometry</th>\n",
       "      <th>categoria</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>185</th>\n",
       "      <td>PA PACHICU. Criado em 23/11/1995, 619.64 km² e...</td>\n",
       "      <td>MULTIPOLYGON (((-39.72711 -2.99860, -39.72305 ...</td>\n",
       "      <td>Assentamento</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>348</th>\n",
       "      <td>PA PATOS. Criado em 12/12/1995, 1598.64 km² e ...</td>\n",
       "      <td>MULTIPOLYGON (((-39.68994 -3.03521, -39.69042 ...</td>\n",
       "      <td>Assentamento</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                             descricao  \\\n",
       "185  PA PACHICU. Criado em 23/11/1995, 619.64 km² e...   \n",
       "348  PA PATOS. Criado em 12/12/1995, 1598.64 km² e ...   \n",
       "\n",
       "                                              geometry     categoria  \n",
       "185  MULTIPOLYGON (((-39.72711 -2.99860, -39.72305 ...  Assentamento  \n",
       "348  MULTIPOLYGON (((-39.68994 -3.03521, -39.69042 ...  Assentamento  "
      ]
     },
     "execution_count": 37,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "area_ass = assentamentos.copy()\n",
    "area_ass['descricao'] = area_ass.apply(lambda x: f'{x[\"nome_proje\"]}. Criado em {x[\"data_de_cr\"]}, {round(float(x[\"area_hecta\"]),2)} km² e {x[\"num_famili\"]} famílias', axis=1)\n",
    "area_ass = area_ass[['descricao','geometry']]\n",
    "area_ass['categoria'] = 'Assentamento'\n",
    "area_ass"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "      <th>descricao</th>\n",
       "      <th>categoria</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.73240 -2.99052, -39.73424 -2.993...</td>\n",
       "      <td>Inscrição Nº 146 - Segundo documento, área de ...</td>\n",
       "      <td>Imóvel</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry  \\\n",
       "0  POLYGON ((-39.73240 -2.99052, -39.73424 -2.993...   \n",
       "\n",
       "                                           descricao categoria  \n",
       "0  Inscrição Nº 146 - Segundo documento, área de ...    Imóvel  "
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "area_146 = gdf_imovel_146.copy()\n",
    "area_146['descricao'] = \"Inscrição Nº 146 - Segundo documento, área de 152 hectares foi negociada pela última vez em 2020 por R$ 160 mil. Terreno abrange a área dentro do manguezal onde foi aberta a estrada sem licenciamento.\"\n",
    "area_146['categoria']=\"Imóvel\"\n",
    "area_146"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>geometry</th>\n",
       "      <th>descricao</th>\n",
       "      <th>categoria</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-39.70824 -3.00480, -39.70732 -3.008...</td>\n",
       "      <td>Inscrição Nº 297 - Área de 126,95 hectares ter...</td>\n",
       "      <td>Imóvel</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry  \\\n",
       "0  POLYGON ((-39.70824 -3.00480, -39.70732 -3.008...   \n",
       "\n",
       "                                           descricao categoria  \n",
       "0  Inscrição Nº 297 - Área de 126,95 hectares ter...    Imóvel  "
      ]
     },
     "execution_count": 39,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "area_297 = gdf_imovel_297.copy()\n",
    "area_297['descricao'] = \"Inscrição Nº 297 - Área de 126,95 hectares teria sido negociada pela última vez em 2010. Avaliada à época em R$ 50 mil. Segundo documento, pertence a um empresário residente em Fortaleza que já foi investigado e responde a processo por fraude em financiamentos - valores de R$ 21 milhões - obtidos junto ao Banco do Nordeste do Brasil.\"\n",
    "area_297['categoria']=\"Imóvel\"\n",
    "area_297"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/alexandre/anaconda3/envs/credos/lib/python3.7/site-packages/pandas/core/dtypes/cast.py:112: ShapelyDeprecationWarning: The array interface is deprecated and will no longer work in Shapely 2.0. Convert the '.coords' to a numpy array instead.\n",
      "  values = construct_1d_object_array_from_listlike(values)\n"
     ]
    }
   ],
   "source": [
    "estradas = pd.read_csv('estradas.csv')\n",
    "\n",
    "estradas['inicio'] = [Point(xy) for xy in zip(estradas.inicio_lng, estradas.inicio_lat)]\n",
    "estradas['fim'] = [Point(xy) for xy in zip(estradas.fim_lng, estradas.fim_lat)]\n",
    "estradas['geometry'] = estradas.apply(lambda row: LineString([row['inicio'], row['fim']]), axis=1)\n",
    "\n",
    "estradas['descricao'] = estradas['nome'] + ': '+ estradas['descricao']\n",
    "estradas['categoria'] = 'Estradas'\n",
    "\n",
    "estradas = estradas[['descricao','categoria','geometry']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/alexandre/anaconda3/envs/credos/lib/python3.7/site-packages/geopandas/geodataframe.py:34: ShapelyDeprecationWarning: The array interface is deprecated and will no longer work in Shapely 2.0. Convert the '.coords' to a numpy array instead.\n",
      "  out = from_shapely(data)\n",
      "/home/alexandre/anaconda3/envs/credos/lib/python3.7/site-packages/geopandas/geoseries.py:128: ShapelyDeprecationWarning: The array interface is deprecated and will no longer work in Shapely 2.0. Convert the '.coords' to a numpy array instead.\n",
      "  data = from_shapely(s.values)\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>descricao</th>\n",
       "      <th>categoria</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>PISTA DO MANGUE: O local de abertura de estrad...</td>\n",
       "      <td>Estradas</td>\n",
       "      <td>LINESTRING (-39.73129 -3.00397, -39.71863 -3.0...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>PISTA ÁREA DA FAIXA DE PRAIA: O outro local da...</td>\n",
       "      <td>Estradas</td>\n",
       "      <td>LINESTRING (-39.73497 -2.99365, -39.73386 -2.9...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                           descricao categoria  \\\n",
       "0  PISTA DO MANGUE: O local de abertura de estrad...  Estradas   \n",
       "1  PISTA ÁREA DA FAIXA DE PRAIA: O outro local da...  Estradas   \n",
       "\n",
       "                                            geometry  \n",
       "0  LINESTRING (-39.73129 -3.00397, -39.71863 -3.0...  \n",
       "1  LINESTRING (-39.73497 -2.99365, -39.73386 -2.9...  "
      ]
     },
     "execution_count": 41,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_estradas = gpd.GeoDataFrame(estradas,crs='EPSG:31984').to_crs('epsg:4326')\n",
    "gdf_estradas"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>index</th>\n",
       "      <th>descricao</th>\n",
       "      <th>geometry</th>\n",
       "      <th>categoria</th>\n",
       "      <th>lng</th>\n",
       "      <th>lat</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>PA PACHICU. Criado em 23/11/1995, 619.64 km² e...</td>\n",
       "      <td>MULTIPOLYGON (((-39.72711 -2.99860, -39.72305 ...</td>\n",
       "      <td>Assentamento</td>\n",
       "      <td>-39.730302</td>\n",
       "      <td>-3.018064</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>PA PATOS. Criado em 12/12/1995, 1598.64 km² e ...</td>\n",
       "      <td>MULTIPOLYGON (((-39.68994 -3.03521, -39.69042 ...</td>\n",
       "      <td>Assentamento</td>\n",
       "      <td>-39.713326</td>\n",
       "      <td>-3.033359</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>Inscrição Nº 146 - Segundo documento, área de ...</td>\n",
       "      <td>POLYGON ((-39.73240 -2.99052, -39.73424 -2.993...</td>\n",
       "      <td>Imóvel</td>\n",
       "      <td>-39.723925</td>\n",
       "      <td>-2.999156</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>Inscrição Nº 297 - Área de 126,95 hectares ter...</td>\n",
       "      <td>POLYGON ((-39.70824 -3.00480, -39.70732 -3.008...</td>\n",
       "      <td>Imóvel</td>\n",
       "      <td>-39.720576</td>\n",
       "      <td>-3.003455</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>Interseção do imóvel 146 com o assentamento PA...</td>\n",
       "      <td>POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...</td>\n",
       "      <td>Área de interseção</td>\n",
       "      <td>-39.722889</td>\n",
       "      <td>-3.002264</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>5</td>\n",
       "      <td>Interseção do imóvel 146 com o assentamento PA...</td>\n",
       "      <td>POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...</td>\n",
       "      <td>Área de interseção</td>\n",
       "      <td>-39.711285</td>\n",
       "      <td>-3.003834</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>6</td>\n",
       "      <td>Interseção do imóvel 297 com o assentamento PA...</td>\n",
       "      <td>POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...</td>\n",
       "      <td>Área de interseção</td>\n",
       "      <td>-39.723087</td>\n",
       "      <td>-3.003863</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>7</td>\n",
       "      <td>Interseção do imóvel 297 com o assentamento PA...</td>\n",
       "      <td>POLYGON ((-39.71787 -3.00739, -39.71595 -3.003...</td>\n",
       "      <td>Área de interseção</td>\n",
       "      <td>-39.712650</td>\n",
       "      <td>-3.007941</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>8</td>\n",
       "      <td>Área de interseção (113.25) entre os imóveis 1...</td>\n",
       "      <td>POLYGON ((-39.72933 -2.99624, -39.73320 -2.993...</td>\n",
       "      <td>Área de interseção</td>\n",
       "      <td>-39.723512</td>\n",
       "      <td>-3.000276</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>9</td>\n",
       "      <td>PISTA DO MANGUE: O local de abertura de estrad...</td>\n",
       "      <td>LINESTRING (-39.73129 -3.00397, -39.71863 -3.0...</td>\n",
       "      <td>Estradas</td>\n",
       "      <td>-39.724962</td>\n",
       "      <td>-3.006519</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>10</td>\n",
       "      <td>PISTA ÁREA DA FAIXA DE PRAIA: O outro local da...</td>\n",
       "      <td>LINESTRING (-39.73497 -2.99365, -39.73386 -2.9...</td>\n",
       "      <td>Estradas</td>\n",
       "      <td>-39.734415</td>\n",
       "      <td>-2.991945</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    index                                          descricao  \\\n",
       "0       0  PA PACHICU. Criado em 23/11/1995, 619.64 km² e...   \n",
       "1       1  PA PATOS. Criado em 12/12/1995, 1598.64 km² e ...   \n",
       "2       2  Inscrição Nº 146 - Segundo documento, área de ...   \n",
       "3       3  Inscrição Nº 297 - Área de 126,95 hectares ter...   \n",
       "4       4  Interseção do imóvel 146 com o assentamento PA...   \n",
       "5       5  Interseção do imóvel 146 com o assentamento PA...   \n",
       "6       6  Interseção do imóvel 297 com o assentamento PA...   \n",
       "7       7  Interseção do imóvel 297 com o assentamento PA...   \n",
       "8       8  Área de interseção (113.25) entre os imóveis 1...   \n",
       "9       9  PISTA DO MANGUE: O local de abertura de estrad...   \n",
       "10     10  PISTA ÁREA DA FAIXA DE PRAIA: O outro local da...   \n",
       "\n",
       "                                             geometry           categoria  \\\n",
       "0   MULTIPOLYGON (((-39.72711 -2.99860, -39.72305 ...        Assentamento   \n",
       "1   MULTIPOLYGON (((-39.68994 -3.03521, -39.69042 ...        Assentamento   \n",
       "2   POLYGON ((-39.73240 -2.99052, -39.73424 -2.993...              Imóvel   \n",
       "3   POLYGON ((-39.70824 -3.00480, -39.70732 -3.008...              Imóvel   \n",
       "4   POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...  Área de interseção   \n",
       "5   POLYGON ((-39.71106 -3.00364, -39.71055 -3.003...  Área de interseção   \n",
       "6   POLYGON ((-39.72305 -3.00010, -39.71838 -3.002...  Área de interseção   \n",
       "7   POLYGON ((-39.71787 -3.00739, -39.71595 -3.003...  Área de interseção   \n",
       "8   POLYGON ((-39.72933 -2.99624, -39.73320 -2.993...  Área de interseção   \n",
       "9   LINESTRING (-39.73129 -3.00397, -39.71863 -3.0...            Estradas   \n",
       "10  LINESTRING (-39.73497 -2.99365, -39.73386 -2.9...            Estradas   \n",
       "\n",
       "           lng        lat  \n",
       "0   -39.730302  -3.018064  \n",
       "1   -39.713326  -3.033359  \n",
       "2   -39.723925  -2.999156  \n",
       "3   -39.720576  -3.003455  \n",
       "4   -39.722889  -3.002264  \n",
       "5   -39.711285  -3.003834  \n",
       "6   -39.723087  -3.003863  \n",
       "7   -39.712650  -3.007941  \n",
       "8   -39.723512  -3.000276  \n",
       "9   -39.724962  -3.006519  \n",
       "10  -39.734415  -2.991945  "
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_completo = pd.concat([area_ass, area_146, area_297, area_int,gdf_estradas])\n",
    "\n",
    "\n",
    "gdf_completo['lng'] = gdf_completo.geometry.apply(lambda p: \"%.6f\" %  p.centroid.x)\n",
    "gdf_completo['lat'] = gdf_completo.geometry.apply(lambda p: \"%.6f\" %  p.centroid.y)\n",
    "\n",
    "gdf_completo.insert(0, 'index', range(0, len(gdf_completo)))\n",
    "gdf_completo = gdf_completo.reset_index( drop=True)\n",
    "gdf_completo"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [],
   "source": [
    "gdf_completo.to_file('resultados/dataset_completo.geojson', driver ='GeoJSON')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>cd_sipra</th>\n",
       "      <th>uf</th>\n",
       "      <th>nome_proje</th>\n",
       "      <th>municipio</th>\n",
       "      <th>area_hecta</th>\n",
       "      <th>capacidade</th>\n",
       "      <th>num_famili</th>\n",
       "      <th>fase</th>\n",
       "      <th>data_de_cr</th>\n",
       "      <th>forma_obte</th>\n",
       "      <th>data_obten</th>\n",
       "      <th>area_calc_</th>\n",
       "      <th>sr</th>\n",
       "      <th>descricao_</th>\n",
       "      <th>geometry</th>\n",
       "      <th>lng</th>\n",
       "      <th>lat</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>34</th>\n",
       "      <td>CE0351000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA CÓRREGO NOVO</td>\n",
       "      <td>ITAREMA</td>\n",
       "      <td>1500.7780</td>\n",
       "      <td>25</td>\n",
       "      <td>22</td>\n",
       "      <td>5</td>\n",
       "      <td>28/08/2007</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>12/12/2005</td>\n",
       "      <td>1534.90</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.94784 -3.01505, -39.94188 ...</td>\n",
       "      <td>-39.969719</td>\n",
       "      <td>-3.021816</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>49</th>\n",
       "      <td>CE0048000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA CAPIM AÇU</td>\n",
       "      <td>ITAREMA</td>\n",
       "      <td>1452.5905</td>\n",
       "      <td>46</td>\n",
       "      <td>52</td>\n",
       "      <td>6</td>\n",
       "      <td>20/10/1989</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>15/08/1988</td>\n",
       "      <td>1451.29</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.97704 -3.09274, -39.96264 ...</td>\n",
       "      <td>-39.962977</td>\n",
       "      <td>-3.111622</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59</th>\n",
       "      <td>CE0088000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA CÓRREGO DAS AROEIRAS</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>1920.1004</td>\n",
       "      <td>64</td>\n",
       "      <td>46</td>\n",
       "      <td>6</td>\n",
       "      <td>28/09/1995</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>12/01/1995</td>\n",
       "      <td>1921.23</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.84980 -3.22224, -39.83318 ...</td>\n",
       "      <td>-39.835365</td>\n",
       "      <td>-3.246047</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>68</th>\n",
       "      <td>CE0090000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA LAGOA DO CACHIMBO</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>1085.7101</td>\n",
       "      <td>35</td>\n",
       "      <td>35</td>\n",
       "      <td>6</td>\n",
       "      <td>28/09/1995</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>12/01/1995</td>\n",
       "      <td>1068.53</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.81329 -3.13211, -39.83820 ...</td>\n",
       "      <td>-39.818697</td>\n",
       "      <td>-3.111558</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>71</th>\n",
       "      <td>CE0025000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA SABIAGUABA</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>718.9365</td>\n",
       "      <td>28</td>\n",
       "      <td>28</td>\n",
       "      <td>6</td>\n",
       "      <td>22/06/1987</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>16/02/1987</td>\n",
       "      <td>718.88</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.55610 -3.09330, -39.55486 ...</td>\n",
       "      <td>-39.564236</td>\n",
       "      <td>-3.098053</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>80</th>\n",
       "      <td>CE0196000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA LAGOA DOS NEGROS/VOLTA</td>\n",
       "      <td>ITAREMA</td>\n",
       "      <td>3193.0697</td>\n",
       "      <td>92</td>\n",
       "      <td>81</td>\n",
       "      <td>5</td>\n",
       "      <td>22/12/1997</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>04/08/1997</td>\n",
       "      <td>3179.36</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.94180 -3.03902, -39.93427 ...</td>\n",
       "      <td>-39.963710</td>\n",
       "      <td>-3.060799</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>81</th>\n",
       "      <td>CE0026000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA LAGOA DO MINEIRO</td>\n",
       "      <td>ITAREMA</td>\n",
       "      <td>5763.8271</td>\n",
       "      <td>135</td>\n",
       "      <td>128</td>\n",
       "      <td>6</td>\n",
       "      <td>13/07/1987</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>26/06/1986</td>\n",
       "      <td>5767.34</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.78835 -3.01933, -39.78673 ...</td>\n",
       "      <td>-39.796962</td>\n",
       "      <td>-3.050461</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>102</th>\n",
       "      <td>CE0091000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA MUFUMBO</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>1305.2419</td>\n",
       "      <td>42</td>\n",
       "      <td>17</td>\n",
       "      <td>6</td>\n",
       "      <td>23/11/1995</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>24/03/1995</td>\n",
       "      <td>1312.74</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.96476 -3.40401, -39.95669 ...</td>\n",
       "      <td>-39.982201</td>\n",
       "      <td>-3.416894</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>148</th>\n",
       "      <td>CE0071000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA LAGOA DO JARDIM</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>1416.7753</td>\n",
       "      <td>55</td>\n",
       "      <td>48</td>\n",
       "      <td>6</td>\n",
       "      <td>31/05/1995</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>19/08/1993</td>\n",
       "      <td>1431.10</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.81329 -3.13211, -39.81288 ...</td>\n",
       "      <td>-39.790212</td>\n",
       "      <td>-3.118689</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>165</th>\n",
       "      <td>CE0239000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA MARRECAS</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>4391.0230</td>\n",
       "      <td>93</td>\n",
       "      <td>78</td>\n",
       "      <td>6</td>\n",
       "      <td>18/12/1998</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>02/09/1998</td>\n",
       "      <td>4393.40</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.95001 -3.37377, -39.94704 ...</td>\n",
       "      <td>-39.924655</td>\n",
       "      <td>-3.405632</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>176</th>\n",
       "      <td>CE0159000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA CANAÃ/MELANCIAS</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>857.8950</td>\n",
       "      <td>30</td>\n",
       "      <td>24</td>\n",
       "      <td>6</td>\n",
       "      <td>26/12/1996</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>24/07/1996</td>\n",
       "      <td>858.42</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.75926 -3.07685, -39.75959 ...</td>\n",
       "      <td>-39.774410</td>\n",
       "      <td>-3.087919</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>252</th>\n",
       "      <td>CE0211000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA MACACO II</td>\n",
       "      <td>ITAREMA</td>\n",
       "      <td>562.2214</td>\n",
       "      <td>18</td>\n",
       "      <td>14</td>\n",
       "      <td>6</td>\n",
       "      <td>30/12/1997</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>13/11/1997</td>\n",
       "      <td>532.33</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.85614 -3.05346, -39.83233 ...</td>\n",
       "      <td>-39.840864</td>\n",
       "      <td>-3.066974</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>261</th>\n",
       "      <td>CE0058000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA SALGADO COMPRIDO</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>1522.4570</td>\n",
       "      <td>60</td>\n",
       "      <td>56</td>\n",
       "      <td>6</td>\n",
       "      <td>31/10/1991</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>13/06/1989</td>\n",
       "      <td>1518.02</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.72569 -3.05441, -39.72572 ...</td>\n",
       "      <td>-39.740110</td>\n",
       "      <td>-3.054477</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>307</th>\n",
       "      <td>CE0203000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA LAGOA DOS VEADOS</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>844.8697</td>\n",
       "      <td>25</td>\n",
       "      <td>22</td>\n",
       "      <td>6</td>\n",
       "      <td>22/12/1997</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>09/07/1997</td>\n",
       "      <td>849.58</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.81429 -3.15654, -39.81479 ...</td>\n",
       "      <td>-39.829888</td>\n",
       "      <td>-3.167009</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>369</th>\n",
       "      <td>CE0202000</td>\n",
       "      <td>CE</td>\n",
       "      <td>PA VEDÓIA</td>\n",
       "      <td>AMONTADA</td>\n",
       "      <td>458.0319</td>\n",
       "      <td>29</td>\n",
       "      <td>28</td>\n",
       "      <td>6</td>\n",
       "      <td>22/12/1997</td>\n",
       "      <td>Desapropriação</td>\n",
       "      <td>16/05/1997</td>\n",
       "      <td>458.33</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>MULTIPOLYGON (((-39.72717 -3.10850, -39.72216 ...</td>\n",
       "      <td>-39.707718</td>\n",
       "      <td>-3.105968</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      cd_sipra  uf                 nome_proje municipio area_hecta  \\\n",
       "34   CE0351000  CE            PA CÓRREGO NOVO   ITAREMA  1500.7780   \n",
       "49   CE0048000  CE               PA CAPIM AÇU   ITAREMA  1452.5905   \n",
       "59   CE0088000  CE    PA CÓRREGO DAS AROEIRAS  AMONTADA  1920.1004   \n",
       "68   CE0090000  CE       PA LAGOA DO CACHIMBO  AMONTADA  1085.7101   \n",
       "71   CE0025000  CE              PA SABIAGUABA  AMONTADA   718.9365   \n",
       "80   CE0196000  CE  PA LAGOA DOS NEGROS/VOLTA   ITAREMA  3193.0697   \n",
       "81   CE0026000  CE        PA LAGOA DO MINEIRO   ITAREMA  5763.8271   \n",
       "102  CE0091000  CE                 PA MUFUMBO  AMONTADA  1305.2419   \n",
       "148  CE0071000  CE         PA LAGOA DO JARDIM  AMONTADA  1416.7753   \n",
       "165  CE0239000  CE                PA MARRECAS  AMONTADA  4391.0230   \n",
       "176  CE0159000  CE         PA CANAÃ/MELANCIAS  AMONTADA   857.8950   \n",
       "252  CE0211000  CE               PA MACACO II   ITAREMA   562.2214   \n",
       "261  CE0058000  CE        PA SALGADO COMPRIDO  AMONTADA  1522.4570   \n",
       "307  CE0203000  CE        PA LAGOA DOS VEADOS  AMONTADA   844.8697   \n",
       "369  CE0202000  CE                  PA VEDÓIA  AMONTADA   458.0319   \n",
       "\n",
       "     capacidade  num_famili  fase  data_de_cr      forma_obte  data_obten  \\\n",
       "34           25          22     5  28/08/2007  Desapropriação  12/12/2005   \n",
       "49           46          52     6  20/10/1989  Desapropriação  15/08/1988   \n",
       "59           64          46     6  28/09/1995  Desapropriação  12/01/1995   \n",
       "68           35          35     6  28/09/1995  Desapropriação  12/01/1995   \n",
       "71           28          28     6  22/06/1987  Desapropriação  16/02/1987   \n",
       "80           92          81     5  22/12/1997  Desapropriação  04/08/1997   \n",
       "81          135         128     6  13/07/1987  Desapropriação  26/06/1986   \n",
       "102          42          17     6  23/11/1995  Desapropriação  24/03/1995   \n",
       "148          55          48     6  31/05/1995  Desapropriação  19/08/1993   \n",
       "165          93          78     6  18/12/1998  Desapropriação  02/09/1998   \n",
       "176          30          24     6  26/12/1996  Desapropriação  24/07/1996   \n",
       "252          18          14     6  30/12/1997  Desapropriação  13/11/1997   \n",
       "261          60          56     6  31/10/1991  Desapropriação  13/06/1989   \n",
       "307          25          22     6  22/12/1997  Desapropriação  09/07/1997   \n",
       "369          29          28     6  22/12/1997  Desapropriação  16/05/1997   \n",
       "\n",
       "     area_calc_    sr descricao_  \\\n",
       "34      1534.90  None       None   \n",
       "49      1451.29  None       None   \n",
       "59      1921.23  None       None   \n",
       "68      1068.53  None       None   \n",
       "71       718.88  None       None   \n",
       "80      3179.36  None       None   \n",
       "81      5767.34  None       None   \n",
       "102     1312.74  None       None   \n",
       "148     1431.10  None       None   \n",
       "165     4393.40  None       None   \n",
       "176      858.42  None       None   \n",
       "252      532.33  None       None   \n",
       "261     1518.02  None       None   \n",
       "307      849.58  None       None   \n",
       "369      458.33  None       None   \n",
       "\n",
       "                                              geometry         lng        lat  \n",
       "34   MULTIPOLYGON (((-39.94784 -3.01505, -39.94188 ...  -39.969719  -3.021816  \n",
       "49   MULTIPOLYGON (((-39.97704 -3.09274, -39.96264 ...  -39.962977  -3.111622  \n",
       "59   MULTIPOLYGON (((-39.84980 -3.22224, -39.83318 ...  -39.835365  -3.246047  \n",
       "68   MULTIPOLYGON (((-39.81329 -3.13211, -39.83820 ...  -39.818697  -3.111558  \n",
       "71   MULTIPOLYGON (((-39.55610 -3.09330, -39.55486 ...  -39.564236  -3.098053  \n",
       "80   MULTIPOLYGON (((-39.94180 -3.03902, -39.93427 ...  -39.963710  -3.060799  \n",
       "81   MULTIPOLYGON (((-39.78835 -3.01933, -39.78673 ...  -39.796962  -3.050461  \n",
       "102  MULTIPOLYGON (((-39.96476 -3.40401, -39.95669 ...  -39.982201  -3.416894  \n",
       "148  MULTIPOLYGON (((-39.81329 -3.13211, -39.81288 ...  -39.790212  -3.118689  \n",
       "165  MULTIPOLYGON (((-39.95001 -3.37377, -39.94704 ...  -39.924655  -3.405632  \n",
       "176  MULTIPOLYGON (((-39.75926 -3.07685, -39.75959 ...  -39.774410  -3.087919  \n",
       "252  MULTIPOLYGON (((-39.85614 -3.05346, -39.83233 ...  -39.840864  -3.066974  \n",
       "261  MULTIPOLYGON (((-39.72569 -3.05441, -39.72572 ...  -39.740110  -3.054477  \n",
       "307  MULTIPOLYGON (((-39.81429 -3.15654, -39.81479 ...  -39.829888  -3.167009  \n",
       "369  MULTIPOLYGON (((-39.72717 -3.10850, -39.72216 ...  -39.707718  -3.105968  "
      ]
     },
     "execution_count": 44,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Selecionando os outros assentamentos na região\n",
    "\n",
    "assentamentos_federais_ce_itarema_amontada = assentamentos_federais_ce.loc[(assentamentos_federais_ce.municipio == 'ITAREMA')|(assentamentos_federais_ce.municipio == 'AMONTADA')]\n",
    "\n",
    "assentamentos_federais_ce_itarema_amontada = assentamentos_federais_ce_itarema_amontada.loc[~(\n",
    "    (assentamentos_federais_ce.nome_proje == 'PA PATOS') | (assentamentos_federais_ce.nome_proje == 'PA PACHICU')\n",
    "    )]\n",
    "\n",
    "assentamentos_federais_ce_itarema_amontada['lng'] = assentamentos_federais_ce_itarema_amontada.geometry.apply(lambda p: \"%.6f\" %  p.centroid.x)\n",
    "assentamentos_federais_ce_itarema_amontada['lat'] = assentamentos_federais_ce_itarema_amontada.geometry.apply(lambda p: \"%.6f\" %  p.centroid.y)\n",
    "\n",
    "assentamentos_federais_ce_itarema_amontada"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [],
   "source": [
    "assentamentos_federais_ce_itarema_amontada.to_file('resultados/outros_assentamentos.geojson', driver ='GeoJSON')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.17"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}