tunned-geobr 0.2.2__py3-none-any.whl → 0.2.4__py3-none-any.whl

tunned_geobr/read_planned_biomass_ute.py

@@ -0,0 +1,166 @@
+import geopandas as gpd
+import os
+import tempfile
+import urllib.parse
+import requests
+import shutil
+from zipfile import ZipFile
+from pathlib import Path
+from io import BytesIO
+import warnings
+import json
+
+def read_planned_biomass_ute(simplified=False, verbose=False):
+    """Download Planned Biomass Thermoelectric Power Plants data from EPE.
+    
+    This function downloads and processes planned biomass thermoelectric power plants data from EPE 
+    (Energy Research Company). The data includes information about planned biomass thermoelectric 
+    power generation projects across Brazil.
+    Original source: EPE (Empresa de Pesquisa Energética)
+    
+    Parameters
+    ----------
+    simplified : boolean, by default False
+        If True, returns a simplified version of the dataset with fewer columns
+    verbose : boolean, by default False
+        If True, prints detailed information about the download process
+
+    Returns
+    -------
+    gpd.GeoDataFrame
+        Geodataframe with planned biomass thermoelectric power plants data
+        
+    Example
+    -------
+    >>> from tunned_geobr import read_planned_biomass_ute
+    
+    # Read planned biomass thermoelectric power plants data
+    >>> planned_biomass_ute = read_planned_biomass_ute()
+    """
+    
+    # URL for the EPE geoserver
+    url = r"https://gisepeprd2.epe.gov.br/arcgis/rest/services/Download_Dados_Webmap_EPE/GPServer/Extract%20Data%20Task/execute?f=json&env%3AoutSR=102100&Layers_to_Clip=%5B%22UTE%20Biomassa%20-%20Expans%C3%A3o%20Planejada%22%5D&Area_of_Interest=%7B%22geometryType%22%3A%22esriGeometryPolygon%22%2C%22features%22%3A%5B%7B%22geometry%22%3A%7B%22rings%22%3A%5B%5B%5B-8655251.47456396%2C-4787514.465591563%5D%2C%5B-8655251.47456396%2C1229608.401015912%5D%2C%5B-3508899.2341809804%2C1229608.401015912%5D%2C%5B-3508899.2341809804%2C-4787514.465591563%5D%2C%5B-8655251.47456396%2C-4787514.465591563%5D%5D%5D%2C%22spatialReference%22%3A%7B%22wkid%22%3A102100%7D%7D%7D%5D%2C%22sr%22%3A%7B%22wkid%22%3A102100%7D%7D&Feature_Format=Shapefile%20-%20SHP%20-%20.shp&Raster_Format=Tagged%20Image%20File%20Format%20-%20TIFF%20-%20.tif"
+ 
+    try:
+        # Disable SSL verification warning
+        warnings.filterwarnings('ignore', message='Unverified HTTPS request')
+        
+        if verbose:
+            print("Requesting data from EPE server...")
+            
+        response = requests.get(url, timeout=60, verify=False)
+        if not response.ok:
+            raise Exception(f"Error getting JSON response: {response.status_code}")
+
+        json_response = response.json()
+        
+        if verbose:
+            print(f"JSON response received: {json.dumps(json_response, indent=2)[:500]}...")
+            
+        if 'results' not in json_response or len(json_response['results']) == 0:
+            raise Exception("Invalid JSON response structure")
+            
+        if 'value' not in json_response['results'][0] or 'url' not in json_response['results'][0]['value']:
+            raise Exception("URL not found in JSON response")
+            
+        file_url = json_response['results'][0]['value']['url']
+        
+        if verbose:
+            print(f"Downloading file from: {file_url}")
+        
+        file_response = requests.get(file_url, stream=True, timeout=60, verify=False)
+        if not file_response.ok:
+            raise Exception(f"Error downloading file: {file_response.status_code}")
+        
+        # Check if content is actually a zip file
+        content = file_response.content
+        if len(content) < 100:
+            if verbose:
+                print(f"Warning: Downloaded content is very small ({len(content)} bytes)")
+                print(f"Content preview: {content[:100]}")
+            
+        # Create a temporary directory to extract the files
+        with tempfile.TemporaryDirectory() as temp_dir:
+            if verbose:
+                print(f"Extracting files to temporary directory: {temp_dir}")
+                
+            try:
+                # Extract the zip file
+                with ZipFile(BytesIO(content)) as zip_ref:
+                    zip_ref.extractall(temp_dir)
+                    
+                    if verbose:
+                        print(f"Files in zip: {zip_ref.namelist()}")
+            except Exception as zip_error:
+                if verbose:
+                    print(f"Error extracting zip: {str(zip_error)}")
+                    print(f"Saving content to debug.zip for inspection")
+                    with open("debug.zip", "wb") as f:
+                        f.write(content)
+                raise Exception(f"Failed to extract zip file: {str(zip_error)}")
+            
+            # Find the shapefile
+            all_files = os.listdir(temp_dir)
+            if verbose:
+                print(f"Files in temp directory: {all_files}")
+                
+            shp_files = [f for f in all_files if f.endswith('.shp')]
+            if not shp_files:
+                # Try looking in subdirectories
+                for root, dirs, files in os.walk(temp_dir):
+                    shp_files.extend([os.path.join(root, f) for f in files if f.endswith('.shp')])
+                
+                if not shp_files:
+                    raise Exception("No shapefile found in the downloaded data")
+            
+            # Read the shapefile
+            shp_path = shp_files[0] if os.path.isabs(shp_files[0]) else os.path.join(temp_dir, shp_files[0])
+            if verbose:
+                print(f"Reading shapefile: {shp_path}")
+                
+            gdf = gpd.read_file(shp_path)
+            
+            # Convert to SIRGAS 2000 (EPSG:4674)
+            gdf = gdf.to_crs(4674)
+            
+            if verbose:
+                print(f"Data loaded successfully with {len(gdf)} records")
+                print(f"Columns: {gdf.columns.tolist()}")
+            
+            if simplified:
+                # Keep only the most relevant columns
+                columns_to_keep = [
+                    'geometry',
+                    'nome',            # Power plant name
+                    'potencia',        # Capacity in MW
+                    'combustivel',     # Fuel type
+                    'leilao',          # Auction
+                    'ceg',             # CEG code
+                    'ano_prev'         # Expected year
+                ]
+                
+                # Filter columns that actually exist in the dataset
+                existing_columns = ['geometry'] + [col for col in columns_to_keep[1:] if col in gdf.columns]
+                if len(existing_columns) <= 1:
+                    if verbose:
+                        print("Warning: No matching columns found for simplified version. Returning all columns.")
+                else:
+                    gdf = gdf[existing_columns]
+    
+    except Exception as e:
+        raise Exception(f"Error downloading or processing planned biomass thermoelectric power plants data: {str(e)}")
+        
+    return gdf
+
+if __name__ == '__main__':
+    try:
+        biomass_ute_data = read_planned_biomass_ute(verbose=True)
+        print(f"Downloaded planned biomass thermoelectric power plants data with {len(biomass_ute_data)} records and {len(biomass_ute_data.columns)} columns")
+        
+        # Test simplified version
+        simplified_data = read_planned_biomass_ute(simplified=True)
+        print(f"Simplified data has {len(simplified_data.columns)} columns: {simplified_data.columns.tolist()}")
+    except Exception as e:
+        print(f"Error: {str(e)}")
+        import traceback
+        traceback.print_exc()

