goesgcp 2.1.0__tar.gz → 3.0.1__tar.gz

{goesgcp-2.1.0/goesgcp.egg-info → goesgcp-3.0.1}/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: goesgcp
-Version: 2.1.0
+Version: 3.0.1
 Summary: A package to download and process GOES-16/17 data
 Home-page: https://github.com/helvecioneto/goesgcp
 Author: Helvecio B. L. Neto
@@ -19,11 +19,12 @@ Classifier: Topic :: Software Development
 Classifier: Topic :: Utilities
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: google-cloud-storage
-Requires-Dist: pyproj
+Requires-Dist: pandas
 Requires-Dist: xarray
+Requires-Dist: google-cloud-storage
 Requires-Dist: netcdf4
-Requires-Dist: rioxarray
+Requires-Dist: tqdm
+Requires-Dist: setuptools
 Dynamic: author
 Dynamic: author-email
 Dynamic: classifier
@@ -31,6 +32,7 @@ Dynamic: description
 Dynamic: description-content-type
 Dynamic: home-page
 Dynamic: license
+Dynamic: license-file
 Dynamic: requires-dist
 Dynamic: summary
 
@@ -62,6 +64,37 @@ pip install goesgcp
 ```
 
 
+Obs: If gdal is not installed, you can install it using the following command:
+
+Linux:
+```bash
+sudo apt-get install gdal-bin
+```
+
+Windows:
+```bash
+conda install -c conda-forge gdal
+```
+
+MacOS:
+```bash
+brew install gdal
+```
+
+Or you can install the wheel file:
+
+```bash
+python -m pip install gdal -f https://girder.github.io/large_image_wheels
+```
+
+and install other dependencies:
+
+```bash
+pip install -r requirements.txt
+```
+
+
+
 ## Usage
 
 ### Available Command-Line Arguments

{goesgcp-2.1.0 → goesgcp-3.0.1}/README.md

@@ -26,6 +26,37 @@ pip install goesgcp
 ```
 
 
+Obs: If gdal is not installed, you can install it using the following command:
+
+Linux:
+```bash
+sudo apt-get install gdal-bin
+```
+
+Windows:
+```bash
+conda install -c conda-forge gdal
+```
+
+MacOS:
+```bash
+brew install gdal
+```
+
+Or you can install the wheel file:
+
+```bash
+python -m pip install gdal -f https://girder.github.io/large_image_wheels
+```
+
+and install other dependencies:
+
+```bash
+pip install -r requirements.txt
+```
+
+
+
 ## Usage
 
 ### Available Command-Line Arguments

{goesgcp-2.1.0 → goesgcp-3.0.1}/goesgcp/main.py

@@ -13,7 +13,8 @@ from google.cloud import storage
 from datetime import datetime, timedelta, timezone
 from pyproj import CRS, Transformer
 from google.api_core.exceptions import GoogleAPIError
-
+import netCDF4
+import numpy as np
 import warnings
 warnings.filterwarnings('ignore')
 
@@ -34,7 +35,7 @@ def get_directory_prefix(year, julian_day, hour):
     return f"{year}/{julian_day}/{str(hour).zfill(2)}/"
 
 
