01_baixarDados_EMFP_ANA.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "53565c72-d29a-490e-9e43-451b51f7ad88",
   "metadata": {},
   "source": [
    "# Notebook para baixar os dados de todas as estações fluviométricas e pluviométricas da rede\n",
    "\n",
    "Os dados estão armazenados um GeoPackage e estão acessíveis a partir do link:\n",
    "\n",
    "- https://metadados.snirh.gov.br/geonetwork/srv/por/catalog.search#/metadata/f85dbf06-a869-414c-afc5-bb01869e9156"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "0f5a4240-abc5-4f55-bb9c-e5c582b1ffaa",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import geopandas as gpd\n",
    "import shapely\n",
    "from pyogrio import read_dataframe"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "5c82dc0c-62e5-4987-b6f0-3469e7bef419",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "url_dataset = \"https://metadados.snirh.gov.br/files/f85dbf06-a869-414c-afc5-bb01869e9156/Estacoes_Fluviometricas_e_Pluviometricas_da_Rede.gpkg\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "4f941fcf-9e2a-4fc7-9493-5c240dfc6b3a",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "dataset_local = \"./dados/ana/Estacoes_Fluviometricas_e_Pluviometricas_da_Rede.gpkg\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "56d80672-bbe9-477f-84ce-389cf0cec76f",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "ename": "DataSourceError",
     "evalue": "'/vsicurl/https://metadados.snirh.gov.br/files/f85dbf06-a869-414c-afc5-bb01869e9156/Estacoes_Fluviometricas_e_Pluviometricas_da_Rede.gpkg' does not exist in the file system, and is not recognized as a supported dataset name.",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mCPLE_OpenFailedError\u001b[0m                      Traceback (most recent call last)",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\_io.pyx:154\u001b[0m, in \u001b[0;36mpyogrio._io.ogr_open\u001b[1;34m()\u001b[0m\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\_err.pyx:179\u001b[0m, in \u001b[0;36mpyogrio._err.exc_wrap_pointer\u001b[1;34m()\u001b[0m\n",
      "\u001b[1;31mCPLE_OpenFailedError\u001b[0m: '/vsicurl/https://metadados.snirh.gov.br/files/f85dbf06-a869-414c-afc5-bb01869e9156/Estacoes_Fluviometricas_e_Pluviometricas_da_Rede.gpkg' does not exist in the file system, and is not recognized as a supported dataset name.",
      "\nDuring handling of the above exception, another exception occurred:\n",
      "\u001b[1;31mDataSourceError\u001b[0m                           Traceback (most recent call last)",
      "Cell \u001b[1;32mIn[13], line 1\u001b[0m\n\u001b[1;32m----> 1\u001b[0m \u001b[43mread_dataframe\u001b[49m\u001b[43m(\u001b[49m\u001b[43murl_dataset\u001b[49m\u001b[43m)\u001b[49m\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\geopandas.py:149\u001b[0m, in \u001b[0;36mread_dataframe\u001b[1;34m(path_or_buffer, layer, encoding, columns, read_geometry, force_2d, skip_features, max_features, where, bbox, fids, sql, sql_dialect, fid_as_index, use_arrow, **kwargs)\u001b[0m\n\u001b[0;32m    146\u001b[0m path_or_buffer \u001b[38;5;241m=\u001b[39m _stringify_path(path_or_buffer)\n\u001b[0;32m    148\u001b[0m read_func \u001b[38;5;241m=\u001b[39m read_arrow \u001b[38;5;28;01mif\u001b[39;00m use_arrow \u001b[38;5;28;01melse\u001b[39;00m read\n\u001b[1;32m--> 149\u001b[0m result \u001b[38;5;241m=\u001b[39m \u001b[43mread_func\u001b[49m\u001b[43m(\u001b[49m\n\u001b[0;32m    150\u001b[0m \u001b[43m    \u001b[49m\u001b[43mpath_or_buffer\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    151\u001b[0m \u001b[43m    \u001b[49m\u001b[43mlayer\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mlayer\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    152\u001b[0m \u001b[43m    \u001b[49m\u001b[43mencoding\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mencoding\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    153\u001b[0m \u001b[43m    \u001b[49m\u001b[43mcolumns\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mcolumns\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    154\u001b[0m \u001b[43m    \u001b[49m\u001b[43mread_geometry\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mread_geometry\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    155\u001b[0m \u001b[43m    \u001b[49m\u001b[43mforce_2d\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mforce_2d\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    156\u001b[0m \u001b[43m    \u001b[49m\u001b[43mskip_features\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mskip_features\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    157\u001b[0m \u001b[43m    \u001b[49m\u001b[43mmax_features\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mmax_features\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    158\u001b[0m \u001b[43m    \u001b[49m\u001b[43mwhere\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mwhere\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    159\u001b[0m \u001b[43m    \u001b[49m\u001b[43mbbox\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mbbox\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    160\u001b[0m \u001b[43m    \u001b[49m\u001b[43mfids\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mfids\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    161\u001b[0m \u001b[43m    \u001b[49m\u001b[43msql\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43msql\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    162\u001b[0m \u001b[43m    \u001b[49m\u001b[43msql_dialect\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43msql_dialect\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    163\u001b[0m \u001b[43m    \u001b[49m\u001b[43mreturn_fids\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mfid_as_index\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    164\u001b[0m \u001b[43m    \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mkwargs\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    165\u001b[0m \u001b[43m\u001b[49m\u001b[43m)\u001b[49m\n\u001b[0;32m    167\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m use_arrow:\n\u001b[0;32m    168\u001b[0m     meta, table \u001b[38;5;241m=\u001b[39m result\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\raw.py:137\u001b[0m, in \u001b[0;36mread\u001b[1;34m(path_or_buffer, layer, encoding, columns, read_geometry, force_2d, skip_features, max_features, where, bbox, fids, sql, sql_dialect, return_fids, **kwargs)\u001b[0m\n\u001b[0;32m    134\u001b[0m dataset_kwargs \u001b[38;5;241m=\u001b[39m _preprocess_options_key_value(kwargs) \u001b[38;5;28;01mif\u001b[39;00m kwargs \u001b[38;5;28;01melse\u001b[39;00m {}\n\u001b[0;32m    136\u001b[0m \u001b[38;5;28;01mtry\u001b[39;00m:\n\u001b[1;32m--> 137\u001b[0m     result \u001b[38;5;241m=\u001b[39m \u001b[43mogr_read\u001b[49m\u001b[43m(\u001b[49m\n\u001b[0;32m    138\u001b[0m \u001b[43m        \u001b[49m\u001b[43mpath\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    139\u001b[0m \u001b[43m        \u001b[49m\u001b[43mlayer\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mlayer\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    140\u001b[0m \u001b[43m        \u001b[49m\u001b[43mencoding\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mencoding\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    141\u001b[0m \u001b[43m        \u001b[49m\u001b[43mcolumns\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mcolumns\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    142\u001b[0m \u001b[43m        \u001b[49m\u001b[43mread_geometry\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mread_geometry\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    143\u001b[0m \u001b[43m        \u001b[49m\u001b[43mforce_2d\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mforce_2d\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    144\u001b[0m \u001b[43m        \u001b[49m\u001b[43mskip_features\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mskip_features\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    145\u001b[0m \u001b[43m        \u001b[49m\u001b[43mmax_features\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mmax_features\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;129;43;01mor\u001b[39;49;00m\u001b[43m \u001b[49m\u001b[38;5;241;43m0\u001b[39;49m\u001b[43m,\u001b[49m\n\u001b[0;32m    146\u001b[0m \u001b[43m        \u001b[49m\u001b[43mwhere\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mwhere\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    147\u001b[0m \u001b[43m        \u001b[49m\u001b[43mbbox\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mbbox\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    148\u001b[0m \u001b[43m        \u001b[49m\u001b[43mfids\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mfids\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    149\u001b[0m \u001b[43m        \u001b[49m\u001b[43msql\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43msql\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    150\u001b[0m \u001b[43m        \u001b[49m\u001b[43msql_dialect\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43msql_dialect\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    151\u001b[0m \u001b[43m        \u001b[49m\u001b[43mreturn_fids\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mreturn_fids\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    152\u001b[0m \u001b[43m        \u001b[49m\u001b[43mdataset_kwargs\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mdataset_kwargs\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    153\u001b[0m \u001b[43m    \u001b[49m\u001b[43m)\u001b[49m\n\u001b[0;32m    154\u001b[0m \u001b[38;5;28;01mfinally\u001b[39;00m:\n\u001b[0;32m    155\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m buffer \u001b[38;5;129;01mis\u001b[39;00m \u001b[38;5;129;01mnot\u001b[39;00m \u001b[38;5;28;01mNone\u001b[39;00m:\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\_io.pyx:971\u001b[0m, in \u001b[0;36mpyogrio._io.ogr_read\u001b[1;34m()\u001b[0m\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\_io.pyx:164\u001b[0m, in \u001b[0;36mpyogrio._io.ogr_open\u001b[1;34m()\u001b[0m\n",
      "\u001b[1;31mDataSourceError\u001b[0m: '/vsicurl/https://metadados.snirh.gov.br/files/f85dbf06-a869-414c-afc5-bb01869e9156/Estacoes_Fluviometricas_e_Pluviometricas_da_Rede.gpkg' does not exist in the file system, and is not recognized as a supported dataset name."
     ]
    }
   ],
   "source": [
    "read_dataframe(url_dataset)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "9bb60a94-fb03-4f61-80ca-c7f867622e5b",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "ename": "DataSourceError",
     "evalue": "'/vsicurl/https://metadados.snirh.gov.br/files/f85dbf06-a869-414c-afc5-bb01869e9156/Estacoes_Fluviometricas_e_Pluviometricas_da_Rede.gpkg' does not exist in the file system, and is not recognized as a supported dataset name.",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mCPLE_OpenFailedError\u001b[0m                      Traceback (most recent call last)",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\_io.pyx:154\u001b[0m, in \u001b[0;36mpyogrio._io.ogr_open\u001b[1;34m()\u001b[0m\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\_err.pyx:179\u001b[0m, in \u001b[0;36mpyogrio._err.exc_wrap_pointer\u001b[1;34m()\u001b[0m\n",
      "\u001b[1;31mCPLE_OpenFailedError\u001b[0m: '/vsicurl/https://metadados.snirh.gov.br/files/f85dbf06-a869-414c-afc5-bb01869e9156/Estacoes_Fluviometricas_e_Pluviometricas_da_Rede.gpkg' does not exist in the file system, and is not recognized as a supported dataset name.",
      "\nDuring handling of the above exception, another exception occurred:\n",
      "\u001b[1;31mDataSourceError\u001b[0m                           Traceback (most recent call last)",
      "Cell \u001b[1;32mIn[8], line 1\u001b[0m\n\u001b[1;32m----> 1\u001b[0m gpkg_gdf \u001b[38;5;241m=\u001b[39m \u001b[43mgpd\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mread_file\u001b[49m\u001b[43m(\u001b[49m\u001b[43murl_dataset\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mlayer\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[38;5;124;43mEstacoes_Fluviometricas_e_Pluviometricas_da_Rede\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mengine\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[38;5;124;43mpyogrio\u001b[39;49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[43m)\u001b[49m\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\geopandas\\io\\file.py:271\u001b[0m, in \u001b[0;36m_read_file\u001b[1;34m(filename, bbox, mask, rows, engine, **kwargs)\u001b[0m\n\u001b[0;32m    268\u001b[0m             from_bytes \u001b[38;5;241m=\u001b[39m \u001b[38;5;28;01mTrue\u001b[39;00m\n\u001b[0;32m    270\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m engine \u001b[38;5;241m==\u001b[39m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mpyogrio\u001b[39m\u001b[38;5;124m\"\u001b[39m:\n\u001b[1;32m--> 271\u001b[0m     \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[43m_read_file_pyogrio\u001b[49m\u001b[43m(\u001b[49m\u001b[43mfilename\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mbbox\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mbbox\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mmask\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mmask\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mrows\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mrows\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mkwargs\u001b[49m\u001b[43m)\u001b[49m\n\u001b[0;32m    273\u001b[0m \u001b[38;5;28;01melif\u001b[39;00m engine \u001b[38;5;241m==\u001b[39m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mfiona\u001b[39m\u001b[38;5;124m\"\u001b[39m:\n\u001b[0;32m    274\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m pd\u001b[38;5;241m.\u001b[39mapi\u001b[38;5;241m.\u001b[39mtypes\u001b[38;5;241m.\u001b[39mis_file_like(filename):\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\geopandas\\io\\file.py:427\u001b[0m, in \u001b[0;36m_read_file_pyogrio\u001b[1;34m(path_or_bytes, bbox, mask, rows, **kwargs)\u001b[0m\n\u001b[0;32m    424\u001b[0m     kwargs[\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mread_geometry\u001b[39m\u001b[38;5;124m\"\u001b[39m] \u001b[38;5;241m=\u001b[39m \u001b[38;5;28;01mFalse\u001b[39;00m\n\u001b[0;32m    426\u001b[0m \u001b[38;5;66;03m# TODO: if bbox is not None, check its CRS vs the CRS of the file\u001b[39;00m\n\u001b[1;32m--> 427\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[43mpyogrio\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mread_dataframe\u001b[49m\u001b[43m(\u001b[49m\u001b[43mpath_or_bytes\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mbbox\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mbbox\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mkwargs\u001b[49m\u001b[43m)\u001b[49m\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\geopandas.py:149\u001b[0m, in \u001b[0;36mread_dataframe\u001b[1;34m(path_or_buffer, layer, encoding, columns, read_geometry, force_2d, skip_features, max_features, where, bbox, fids, sql, sql_dialect, fid_as_index, use_arrow, **kwargs)\u001b[0m\n\u001b[0;32m    146\u001b[0m path_or_buffer \u001b[38;5;241m=\u001b[39m _stringify_path(path_or_buffer)\n\u001b[0;32m    148\u001b[0m read_func \u001b[38;5;241m=\u001b[39m read_arrow \u001b[38;5;28;01mif\u001b[39;00m use_arrow \u001b[38;5;28;01melse\u001b[39;00m read\n\u001b[1;32m--> 149\u001b[0m result \u001b[38;5;241m=\u001b[39m \u001b[43mread_func\u001b[49m\u001b[43m(\u001b[49m\n\u001b[0;32m    150\u001b[0m \u001b[43m    \u001b[49m\u001b[43mpath_or_buffer\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    151\u001b[0m \u001b[43m    \u001b[49m\u001b[43mlayer\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mlayer\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    152\u001b[0m \u001b[43m    \u001b[49m\u001b[43mencoding\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mencoding\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    153\u001b[0m \u001b[43m    \u001b[49m\u001b[43mcolumns\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mcolumns\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    154\u001b[0m \u001b[43m    \u001b[49m\u001b[43mread_geometry\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mread_geometry\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    155\u001b[0m \u001b[43m    \u001b[49m\u001b[43mforce_2d\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mforce_2d\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    156\u001b[0m \u001b[43m    \u001b[49m\u001b[43mskip_features\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mskip_features\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    157\u001b[0m \u001b[43m    \u001b[49m\u001b[43mmax_features\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mmax_features\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    158\u001b[0m \u001b[43m    \u001b[49m\u001b[43mwhere\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mwhere\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    159\u001b[0m \u001b[43m    \u001b[49m\u001b[43mbbox\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mbbox\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    160\u001b[0m \u001b[43m    \u001b[49m\u001b[43mfids\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mfids\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    161\u001b[0m \u001b[43m    \u001b[49m\u001b[43msql\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43msql\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    162\u001b[0m \u001b[43m    \u001b[49m\u001b[43msql_dialect\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43msql_dialect\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    163\u001b[0m \u001b[43m    \u001b[49m\u001b[43mreturn_fids\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mfid_as_index\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    164\u001b[0m \u001b[43m    \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mkwargs\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    165\u001b[0m \u001b[43m\u001b[49m\u001b[43m)\u001b[49m\n\u001b[0;32m    167\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m use_arrow:\n\u001b[0;32m    168\u001b[0m     meta, table \u001b[38;5;241m=\u001b[39m result\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\raw.py:137\u001b[0m, in \u001b[0;36mread\u001b[1;34m(path_or_buffer, layer, encoding, columns, read_geometry, force_2d, skip_features, max_features, where, bbox, fids, sql, sql_dialect, return_fids, **kwargs)\u001b[0m\n\u001b[0;32m    134\u001b[0m dataset_kwargs \u001b[38;5;241m=\u001b[39m _preprocess_options_key_value(kwargs) \u001b[38;5;28;01mif\u001b[39;00m kwargs \u001b[38;5;28;01melse\u001b[39;00m {}\n\u001b[0;32m    136\u001b[0m \u001b[38;5;28;01mtry\u001b[39;00m:\n\u001b[1;32m--> 137\u001b[0m     result \u001b[38;5;241m=\u001b[39m \u001b[43mogr_read\u001b[49m\u001b[43m(\u001b[49m\n\u001b[0;32m    138\u001b[0m \u001b[43m        \u001b[49m\u001b[43mpath\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    139\u001b[0m \u001b[43m        \u001b[49m\u001b[43mlayer\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mlayer\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    140\u001b[0m \u001b[43m        \u001b[49m\u001b[43mencoding\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mencoding\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    141\u001b[0m \u001b[43m        \u001b[49m\u001b[43mcolumns\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mcolumns\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    142\u001b[0m \u001b[43m        \u001b[49m\u001b[43mread_geometry\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mread_geometry\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    143\u001b[0m \u001b[43m        \u001b[49m\u001b[43mforce_2d\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mforce_2d\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    144\u001b[0m \u001b[43m        \u001b[49m\u001b[43mskip_features\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mskip_features\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    145\u001b[0m \u001b[43m        \u001b[49m\u001b[43mmax_features\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mmax_features\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;129;43;01mor\u001b[39;49;00m\u001b[43m \u001b[49m\u001b[38;5;241;43m0\u001b[39;49m\u001b[43m,\u001b[49m\n\u001b[0;32m    146\u001b[0m \u001b[43m        \u001b[49m\u001b[43mwhere\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mwhere\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    147\u001b[0m \u001b[43m        \u001b[49m\u001b[43mbbox\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mbbox\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    148\u001b[0m \u001b[43m        \u001b[49m\u001b[43mfids\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mfids\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    149\u001b[0m \u001b[43m        \u001b[49m\u001b[43msql\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43msql\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    150\u001b[0m \u001b[43m        \u001b[49m\u001b[43msql_dialect\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43msql_dialect\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    151\u001b[0m \u001b[43m        \u001b[49m\u001b[43mreturn_fids\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mreturn_fids\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    152\u001b[0m \u001b[43m        \u001b[49m\u001b[43mdataset_kwargs\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mdataset_kwargs\u001b[49m\u001b[43m,\u001b[49m\n\u001b[0;32m    153\u001b[0m \u001b[43m    \u001b[49m\u001b[43m)\u001b[49m\n\u001b[0;32m    154\u001b[0m \u001b[38;5;28;01mfinally\u001b[39;00m:\n\u001b[0;32m    155\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m buffer \u001b[38;5;129;01mis\u001b[39;00m \u001b[38;5;129;01mnot\u001b[39;00m \u001b[38;5;28;01mNone\u001b[39;00m:\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\_io.pyx:971\u001b[0m, in \u001b[0;36mpyogrio._io.ogr_read\u001b[1;34m()\u001b[0m\n",
      "File \u001b[1;32m~\\workspace\\github\\estacoesMetereoligicas-piaui\\.venv\\Lib\\site-packages\\pyogrio\\_io.pyx:164\u001b[0m, in \u001b[0;36mpyogrio._io.ogr_open\u001b[1;34m()\u001b[0m\n",
      "\u001b[1;31mDataSourceError\u001b[0m: '/vsicurl/https://metadados.snirh.gov.br/files/f85dbf06-a869-414c-afc5-bb01869e9156/Estacoes_Fluviometricas_e_Pluviometricas_da_Rede.gpkg' does not exist in the file system, and is not recognized as a supported dataset name."
     ]
    }
   ],
   "source": [
    "\n",
    "gpkg_gdf = gpd.read_file(url_dataset, layer=\"Estacoes_Fluviometricas_e_Pluviometricas_da_Rede\", engine='pyogrio')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "id": "738b4f65-a12d-4368-8897-6fded41df70a",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "gdf_br = gpd.read_file(path_gpkg_br)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "id": "531e3a5b-75c2-4631-ba3d-6f699d8fc445",
   "metadata": {
    "tags": []
   },
   "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>OBJECTID</th>\n",
       "      <th>Nome</th>\n",
       "      <th>TipoEstacao</th>\n",
       "      <th>TipoEstacaoCodigo</th>\n",
       "      <th>Operando</th>\n",
       "      <th>CodigoAdicional</th>\n",
       "      <th>Latitude</th>\n",
       "      <th>Longitude</th>\n",
       "      <th>Altitude</th>\n",
       "      <th>AreaDrenagem</th>\n",
       "      <th>...</th>\n",
       "      <th>RedeNavegacao</th>\n",
       "      <th>RedeCursoDagua</th>\n",
       "      <th>RedeEstrategica</th>\n",
       "      <th>RedeCaptacao</th>\n",
       "      <th>RedeRHNR</th>\n",
       "      <th>RedeQA</th>\n",
       "      <th>RedeClasseVazao</th>\n",
       "      <th>Descricao</th>\n",
       "      <th>DataAlteracao</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>PATATE DJ AMBATO</td>\n",
       "      <td>Fluviométrica</td>\n",
       "      <td>1</td>\n",
       "      <td>Sim</td>\n",
       "      <td>H-817</td>\n",
       "      <td>-1.2500</td>\n",
       "      <td>-78.5100</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2010-09-09</td>\n",
       "      <td>POINT (-78.51000 -1.25000)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>NEGRO AJ PUMANCUCHI</td>\n",
       "      <td>Fluviométrica</td>\n",
       "      <td>1</td>\n",
       "      <td>Sim</td>\n",
       "      <td>H-795</td>\n",
       "      <td>-0.7758</td>\n",
       "      <td>-78.7731</td>\n",
       "      <td>3060.0</td>\n",
       "      <td>91.0</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaT</td>\n",
       "      <td>POINT (-78.77310 -0.77580)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>NAGSICHE P E CUSUBAMBA</td>\n",
       "      <td>Fluviométrica</td>\n",
       "      <td>1</td>\n",
       "      <td>Sim</td>\n",
       "      <td>H-793</td>\n",
       "      <td>-1.0617</td>\n",
       "      <td>-78.6903</td>\n",
       "      <td>2960.0</td>\n",
       "      <td>275.0</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaT</td>\n",
       "      <td>POINT (-78.69030 -1.06170)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4</td>\n",
       "      <td>AMBATO EN AMBATO</td>\n",
       "      <td>Fluviométrica</td>\n",
       "      <td>1</td>\n",
       "      <td>Sim</td>\n",
       "      <td>H-801</td>\n",
       "      <td>-1.2589</td>\n",
       "      <td>-78.6628</td>\n",
       "      <td>2680.0</td>\n",
       "      <td>763.0</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaT</td>\n",
       "      <td>POINT (-78.66280 -1.25890)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>5</td>\n",
       "      <td>PACHANLICA AJ AMBATO</td>\n",
       "      <td>Fluviométrica</td>\n",
       "      <td>1</td>\n",
       "      <td>Sim</td>\n",
       "      <td>H-803</td>\n",
       "      <td>-1.2364</td>\n",
       "      <td>-78.5494</td>\n",
       "      <td>2310.0</td>\n",
       "      <td>407.0</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaT</td>\n",
       "      <td>POINT (-78.54940 -1.23640)</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 74 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "   OBJECTID                    Nome    TipoEstacao  TipoEstacaoCodigo  \\\n",
       "0         1        PATATE DJ AMBATO  Fluviométrica                  1   \n",
       "1         2     NEGRO AJ PUMANCUCHI  Fluviométrica                  1   \n",
       "2         3  NAGSICHE P E CUSUBAMBA  Fluviométrica                  1   \n",
       "3         4        AMBATO EN AMBATO  Fluviométrica                  1   \n",
       "4         5    PACHANLICA AJ AMBATO  Fluviométrica                  1   \n",
       "\n",
       "  Operando CodigoAdicional  Latitude  Longitude  Altitude  AreaDrenagem  ...  \\\n",
       "0      Sim           H-817   -1.2500   -78.5100       NaN           NaN  ...   \n",
       "1      Sim           H-795   -0.7758   -78.7731    3060.0          91.0  ...   \n",
       "2      Sim           H-793   -1.0617   -78.6903    2960.0         275.0  ...   \n",
       "3      Sim           H-801   -1.2589   -78.6628    2680.0         763.0  ...   \n",
       "4      Sim           H-803   -1.2364   -78.5494    2310.0         407.0  ...   \n",
       "\n",
       "  RedeNavegacao  RedeCursoDagua RedeEstrategica  RedeCaptacao RedeRHNR  \\\n",
       "0           Não            Nulo             Não          Nulo      Não   \n",
       "1           Não            Nulo             Não          Nulo      Não   \n",
       "2           Não            Nulo             Não          Nulo      Não   \n",
       "3           Não            Nulo             Não          Nulo      Não   \n",
       "4           Não            Nulo             Não          Nulo      Não   \n",
       "\n",
       "   RedeQA RedeClasseVazao  Descricao DataAlteracao                    geometry  \n",
       "0     Não             Sim        NaN    2010-09-09  POINT (-78.51000 -1.25000)  \n",
       "1     Não             Sim        NaN           NaT  POINT (-78.77310 -0.77580)  \n",
       "2     Não             Sim        NaN           NaT  POINT (-78.69030 -1.06170)  \n",
       "3     Não             Sim        NaN           NaT  POINT (-78.66280 -1.25890)  \n",
       "4     Não             Sim        NaN           NaT  POINT (-78.54940 -1.23640)  \n",
       "\n",
       "[5 rows x 74 columns]"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_br.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "id": "d67c1efc-5ad4-4959-9dba-5f187b45ad96",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['OBJECTID', 'Nome', 'TipoEstacao', 'TipoEstacaoCodigo', 'Operando',\n",
       "       'CodigoAdicional', 'Latitude', 'Longitude', 'Altitude', 'AreaDrenagem',\n",
       "       'Bacia', 'BaciaCodigo', 'SubBacia', 'SubBaciaCodigo', 'Rio',\n",
       "       'RioCodigo', 'UF', 'UFCodigo', 'Municipio', 'MunicipioCodigo',\n",
       "       'Responsavel', 'ResponsavelSigla', 'ResponsavelCodigo',\n",
       "       'ResponsavelUnidade', 'Operadora', 'OperadoraSigla', 'OperadoraCodigo',\n",
       "       'OperadoraUnidade', 'OperadoraRoteiro', 'EscalaNivel',\n",
       "       'EscalaNivelInicio', 'EscalaNivelFim', 'RegistradorNivel',\n",
       "       'RegistradorNivelInicio', 'RegistradorNivelFim',\n",
       "       'MedicaoDescargaLiquida', 'MedicaoDescargaLiquidaInicio',\n",
       "       'MedicaoDescargaLiquidaInicioFim', 'MedicaoDescargaSolida',\n",
       "       'MedicaoDescargaSolidaInicio', 'MedicaoDescargaSolidaFim', 'MedicaoQA',\n",
       "       'MedicaoQAInicio', 'MedicaoQAFim', 'PluviometroConvencional',\n",
       "       'PluviometroConvencionalInicio', 'PluviometroConvencionalFim',\n",
       "       'RegistradorChuva', 'RegistradorChuvaInicio', 'RegistradorChuvaFim',\n",
       "       'TanqueEvaporimetrico', 'TanqueEvaporimetricoInicio',\n",
       "       'TanqueEvaporimetricoFim', 'EstacaoClimatologica',\n",
       "       'EstacaoClimatologicaInicio', 'EstacaoClimatologicaFim',\n",
       "       'EstacaoPiezometria', 'EstacaoPiezometriaInicio',\n",
       "       'EstacaoPiezometriaFim', 'EstacaoTelemetrica',\n",
       "       'EstacaoTelemetricaInicio', 'EstacaoTelemetricaFim', 'RedeBasica',\n",
       "       'RedeEnergetica', 'RedeNavegacao', 'RedeCursoDagua', 'RedeEstrategica',\n",
       "       'RedeCaptacao', 'RedeRHNR', 'RedeQA', 'RedeClasseVazao', 'Descricao',\n",
       "       'DataAlteracao', 'geometry'],\n",
       "      dtype='object')"
      ]
     },
     "execution_count": 35,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_br.columns"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3f812bd2-71ec-4406-ade1-310eec171148",
   "metadata": {},
   "source": [
    "## Filtrando todas as estações contidas no PIAUÍ"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "id": "ead0ee88-0723-4809-981f-07dba7fd2caf",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "UF = 'PIUAÍ'\n",
    "gdf_pi = gdf_br.query(\"UF == @UF\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "id": "fd8eb05d-4a18-41e2-ba58-d6ec00078c06",
   "metadata": {
    "tags": []
   },
   "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>OBJECTID</th>\n",
       "      <th>Nome</th>\n",
       "      <th>TipoEstacao</th>\n",
       "      <th>TipoEstacaoCodigo</th>\n",
       "      <th>Operando</th>\n",
       "      <th>CodigoAdicional</th>\n",
       "      <th>Latitude</th>\n",
       "      <th>Longitude</th>\n",
       "      <th>Altitude</th>\n",
       "      <th>AreaDrenagem</th>\n",
       "      <th>...</th>\n",
       "      <th>RedeNavegacao</th>\n",
       "      <th>RedeCursoDagua</th>\n",
       "      <th>RedeEstrategica</th>\n",
       "      <th>RedeCaptacao</th>\n",
       "      <th>RedeRHNR</th>\n",
       "      <th>RedeQA</th>\n",
       "      <th>RedeClasseVazao</th>\n",
       "      <th>Descricao</th>\n",
       "      <th>DataAlteracao</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1596</th>\n",
       "      <td>1597</td>\n",
       "      <td>BOA ESPERANÇA</td>\n",
       "      <td>Fluviométrica</td>\n",
       "      <td>1</td>\n",
       "      <td>Não</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-6.750000</td>\n",
       "      <td>-43.566667</td>\n",
       "      <td>85.0</td>\n",
       "      <td>84900.0</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaT</td>\n",
       "      <td>POINT (-43.56667 -6.75000)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1597</th>\n",
       "      <td>1598</td>\n",
       "      <td>BOA ESPERANÇA II</td>\n",
       "      <td>Fluviométrica</td>\n",
       "      <td>1</td>\n",
       "      <td>Não</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-6.800000</td>\n",
       "      <td>-43.516667</td>\n",
       "      <td>85.0</td>\n",
       "      <td>85100.0</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaT</td>\n",
       "      <td>POINT (-43.51667 -6.80000)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1598</th>\n",
       "      <td>1599</td>\n",
       "      <td>VEADOS</td>\n",
       "      <td>Fluviométrica</td>\n",
       "      <td>1</td>\n",
       "      <td>Não</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-6.833333</td>\n",
       "      <td>-43.483333</td>\n",
       "      <td>80.0</td>\n",
       "      <td>88700.0</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaT</td>\n",
       "      <td>POINT (-43.48333 -6.83333)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1599</th>\n",
       "      <td>1600</td>\n",
       "      <td>FAZENDA BANANEIRA</td>\n",
       "      <td>Fluviométrica</td>\n",
       "      <td>1</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-10.016400</td>\n",
       "      <td>-45.024700</td>\n",
       "      <td>420.0</td>\n",
       "      <td>681.0</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaT</td>\n",
       "      <td>POINT (-45.02470 -10.01640)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1600</th>\n",
       "      <td>1601</td>\n",
       "      <td>SÃO GONÇALO</td>\n",
       "      <td>Fluviométrica</td>\n",
       "      <td>1</td>\n",
       "      <td>Não</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-10.033600</td>\n",
       "      <td>-45.303600</td>\n",
       "      <td>NaN</td>\n",
       "      <td>699.0</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2018-07-20</td>\n",
       "      <td>POINT (-45.30360 -10.03360)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29952</th>\n",
       "      <td>29953</td>\n",
       "      <td>BARRA</td>\n",
       "      <td>Pluvometrica</td>\n",
       "      <td>2</td>\n",
       "      <td>Não</td>\n",
       "      <td>4609187</td>\n",
       "      <td>-10.050000</td>\n",
       "      <td>-45.066700</td>\n",
       "      <td>490.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2010-04-28</td>\n",
       "      <td>POINT (-45.06670 -10.05000)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29953</th>\n",
       "      <td>29954</td>\n",
       "      <td>GALHEIRO</td>\n",
       "      <td>Pluvometrica</td>\n",
       "      <td>2</td>\n",
       "      <td>Não</td>\n",
       "      <td>4609127</td>\n",
       "      <td>-10.050000</td>\n",
       "      <td>-45.366700</td>\n",
       "      <td>630.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2010-04-28</td>\n",
       "      <td>POINT (-45.36670 -10.05000)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29954</th>\n",
       "      <td>29955</td>\n",
       "      <td>JENIPAPEIRO</td>\n",
       "      <td>Pluvometrica</td>\n",
       "      <td>2</td>\n",
       "      <td>Não</td>\n",
       "      <td>4619364</td>\n",
       "      <td>-10.650000</td>\n",
       "      <td>-45.183300</td>\n",
       "      <td>510.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2010-04-28</td>\n",
       "      <td>POINT (-45.18330 -10.65000)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29955</th>\n",
       "      <td>29956</td>\n",
       "      <td>FAZENDA DO BENÉ</td>\n",
       "      <td>Pluvometrica</td>\n",
       "      <td>2</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-10.153300</td>\n",
       "      <td>-45.192200</td>\n",
       "      <td>456.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaT</td>\n",
       "      <td>POINT (-45.19220 -10.15330)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29956</th>\n",
       "      <td>29957</td>\n",
       "      <td>SÍTIO BARRA DO BREJO</td>\n",
       "      <td>Pluvometrica</td>\n",
       "      <td>2</td>\n",
       "      <td>Sim</td>\n",
       "      <td>NaN</td>\n",
       "      <td>-10.044200</td>\n",
       "      <td>-45.325000</td>\n",
       "      <td>440.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Nulo</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>Não</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaT</td>\n",
       "      <td>POINT (-45.32500 -10.04420)</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>544 rows × 74 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       OBJECTID                  Nome    TipoEstacao  TipoEstacaoCodigo  \\\n",
       "1596       1597         BOA ESPERANÇA  Fluviométrica                  1   \n",
       "1597       1598      BOA ESPERANÇA II  Fluviométrica                  1   \n",
       "1598       1599                VEADOS  Fluviométrica                  1   \n",
       "1599       1600     FAZENDA BANANEIRA  Fluviométrica                  1   \n",
       "1600       1601           SÃO GONÇALO  Fluviométrica                  1   \n",
       "...         ...                   ...            ...                ...   \n",
       "29952     29953                 BARRA   Pluvometrica                  2   \n",
       "29953     29954              GALHEIRO   Pluvometrica                  2   \n",
       "29954     29955           JENIPAPEIRO   Pluvometrica                  2   \n",
       "29955     29956       FAZENDA DO BENÉ   Pluvometrica                  2   \n",
       "29956     29957  SÍTIO BARRA DO BREJO   Pluvometrica                  2   \n",
       "\n",
       "      Operando CodigoAdicional   Latitude  Longitude  Altitude  AreaDrenagem  \\\n",
       "1596       Não             NaN  -6.750000 -43.566667      85.0       84900.0   \n",
       "1597       Não             NaN  -6.800000 -43.516667      85.0       85100.0   \n",
       "1598       Não             NaN  -6.833333 -43.483333      80.0       88700.0   \n",
       "1599       Sim             NaN -10.016400 -45.024700     420.0         681.0   \n",
       "1600       Não             NaN -10.033600 -45.303600       NaN         699.0   \n",
       "...        ...             ...        ...        ...       ...           ...   \n",
       "29952      Não         4609187 -10.050000 -45.066700     490.0           NaN   \n",
       "29953      Não         4609127 -10.050000 -45.366700     630.0           NaN   \n",
       "29954      Não         4619364 -10.650000 -45.183300     510.0           NaN   \n",
       "29955      Sim             NaN -10.153300 -45.192200     456.0           NaN   \n",
       "29956      Sim             NaN -10.044200 -45.325000     440.0           NaN   \n",
       "\n",
       "       ... RedeNavegacao  RedeCursoDagua RedeEstrategica  RedeCaptacao  \\\n",
       "1596   ...           Não            Nulo             Não          Nulo   \n",
       "1597   ...           Não            Nulo             Não          Nulo   \n",
       "1598   ...           Não            Nulo             Não          Nulo   \n",
       "1599   ...           Não            Nulo             Não          Nulo   \n",
       "1600   ...           Não            Nulo             Não          Nulo   \n",
       "...    ...           ...             ...             ...           ...   \n",
       "29952  ...           Não             Não             Não          Nulo   \n",
       "29953  ...           Não            Nulo             Não          Nulo   \n",
       "29954  ...           Não             Não             Não          Nulo   \n",
       "29955  ...           Não             Não             Não          Nulo   \n",
       "29956  ...           Não             Não             Não          Nulo   \n",
       "\n",
       "      RedeRHNR  RedeQA RedeClasseVazao  Descricao DataAlteracao  \\\n",
       "1596       Não    Nulo             Sim        NaN           NaT   \n",
       "1597       Não    Nulo             Sim        NaN           NaT   \n",
       "1598       Não     Não             Não        NaN           NaT   \n",
       "1599       Não     Não             Sim        NaN           NaT   \n",
       "1600       Não    Nulo             Sim        NaN    2018-07-20   \n",
       "...        ...     ...             ...        ...           ...   \n",
       "29952      Não     Não             Não        NaN    2010-04-28   \n",
       "29953      Não    Nulo             Não        NaN    2010-04-28   \n",
       "29954      Não     Não             Não        NaN    2010-04-28   \n",
       "29955      Não     Não             Não        NaN           NaT   \n",
       "29956      Não     Não             Não        NaN           NaT   \n",
       "\n",
       "                          geometry  \n",
       "1596    POINT (-43.56667 -6.75000)  \n",
       "1597    POINT (-43.51667 -6.80000)  \n",
       "1598    POINT (-43.48333 -6.83333)  \n",
       "1599   POINT (-45.02470 -10.01640)  \n",
       "1600   POINT (-45.30360 -10.03360)  \n",
       "...                            ...  \n",
       "29952  POINT (-45.06670 -10.05000)  \n",
       "29953  POINT (-45.36670 -10.05000)  \n",
       "29954  POINT (-45.18330 -10.65000)  \n",
       "29955  POINT (-45.19220 -10.15330)  \n",
       "29956  POINT (-45.32500 -10.04420)  \n",
       "\n",
       "[544 rows x 74 columns]"
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_pi"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c2375194-8cbe-4c52-b055-a41b99c023f3",
   "metadata": {},
   "source": [
    "## Filtrando todas as estações contidas no Piauí e em Operação"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "id": "46f9c108-cb9f-4e53-b9cd-1f92bd6c2418",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(257, 74)"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "query_on = 'Sim'\n",
    "gdf_pi_on = gdf_pi.query(\"Operando == @query_on\") # retorna GeoDataFrame\n",
    "\n",
    "gdf_pi_on.shape #Retorna um tupla"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "id": "96f66c97-cd0d-416e-b2c6-109891eaf1bf",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "gdf_pi_on.to_parquet('./dados/ana/EstacoesFluvioPluvio_PiauiOperando.parquet')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "id": "b069a64e-27af-4354-a941-642f3deeb852",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "gdf_pi_on.to_file('./dados/ana/EstacoesFluvioPluvio_PiauiOperando.geojson', Driver='GeoJson')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "id": "9a45ce4a-c619-4376-9ece-b225c775fd56",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['CPRM', 'ANA', 'DNOCS', 'SEMARH-PI', 'SUDENE', 'SEMA-MA', 'CHESF',\n",
       "       'CEMADEN', 'INMET', 'DEPV', 'EMATER', 'SNIRH-MA'], dtype=object)"
      ]
     },
     "execution_count": 56,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_pi_on.ResponsavelSigla.unique()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "id": "aebd30ba-2413-4daf-be78-4a10d1744413",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Responsavel\n",
       "Agência Nacional de Águas                                                 108\n",
       "Centro Nacional de Monitoramento e Alertas de Desastres Naturais.          28\n",
       "Cia de Pesquisa de Recursos Minerais                                       10\n",
       "Companhia Hidro Elétrica do São Francisco-CHESF                             8\n",
       "Departamento Nacional de Obras Contra as Secas                             46\n",
       "Diretoria de Eletrônica e Proteção ao Voo do Ministério da Aeronáutica      1\n",
       "Empresa Estadual de Assistência Tec. e Ext. Rural                           2\n",
       "Instituto Nacional de Meteorologia                                         15\n",
       "Secretaria de Estado de Meio Ambiente e Recursos Naturais - Maranhão        1\n",
       "Secretaria do Meio Ambiente e Recursos Hídricos - PI                       36\n",
       "Sistema Nacional de Informações sobre Recursos Hídricos - Maranhão          1\n",
       "Superintendência do Desenvolvimento do Nordeste                             1\n",
       "Name: OperadoraSigla, dtype: int64"
      ]
     },
     "execution_count": 63,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_pi_on.OperadoraSigla.groupby(gdf_pi_on['Responsavel']).count()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "df7d2890-efe1-4e0a-9a90-d15f0055a5c2",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}