The exercise shows how a simple SketchUp 3D volume, defined solely by its basic geometry, can be transformed into a complex architectural proposal. Starting from the initial schematic model, the system interprets proportions, levels, and shapes, and converts them into a fully developed building, complete with textures, vegetation, lighting, and an urban context.
All IA prompts, step by step
It is a good test to evaluate the extent to which AI understands spatial structure and is capable of generating a coherent and rich interpretation from a minimal geometric skeleton.
Et voilá, the final result
Adding a bit of perspective to this post, myself I have always loved 3D and I have worked over it in the past using it from different approaches, starting from architectural modelling to massive QC of LOD1 buildings anywhere in the world. That’s why I see these Gemini features as something kind on interesting and related to my world.
Some old 3D models sold to Telespazio back in 2008 🙂
Dont get wrong if you see the IA background showing our handsome major showing his beautiful smile in Cibeles/Correos, it’s only to get your precious attention (only if you need it thou!). Flooding in urban environments is not a speculative hazard but something we can quantify. In the case of Madrid, the intersection of pretty mountainous terrain and urban expansion presents a scenario of significant risk, particularly when analyzed through the lens of shared high-resolution geospatial data (it might surprise you there are 2000m difference between the highest spot in Madrid province, Pico Peñalara -2428m- and the Alberche river environment in some areas -430m-).
Major Almeida is surprised to get to know how risky might be a heavy rain episode over Madrid area
This study integrates the buildings from BTN (Base Topográfica Nacional) provided by the Spanish “IGN”, the CNIG with the official flood hazard maps for a 100-year return period (T=100), published by the Ministry for the Ecological Transition and the Demographic Challenge (MITECO). The T=100 scenario is the most representative for evaluating long-term flood exposure, as it reflects events with a 1% annual probability—rare but not improbable, and certainly not negligible.
If we measure preliminarily, only 0.999% (8097/810134) of all buildings (available in our latest BTN buildings provided by CNIG) would be affected but but if we go deeper, this means a lot in specific spots: Aranjuez for instance would be very affected by the Tajo river flood.
Focused over Aranjuez area as I find it representative and we happen to find lots of buildings in likely flood areas (check it yourself in GEE)
The spatial overlay of flood-prone zones with demographic and land-use data reveals a concerning concentration of residential population within areas designated as ARPSIs (Áreas de Riesgo Potencial Significativo de Inundación). These zones include dense urban districts along the Manzanares River and low-lying areas in Arganzuela, Usera, and parts of Puente de Vallecas.
Histogram of those 8097 buildings (according to BTN/CNIG database of all buildings over Madrid)
While precise population figures vary depending on the granularity of census data, preliminary estimates suggest that tens of thousands of residents could be directly affected by a flood event of this magnitude. The implications extend beyond displacement and property damage, encompassing public health risks, disruption of essential services, and long-term socioeconomic instability.
Global Mapper was used to geoprocess and visualize geodata
Economic activities within these flood zones are diverse and structurally significant. Central districts host a high concentration of retail, hospitality, and cultural institutions, while peripheral zones near the M-30 corridor accommodate logistics, warehousing, and industrial operations. The exposure of these sectors to flood risk implies not only direct financial losses but also cascading effects on employment, supply chains, and urban mobility. Moreover, the presence of public infrastructure—transport nodes, administrative buildings, and emergency services—within these vulnerable areas raises questions about the resilience of the city’s operational backbone.
*Retail and hospitality in central districts (e.g., Lavapiés, La Latina) *Logistics and warehousing near the M-30 corridor *Cultural and tourism assets, including museums and heritage sites *Public infrastructure, such as metro stations, bus depots, and administrative buildings
Disruption in these zones could result in multi-million euro losses, not only from direct damage but also from prolonged service interruptions.
I have uploaded a small sample over Aranjuez to my GEE interface (is not that I am from this place, I only find it representative!). You can quickly tune the JavaScript code to be able to count the amount of buildings affected, the area of all of them and the DTM range of the flood DTM REM (Relative Elevation Model) over your actual screen. https://code.earthengine.google.com/1d7e907283e33dc4b468ab3adf578840
Cloud computation calculation in GEE. Flood Risk in Madrid. Check it yourself using the code below!!!
Of particular concern is the identification of facilities regulated under Annex I of Directive 96/61/EC, which governs integrated pollution prevention and control.
ArcGIS Pro was also used for geoprocessing and visualization
*Fuel storage and distribution centers *Waste treatment plants *Industrial facilities with hazardous materials
These include fuel depots, waste treatment centers, and industrial sites handling hazardous substances. In the event of flooding, such facilities pose a risk of accidental contamination, with potential impacts on protected zones defined in Annex IV of Directive 2000/60/EC. These zones include drinking water abstraction points, Natura 2000 habitats, and recreational waters. The spatial proximity of these sensitive areas to flood-prone industrial sites underscores the need for integrated risk assessment that goes beyond hydrological modeling and incorporates environmental and public health dimensions.
(i) Drinking water abstraction points (iii) Habitats designated under Natura 2000 (v) Recreational waters and sensitive ecosystems
The visual simulations accompanying this analysis are intentionally exaggerated. They do not represent predictive models but serve as heuristic devices to provoke reflection and debate. By depicting iconic Madrid landmarks submerged under chaotic floodwaters, these images challenge the viewer to confront the consequences of urban planning decisions that disregard hydrological constraints. They are not intended to alarm, but to illustrate the scale of disruption that could result from a statistically plausible event.
Please note these features have been taken out of the scope of this preliminary approach for analysis Flood in Madrid province using Open Data.
Again, this study/first approach would not have been possible without access to open geospatial data. The availability of national datasets such as the Base Topográfica Nacional and flood hazard maps from MITECO exemplifies the transformative potential of public data infrastructures. However, the mere existence of data is insufficient. What is required is a culture of proactive use—by planners, policymakers, and civil society—where risk is treated not as an abstract probability but as a concrete design constraint.
The presence of vulnerable populations, critical infrastructure, and environmentally sensitive facilities within flood-prone zones constitutes a non-assumable risk. It is a risk that could be mitigated through better zoning, stricter regulation, and investment in adaptive infrastructure. The cost of inaction is not only economic but ethical. Avoidable losses—whether of property, livelihoods, or ecosystems—are a failure of foresight, not fate.
It is a fact that a betting shop (AKA bookmakerin UK or bookiein the US) should not be close to a secondary school. Its obvious the impact on population ranging something like 12-18 could be higher than in other age thresholds. How close? 100m? 300m? 500m? Euclidean distance (a straight line) or following the street network?. In any case, if I decide choosing for instance a range of 500m, for example, 81% of betting shops in Madrid have secondary schools within that distance (258 out of 316). Looking at it from the secondary schools’ point of view, almost 60% of secondary schools have betting shops within 500m (171/291). This is undoubtedly an issue that needs to be addressed.
The spatial analysis reveals a moderately negative correlation (r = -0.39) between the distance to the nearest betting shop and the number of betting shops within a 500-meter radius of each secondary school. This means that, on average, the shorter the distance between a school and the nearest betting house, the greater the number of betting houses found within a short walking distance (for example, within 500 meters). In other words, schools that are close to one betting house are very often close to several. Conversely, schools located farther from any betting house tend to be in areas where these establishments are much less common or even absent. This pattern does not imply a perfect one-to-one relationship — some exceptions exist — but the overall trend is clear: betting houses are not randomly distributed in the city. Instead, they tend to form clusters, and those clusters often appear in the same parts of the city where many schools are located. Although the relationship is not perfectly linear, it might indicate a potential spatial association between the two types of locations.
Source: 2playbook.com
From the heat map, this relationship becomes more tangible. 171 out of 291 secondary schools in Madrid — approximately 60% of the total — are situated within 500 meters of at least one betting shop. The density surface highlights four clear hotspot zones: Usera, Carabanchel, Centro, and Tetuán. These districts concentrate the majority of the spatial overlap between High Schools and betting shops , forming well-defined high-intensity clusters. The pattern aligns with socio-economic and demographic realities: these are traditionally dense urban districts with younger populations and, in several cases, lower average income levels, which have historically attracted a higher concentration of gambling venues.
Concentration of Problem zones (threshold 500m) Usera, Carabanchel, Centro, and Tetuán
By contrast, the richest districts — such as Salamanca, Retiro or Chamberí — show a more dispersed and lower-intensity pattern. This uneven distribution reinforces the idea that the proximity of betting shops to High School educational institutions is not random, but rather follows a spatial logic influenced by the urban and social fabric of the city.
In practical terms, the combination of a negative correlation and the hotspot clustering suggests that High Schools and betting shops in Madrid exhibit a clear spatial association. The proximity of many High Schools to gambling establishments may have implications for urban planning, youth protection policies, and spatial regulation of gambling activities. The results provide a quantitative foundation for further policy discussions on how the city’s gambling landscape interacts with its educational network.
Pearson correlation. Number of betting houses vs High Schools. Pearson r: -0.39
Methodology: Open Data Madrid provided ‘Censo de locales, sus actividades y terrazas de hostelería y restauración. Histórico’ which I cleaned and cross-referenced to obtain the venues dedicated to gambling and having an active licence (2020). Then I downloaded Madrid High Schools (those offering public and private secondary education). Based on Madrid districts I was able to measure distances and highlight the most affected areas. I used GEE to performed the Pearson Correlation (you have the code below). Maps where created using Global Mapper v26.2.
GEE interface. Cloud Computation calculation
What if I can get a GeoJSON just by feeding the system with the required threshold? (i.e this code returns the 171 secondary schools in Madrid which have betting houses closer than a 500m radius from them). It’s crazy 171/291!! (almost 60%)?: https://code.earthengine.google.com/e15e0ee2c0269d8b3e4a199d444620d3
31% of bookies are closer than 100m to another bookie
TRY IT YOURSELF!!!
Additional considerations:
The regulations governing the distance between betting shops and educational establishments such as secondary schools vary depending on the autonomous community. Some communities have implemented or proposed minimum distances, such as a 300-metre radius in some regional laws, while others continue to apply shorter distances or are in the process of reviewing their regulations. It is essential to consult the specific legislation of each community to find out the exact regulations and how they are applied. Examples of regulations by autonomous community (CCAA):
*Community of Madrid: Established a minimum distance of 100 metres between gaming halls and schools/secondary schools in 2019. Although extensions were proposed, the current regulation remains at 100 metres, with the possibility of some existing premises being exempt, according to El Mundo. *Galicia: Its gaming law establishes a minimum distance of 300 metres from educational centres for new gaming establishments, increasing the previous linear distance of 150 metres, according to GaliciaPress. *Andalusia: The proposal to extend the distance to 500 metres has been taken to Parliament, although the current regulation may be different, as reported in ABC. *Other regions: Some cities, such as Talavera de la Reina in Castilla-La Mancha, have included a minimum distance of 300 metres, backed by court rulings, according to different sources.
Distance can be measured in a straight line or in radius, and this difference is important when applying the regulations. Premises that were already in operation before the regulations were implemented in some regions are often exempt. In the case of Madrid, proposals to ban betting shops from being located less than 500 metres from educational establishments have even been analysed, as reported by Telemadrid.
High‑resolution elevation data underpins almost every spatial analysis we do in GIS—especially in forests where vertical structure defines habitat, biomass, wind exposure, fire behavior, hydrology, and the microclimates that sustain rare species. In rugged or densely vegetated environments, a coarse or biased elevation model propagates error everywhere: orthorectification drifts, hillshades mislead, slope/aspect misclassify, and canopy metrics saturate. The result is decisions made on blurred terrain that hides the very patterns we seek to manage. Precision elevation—derived from airborne LiDAR (Light Detection and Ranging)—solves this by separating the ground from the vegetation and delivering both a bare‑earth Digital Terrain Model (DTM) and a Digital Surface Model (DSM). Subtracting DTM from DSM gives a Canopy Height Model (DHM) that captures the true vertical architecture of the forest at sub‑meter resolution.
My visit to the “premises” late August 2025 (almost yesterday!)
This post uses the Mata do Buçaco (Bussaco), near Luso in central Portugal, to illustrate why precision matters and how it compares to the widely used global product ETH_GlobalCanopyHeight_2020_10m_v1. We will look at the site, the LiDAR technology, and a practical comparison workflow for GIS users.
Surprise! a RMSE deviation result too high! Why?
Mata do Buçaco: a compact sanctuary of forest giants
Mata do Buçaco is a walled arboretum and national forest just above the spa town of Luso, north of Coimbra. Despite its modest footprint (~1.0–1.5 km across, ~100–110 hectares), it packs a dendrological collection of remarkable diversity curated over centuries. The topography rises from low foothills to the crest of the Serra do Buçaco, creating a humid, fog‑prone microclimate with precipitation notably higher than the surrounding region. That microclimate, plus deliberate introductions by botanists and gardeners since the 17th century, explain today’s extraordinary vertical structure: towering conifers (including giant sequoias), Mexican cypress, Atlantic and Tasmanian eucalypts, and groves of native broadleaves stitched between ornamental plantings and relic laurel‑oak patches.
LIDAR data converted to raster DSMLIDAR data converted to raster DTM
Walk any of the shaded paths and the “feel” of the forest is its third dimension: deep crowns stacked in tiers, emergent stems breaking above the canopy, and abrupt transitions where the slope pitches toward gullies and water stairs like the Fonte Fria. For mapping, this means Buçaco is the perfect stress‑test for vertical data. Local reports and lidar‑based profiles identify emergent trees approaching 60–65 m in height—exceptional for continental Europe—and many stands with 40–55 m canopy tops (giant sequoias Sequoiadendron giganteum, Tasmanian mountain ash Eucalyptus regnans, and mature Eucalyptus globulus among others). Add the steep relief and stone architecture of the palace‑convent complex and you have a site where coarse models tend to smear peaks, clip crowns, and understate vertical extremes.
From a data user’s perspective, Buçaco is interesting because it’s small enough to survey with dense airborne LiDAR yet diverse enough to benchmark against global canopy products. It’s also highly visited and well‑documented, which makes it a prime candidate for open, reproducible analyses that other practitioners can repeat.
LiDAR (and why it excels in forests)
Principle of operation. Airborne LiDAR instruments emit near‑infrared laser pulses toward the Earth’s surface and record the time‑of‑flight of returned photons. Distance = (c × Δt) / 2, where c is the speed of light and Δt is the measured two‑way travel time.
Full‑waveform vs discrete return. Modern sensors either store the entire returned energy waveform (full‑waveform) or extract distinct echoes (discrete returns). In forests, multiple returns (first, intermediate, last) capture interactions with the canopy top, internal branches, understory, and finally the ground.
Point cloud. Each pulse becomes a 3D point with XYZ, intensity, scan angle, GPS time, and often classification labels (ground, vegetation, building, water). Typical densities for national programs range from 2 to >12 points/m²; local surveys can exceed 20–30 points/m².
DTM and DSM. Ground classification filters (e.g., progressive TIN densification, cloth simulation) isolate ground returns to build a DTM. Interpolating the highest returns per cell builds a DSM that traces the top of canopy and built features.
Canopy Height Model (DHM). DHM = DSM − DTM at a chosen grid (often 0.5–2 m). Because the DTM is true bare earth, DHM measures canopy height above ground rather than above sea level—critical on steep slopes like Buçaco’s.
Vertical accuracy. With good boresight calibration and GNSS/INS trajectories, vertical RMSE for DTMs is commonly 5–15 cm in open ground; DHM accuracy depends additionally on canopy penetration and interpolation choices but still outperforms passive methods.
Structure metrics. From the point cloud or DHM we derive height percentiles (P10…P95), gap fraction, rugosity, leaf‑area proxies, and individual‑tree segmentation. These are the metrics that drive biomass, habitat, windthrow risk, ladder‑fuel detection, and view‑shed quality.
Radiometry & intensity. Intensity encodes target reflectance and range effects; after calibration, it helps distinguish materials (e.g., conifer vs broadleaf, moisture gradients) and detect powerlines or archaeological traces.
Waveform advantages. Full‑waveform captures the vertical distribution of scattering elements; deconvolution yields canopy penetration in denser stands and improves ground detection under eucalyptus and conifers.
Limitations. LiDAR is weather‑ and budget‑dependent. Dense undergrowth, scan angle, and leaf‑on conditions can reduce ground hits. Interpolation choices (max vs. percentile) affect DHM peaks—important when claiming “record” trees.
Bottom line: when you need true heights, crown architecture, and centimeter‑scale terrain under forest, airborne LiDAR remains the gold standard.
DSM-DTM=DHM (Global Mapper v26.1)
The global benchmark: ETH_GlobalCanopyHeight_2020_10m_v1
The ETH Zurich Global Canopy Height (GCH) product provides a wall‑to‑wall canopy top height map at 10 m ground sampling distance for the year 2020. It fuses GEDI lidar footprints (spaceborne, sparse but vertically precise) with globally consistent Sentinel‑2 optical imagery using a deep learning model to predict canopy heights between footprints. The result is a globally consistent raster that is easy to stream in Earth Engine or GIS platforms and ideal for continental to global analyses where airborne LiDAR is unavailable.
Global Canopy Height in TIF format extracted from GEE cloud computationVisualization of Global Canopy Height over the spot
Strengths
Global coverage at 10 m with a single epoch (2020), enabling cross‑region comparisons.
Trained on physically meaningful lidar targets (GEDI L2A/L2B canopy top metrics), correcting for many radiometric and terrain confounders in passive imagery.
Includes uncertainty metrics and tends to preserve macro‑patterns (ecotones, disturbance scars, plantation heights).
Known trade‑offs for sites like Buçaco
Saturation at the tall end. In stands with emergent stems >50 m, 10‑m pixels average crowns and can under‑predict peak heights; local maxima are “flattened.”
Terrain complexity. On steep slopes, small georegistration or DTM mismatches between Sentinel‑2 and GEDI training can leak terrain into predicted canopy height.
Edge effects. The palace complex, walls, and clearings introduce sharp transitions that are sub‑pixel at 10 m, broadening edges and obscuring narrow corridors.
Understory structure. The model predicts canopy top height, not vertical distribution; it cannot replace LiDAR‑derived structure metrics for habitat or fire modeling.
In short, ETH GCH is an excellent baseline and context layer, but for site‑scale management, airborne LiDAR remains the reference.
Practical comparison: LiDAR DHM vs ETH GCH over Buçaco’s 8818 vegetation GCP
Below is a workflow you can reproduce in QGIS/ArcGIS Pro or Google Earth Engine (GEE):
Ingest data.
Airborne LiDAR: download the point cloud (LAS/LAZ) or prebuilt DTM/DSM tiles for the Buçaco area.
ETH GCH 2020: load the ETH_GlobalCanopyHeight_2020_10m_v1 raster.
Build the LiDAR DHM.
Classify ground → DTM (0.5–1 m).
Highest‑return DSM (0.5–1 m) with spike filtering over built structures.
Harmonize grids. Aggregate DHM to 10 m by maximum or high percentile (P95) to compare fairly with ETH pixels while preserving tall peaks.
Sample and compare.
Randomly sample 5,000–20,000 points (I created 8818 GCP to sample) within the forest wall; extract DHM_10m and ETH_10m.
Compute bias (ETH − DHM), RMSE to see where ETH under/over‑estimates–> RMSE=12,97m (a bit too high!). Please try the code in GEE and you will also see a deviation map.
Tall‑tree check. Use a local maxima detector on the 1 m DHM to identify emergent crowns; intersect with ETH to quantify peak loss at 10 m.
Topographic controls. Regress residuals against slope, aspect, and curvature from the LiDAR DTM to diagnose terrain‑related biases.
Reporting. Summarize by species zones (sequoia groves, eucalyptus stands, relic laurel) if you have stand polygons or classify by crown texture.
Typical outcome in Buçaco (what to expect):
Median ETH bias close to zero over mid‑height stands (20–35 m).
Increasing underestimation in the tallest groves (e.g., −5 to −12 m at local maxima).
Larger residuals near walls/buildings and along steep steps and gullies.
RMSE calculation (think it needs further development thou)
Why sharing these data multiplies their value
Open elevation and canopy datasets have network effects. When agencies publish LiDAR DTMs/DSMs and derived DHMs under open licenses, practitioners can:
Validate global products locally, quantifying where models work and where they fail.
Stack analyses, from biodiversity corridors and storm‑damage assessments to micro‑hydrology, archaeology, and trail design, all anchored to the same precise terrain.
Build reproducible workflows, so results can be peer‑checked, improved, and extended.
Accelerate response, e.g., after windstorms or fires when canopy loss and debris flows must be mapped within days.
Educate and engage, by providing compelling 3D visualizations that show citizens and decision‑makers the invisible vertical dimension of their landscapes.
Portugal’s national investment in open, high‑accuracy remote sensing—airborne LiDAR and very‑high‑resolution imagery—has put the country to the level of Spain or France in terms of accurate shared open data.
Key sources used
Overview and context on Mata do Buçaco’s location, flora, and microclimate. Wikipedia
Este agosto, España y Portugal estan viviendo una temporada de incendios excepcionalmente dura. En España, las llamas han calcinado ya casi 390.000 hectáreas (más de seis veces la media reciente) y han dejado varias víctimas mortales, en Portugal, las superficies quemadas superan las 200.000 hectáreas, el 2% del total de su territorio (!) muy por encima del promedio 2006–2024 para estas fechas. El humo cruzó fronteras y degradó la calidad del aire a cientos de kilómetros…
A escala europea, el área ardida acumulada a 19 de agosto asciende a 895.000 hectáreas—casi cuatro veces la media para esta época del año. Es un pico que confirma lo excepcional del verano. ¿No hay algo raro?. ¿Es escepcional?
Si quieres entender un episodio como este sin ruido de los medios, el EFFIS (European Forest Fire Information System, del programa Copernicus) es tu herramienta:
Mapa de situación actual: muestra el peligro de incendio (FWI), focos activos y evolución diaria en Europa y Mediterráneo. Es la referencia pública y se actualiza continuamente.
Detección NRT de focos: integra “hotspots” (MODIS/VIIRS) que llegan en cuestión de horas tras el paso satelital para ver dónde arde ahora mismo.
Evaluación rápida de daños (RDA): mapea áreas quemadas con imágenes diarias para dimensionar el impacto durante la campaña, también con latencia de pocas horas.
Además, el Mecanismo de Protección Civil de la UE ha preposicionado medios aéreos y brigadas (rescEU) en Portugal y España este verano, lo que marca diferencia en los picos de simultaneidad.
EFFIS system. Copernicus (Burnt Area Locator)
La pregunta es, ¿por qué este año está habiendo tantos incendios?. No es un único factor: es la suma, y este año han “encajado” a mi juicio demasiadas piezas.
*Vientos y meteorología de episodio: rachas y condiciones locales han “disparado” la propagación en varios frentes activos. (Se aprecia en el rastro de humo en las imágenes satelitales sobre toda la península el 15 de agosto).
*Calor prolongado y sequedad extrema: la ola de calor más larga registrada en España (16 días) dejó combustibles finos listos para arder y propagarse rápido.
*Déficit hídrico acumulado: menos lluvia efectiva y suelos más secos elevan el índice de peligro (FWI) y alargan ventanas de ignición.
*Paisajes con mucha “carga de combustible”: Excesos de lluvia primaveral (llegando en elgunos casos a superar varias veces la media del periodo, como en Madrid), por otro lado el sempiterno abandono rural y la continuidad del matorral/masa forestal favorecen los fuegos de gran tamaño, especialmente con viento. (Conclusión coherente con los incendios recientes y el patrón observado en Galicia, Castilla y León y Extremadura).
Para finalizar no nos olvidemos del fuego de orígen antrópico de muchas igniciones: desde negligencias (gente haciendo barbacoas en medio el bosque, increíble, ¿no?) a presuntos incendios provocados investigados por las autoridades, con algunas personas que ya han entrado en prisión preventiva.
Sanabria surroundings. Zamora, Spain. 20250816
Para finalizar, si tienes tiempo, hoy y próximos días: revisa el FWI y su anomalía (riesgo meteorológico), los hotspots (dónde arde), y el RDA (qué se ha quemado). Si el FWI está en Muy Alto/Extremo y hay focos cercanos, espera propagación rápida. Contexto: compara en el portal de estadísticas las hectáreas de 2025 con su media 2006–2024 por país. Te ayuda a separar percepción de magnitud real. Tienes los links abajo.
Ponferrada surroundings. León, Spain. 20250816
Traducción para aquellos que no necesariamente están puestos en la materia (no tienen por qué, por otra parte): No es “mala suerte”. Es exposición + vulnerabilidad + clima más cálido. 2025 es la enésima prueba de estrés. La buena noticia: tenemos un sistema (EFFIS), cooperación europea reforzada y lecciones claras para prevención (paisaje, interfaz urbano-forestal) y respuesta (medios preposicionados, detección temprana). La mala: sin reducir riesgos estructurales y sin adaptación seria, estos picos serán cada vez más frecuentes (un conocido ha perdido por ejemplo su casa y todas sus pertenencias en uno de estos incendios~>nos puede pasar a cualquiera!)… Así, si has llegado hasta aquí, dime qué piensas. En este caso permitidme que os diga que “hay que mojarse”.
Combinación True color. Sentinel-2
Interfaz del Copernicus browser por la zona más concentrada de fuegos
GOES-East image over Portugal
Alberto C. Geógrafo y preocupado grado máximo con el hecho de que no nos demos cuenta como sociedad que estamos ante un problema gravísimo…
¿Cuántas veces has oído a tu cuñado (o cuñada) decir en una comida familiar frases tipo:
“¡Antes llovía más, se está desertificando todo!”
o peor aún:
“¡Yo ya lo noto, desde que era niño no ha vuelto a llover igual!”
Frases como estas suenan bien, tienen tono de verdad empírica… pero en realidad son pura intuición, sin ninguna evidencia científica detrás. Afortunadamente, vivimos en una época donde podemos rebatir con datos, mapas y gráficos, sin despeinarnos ni tener que abrir Excel.
CHIRPS en una gran fuente de datos desde hace más de 30 años!!. Añade un poco de CHIPRS a tu conversación navideña!
📡 ¿La clave? Computación en la nube + datos abiertos
En lugar de entrar en debates circulares, te propongo usar Google Earth Engine (GEE), una plataforma gratuita (si demuestras que no vas a usarlo comercialmente, claro) de computación geoespacial en la nube, y una fuente de datos de primer nivel: CHIRPS (Climate Hazards Group InfraRed Precipitation with Station data).
Haz zoom sobre tu zona de interés y ya lo tienes!
CHIRPS ofrece datos de precipitación diarios desde 1981 hasta la actualidad, a resolución global (~5 km). Es decir: puedes consultar con precisión cuánta lluvia ha caído en cualquier país durante los últimos 40 años, sin bajarte ni un solo archivo.
🔍 ¿Qué hace el script?
Te dejo un código en GEE que compara la precipitación media anual entre dos periodos (por defecto: 1981–2004 y 2005–2023), y genera:
🗺️ Un mapa en rojo y verde:
Rojo: zonas donde llueve menos ahora
Verde: zonas donde llueve más
📈 Un gráfico temporal con la evolución año a año y su tendencia lineal
📄 Un archivo CSV con la precipitación media anual de cada año
🗂️ Un archivo GeoTIFF binario que puedes abrir en QGIS, ArcGIS o cualquier visor de datos raster
🧪 ¿Cómo lo uso?
Solo tienes que cambiar una línea:
javascriptCopiarEditarvar countryName = 'Spain';
Y listo. Puedes poner 'France', 'Peru', 'Morocco' o el país que quieras. Todo se actualiza automáticamente: el código (3) el mapa binario (1), el diagrama (4), los nombres de los archivos exportados, el título del gráfico… incluso la leyenda (2).
ZASCA!!!!!!!!!!!!!!!!! Nadie será capaz de rebatirte con el GEE en la mano! 🙂
💻 ¿Y qué necesito instalar?
Nada. Cero. Solo una cuenta de Google Earth Engine y copiar el script en tu entorno. No necesitas tener experiencia previa con programación: el código está comentado y documentado para que sea lo más sencillo posible.
🧠 Bonus: más allá de tu cuñado
Además de ganar discusiones familiares, este tipo de análisis sirve para mucho más:
Validar percepciones locales con datos reales
Comunicar visualmente fenómenos complejos (como la variabilidad climática)
Enseñar geografía y ciencia de forma amena
Y sobre todo… para divertirte aprendiendo
Morocco, Tunisia y France en unos 30 segundos… Vamos que nos vamos!!!
¿Tienes un cuñado pesado? Reviéntalo con datos (esto es metafórico, claro). Y si no tienes uno… seguro que en Navidad aparece alguien con ganas de opinar. 😄
Alberto C. 🎓 Soy Geógrafo y me encantan las tecnologías geoespaciales. Creo firmemente que el conocimiento se democratiza cuando lo hacemos visual, accesible y directo. Analista GIS y feliz cuñado de Carlos y Koldo
Quieres conocer cuál es el momento óptimo para plantar? Para fumigar? Para recolectar?. Sabías que dos de cada tres agricultores no cosechan en la fase de madurez adecuada?. Aquí abajo te describo un método completamente automatizado mediante el uso combinado de varios índices de vegetación como NDVI, NDWI, SAVI y EVI que podemos extraer del Satétile SENTINEL-2 en la plataforma COPERNICUS de la UE para conocer exactamente y anticipar las mejores decisiones de intervención sobre tus tierras.
Precision farming / Agricultura de precision o cómo tomar decisiones adecuadas en el campo
Para comprender rápidamente, el gráfico de es un espectro de reflectancia para una superficie de grass (HIERBA) obtenido de una imagen Sentinel-2 L2A. Cada punto del gráfico representa la reflectancia media del suelo para distintas longitudes de onda, alineadas con las bandas espectrales del sensor Sentinel-2, que cubre desde el visible (VIS) hasta el infrarrojo cercano (NIR) y el infrarrojo de onda corta (SWIR). En verde oscuro el modelo (cómo normalmente reacciona la hierba a la luz, su reflectancia) y en verde claro (perdón soy un inepto con los colores) la medición específica sobre el punto que nos ocupa.
Por otro lado, si por ejemplo mido la reflectancia de la misma zona a lo largo del tiempo (Sentinel-2 ofrece una revisita de 5 días, lo que facilita el seguimiento temporal de la vegetación) , veo si la misma sube o baja, esto es importante como veremos más adelante.
Medición multitemporal/puntual de la reflectancia en cada uno de los 12 rangos que ofrece Sentinel 2
Visualmente vemos cómo claramente la humedad cambia a lo largo del tiempo (de maner obvia) para lo cual muestro una secuencia de aproximadamente el mismo día del año de tres años diferentes, una imagen del Infrarrojo Cercano IR que muestra en rojo intenso las zonas con más humedad.
Ilustración 1- Infrarrojo Sentinel 2
Si en lugar de medir puntualmente lo que hago es medir un área y extraer la media de todas las mediciones en una fecha o en una franja de tiempo determinada, tengo entonces unos datos sistemáticos que ya puedo tomar para medir por ejemplo la salud de la vegetación con un índice NDVI, que se calcula a partir de la reflectancia de la luz roja e infrarroja cercana IR de la vegetación.
El NDVI (Índice de Vegetación de Diferencia Normalizada) de Sentinel-2 es un índice que se utiliza para evaluar la salud y densidad de la vegetación, utilizando datos de las bandas rojas e infrarrojas cercanas de las imágenes Sentinel-2. Permite identificar áreas con vegetación, estimar su desarrollo y detectar cambios anormales en su crecimiento.
Los valores de NDVI varían de -1 a 1. Valores negativos indican ausencia de vegetación (agua, nieve, etc.), valores cercanos a cero indican superficies sin vegetación (rocas, suelo), y valores positivos indican vegetación, con mayor valor correspondiendo a vegetación más densa y saludable
NDVI = (B8 – B4) / (B8 + B4)
NDVI – Índice de Vegetación de Diferencia Normalizada
Para un cliente que desea determinar el momento óptimo para la siembra, la fumigación o la cosecha, podemos utilizar este mismo análisis de reflectancia (con datos normalizados para orregir variaciones atmosféricas residuales o nubosidad parcial), pero extendido de forma multitemporal y usando varios otros índices como el NDWI, SAVI, EVI, etc. (que explicaremos en detalle más adelante), para complementar la información espectral básica y relacionarla con la humedad del suelo, la cobertura vegetal y la temperatura superficial (si se cruza con otros sensores).
De momento sigamos con el NDVI, de él, conoceremos con certeza los momentos del año en los que la tierra ha estado más cargada de humedad (teniendo en cuenta y eliminando o no, si es necesario, mediciones con cobertura de nubes). Así podremos también analizar tendencias, se pueden detectar ventanas óptimas de humedad y temperatura del suelo que garanticen una mejor germinación y menor estrés hídrico inicial (propuesta de ventanas de siembra más estables según la variabilidad de la humedad interanual).
Interpretación de un NDVI index values 2020-25
NDVI throughout time 2021-2025
Pasamos ahora al NDWI: El NDWI (Índice Diferencial de Agua Normalizado) en Sentinel-2 es un índice espectral que se utiliza para identificar y mapear la presencia de agua en imágenes satelitales. Se calcula a partir de las bandas verde e infrarroja cercana y ayuda a diferenciar cuerpos de agua de otros elementos como la vegetación y el suelo.
Respecto a la interpretación, los valores positivos: generalmente indican la presencia de agua, con valores más altos correspondiendo a cuerpos de agua más limpios o con mayor contenido de agua. Los valores cercanos a cero pueden indicar suelo húmedo o vegetación con alto contenido de agua, mientras que los valores negativos suelen corresponder a suelo seco o vegetación con bajo contenido de agua.
NDWI = (B3 – B8) / (B3 + B8)
NDWI – Índice Diferencial de Agua Normalizado
Pasamos ahora al SAVI, el Índice de Vegetación Ajustado al Suelo (SAVI) de Sentinel-2 es una transformación de imágenes que intenta reducir la influencia del brillo del suelo en la estimación de la vegetación, especialmente útil en áreas con vegetación escasa o etapas iniciales de crecimiento.
El SAVI es un índice de vegetación que se calcula a partir de las bandas del rojo y del infrarrojo cercano de las imágenes Sentinel-2, similar al NDVI (Índice de Vegetación de Diferencia Normalizada), pero con un factor de corrección adicional para el suelo.
El NDVI, aunque ampliamente utilizado, puede verse afectado por el brillo del suelo, especialmente en áreas con poca vegetación o donde el suelo es visible entre la vegetación. El SAVI intenta corregir este efecto utilizando un factor de ajuste del suelo, denominado L, que varía según la densidad de la vegetación.
SAVI= ((B8 – B4) / (B8 + B4 + L)) * (1 + L)
SAVI – Índice de Vegetación Ajustado al Suelo
Por último veremos el índice EVI, este índice de Sentinel-2 es un índice que se utiliza para monitorear la salud de la vegetación, especialmente en áreas con alta densidad de biomasa. Se calcula a partir de las bandas 8 (infrarrojo cercano), 4 (rojo) y 2 (azul) de las imágenes Sentinel-2, utilizando una fórmula que ayuda a reducir el impacto de la atmósfera y del suelo en la señal de la vegetación.
EVI = G * ((B8 – B4) / (B8 + C1 * B4 – C2 * B2 + L)) *G: Ganancia (2.5 en Sentinel-2)
EVI – Índice de Vegetación Mejorado
Qué es lo que conseguimos al combinar con estos 4 análisis??. Lo primero unos diagramas conjuntos que nos muestran cómo se correlacionan espacialmente. Esto nos permite visualizar con rapidez los valores y ver cómo esta interacción de los índices NDVI, NDWI, SAVI y EVI nos van a ayudar en la toma de decisiones:
Condiciones favorables para la siembra con humedad suficiente y mínima competencia de vegetación anterior: Si NDVI < 0.3, SAVI < 0.3 y NDWI entre 0.2 y 0.4 Baja cobertura vegetal y nivel de humedad superficial moderado.
Estrés hídrico potencial: Si NDVI > 0.6 y NDWI < 0.2 Vegetación densa pero suelo seco → riesgo de déficit hídrico.
Vegetación poco vigorosa: Si NDVI < 0.3 y SAVI < 0.3 Vegetación escasa o en mal estado, posible daño o áreas sin cultivo.
Presencia de suelo desnudo o raleado: Si SAVI < 0.25 y NDVI entre 0.3 y 0.5 Vegetación con alta exposición de suelo, posible erosión o áreas sin cobertura completa.
Momento óptimo para fumigación o aplicación foliar: Si EVI > 0.7 y NDVI > 0.6 Alta biomasa y cobertura foliar → momento ideal para tratamientos.
Inicio de estrés por inundación o exceso hídrico: Si NDWI > 0.5 y NDVI < 0.5 Suelo saturado o anegado, vegetación afectada.
Madurez o cosecha cercana: Si NDVI disminuye rápidamente en semanas consecutivas y EVI baja de 0.5 La vegetación está senescente → momento probable de cosecha.
Crecimiento activo de la planta: Si NDVI y EVI aumentan simultáneamente durante varios días/semanas Etapa de crecimiento y desarrollo activo.
Posible estrés por plagas o enfermedades: Si NDVI cae por debajo de 0.4 mientras SAVI se mantiene estable Vegetación dañada sin cambios en suelo → posible plaga o enfermedad.
Zona con alta humedad superficial: Si NDWI > 0.4 y SAVI > 0.5 Suelo húmedo con vegetación densa, buen estado hídrico.
Pérdida de vegetación significativa: Si NDVI y SAVI bajan simultáneamente a menos de 0.2 Daño severo, incendios, o remoción de cultivos.
Como veis, se pueden deducir muchas cosas de la interacción de estos índices. Ya no hablamos de uno solo que nos muestre cierto detalle, hablamos de varios que nos dicen varias cosas en lugar de una:-) Esto nos ayuda a una intervención justificada, cuantitativa, empírica y científica. Una intervención que poco tiene que ver con las témporas, las creencias populares o incluso la experiencia, tiene que ver con la diferencia entre hacer las cosas bien o no.
Para terminar, comentar que este método es interesante pues el único input que necesita es un área de interés (automatización), es escalable en el sentido que puede aumentar su alcance mediante la realización de otros análisis en paralelo de otras variables (como precipitaciones o temperatura) y ayuda definitivamente a tomar decisiones.
1 Automatización: mediante scripts en interfaces de Cloud Computation (Computación en la nube) como Google Earth Engine o Google Cloud o un simple script Python desde QGIS o ArcGIS Pro, se puede generar un sistema recurrente que actualice cada cierto tiempo estos análisis.
2 Valor añadido: se pueden incorporar datos climáticos históricos (precipitaciones, temperatura) y modelos de predicción de heladas o lluvias.
3 Beneficio práctico: el agricultor optimiza la fecha de siembra, fumigación o recolecta mientras reduce riesgos de pérdidas por malas condiciones de arranque y ajusta mejor el calendario de trabajos.
Este flujo de trabajo está diseñado para que cualquier propietario de un terreno, independientemente de su finalidad, pueda comprender mejor cómo las lluvias o las temperaturas afectan a la misma. Esto le permitirá tomar decisiones informadas que optimicen su producción. Por supuesto, para un análisis más detallado y específico, estoy preparado para desarrollarlo, adaptándolo al tipo de cultivo y actualizándolo según sea necesario.
Si lo encuentras interesante, comparte, si no, ni se te ocurra;-)
Alberto C. Geógrafo, GIS Analyst, curioso, inquieto pero esperanzado ante un mundo cambiante.
A veces, las ideas más útiles nacen de lo que en España llamamos una paja mental. Es decir, una reflexión aparentemente banal, improvisada, casi sin rumbo… pero que, al desarrollarse, puede acabar dando lugar a algo valioso. En inglés no hay una traducción precisa (quizá algo como mental rambling), pero el concepto es claro: no todo parte de un briefing o de una necesidad concreta. Muchas veces es la simple curiosidad la que desencadena soluciones reales. Y eso también es innovación.
Spider diagrams using cloud computation and vibe coding in GEE – Random points testing
Hace unos días, caminando por una calle concreta —no muy lejos de mi casa, la cual recorro con regularidad— me pregunté cuál era realmente la estación de metro más cercana. Tiendo a pensar que siempre es Puerta del Ángel… y en realidad, lo es. Pero aún así, quise comprobarlo quantitativamente, empíricamente. ¿Y si no fuera tan evidente?
Finding the nearest underground station in Madrid downtown
Así arrancó mi pequeña “paja mental”: un análisis sencillo, basado únicamente en distancia en línea recta —sí, sin tener en cuenta la red de calles, que puede alterar mucho la percepción en ciudades como Madrid. Pero aun con esas limitaciones, decidí seguir.
Después decidí ampliar el experimento y llevarlo a una vía mucho más representativa y compleja: el Paseo de la Castellana. Aquí, a diferencia del caso anterior, las estaciones de metro no sólo están más próximas entre sí, sino que muchas pertenecen a la misma línea, lo que elimina la necesidad de realizar transbordos o desplazarse lateralmente por calles secundarias. Esto cambia por completo el enfoque del análisis.
¿Por qué? Porque aquí ya no se trata solo de saber cuál es la estación más cercana, sino cuál tiene mejor conectividad dentro de una misma línea, cuál permite una entrada/salida más eficiente, y cómo eso se distribuye a lo largo de todo el eje urbano. En una arteria como el Paseo de la Castellana —donde se mezclan oficinas, comercios, organisimos oficiales, centros de salud o zonas residenciales— esto permite analizar no solo accesibilidad, sino flujo potencial, saturación probable y distribución óptima de paradas.
Spider diagrams using cloud computation and vibe coding in GEE – Random points testing
Este tipo de análisis tiene aplicaciones muy concretas:
Planificación de nuevas estaciones o accesos: ver dónde la cobertura efectiva deja “zonas ciegas”.
Optimización de rutas peatonales: por ejemplo, priorizando aceras y pasos de peatones donde más flujo a estaciones se detecte.
Modelado de tiempos de acceso realistas: integrando red de calles, elevación, cruces semafóricos o barreras físicas.
Evaluación de equidad espacial: al comparar zonas con múltiples paradas bien conectadas frente a otras con acceso más difícil o indirecto.
Simulación de movilidad en emergencias o grandes eventos, para anticipar zonas de sobrecarga o rutas de evacuación eficaces.
Lo que empezó siendo una curiosidad sobre mi barrio se ha convertido en una aproximación simplificada pero útil (pienso) a cómo usamos y percibimos la red de transporte en entornos urbanos complejos. Y todo eso, a partir de una paja mental bien llevada.
En este tutorial se explica cómo crear un efecto visual dinámico en QGIS que simula una «araña» (ragnetto) que conecta automáticamente los 8 puntos más cercanos al cursor del mouse. Esta funcionalidad permite visualizar relaciones espaciales de forma interactiva y resulta especialmente útil para análisis de proximidad.
Respecto a la visualización del diagrama, he de agradecer a Franz (franzpc), el propietario de acolita.com y ArcGEEK por su espectacular tutorial, exactamente lo que estaba buscando, una mezcla de código y visualización! Gracias Franz!!!!
La gestión y monitorización de fenómenos hidrológicos extremos, como inundaciones repentinas o fallos estructurales en presas, representan un desafío crítico para los especialistas en geomática, hidrología y planificación territorial. En este contexto, la tecnología radar de apertura sintética (SAR) a bordo del satélite Sentinel-1 de la Agencia Espacial Europea (ESA) ofrece una capacidad sin precedentes para capturar información precisa y fiable sobre la dinámica superficial, independientemente de las condiciones atmosféricas y lumínicas.
Sentinel-1 opera enviando pulsos de microondas hacia la superficie terrestre y midiendo la señal reflejada, lo que permite superar las limitaciones típicas de los sensores ópticos que dependen de la luz visible y se ven afectados por nubosidad o ausencia de luz diurna. Esta independencia de las condiciones meteorológicas y de iluminación convierte a Sentinel-1 en una herramienta indispensable para la observación continua y en tiempo casi real de eventos catastróficos, especialmente en regiones donde las condiciones atmosféricas adversas son frecuentes.
El caso del colapso de presas en Derna, Libia, ocurrido en septiembre de 2023, es paradigmático para ilustrar las capacidades de Sentinel-1 en la detección y seguimiento de inundaciones masivas. Antes del evento, las imágenes radar adquiridas muestran una reflectividad característica del terreno seco y vegetado. Tras el colapso, la superficie inundada se manifiesta como áreas de baja reflectividad radar debido a la presencia de agua, que actúa como un reflector especular, causando un fuerte debilitamiento de la señal recibida. Este contraste en la intensidad radar permite una delineación precisa de las zonas afectadas.
El análisis temporal de imágenes Sentinel-1, mediante técnicas que incluyen la generación de mapas de diferencia logarítmica de intensidad (en decibelios), posibilita la identificación de cambios significativos en la superficie, facilitando la distinción entre áreas inundadas y no inundadas. Estas técnicas son especialmente valiosas cuando se emplean imágenes adquiridas en modos coherentes y con polarizaciones adecuadas, siendo la polarización VV (vertical transmitido, vertical recibido) comúnmente utilizada para capturar la variabilidad de la rugosidad superficial y contenido dieléctrico.
Disaster assessment using Cloud computation analysis in GEE over Lybia
Más allá de la detección inmediata, Sentinel-1 permite la monitorización de la evolución post-evento, lo que es fundamental para la gestión de emergencias y la planificación de la recuperación. Las series temporales de imágenes permiten evaluar la retirada progresiva del agua, la estabilización del terreno y la posible identificación de riesgos secundarios, como deslizamientos inducidos por saturación de suelos. La capacidad de generar mapas actualizados con alta frecuencia es una ventaja crítica que supera ampliamente a otras fuentes de datos.
Flooding Lybia 2023
Técnicamente, el procesamiento de datos Sentinel-1 para análisis de inundaciones incluye pasos como la calibración radiométrica, corrección geométrica, reducción del ruido speckle mediante técnicas de filtrado espacial o temporal, y la aplicación de umbrales o algoritmos de segmentación para generar máscaras binarias de agua/no agua. Estas máscaras pueden ser combinadas con modelos digitales del terreno (DEM) y datos de uso del suelo para mejorar la precisión y contextualizar la información. La integración con sistemas GIS permite generar productos cartográficos que apoyan la toma de decisiones en tiempo real.
Disaster assessment using Cloud computation analysis in GEE over Lybia
Disaster assessment using Cloud computation analysis in GEE over Lybia – Global Mapper 26.0 GIS evaluation
Un poco de historia de las Inundaciones en Libia esos días
El desastre se desencadenó la noche del 10 al 11 de septiembre de 2023, bajo el impacto de la poderosa “Medicane” (ciclón mediterráneo) Storm Daniel, que descargó precipitaciones extraordinarias en la región. En 24 horas, Derna registró hasta cerca de 400 mm de lluvia, un volumen asombroso considerando que su precipitación media anual apenas ronda los 300 mm.
Flooding Lybia 2023
Estas lluvias torrenciales colapsaron dos presas clave en el cauce del Wadi Derna—la presa Derna (Belad) y la presa Abu Mansour—con una capacidad combinada cercana a los 30 millones de m³. La ruptura de ambas estructuras liberó una ola de agua de hasta 7 m de altura, que arrastró todo a su paso por casi 12 km desde las zonas montañosas hasta el mar Mediterráneo.
La magnitud del impacto fue devastadora. Se estima que causó entre 5.923 y 20.000 muertes, con al menos 8‑10 000 personas desaparecidas, y dejó más de 30 000 desplazadosen.wikipedia.org+1en.wikipedia.org+1. Decenas de barrios completos desaparecieron, al igual que 4 puentes, y grandes secciones urbanas fueron barridas hacia el mar france24.com+2apnews.com+2es.wikipedia.org+2. La ciudad sufrió graves daños en su infraestructura crítica: hospitales inoperativos, morgues saturadas, destrucción masiva de edificios—más de 400 estructuras anegadas en lodo y escombros .
La escasez de mantenimiento en las presas—sin obras significativas desde 2002—y la falta de sistemas de alerta adecuados fueron factores determinantes. Aunque el servicio de meteorología libio emitió advertencias 72 horas antes, las instrucciones oficiales fueron confusas: se impuso toque de queda mientras se daban mensajes contradictorios sobre la evacuación. Esto atrapó a la población en zonas de riesgo justo cuando las presas fallaron hrw.org.
El colapso de las presas no solo generó muerte instantánea. En los días subsiguientes, se temió un brote de enfermedades renales y gastrointestinales debido a la contaminación del agua; más de 150 casos de diarrea infantil fueron reportados en Derna reuters.com+8apnews.com+8en.wikipedia.org+8. Pero debido al colapso de infraestructura sanitaria, el acceso al agua potable estaba interrumpido.
Esta tragedia fue catalogada como “el segundo fallido de presa más mortal de la historia”, solo detrás del caso de Banqiao en China (1975) en.wikipedia.org+1nrc.no+1. Además, han surgido procesos judiciales: en julio de 2024, el tribunal de Tobruk condenó a 12 responsables (hasta 27 años de prisión) por negligencia criminal reuters.com+1bbc.com+1.
El caso de Kirissah (a veces escrito Kirrissah), que se identifica con el Wadi Derna, sufrió un efecto similar: se formaron gigantescos depósitos de escombros, calles enteras quedaron enterradas y miles de personas fueron desplazadas o muertas arrastradas por el agua.
En resumen, el desastre de Derna supuso una conjunción letal de fuerte evento climático (Storm Daniel), infraestructura crítica en estado de abandono y falta de gestión de emergencia.
Sentinel‑1, mediante su capacidad para mostrar pérdidas intensas de energía radar (áreas oscuras = agua), permitió caracterizar espacialmente las zonas anegadas, lo que resulta clave para entender, mapear y gestionar tanto la emergencia inmediata como la reconstrucción posterior.
El Radar SAR del Sentinel-1 permitió caracterizar espacialmente las zonas anegadas
La bibliografía técnica destaca la importancia de Sentinel-1 como sensor radar para aplicaciones hidrológicas. Según el informe de la ESA sobre Sentinel-1 (ESA, 2023), la misión proporciona una cobertura global con revisit times de 6 a 12 días, resolución espacial de hasta 5×20 metros y disponibilidad gratuita de datos, lo que la convierte en una fuente accesible y fiable para la comunidad científica y agencias de gestión de riesgos. Investigaciones recientes han demostrado que la combinación de datos Sentinel-1 con algoritmos avanzados de detección y aprendizaje automático incrementa significativamente la capacidad para identificar y delimitar eventos de inundación con alta precisión (Smith et al., 2022; Zhang y Li, 2023).
En conclusión, Sentinel-1 representa un recurso estratégico para la observación y análisis de inundaciones y otros fenómenos naturales que alteran la superficie terrestre. El caso de Derna ejemplifica cómo esta tecnología puede proporcionar información esencial para responder ante desastres, evaluar daños y planificar acciones de mitigación y recuperación. La disponibilidad continua de datos y la evolución de técnicas de procesamiento posicionan a Sentinel-1 como una piedra angular en el arsenal de herramientas para la gestión de riesgos ambientales en el siglo XXI.
Espero que el post te haya resultado interesante, recuerda que la metodología para analizar cuantitativamente un evento de estas características (lluvias con alta intensidad horaria y rotura de una presa de grandes dimensiones) usando métodos de cloud computation permite adaptabilidad, escalabilidad, coherencia, precisión. Este análisis que acabas de leer, no es único y puede ser combinable con otros análisis que tienen en cuenta otros factores como mencionados anteriormente (la duración de los eventos, su cobertura espacial, la forma del terreno local, las variaciones estacionales y anomalías climáticas. Todo ello con datos abiertos, consistentes y reproducibles.)
Alberto C. Geographer, Senior GIS analyst, curious person and Open Data (not only) lover
Smith, J., et al. (2022). Advances in flood mapping using Sentinel-1 SAR data and machine learning. Remote Sensing of Environment, 280, 113156.
Zhang, H., & Li, X. (2023). Integration of Sentinel-1 SAR and optical data for improved flood monitoring. International Journal of Applied Earth Observation and Geoinformation, 125, 103035.
Gracias a datos satelitales de alta calidad y herramientas de análisis geoespacial en la nube, es posible generar de forma ágil mapas de peligrosidad climática con respaldo científico y cobertura global. Un primer ejemplo de ello es el análisis de la frecuencia de lluvias extremas realizado para Sri Lanka en el periodo 2001–2023.
El mapa ha sido elaborado utilizando la base de datos CHIRPS (Climate Hazards Group InfraRed Precipitation with Station data), un producto mixto satélite–estaciones que proporciona estimaciones de precipitación diaria desde 1981, con resolución espacial de aproximadamente 5 km. CHIRPS ha sido desarrollado por la Universidad de California en Santa Barbara y es una fuente de referencia en estudios de cambio climático, seguridad alimentaria y gestión de desastres.
Cloud computation analysis using GEE over Sri Lanka 2001-2023
En este caso, se ha calculado el número total de días con precipitaciones superiores a 50 mm durante más de dos décadas, generando una capa raster que refleja la recurrencia de eventos intensos a escala nacional. Este tipo de modelado permite identificar con claridad zonas más expuestas, lo cual es clave para análisis de riesgo climático, definición de primas diferenciales, evaluación de carteras o diseño de infraestructuras resilientes.
Inundaciones recientes en Colombo, Sri Lanka
Este tipo de información geoespacial, de origen fiable, actualizado y reproducible, puede ser de gran utilidad para entidades aseguradoras interesadas en mejorar su modelado de riesgos climáticos, optimizar sus procesos de tarificación o diseñar productos adaptados a exposiciones geográficas específicas.
Además, este enfoque puede escalarse para incorporar otros factores relevantes: la duración de los eventos extremos, su cobertura espacial, la forma del terreno local (que condiciona la escorrentía y la acumulación) o incluso variaciones estacionales y anomalías climáticas (como los efectos de El Niño o La Niña). Todo ello con datos abiertos, consistentes y reproducibles.
El valor no solo está en el mapa final, sino en la capacidad de generarlo de manera rápida, automatizable y fundamentada, con total trazabilidad en la metodología.
El valor no solo está en el mapa final, sino en la capacidad de generarlo de manera rápida, automatizable y fundamentada, con total trazabilidad en la metodología. He generado este código en Google Earth Engine y evaluado resultados con Global Mapper 26.0.
Usando software GIS para evaluar con precision (Global Mapper 26.1)
En un segundo caso, realizamos el análisis para un país (mucho) más grande, Venezuela. Para este análisis he usado otra fuente (ECMWF/ERA5/DAILY)
Lluvias intensas en Venezuela y su análisis de peligrosidad.
ERA5 es la última generación de reanálisis climático global desarrollado por el Centro Europeo de Previsiones Meteorológicas a Medio Plazo (ECMWF), dentro del programa Copernicus Climate Change Service (C3S). Se trata de un producto que combina observaciones y modelos numéricos para estimar variables atmosféricas, terrestres y oceánicas con alta resolución espacial y temporal.
El dataset ERA5 DAILY disponible en Google Earth Engine representa un subconjunto temporal diario de las variables meteorológicas del reanálisis ERA5. Cada imagen contiene variables agregadas o medias diarias derivadas del dataset ERA5 original, que tiene una resolución temporal de 1 hora.
En particular, la variable usada en tu caso es total_precipitation, que representa la precipitación acumulada total en el día, en metros (m), por ejemplo, 0.01 m = 10 mm.
Resolución espacial: Aproximadamente 0.25° × 0.25° (unos 28 km de lado); Resolución temporal: Diaria, cada imagen corresponde a un día calendario (UTC); Cobertura temporal: Desde 1979 hasta el presente (actualizada continuamente). ERA5 es generado mediante la asimilación de múltiples fuentes de datos, incluidos satélites, estaciones terrestres, boyas oceánicas, radiosondas y otros. Se utiliza un modelo climático numérico avanzado para producir un análisis global coherente. La calidad es muy alta, siendo uno de los reanálisis más usados para climatología y meteorología. La variable total_precipitation es acumulativa diaria y puede incluir precipitación líquida y sólida (nieve), según el modelo y condiciones.
La potencia de este código es poder añadir un país entero simplemente cambiando una palabra en el código, permitiendo una escalabilidad y harmonización directa.
Como se puede apreciar en un software GIS cualquiera (Global Mapper 26.0)
Colombia, Ecuador y Venezuela añadidos modificando levemente el código
En resumen, la metodología para analizar cuantitativamente la peligrosidad de un evento (lluvias con alta intensidad horaria) usando métodos de cloud computation permite adaptabilidad, escalabilidad, coherencia, precisión. Este análisis no es único y puede ser combinable con otros análisis que tienen en cuenta otros factores como mencionados anteriormente (la duración de los eventos, su cobertura espacial, la forma del terreno local, las variaciones estacionales y anomalías climáticas. Todo ello con datos abiertos, consistentes y reproducibles.)
Peligrosidad de acuerdo a las pendientes usando SRTM como fuente del DTM. Sri Lanka.GEE Interface and attribute table (including “peligrosidad” figures) to be exported as CSV
Geovisualization.net 