Skip to content

Instantly share code, notes, and snippets.

View AdrianoPereira's full-sized avatar
🎯
Focusing

Adriano Almeida AdrianoPereira

🎯
Focusing
View GitHub Profile
@AdrianoPereira
AdrianoPereira / logo.md
Created August 25, 2026 00:04
Solar wind forecast logo
aquiry-forecast-logo
@AdrianoPereira
AdrianoPereira / download_imerg.py
Created July 24, 2025 20:13
Download IMERG product, just change DOI
import os
import earthaccess
# from dotenv import load_dotenv
import argparse
# load_dotenv()
os.environ['EARTHDATA_USERNAME'] = '<your-user-here>'
os.environ['EARTHDATA_PASSWORD'] = '<your-password-here>'
os.environ['OUTPUT_DIR'] = 'path/to/save'
@AdrianoPereira
AdrianoPereira / download_gsmap_data.py
Created July 24, 2025 20:02
AINPP products: MVK - /standard/v8/hourly, NOW - /now/latest, NRT - /realtime_ver/v8/latest
import ftplib
import os
from datetime import datetime, timedelta
import requests
import pathlib
import tempfile
import gzip
import numpy as np
from dotenv import load_dotenv
load_dotenv()
@AdrianoPereira
AdrianoPereira / Extract_GSMaP_Region_Data.py
Created May 13, 2025 16:08
Single process code for load, extract and save cropped data from GSMaP from global to specigic region bound box
import os
import numpy as np
import gzip
from pathlib import Path
input_base_dir = "/storage/hokusai/GSMaP_NOW/now/half_hour/2024"
output_base_dir = "/storage/hokusai/GSMaP_NOW/SA_2024"
# global coordinates
glat_min = -60.0
@AdrianoPereira
AdrianoPereira / Parallel_Extract_GSMaP_Region_Data.py
Created May 13, 2025 16:07
Parallel code for load, extract and save cropped data from GSMaP from global to specigic region bound box
import os
import numpy as np
import gzip
from pathlib import Path
from concurrent.futures import ProcessPoolExecutor, as_completed
from tqdm import tqdm
input_base_dir = "/storage/hokusai/GSMaP_NOW/now/half_hour/2024"
output_base_dir = "/storage/hokusai/GSMaP_NOW/SA_2024"
import ee
ee.Initialize()
# Coleção de departamentos do Peru
regions = ee.FeatureCollection("FAO/GAUL/2015/level1").filter(
ee.Filter.eq("ADM0_NAME", "Peru")
)
# Imagem SRTM de elevação
srtm = ee.Image("USGS/SRTMGL1_003")
# mosaic
import rasterio
from rasterio.merge import merge
from pathlib import Path
import glob
# Caminho para a pasta onde estão os GeoTIFFs baixados
input_dir = Path("/home/adriano/Doutorado/landcover/EarthEngine_Exports/") # ajuste o nome da pasta se necessário
import ee
ee.Initialize()
# Peru information
regions = ee.FeatureCollection("FAO/GAUL/2015/level1").filter(
ee.Filter.eq("ADM0_NAME", "Peru")
)
# Landcover product
landcover = ee.Image("COPERNICUS/Landcover/100m/Proba-V-C3/Global/2019")
@AdrianoPereira
AdrianoPereira / pyForTraCC_GOES13_Track.py
Created April 25, 2025 13:59
Python script for tracking GOES-13 Infrared (B04) data in a specific interest area.
import gzip
import pyfortracc
import numpy as np
import xarray as xr
scale = 0.01
lon_min, lon_max = -100.0, -25.24
lat_min, lat_max = -56.0, 12.52
res_x, res_y = 0.04, 0.04
@AdrianoPereira
AdrianoPereira / STAC_GOES13_Download.py
Created April 24, 2025 17:14
Script for download GOES13 data from STAC
import pystac_client
from tqdm import tqdm
import os
from datetime import datetime
import requests
STAC_API_URL = 'https://data.inpe.br/bdc/stac/v1/'
BASE_DIR = '/storage/goes/goes13'
CHANNEL = 'B04'