Download the exact captured file
---
name: academic-geosciences-remote-sensing-gis
description: "Specializes in Geosciences, Remote Sensing, Geoprocessing (GIS), Cartography, and Land-Use Planning building on Remote Sensing and Image Interpretation (Lillesand, Kiefer, Chipman), Geographic Information Systems and Science (Longley et al.), and Decifrando a Terra (Teixeira et al.). Covers Satellite Image Spectrometry (Sentinel-2, Landsat-8/9, CBERS-4A), Spectral Indices (NDVI, NDWI, NBR, EVI), Land Use and Land Cover Classification (Random Forest, Google Earth Engine), Digital Elevation Modeling (DEM/SRTM, Topographic Wetness Index TWI, Watershed Delineation), Spatial Analysis in GIS and Geospatial Databases (QGIS, PostGIS, GDAL/GeoPandas), Geodesy and Surveying (GNSS RTK, UTM Projections, SIRGAS2000 Datum), Thematic Cartography (Bertin's Visual Variables), Pedology (SiBCS), Hazard Geomorphology (CPRM/IPT Landslides), and Territorial/Environmental Legislation (CAR/SICAR, EIA/RIMA, City Statute and Master Plans)."
---
# Geosciences, Remote Sensing, and Spatial Analysis (Lillesand & Longley)
This skill establishes the engineering of geospatial data processing, the physics of electromagnetic radiation applied to remote sensing, relief and watershed modeling, spatial geostatistics, and compliance with legal frameworks for territorial and environmental planning.
---
## 🛰️ 1. Remote Sensing and Satellite Spectral Indices
```mermaid
flowchart LR
A["Incident Solar Radiation"] --> B["Terrestrial Surface (Soil, Vegetation, Water)"]
B -->|Spectral Reflectance| C["Orbital Sensors (Sentinel-2 MSI / Landsat-8 OLI)"]
C --> D["Pre-Processing (DOS1 Atmospheric Correction / Sen2Cor / TOA to BOA)"]
D --> E["Biophysical Index Calculation & Supervised Classification"]
```
### 1.1 Spectral Signatures and Normalized Difference Biophysical Indices
Before combining bands or dates, verify coordinate systems, grid alignment, resolution, radiometric scale and processing level. Preserve nodata, cloud and quality masks; select resampling appropriate to continuous or categorical data. For normalized indices, handle zero denominators explicitly rather than assigning an invented value.
- **NDVI (Normalized Difference Vegetation Index)**: Vigor of photosynthetically active biomass:
$$NDVI = \frac{\rho_{NIR} - \rho_{Red}}{\rho_{NIR} + \rho_{Red}}$$
- **EVI (Enhanced Vegetation Index)**: Optimized for high biomass densities with soil and atmospheric aerosol correction:
$$EVI = 2.5 \times \frac{\rho_{NIR} - \rho_{Red}}{\rho_{NIR} + 6 \rho_{Red} - 7.5 \rho_{Blue} + 1}$$
- **NDWI (Normalized Difference Water Index - Gao / McFeeters)**:
$$NDWI_{McFeeters} = \frac{\rho_{Green} - \rho_{NIR}}{\rho_{Green} + \rho_{NIR}}, \quad NDWI_{Gao} = \frac{\rho_{NIR} - \rho_{SWIR}}{\rho_{NIR} + \rho_{SWIR}}$$
- **NBR (Normalized Burn Ratio)**: Mapping of fire scars and burn severity:
$$NBR = \frac{\rho_{NIR} - \rho_{SWIR2}}{\rho_{NIR} + \rho_{SWIR2}}$$
### 1.2 Digital Land Use and Land Cover (LULC) Classification
- **Machine Learning Supervised**: Training Random Forest and Support Vector Machine (SVM) classifiers on band and spectral index sets in cloud environments (Google Earth Engine).
- **Confusion Matrix and Validation**: Kappa Index ($\hat{K}$) and Overall Accuracy:
$$\hat{K} = \frac{N \sum_{i=1}^k x_{ii} - \sum_{i=1}^k (x_{i+} x_{+i})}{N^2 - \sum_{i=1}^k (x_{i+} x_{+i})}$$
---
## 🗺️ 2. Digital Cartography, Geodesy, and Spatial Analysis (GIS)
### 2.1 Coordinate Systems and Geodesy
- **Geodetic Model**: Distinction between the Geoid (equipotential surface of the gravity field) and the Ellipsoid of Revolution (mathematical model).
- **Official Geodetic Datum in Brazil**: **SIRGAS2000** (Geocentric Reference System for the Americas, GRS80 Ellipsoid), compatible with **WGS84**.
- **Universal Transverse Mercator (UTM) Projection**: Conformal projection in 60 zones of $6^\circ$ longitude with central scale factor $k_0 = 0.9996$, central meridian with false easting $X_0 = 500,000\text{ m}$ and false northing $Y_0 = 10,000,000\text{ m}$ (southern hemisphere).
### 2.2 Geospatial Analysis with PostGIS and GeoPandas
Verify the source coordinate reference system before use: assigning a label does not transform coordinates. Transform to a suitable system using the correct datum operation and units before measuring distance or area. Check geometry validity and spatial-join cardinality; overlapping features can multiply rows and area totals.
The example below assumes both geometry columns already use the same suitable projected coordinate system with metre units. Without that precondition, dividing the planar area by 10,000 does not produce reliable hectares.
```sql
-- Example: Spatial query with buffer and intersection in PostGIS
SELECT
l.id_imovel,
ST_Area(ST_Intersection(l.geom, a.geom)) / 10000.0 AS area_sobreposta_ha
FROM
imoveis_rurais_car l
JOIN
areas_preservacao_permanente a
ON
ST_Intersects(l.geom, a.geom)
WHERE
l.municipio = 'Ribeirão Preto';
```
---
## 🏔️ 3. Geomorphology, Pedology, and Hydrological Modeling
### Marine carbonate observations
For seawater carbonate calculations, record the measured parameters and units, sample identity and quality flags, pH scale, salinity, nutrients, temperature and pressure at measurement and target conditions. Use an appropriate pair of independent carbonate parameters and document the selected equilibrium assumptions. Compare additional independent measurements against the calculation with their uncertainty rather than hiding disagreements. Equilibrium speciation alone does not establish air-sea flux or atmospheric carbon removal. See [NOAA's carbonate-system overview](https://www.pmel.noaa.gov/co2/story/Laboratory%2BAnalysis) for the measurement context.
### 3.1 Digital Terrain Modeling (DTM / DEM)
- **Topographic Wetness Index (TWI)**:
$$TWI = \ln\left( \frac{\alpha}{\tan\beta} \right)$$
where $\alpha$ is the specific upslope contributing area per unit contour ($m^2/m$) and $\beta$ is the local slope in radians.
- **Strahler Stream Ordering**: Automatic delineation of drainage networks where two streams of order $u$ converge to form an order $u+1$.
### 3.2 Pedology and Soil Classification in Brazil (SiBCS - Embrapa)
- **Order of Diagnostic Soils**:
- *Latossolos*: Latosolic B horizon (Bw), advanced weathering, deep, well-drained, and alic/dystrophic.
- *Argissolos*: Textural B horizon (Bt), strong textural gradient with susceptibility to sheet erosion and gullies.
- *Neossolos*: Shallow, poorly developed soils over parent rock (A horizon over R or C).
---
## 📋 4. Territorial Management, Environmental Forensics, and Legal Frameworks
| Instrument / Legislation | Legal Framework | Technical Application by the Geographer |
| :--- | :--- | :--- |
| **Brazilian Forest Code** | Federal Law No. 12,651/2012 | Georeferenced delineation of Permanent Preservation Areas (PPAs of watercourses, hilltops, and slopes $> 45^\circ$) and Legal Reserve (LR - $80\%$ Amazon, $35\%$ Cerrado in the Legal Amazon, $20\%$ other regions) via CAR/SICAR. |
| **City Statute** | Federal Law No. 10,257/2001 | Preparation of Participatory Master Plans, Urban Zoning (ZEIS), Progressive Property Tax over Time, and Onerous Grant of the Right to Build. |
| **Geological Hazard Zoning** | CPRM / IPT / Civil Defense Methodology | Mapping of hazard and vulnerability to planar/rotational landslides, debris flows, and gradual/sudden floods. |
| **Environmental Licensing** | CONAMA Resolutions 001/86 and 237/97 | Coordination of Environmental Impact Studies and Environmental Impact Reports (EIA/RIMA) and Degraded Area Recovery Plans (PRAD). |