tunned_geobr/read_planned_cgh.py

@@ -0,0 +1,166 @@
+import geopandas as gpd
+import os
+import tempfile
+import urllib.parse
+import requests
+import shutil
+from zipfile import ZipFile
+from pathlib import Path
+from io import BytesIO
+import warnings
+import json
+
+def read_planned_cgh(simplified=False, verbose=False):
+    """Download Planned Mini Hydroelectric Power Plants (CGH) data from EPE.
+    
+    This function downloads and processes planned mini hydroelectric power plants (CGH) data from EPE 
+    (Energy Research Company). The data includes information about planned mini hydroelectric 
+    power generation projects across Brazil.
+    Original source: EPE (Empresa de Pesquisa Energética)
+    
+    Parameters
+    ----------
+    simplified : boolean, by default False
+        If True, returns a simplified version of the dataset with fewer columns
+    verbose : boolean, by default False
+        If True, prints detailed information about the download process
+
+    Returns
+    -------
+    gpd.GeoDataFrame
+        Geodataframe with planned mini hydroelectric power plants data
+        
+    Example
+    -------
+    >>> from tunned_geobr import read_planned_cgh
+    
+    # Read planned mini hydroelectric power plants data
+    >>> planned_cgh = read_planned_cgh()
+    """
+    
+    # URL for the EPE geoserver
+    url = r"https://gisepeprd2.epe.gov.br/arcgis/rest/services/Download_Dados_Webmap_EPE/GPServer/Extract%20Data%20Task/execute?f=json&env%3AoutSR=102100&Layers_to_Clip=%5B%22CGH%20-%20Expans%C3%A3o%20Planejada%22%5D&Area_of_Interest=%7B%22geometryType%22%3A%22esriGeometryPolygon%22%2C%22features%22%3A%5B%7B%22geometry%22%3A%7B%22rings%22%3A%5B%5B%5B-8655251.47456396%2C-4787514.465591563%5D%2C%5B-8655251.47456396%2C1229608.401015912%5D%2C%5B-3508899.2341809804%2C1229608.401015912%5D%2C%5B-3508899.2341809804%2C-4787514.465591563%5D%2C%5B-8655251.47456396%2C-4787514.465591563%5D%5D%5D%2C%22spatialReference%22%3A%7B%22wkid%22%3A102100%7D%7D%7D%5D%2C%22sr%22%3A%7B%22wkid%22%3A102100%7D%7D&Feature_Format=Shapefile%20-%20SHP%20-%20.shp&Raster_Format=Tagged%20Image%20File%20Format%20-%20TIFF%20-%20.tif"
+ 
+    try:
+        # Disable SSL verification warning
+        warnings.filterwarnings('ignore', message='Unverified HTTPS request')
+        
+        if verbose:
+            print("Requesting data from EPE server...")
+            
+        response = requests.get(url, timeout=60, verify=False)
+        if not response.ok:
+            raise Exception(f"Error getting JSON response: {response.status_code}")
+
+        json_response = response.json()
+        
+        if verbose:
+            print(f"JSON response received: {json.dumps(json_response, indent=2)[:500]}...")
+            
+        if 'results' not in json_response or len(json_response['results']) == 0:
+            raise Exception("Invalid JSON response structure")
+            
+        if 'value' not in json_response['results'][0] or 'url' not in json_response['results'][0]['value']:
+            raise Exception("URL not found in JSON response")
+            
+        file_url = json_response['results'][0]['value']['url']
+        
+        if verbose:
+            print(f"Downloading file from: {file_url}")
+        
+        file_response = requests.get(file_url, stream=True, timeout=60, verify=False)
+        if not file_response.ok:
+            raise Exception(f"Error downloading file: {file_response.status_code}")
+        
+        # Check if content is actually a zip file
+        content = file_response.content
+        if len(content) < 100:
+            if verbose:
+                print(f"Warning: Downloaded content is very small ({len(content)} bytes)")
+                print(f"Content preview: {content[:100]}")
+            
+        # Create a temporary directory to extract the files
+        with tempfile.TemporaryDirectory() as temp_dir:
+            if verbose:
+                print(f"Extracting files to temporary directory: {temp_dir}")
+                
+            try:
+                # Extract the zip file
+                with ZipFile(BytesIO(content)) as zip_ref:
+                    zip_ref.extractall(temp_dir)
+                    
+                    if verbose:
+                        print(f"Files in zip: {zip_ref.namelist()}")
+            except Exception as zip_error:
+                if verbose:
+                    print(f"Error extracting zip: {str(zip_error)}")
+                    print(f"Saving content to debug.zip for inspection")
+                    with open("debug.zip", "wb") as f:
+                        f.write(content)
+                raise Exception(f"Failed to extract zip file: {str(zip_error)}")
+            
+            # Find the shapefile
+            all_files = os.listdir(temp_dir)
+            if verbose:
+                print(f"Files in temp directory: {all_files}")
+                
+            shp_files = [f for f in all_files if f.endswith('.shp')]
+            if not shp_files:
+                # Try looking in subdirectories
+                for root, dirs, files in os.walk(temp_dir):
+                    shp_files.extend([os.path.join(root, f) for f in files if f.endswith('.shp')])
+                
+                if not shp_files:
+                    raise Exception("No shapefile found in the downloaded data")
+            
+            # Read the shapefile
+            shp_path = shp_files[0] if os.path.isabs(shp_files[0]) else os.path.join(temp_dir, shp_files[0])
+            if verbose:
+                print(f"Reading shapefile: {shp_path}")
+                
+            gdf = gpd.read_file(shp_path)
+            
+            # Convert to SIRGAS 2000 (EPSG:4674)
+            gdf = gdf.to_crs(4674)
+            
+            if verbose:
+                print(f"Data loaded successfully with {len(gdf)} records")
+                print(f"Columns: {gdf.columns.tolist()}")
+            
+            if simplified:
+                # Keep only the most relevant columns
+                columns_to_keep = [
+                    'geometry',
+                    'nome',            # Power plant name
+                    'potencia',        # Capacity in MW
+                    'rio',             # River name
+                    'bacia',           # Basin
+                    'sub_bacia',       # Sub-basin
+                    'ano_prev'         # Expected year
+                ]
+                
+                # Filter columns that actually exist in the dataset
+                existing_columns = ['geometry'] + [col for col in columns_to_keep[1:] if col in gdf.columns]
+                if len(existing_columns) <= 1:
+                    if verbose:
+                        print("Warning: No matching columns found for simplified version. Returning all columns.")
+                else:
+                    gdf = gdf[existing_columns]
+    
+    except Exception as e:
+        raise Exception(f"Error downloading or processing planned mini hydroelectric power plants data: {str(e)}")
+        
+    return gdf
+
+if __name__ == '__main__':
+    try:
+        cgh_data = read_planned_cgh(verbose=True)
+        print(f"Downloaded planned mini hydroelectric power plants data with {len(cgh_data)} records and {len(cgh_data.columns)} columns")
+        
+        # Test simplified version
+        simplified_data = read_planned_cgh(simplified=True)
+        print(f"Simplified data has {len(simplified_data.columns)} columns: {simplified_data.columns.tolist()}")
+    except Exception as e:
+        print(f"Error: {str(e)}")
+        import traceback
+        traceback.print_exc()