-def get_files_period(connection, bucket_name, base_prefix, pattern, 
+def get_files_period(connection, bucket_name, base_prefix, pattern,
                      start, end, bt_hour=[], bt_min=[], freq=None):
     """
     Fetches files from a GCP bucket within a specified time period and returns them as a DataFrame.
@@ -87,7 +88,7 @@ def get_files_period(connection, bucket_name, base_prefix, pattern,
         print("No files found matching the pattern and time range.")
         print(prefix)
         sys.exit(1)
-    
+
     # Transform file_name to datetime
     df['last_modified'] = pd.to_datetime(df['file_name'].str.extract(r'(\d{4}\d{3}\d{2}\d{2})').squeeze(), format='%Y%j%H%M')
 
@@ -161,7 +162,24 @@ def crop_reproject(args):
     Crops and reprojects a GOES-16 file to EPSG:4326.
     """
 
-    file, output, var_name, lat_min, lat_max, lon_min, lon_max, resolution, save_format = args
+    file, output, var_name, lat_min, lat_max, lon_min, lon_max, resolution, save_format, \
+    more_info, file_pattern, classic_format, remap, method = args
+
+    if more_info:
+        # Open file using netCDF4
+        ds_s = xr.open_dataset(file, engine="netcdf4", decode_cf=False)
+        if var_name is None:
+            var_names = [var for var in ds_s.data_vars if len(ds_s[var].dims) == 2]
+            var_names = [var for var in var_names if 'DQF' not in var]
+        else:
+            var_names = [var_name]
+        scale_factors = [ds_s[var].attrs["scale_factor"] for var in var_names]
+        add_offsets = [ds_s[var].attrs["add_offset"] for var in var_names]
+        fill_values = [ds_s[var].attrs["_FillValue"] for var in var_names]
+        units = [ds_s[var].attrs["units"] for var in var_names]
+        sat_lat = ds_s["goes_imager_projection"].attrs["latitude_of_projection_origin"]
+        sat_lon = ds_s["goes_imager_projection"].attrs["longitude_of_projection_origin"]
+        ds_s.close()
 
     # Open the file
     ds = xr.open_dataset(file, engine="netcdf4")
@@ -169,8 +187,6 @@ def crop_reproject(args):
     if var_name is None:
         # Get all variables are 2D
         var_names = [var for var in ds.data_vars if len(ds[var].dims) == 2]
-
-        # Remove DQF variables
         var_names = [var for var in var_names if 'DQF' not in var]
     else:
         var_names = [var_name]
@@ -188,118 +204,151 @@ def crop_reproject(args):
     crs = CRS.from_cf(ds["goes_imager_projection"].attrs)
     ds = ds.rio.write_crs(crs)
 
-    # Try to reduce the size of the dataset
-    try:
-        # Create a transformer
-        transformer = Transformer.from_crs(CRS.from_epsg(4326), crs)
-        # Calculate the margin
-        margin_ratio = 0.40  # 40% margin
-
-        # Get the bounding box
-        min_x, min_y = transformer.transform(lat_min, lon_min)
-        max_x, max_y = transformer.transform(lat_max, lon_max)
-
-        # Calculate the range
-        x_range = abs(max_x - min_x)
-        y_range = abs(max_y - min_y)
-
-        margin_x = x_range * margin_ratio
-        margin_y = y_range * margin_ratio
-
-        # Expand the bounding box
-        min_x -= margin_x
-        max_x += margin_x
-        min_y -= margin_y
-        max_y += margin_y
-
-        # Select the region
-        if ds["y"].values[0] > ds["y"].values[-1]:  # Eixo y decrescente
-            ds_ = ds.sel(x=slice(min_x, max_x), y=slice(max_y, min_y))
-        else:  # Eixo y crescente
-            ds_ = ds.sel(x=slice(min_x, max_x), y=slice(min_y, max_y))
-        # Sort by y
-        if ds_["y"].values[0] > ds_["y"].values[-1]:
-            ds_ = ds_.sortby("y")
-        # Assign to ds
-        ds = ds_
-    except:
-        pass
+    # Create a transformer
+    transformer = Transformer.from_crs(CRS.from_epsg(4326), crs)
+    # Calculate the margin
+    margin_ratio = 0.40  # 40% margin
 
-    try:
-        # Reproject to EPSG:4326
-        ds = ds.rio.reproject("EPSG:4326", resolution=resolution)
+    # Get the bounding box
+    min_x, min_y = transformer.transform(lat_min, lon_min)
+    max_x, max_y = transformer.transform(lat_max, lon_max)
+
+    # Calculate the range
+    x_range = abs(max_x - min_x)
+    y_range = abs(max_y - min_y)
 
-        # Rename lat/lon coordinates
-        ds = ds.rename({"x": "lon", "y": "lat"})
+    margin_x = x_range * margin_ratio
+    margin_y = y_range * margin_ratio
 
-        # Add resolution to attributes
+    # Expand the bounding box
+    min_x -= margin_x
+    max_x += margin_x
+    min_y -= margin_y
+    max_y += margin_y
+
+    # Sort the values
+    y_min, y_max = sorted([min_y, max_y], reverse=ds["y"].values[0] > ds["y"].values[-1])
+
+    # Crop the dataset based on the bounding box
+    ds = ds.sel(x=slice(min_x, max_x), y=slice(y_min, y_max))
+
+    # Sort the values
+    if ds["y"].values[0] > ds["y"].values[-1]:
+        ds = ds.sortby("y")
+
+    # Reproject the dataset in serial and convert values to short
+    ds = ds.rio.reproject("EPSG:4326", resolution=resolution)
+
+    # Rename lat/lon coordinates
+    ds = ds.rename({"x": "lon", "y": "lat"})
+
+    # Check if remap is not a string
+    if type(remap) != str:
         for var in var_names:
             ds[var].attrs['resolution'] = "x={:.2f} y={:.2f} degree".format(resolution, resolution)
             ds[var].attrs['comments'] = 'Cropped and reprojected to EPSG:4326 by goesgcp'
-
-        # Crop using lat/lon coordinates, in parallel
+        # Crop using lat/lon coordinates
         ds = ds.rio.clip_box(minx=lon_min, miny=lat_min, maxx=lon_max, maxy=lat_max)
+    else:
+        # Add _FillValue to the variables
+        for var in var_names:
+            ds[var].attrs['_FillValue'] = float(fill_values[var_names.index(var)])
 
-        # Add global metadata comments
-        ds.attrs['comments'] = "Data processed by goesgcp, author: Helvecio B. L. Neto (helvecioblneto@gmail.com)"
-    except Exception as e:
-        print(f"Error processing file {file}: {e}")
-        pass    
-    
-    if save_format == 'by_date':
-        file_datetime = datetime.strptime(ds.time_coverage_start, 
+    # Add global metadata comments
+    ds.attrs['comments'] = "Data processed by goesgcp, author: Helvecio B. L. Neto (helvecioblneto@gmail.com)"
+
+    # Get the file datetime
+    file_datetime = datetime.strptime(ds.time_coverage_start,
                                           "%Y-%m-%dT%H:%M:%S.%fZ")
+    if save_format == 'by_date':
         year = file_datetime.strftime("%Y")
         month = file_datetime.strftime("%m")
         day = file_datetime.strftime("%d")
-        output_directory = f"{output}{year}/{month}/{day}/"
+        output_directory = f"{output}{year}/{month}/"
     elif save_format == 'julian':
-        file_datetime = datetime.strptime(ds.time_coverage_start, 
+        file_datetime = datetime.strptime(ds.time_coverage_start,
                                           "%Y-%m-%dT%H:%M:%S.%fZ")
         year = file_datetime.strftime("%Y")
         julian_day = file_datetime.timetuple().tm_yday
         output_directory = f"{output}{year}/{julian_day}/"
     else:
         output_directory = output
-        
 
     # Create the output directory
     pathlib.Path(output_directory).mkdir(parents=True, exist_ok=True)
 
-    # Save the file
-    output_file = f"{output_directory}{file.split('/')[-1]}"
-    ds.to_netcdf(output_file, mode='w')
+    # Apply file pattern
+    if file_pattern is not None:
+        # Get the timestamp
+        file_datetime = datetime.strptime(ds.time_coverage_start,
+                                          "%Y-%m-%dT%H:%M:%S.%fZ")
+        file_pattern = file_datetime.strftime(file_pattern)
+        output_file = f"{output_directory}{file_pattern}.nc"
+    else:
+        output_file = f"{output_directory}{file.split('/')[-1]}"
 
-    # Close the file
+    # Write the file
+    if classic_format:
+        # Change the data type to round values to int16
+        ds.to_netcdf(output_file, mode='w', format='NETCDF3_CLASSIC', encoding={var: {'dtype': np.int16} for var in var_names
+        })
+    else:
+        ds.to_netcdf(output_file, mode='w', encoding={var: {'zlib': True} for var in var_names})
     ds.close()
 
-    return output_file
-
+    # Remap the file
+    if remap:
+        remap_file((remap, output_file, output, method))
+
+    if more_info:
+        with netCDF4.Dataset(output_file, 'r+') as ds:
+            # Clear old attributes
+            for var_name in var_names:
+                var = ds.variables[var_name]
+                for attr in var.ncattrs():
+                    if attr == 'long_name' or attr == '_FillValue':
+                        continue
+                    var.delncattr(attr)
+            # Add new attributes
+            for var in range(len(var_names)):
+                ds[var_names[var]].setncattr('scale_factor', scale_factors[var])
+                ds[var_names[var]].setncattr('add_offset', add_offsets[var])
+                ds[var_names[var]].setncattr('missing_value', fill_values[var])
+                ds[var_names[var]].setncattr('units', np.float32(units[var]))
+            # Add variable satlat
+            ds.createDimension('satlat', 1)
+            ds.createVariable('satlat', 'f4', ('satlat',))
+            ds.variables['satlat'][:] = sat_lat
+            ds.variables['satlat'].long_name = 'Satellite Latitude'
+            ds.variables['satlat'].units = 'degrees_north'
+            ds.createDimension('satlon', 1)
+            ds.createVariable('satlon', 'f4', ('satlon',))
+            ds.variables['satlon'][:] = sat_lon
+            ds.variables['satlon'].long_name = 'Satellite Longitude'
+            ds.variables['satlon'].units = 'degrees_east'
+            ds.createDimension('julian_day', 1)
+            ds.createVariable('julian_day', 'i2', ('julian_day',))
+            ds.variables['julian_day'][:] = int(file_datetime.timetuple().tm_yday)
+            ds.variables['julian_day'].long_name = 'Julian day'
+            ds.variables['julian_day'].units = 'day'
+            time_of_day_char = netCDF4.stringtochar(np.array([str(file_datetime.strftime("%H%M"))], 'S4'))
+             # Add variable time_of_day
+            ds.createDimension('time_of_day', 4)
+            ds.createVariable('time_of_day', 'S1', ('time_of_day',))
+            ds.variables['time_of_day'][:] = time_of_day_char
+            ds.variables['time_of_day'].long_name = 'Time of day'
+            ds.variables['time_of_day'].units = 'hour and minute'
+            ds.variables['time_of_day'].comment = str(file_datetime.strftime("%H%M"))
 
 
 def remap_file(args):
     """ Remap the download file based on the input file. """
 
-    base_file, target_file, var_name, save_format, output = args
-
-    # Open the files
-    base_ds = xr.open_dataset(base_file, engine="netcdf4")
+    base_file, target_file, output, method = args
 
-    if save_format == 'by_date':
-        file_datetime = datetime.strptime(ds.time_coverage_start, 
-                                          "%Y-%m-%dT%H:%M:%S.%fZ")
-        year = file_datetime.strftime("%Y")
-        month = file_datetime.strftime("%m")
-        day = file_datetime.strftime("%d")
-        output_directory = f"{output}{year}/{month}/{day}/"
-    elif save_format == 'julian':
-        file_datetime = datetime.strptime(ds.time_coverage_start, 
-                                          "%Y-%m-%dT%H:%M:%S.%fZ")
-        year = file_datetime.strftime("%Y")
-        julian_day = file_datetime.timetuple().tm_yday
-        output_directory = f"{output}{year}/{julian_day}/"
-    else:
-        output_directory = output
+    # Get output directory based on target_file
+    output_file = f"{output}{target_file.split('/')[-1]}"
+    output_directory = output_file.replace(target_file.split('/')[-1], "")
 
     # Create the output directory
     pathlib.Path(output_directory).mkdir(parents=True, exist_ok=True)
@@ -311,7 +360,7 @@ def remap_file(args):
 
     # Run the cdo command
     cdo_command = [
-        "cdo", "remapnn," + base_file, target_file, output_file
+        "cdo", method+"," + base_file, target_file, output_file
     ]
 
     try:
@@ -320,9 +369,6 @@ def remap_file(args):
         print(f"Error remapping file {target_file}: {e}")
         pass
 
-    # Close the files
-    base_ds.close()
-
     # Delete the target file
     pathlib.Path(target_file).unlink()
 
@@ -334,18 +380,17 @@ def process_file(args):
     """
     Downloads and processes a GOES-16 file.
     """
-    
+
     bucket_name, blob_name, local_path, output_path, var_name, lat_min, lat_max, lon_min, lon_max, resolution, \
-    save_format, retries, remap = args
+    save_format, retries, remap, met, more_info, file_pattern, classic_format = args
 
+    # Download the file
     attempt = 0
     while attempt < retries:
         try:
             # Connect to the bucket
             bucket = storage_client.bucket(bucket_name)
             blob = bucket.blob(blob_name)
-
-            # Download the file
             blob.download_to_filename(local_path, timeout=120)
             break  # Exit the loop if the download is successful
         except (GoogleAPIError, Exception) as e:  # Catch any exception
@@ -353,21 +398,20 @@ def process_file(args):
             if attempt < retries:
                 time.sleep(2 ** attempt)  # Backoff exponencial
             else:
-                # Log the error to a file
                 with open('fail.log', 'a') as log_file:
                     log_file.write(f"Failed to download {blob_name} after {retries} attempts. Error: {e}\n")
-
     # Crop the file
-    output_file = crop_reproject((local_path, output_path, var_name,
-                                   lat_min, lat_max, lon_min, lon_max, resolution, save_format))
-
-    # Remap the file
-    if remap is not None:
-        # Remap the file
-        remap_file((remap, output_file, var_name, save_format, output_path))
-
-    # Remove the local file
-    pathlib.Path(local_path).unlink()
+    try:
+        crop_reproject((local_path, output_path,
+                    var_name, lat_min, lat_max, lon_min, lon_max,
+                    resolution, save_format,
+                    more_info, file_pattern, classic_format, remap, met))
+        # Remove the local file
+        pathlib.Path(local_path).unlink()
+    except Exception as e:
+        with open('fail.log', 'a') as log_file:
+            log_file.write(f"Failed to process {blob_name}. Error: {e}\n")
+        pass
 
 # Create connection
 storage_client = storage.Client.create_anonymous_client()
@@ -377,7 +421,7 @@ def main():
 
     epilog = """
     Example usage:
-    
+
     - To download recent 3 files from the GOES-16 satellite for the ABI-L2-CMIPF product,
     change resolution to 0.045, and crop the files between latitudes -35 and 5 and longitudes -80 and -30:
 
@@ -411,9 +455,9 @@ def main():
     parser = argparse.ArgumentParser(description='Download and process GOES Satellite data files from GCP.',
                                     epilog=epilog,
                                     formatter_class=argparse.RawDescriptionHelpFormatter)
-    
+
     # Satellite and product settings
-    parser.add_argument('--satellite', type=str, default='goes-16', choices=['goes-16', 'goes-18'], help='Name of the satellite (e.g., goes16)')
+    parser.add_argument('--satellite', type=str, default='goes-19', choices=['goes-16', 'goes-18', 'goes-19'], help='Name of the satellite (e.g., goes16)')
     parser.add_argument('--product', type=str, default='ABI-L2-CMIPF', help='Name of the satellite product', choices=product_names)
     parser.add_argument('--var_name', type=str, default=None, help='Variable name to extract (e.g., CMI)')
     parser.add_argument('--channel', type=int, default=13, help='Channel to use (e.g., 13)')
@@ -430,24 +474,26 @@ def main():
     parser.add_argument('--bt_min', nargs=2, type=int, default=[0, 60], help='Filter data between these minutes (e.g., 0 60)')
 
     # Geographic bounding box
-    parser.add_argument('--lat_min', type=float, default=-81.3282, help='Minimum latitude of the bounding box')
-    parser.add_argument('--lat_max', type=float, default=81.3282, help='Maximum latitude of the bounding box')
-    parser.add_argument('--lon_min', type=float, default=-156.2995, help='Minimum longitude of the bounding box')
-    parser.add_argument('--lon_max', type=float, default=6.2995, help='Maximum longitude of the bounding box')
-    parser.add_argument('--resolution', type=float, default=0.03208, help='Resolution of the output file')
-    parser.add_argument('--output', type=str, default='output/', help='Path for saving output files')
+    parser.add_argument('--lat_min', type=float, default=-56, help='Minimum latitude of the bounding box')
+    parser.add_argument('--lat_max', type=float, default=35, help='Maximum latitude of the bounding box')
+    parser.add_argument('--lon_min', type=float, default=-116, help='Minimum longitude of the bounding box')
+    parser.add_argument('--lon_max', type=float, default=-25, help='Maximum longitude of the bounding box')
+    parser.add_argument('--resolution', type=float, default=0.03, help='Resolution of the output file')
+    parser.add_argument('--output', type=str, default='./output/', help='Path for saving output files')
 
     # Remap
     parser.add_argument('--remap', type=str, default=None, help='Give a input file to remap the output')
+    parser.add_argument('--method', type=str, default='remapnn', help='Remap method to use (e.g., remapnn)')
 
     # Other settings
     parser.add_argument('--parallel', type=lambda x: bool(strtobool(x)), default=True, help='Use parallel processing')
     parser.add_argument('--processes', type=int, default=4, help='Number of processes for parallel execution')
     parser.add_argument('--max_attempts', type=int, default=3, help='Number of attempts to download a file')
+    parser.add_argument('--info', type=lambda x: bool(strtobool(x)), default=False, help='Show information messages')
     parser.add_argument('--save_format', type=str, default='flat', choices=['flat', 'by_date','julian'],
                     help="Save the files in a flat structure or by date")
-
-
+    parser.add_argument('--file_pattern', type=str, default=None, help='Pattern for the files')
+    parser.add_argument('--netcdf_classic', type=lambda x: bool(strtobool(x)), default=False, help='Save the files in netCDF classic format')
     # Parse arguments
     args = parser.parse_args()
 
@@ -477,7 +523,10 @@ def main():
     bt_min = args.bt_min
     save_format = args.save_format
     remap = args.remap
-
+    method = args.method
+    more_info = args.info
+    file_pattern = args.file_pattern
+    classic_format = args.netcdf_classic
 
     # Check mandatory arguments
     if not args.recent and not (args.start and args.end):
@@ -510,8 +559,8 @@ def main():
     # Check operational mode if is recent or specific date
     if start and end:
         files_list = get_files_period(storage_client, bucket_name,
-                                       product, pattern, start, end,
-                                        bt_hour, bt_min, freq)
+                                    product, pattern, start, end,
+                                    bt_hour, bt_min, freq)
     else:
         # Get recent files
         files_list = get_recent_files(storage_client, bucket_name, product, pattern, recent)
@@ -529,12 +578,13 @@ def main():
     loading_bar = tqdm.tqdm(total=len(files_list), ncols=100, position=0, leave=True,
                         bar_format='{l_bar}{bar}| {n_fmt}/{total_fmt} + \
                         [Elapsed:{elapsed} Remaining:<{remaining}]')
-    
+
     if parallel: # Run in parallel
         # Create a list of tasks
-        tasks = [(bucket_name, file, f"tmp/{file.split('/')[-1]}", output_path, var_name, 
+        tasks = [(bucket_name, file, f"tmp/{file.split('/')[-1]}", output_path, var_name,
         lat_min, lat_max, lon_min, lon_max, resolution,
-        save_format, max_attempts, remap) for file in files_list]
+        save_format, max_attempts, remap, method,
+        more_info, file_pattern, classic_format) for file in files_list]
 
         # Download files in parallel
         with Pool(processes=args.processes) as pool:
@@ -546,10 +596,12 @@ def main():
             local_path = f"tmp/{file.split('/')[-1]}"
             process_file((bucket_name, file, local_path, output_path, var_name,
             lat_min, lat_max, lon_min, lon_max, resolution,
-            save_format, max_attempts, remap))
+            save_format, max_attempts, remap, method, more_info,
+              file_pattern, classic_format))
             loading_bar.update(1)
         loading_bar.close()
 
+    # Clean up the temporary directory
     shutil.rmtree('tmp/')
 
 if __name__ == '__main__':

{goesgcp-2.1.0 → goesgcp-3.0.1/goesgcp.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: goesgcp
-Version: 2.1.0
+Version: 3.0.1
 Summary: A package to download and process GOES-16/17 data
 Home-page: https://github.com/helvecioneto/goesgcp
 Author: Helvecio B. L. Neto
@@ -19,11 +19,12 @@ Classifier: Topic :: Software Development
 Classifier: Topic :: Utilities
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: google-cloud-storage
-Requires-Dist: pyproj
+Requires-Dist: pandas
 Requires-Dist: xarray
+Requires-Dist: google-cloud-storage
 Requires-Dist: netcdf4
-Requires-Dist: rioxarray
+Requires-Dist: tqdm
+Requires-Dist: setuptools
 Dynamic: author
 Dynamic: author-email
 Dynamic: classifier
@@ -31,6 +32,7 @@ Dynamic: description
 Dynamic: description-content-type
 Dynamic: home-page
 Dynamic: license
+Dynamic: license-file
 Dynamic: requires-dist
 Dynamic: summary
 
@@ -62,6 +64,37 @@ pip install goesgcp
 ```
 
 
+Obs: If gdal is not installed, you can install it using the following command:
+
+Linux:
+```bash
+sudo apt-get install gdal-bin
+```
+
+Windows:
+```bash
+conda install -c conda-forge gdal
+```
+
+MacOS:
+```bash
+brew install gdal
+```
+
+Or you can install the wheel file:
+
+```bash
+python -m pip install gdal -f https://girder.github.io/large_image_wheels
+```
+
+and install other dependencies:
+
+```bash
+pip install -r requirements.txt
+```
+
+
+
 ## Usage
 
 ### Available Command-Line Arguments

{goesgcp-2.1.0 → goesgcp-3.0.1}/goesgcp.egg-info/requires.txt

@@ -1,5 +1,6 @@
-google-cloud-storage
-pyproj
+pandas
 xarray
+google-cloud-storage
 netcdf4
-rioxarray
+tqdm
+setuptools

{goesgcp-2.1.0 → goesgcp-3.0.1}/requirements.txt

@@ -1,5 +1,6 @@
-google-cloud-storage
-pyproj
+pandas
 xarray
+google-cloud-storage
 netcdf4
-rioxarray
+tqdm
+setuptools

goesgcp-3.0.1/setup.py

@@ -0,0 +1,101 @@
+import subprocess
+
+# Upgrade pip and setuptools before installing dependencies
+try:
+    subprocess.run(["python", "-m", "ensurepip", "--upgrade"], check=True)
+    subprocess.run(["python", "-m", "pip", "install", "--upgrade", "pip"], check=True)
+    subprocess.run(["python", "-m", "pip", "install", "--upgrade", "setuptools"], check=True)
+except subprocess.CalledProcessError:
+    print("Error upgrading pip and setuptools.")
+
+import os
+import sys
+import platform
+from setuptools import setup, find_packages
+
+
+def install_gdal():
+    """Installs GDAL based on the operating system."""
+    system = platform.system()
+
+    try:
+        if system == "Linux":
+            print("Installing GDAL on Linux...")
+            subprocess.run(
+                ["python", "-m", "pip", "install", "gdal", "-f", "https://girder.github.io/large_image_wheels"],
+                check=True
+            )
+
+        elif system == "Windows":
+            python_version = sys.version_info
+            if python_version[0] == 3:
+                version_map = {
+                    10: "cp310",
+                    11: "cp311",
+                    12: "cp312",
+                    13: "cp313"
+                }
+                if python_version[1] in version_map:
+                    gdal_url = f"https://github.com/cgohlke/geospatial-wheels/releases/download/v2025.1.20/GDAL-3.10.1-{version_map[python_version[1]]}-win_amd64.whl"
+                    print(f"Downloading GDAL for Python {python_version[0]}.{python_version[1]}...")
+                    subprocess.run(["python", "-m", "pip", "install", gdal_url], check=True)
+                else:
+                    print("Unsupported Python version for GDAL.")
+                    sys.exit(1)
+            else:
+                print("Unsupported Python version for GDAL.")
+                sys.exit(1)
+
+        elif system == "Darwin":
+            print("GDAL is not available for macOS via this installer.")
+            sys.exit(1)
+
+        else:
+            print("Unsupported operating system.")
+            sys.exit(1)
+    
+    except subprocess.CalledProcessError:
+        print("Error installing GDAL.")
+
+# Install GDAL before running setup
+install_gdal()
+
+# Read dependencies from requirements.txt
+req_file = os.path.join(os.path.dirname(__file__), "requirements.txt")
+if os.path.exists(req_file):
+    with open(req_file) as f:
+        requirements = f.read().splitlines()
+else:
+    requirements = []
+
+setup(
+    name="goesgcp",
+    version="3.0.1",
+    author="Helvecio B. L. Neto",
+    author_email="helvecioblneto@gmail.com",
+    description="A package to download and process GOES-16/17 data",
+    long_description=open("README.md").read(),
+    long_description_content_type="text/markdown",
+    url="https://github.com/helvecioneto/goesgcp",
+    packages=find_packages(),
+    install_requires=requirements,
+    license="LICENSE",
+    classifiers=[
+        "Programming Language :: Python",
+        "Development Status :: 5 - Production/Stable",
+        "Operating System :: OS Independent",
+        "Programming Language :: Python :: 3.10",
+        "Programming Language :: Python :: 3.11",
+        "Programming Language :: Python :: 3.12",
+        "Topic :: Scientific/Engineering :: Atmospheric Science",
+        "Topic :: Scientific/Engineering :: GIS",
+        "Topic :: Scientific/Engineering",
+        "Topic :: Software Development",
+        "Topic :: Utilities",
+    ],
+    entry_points={
+        'console_scripts': [
+            'goesgcp=goesgcp.main:main',
+        ],
+    },
+)

goesgcp-2.1.0/setup.py

@@ -1,49 +0,0 @@
-
-import os
-
-try:
-    os.system("python -m ensurepip --upgrade")
-    os.system("python -m pip install --upgrade setuptools")
-except:
-    pass
-
-from setuptools import setup, find_packages
-
-req_file = os.path.join(os.path.dirname(__file__), "requirements.txt")
-if os.path.exists(req_file):
-    with open(req_file) as f:
-        requirements = f.read().splitlines()
-else:
-    requirements = []
-
-setup(
-    name="goesgcp",
-    version="2.1.0",
-    author="Helvecio B. L. Neto",
-    author_email="helvecioblneto@gmail.com",
-    description="A package to download and process GOES-16/17 data",
-    long_description=open("README.md").read(),
-    long_description_content_type="text/markdown",
-    url="https://github.com/helvecioneto/goesgcp",
-    packages=find_packages(),
-    install_requires=requirements,
-    license="LICENSE",
-    classifiers=[
-        "Programming Language :: Python",
-        "Development Status :: 5 - Production/Stable",
-        "Operating System :: OS Independent",
-        "Programming Language :: Python :: 3.10",
-        "Programming Language :: Python :: 3.11",
-        "Programming Language :: Python :: 3.12",
-	    "Topic :: Scientific/Engineering :: Atmospheric Science",
-        "Topic :: Scientific/Engineering :: GIS",
-        "Topic :: Scientific/Engineering",
-        "Topic :: Software Development",
-        "Topic :: Utilities",
-    ],
-    entry_points={
-        'console_scripts': [
-            'goesgcp=goesgcp.main:main',
-        ],
-    },
-)