tunned_geobr/read_planned_eolic.py

@@ -0,0 +1,165 @@
+import geopandas as gpd
+import os
+import tempfile
+import urllib.parse
+import requests
+import shutil
+from zipfile import ZipFile
+from pathlib import Path
+from io import BytesIO
+import warnings
+import json
+
+def read_planned_eolic(simplified=False, verbose=False):
+    """Download Planned Eolic Power Plants data from EPE.
+    
+    This function downloads and processes planned eolic (wind) power plants data from EPE 
+    (Energy Research Company). The data includes information about planned wind power 
+    generation projects across Brazil.
+    Original source: EPE (Empresa de Pesquisa Energética)
+    
+    Parameters
+    ----------
+    simplified : boolean, by default False
+        If True, returns a simplified version of the dataset with fewer columns
+    verbose : boolean, by default False
+        If True, prints detailed information about the download process
+
+    Returns
+    -------
+    gpd.GeoDataFrame
+        Geodataframe with planned eolic power plants data
+        
+    Example
+    -------
+    >>> from tunned_geobr import read_planned_eolic
+    
+    # Read planned eolic power plants data
+    >>> planned_eolic = read_planned_eolic()
+    """
+    
+    # URL for the EPE geoserver WFS service
+    url = r"https://gisepeprd2.epe.gov.br/arcgis/rest/services/Download_Dados_Webmap_EPE/GPServer/Extract%20Data%20Task/execute?f=json&env%3AoutSR=102100&Layers_to_Clip=%5B%22EOL%20-%20Expans%C3%A3o%20Planejada%22%5D&Area_of_Interest=%7B%22geometryType%22%3A%22esriGeometryPolygon%22%2C%22features%22%3A%5B%7B%22geometry%22%3A%7B%22rings%22%3A%5B%5B%5B-8202935.359686549%2C-4447637.219321614%5D%2C%5B-8202935.359686549%2C1119424.4247428626%5D%2C%5B-3487076.4626055676%2C1119424.4247428626%5D%2C%5B-3487076.4626055676%2C-4447637.219321614%5D%2C%5B-8202935.359686549%2C-4447637.219321614%5D%5D%5D%2C%22spatialReference%22%3A%7B%22wkid%22%3A102100%7D%7D%7D%5D%2C%22sr%22%3A%7B%22wkid%22%3A102100%7D%7D&Feature_Format=Shapefile%20-%20SHP%20-%20.shp&Raster_Format=Tagged%20Image%20File%20Format%20-%20TIFF%20-%20.tif"
+ 
+    try:
+        # Disable SSL verification warning
+        warnings.filterwarnings('ignore', message='Unverified HTTPS request')
+        
+        if verbose:
+            print("Requesting data from EPE server...")
+            
+        response = requests.get(url, timeout=60, verify=False)
+        if not response.ok:
+            raise Exception(f"Error getting JSON response: {response.status_code}")
+
+        json_response = response.json()
+        
+        if verbose:
+            print(f"JSON response received: {json.dumps(json_response, indent=2)[:500]}...")
+            
+        if 'results' not in json_response or len(json_response['results']) == 0:
+            raise Exception("Invalid JSON response structure")
+            
+        if 'value' not in json_response['results'][0] or 'url' not in json_response['results'][0]['value']:
+            raise Exception("URL not found in JSON response")
+            
+        file_url = json_response['results'][0]['value']['url']
+        
+        if verbose:
+            print(f"Downloading file from: {file_url}")
+        
+        file_response = requests.get(file_url, stream=True, timeout=60, verify=False)
+        if not file_response.ok:
+            raise Exception(f"Error downloading file: {file_response.status_code}")
+        
+        # Check if content is actually a zip file
+        content = file_response.content
+        if len(content) < 100:
+            if verbose:
+                print(f"Warning: Downloaded content is very small ({len(content)} bytes)")
+                print(f"Content preview: {content[:100]}")
+            
+        # Create a temporary directory to extract the files
+        with tempfile.TemporaryDirectory() as temp_dir:
+            if verbose:
+                print(f"Extracting files to temporary directory: {temp_dir}")
+                
+            try:
+                # Extract the zip file
+                with ZipFile(BytesIO(content)) as zip_ref:
+                    zip_ref.extractall(temp_dir)
+                    
+                    if verbose:
+                        print(f"Files in zip: {zip_ref.namelist()}")
+            except Exception as zip_error:
+                if verbose:
+                    print(f"Error extracting zip: {str(zip_error)}")
+                    print(f"Saving content to debug.zip for inspection")
+                    with open("debug.zip", "wb") as f:
+                        f.write(content)
+                raise Exception(f"Failed to extract zip file: {str(zip_error)}")
+            
+            # Find the shapefile
+            all_files = os.listdir(temp_dir)
+            if verbose:
+                print(f"Files in temp directory: {all_files}")
+                
+            shp_files = [f for f in all_files if f.endswith('.shp')]
+            if not shp_files:
+                # Try looking in subdirectories
+                for root, dirs, files in os.walk(temp_dir):
+                    shp_files.extend([os.path.join(root, f) for f in files if f.endswith('.shp')])
+                
+                if not shp_files:
+                    raise Exception("No shapefile found in the downloaded data")
+            
+            # Read the shapefile
+            shp_path = shp_files[0] if os.path.isabs(shp_files[0]) else os.path.join(temp_dir, shp_files[0])
+            if verbose:
+                print(f"Reading shapefile: {shp_path}")
+                
+            gdf = gpd.read_file(shp_path)
+            
+            # Convert to SIRGAS 2000 (EPSG:4674)
+            gdf = gdf.to_crs(4674)
+            
+            if verbose:
+                print(f"Data loaded successfully with {len(gdf)} records")
+                print(f"Columns: {gdf.columns.tolist()}")
+            
+            if simplified:
+                # Keep only the most relevant columns
+                columns_to_keep = [
+                    'geometry',
+                    'nome',            # Power plant name
+                    'potencia',        # Capacity in MW
+                    'leilao',          # Auction
+                    'ceg',             # CEG code
+                    'ano_prev'         # Expected year
+                ]
+                
+                # Filter columns that actually exist in the dataset
+                existing_columns = ['geometry'] + [col for col in columns_to_keep[1:] if col in gdf.columns]
+                if len(existing_columns) <= 1:
+                    if verbose:
+                        print("Warning: No matching columns found for simplified version. Returning all columns.")
+                else:
+                    gdf = gdf[existing_columns]
+    
+    except Exception as e:
+        raise Exception(f"Error downloading or processing planned eolic power plants data: {str(e)}")
+        
+    return gdf
+
+if __name__ == '__main__':
+    try:
+        eolic_data = read_planned_eolic(verbose=True)
+        print(f"Downloaded planned eolic power plants data with {len(eolic_data)} records and {len(eolic_data.columns)} columns")
+        
+        # Test simplified version
+        simplified_data = read_planned_eolic(simplified=True)
+        print(f"Simplified data has {len(simplified_data.columns)} columns: {simplified_data.columns.tolist()}")
+    except Exception as e:
+        print(f"Error: {str(e)}")
+        import traceback
+        traceback.print_exc()

tunned_geobr/read_planned_fossile_ute.py

@@ -0,0 +1,166 @@
+import geopandas as gpd
+import os
+import tempfile
+import urllib.parse
+import requests
+import shutil
+from zipfile import ZipFile
+from pathlib import Path
+from io import BytesIO
+import warnings
+import json
+
+def read_planned_fossile_ute(simplified=False, verbose=False):
+    """Download Planned Fossil Thermoelectric Power Plants data from EPE.
+    
+    This function downloads and processes planned fossil thermoelectric power plants data from EPE 
+    (Energy Research Company). The data includes information about planned fossil thermoelectric 
+    power generation projects across Brazil.
+    Original source: EPE (Empresa de Pesquisa Energética)
+    
+    Parameters
+    ----------
+    simplified : boolean, by default False
+        If True, returns a simplified version of the dataset with fewer columns
+    verbose : boolean, by default False
+        If True, prints detailed information about the download process
+
+    Returns
+    -------
+    gpd.GeoDataFrame
+        Geodataframe with planned fossil thermoelectric power plants data
+        
+    Example
+    -------
+    >>> from tunned_geobr import read_planned_fossile_ute
+    
+    # Read planned fossil thermoelectric power plants data
+    >>> planned_fossile_ute = read_planned_fossile_ute()
+    """
+    
+    # URL for the EPE geoserver
+    url = r"https://gisepeprd2.epe.gov.br/arcgis/rest/services/Download_Dados_Webmap_EPE/GPServer/Extract%20Data%20Task/execute?f=json&env%3AoutSR=102100&Layers_to_Clip=%5B%22UTE%20F%C3%B3ssil%20-%20Expans%C3%A3o%20Planejada%22%5D&Area_of_Interest=%7B%22geometryType%22%3A%22esriGeometryPolygon%22%2C%22features%22%3A%5B%7B%22geometry%22%3A%7B%22rings%22%3A%5B%5B%5B-8655251.47456396%2C-4787514.465591563%5D%2C%5B-8655251.47456396%2C1229608.401015912%5D%2C%5B-3508899.2341809804%2C1229608.401015912%5D%2C%5B-3508899.2341809804%2C-4787514.465591563%5D%2C%5B-8655251.47456396%2C-4787514.465591563%5D%5D%5D%2C%22spatialReference%22%3A%7B%22wkid%22%3A102100%7D%7D%7D%5D%2C%22sr%22%3A%7B%22wkid%22%3A102100%7D%7D&Feature_Format=Shapefile%20-%20SHP%20-%20.shp&Raster_Format=Tagged%20Image%20File%20Format%20-%20TIFF%20-%20.tif"
+ 
+    try:
+        # Disable SSL verification warning
+        warnings.filterwarnings('ignore', message='Unverified HTTPS request')
+        
+        if verbose:
+            print("Requesting data from EPE server...")
+            
+        response = requests.get(url, timeout=60, verify=False)
+        if not response.ok:
+            raise Exception(f"Error getting JSON response: {response.status_code}")
+
+        json_response = response.json()
+        
+        if verbose:
+            print(f"JSON response received: {json.dumps(json_response, indent=2)[:500]}...")
+            
+        if 'results' not in json_response or len(json_response['results']) == 0:
+            raise Exception("Invalid JSON response structure")
+            
+        if 'value' not in json_response['results'][0] or 'url' not in json_response['results'][0]['value']:
+            raise Exception("URL not found in JSON response")
+            
+        file_url = json_response['results'][0]['value']['url']
+        
+        if verbose:
+            print(f"Downloading file from: {file_url}")
+        
+        file_response = requests.get(file_url, stream=True, timeout=60, verify=False)
+        if not file_response.ok:
+            raise Exception(f"Error downloading file: {file_response.status_code}")
+        
+        # Check if content is actually a zip file
+        content = file_response.content
+        if len(content) < 100:
+            if verbose:
+                print(f"Warning: Downloaded content is very small ({len(content)} bytes)")
+                print(f"Content preview: {content[:100]}")
+            
+        # Create a temporary directory to extract the files
+        with tempfile.TemporaryDirectory() as temp_dir:
+            if verbose:
+                print(f"Extracting files to temporary directory: {temp_dir}")
+                
+            try:
+                # Extract the zip file
+                with ZipFile(BytesIO(content)) as zip_ref:
+                    zip_ref.extractall(temp_dir)
+                    
+                    if verbose:
+                        print(f"Files in zip: {zip_ref.namelist()}")
+            except Exception as zip_error:
+                if verbose:
+                    print(f"Error extracting zip: {str(zip_error)}")
+                    print(f"Saving content to debug.zip for inspection")
+                    with open("debug.zip", "wb") as f:
+                        f.write(content)
+                raise Exception(f"Failed to extract zip file: {str(zip_error)}")
+            
+            # Find the shapefile
+            all_files = os.listdir(temp_dir)
+            if verbose:
+                print(f"Files in temp directory: {all_files}")
+                
+            shp_files = [f for f in all_files if f.endswith('.shp')]
+            if not shp_files:
+                # Try looking in subdirectories
+                for root, dirs, files in os.walk(temp_dir):
+                    shp_files.extend([os.path.join(root, f) for f in files if f.endswith('.shp')])
+                
+                if not shp_files:
+                    raise Exception("No shapefile found in the downloaded data")
+            
+            # Read the shapefile
+            shp_path = shp_files[0] if os.path.isabs(shp_files[0]) else os.path.join(temp_dir, shp_files[0])
+            if verbose:
+                print(f"Reading shapefile: {shp_path}")
+                
+            gdf = gpd.read_file(shp_path)
+            
+            # Convert to SIRGAS 2000 (EPSG:4674)
+            gdf = gdf.to_crs(4674)
+            
+            if verbose:
+                print(f"Data loaded successfully with {len(gdf)} records")
+                print(f"Columns: {gdf.columns.tolist()}")
+            
+            if simplified:
+                # Keep only the most relevant columns
+                columns_to_keep = [
+                    'geometry',
+                    'nome',            # Power plant name
+                    'potencia',        # Capacity in MW
+                    'combustivel',     # Fuel type
+                    'leilao',          # Auction
+                    'ceg',             # CEG code
+                    'ano_prev'         # Expected year
+                ]
+                
+                # Filter columns that actually exist in the dataset
+                existing_columns = ['geometry'] + [col for col in columns_to_keep[1:] if col in gdf.columns]
+                if len(existing_columns) <= 1:
+                    if verbose:
+                        print("Warning: No matching columns found for simplified version. Returning all columns.")
+                else:
+                    gdf = gdf[existing_columns]
+    
+    except Exception as e:
+        raise Exception(f"Error downloading or processing planned fossil thermoelectric power plants data: {str(e)}")
+        
+    return gdf
+
+if __name__ == '__main__':
+    try:
+        fossile_ute_data = read_planned_fossile_ute(verbose=True)
+        print(f"Downloaded planned fossil thermoelectric power plants data with {len(fossile_ute_data)} records and {len(fossile_ute_data.columns)} columns")
+        
+        # Test simplified version
+        simplified_data = read_planned_fossile_ute(simplified=True)
+        print(f"Simplified data has {len(simplified_data.columns)} columns: {simplified_data.columns.tolist()}")
+    except Exception as e:
+        print(f"Error: {str(e)}")
+        import traceback
+        traceback.print_exc()