Ders 16: Uzaktan Algılama ve Uydu Verileri
1. Giriş: Gökyüzünden Tarlaya Bakmak
Bir çiftçi tarlasını yürüyerek kontrol edebilir, ancak 500 hektarlık bir araziyi gözle kontrol etmek günler alır. Uzaktan algılama (remote sensing), uydu, uçak veya drone üzerindeki sensörlerle yeryüzeyini görüntüleme ve analiz etme bilimidir. Makine öğrenimi ile birleştiğinde, milyonlarca hektarlık tarım arazisini otomatik olarak izleme, mahsul sınıflandırma, verim tahmin etme ve çevresel stres faktörlerini tespit etme mümkün hale gelir.
Uzaktan algılamanın tarım için önemi şu verilerle somutlaşır: Dünya genelinde yaklaşık 1.5 milyar hektar tarım arazisi vardır ve bu arazilerin düzenli izlenmesi gıda güvenliği, su yönetimi ve iklim değişikliği adaptasyonu açısından önemli kabul edilir. Geleneksel tarla ölçümleri (in-situ) yalnızca belirli noktalarda doğru bilgi verir; uzaktan algılama ise geniş alanlı, tekrarlanabilir ve zamansal süreklilikte veri sağlar.
2. Uydu Platformları
2.1 Sentinel Serisi (ESA - Avrupa Uzay Ajansı)
Sentinel-1 (SAR Radarı): - Sentetik Açıklıklı Radar (SAR) taşıyan aktif uydu - C-bantı (5.405 GHz) ile çalışan, bulutlu havada bile görüntüleme yapabilir - Polarizasyonlar: VV (dikey-dikey) ve VH (dikey-yatay) - Mekansal çözünürlük: 5-20m - Tekrar süresi: 6 gün - Tarımsal kullanım: Toprak nemi tahmini, mahsul büyüme izleme, sel haritalama
Sentinel-2 (Çok Spektrumlu Optik): - 13 spektral bant ile görüntüleme - Mekansal çözünürlük: 10m (görünür ve yakın kızılötesi), 20m (kırmızı kenar ve kısa dalga kızılötesi), 60m (atmosferik bantlar) - Tekrar süresi: 5 gün - Tarımsal kullanım: Bitki indeksleri hesaplama, mahsul sınıflandırma, hastalık tespiti
Sentinel-2 Bantları ve Tarımsal Önemleri:
| Bant | Dalga Boyu (nm) | Çözünürlük (m) | Tarımsal Kullanım |
|---|---|---|---|
| B2 (Mavi) | 490 | 10 | Atmosferik saçınım düzeltmesi |
| B3 (Yeşil) | 560 | 10 | Bitki canlılık değerlendirmesi |
| B4 (Kırmızı) | 665 | 10 | Klorofil absorpsiyonu, NDVI hesabı |
| B5 (Kırmızı Kenar 1) | 705 | 20 | Bitki stres tespiti |
| B6 (Kırmızı Kenar 2) | 740 | 20 | Yaprak alan indeksi (LAI) |
| B7 (Kırmızı Kenar 3) | 783 | 20 | Kanopi yapısal bilgi |
| B8 (NIR) | 842 | 10 | NDVI hesabı, biyokütle tahmini |
| B8A (Dar NIR) | 865 | 20 | Su buharı düzeltmesi |
| B11 (SWIR 1) | 1610 | 20 | Nem içeriği, toprak nemi |
| B12 (SWIR 2) | 2190 | 20 | Bitki su stresi |
2.2 Landsat Serisi (NASA/USGS)
- 1972'den beri sürekli veri sağlayan en uzun ömürlü uydu programı
- Landsat-8/9: 30m çözünürlük, 16 gün tekrar süresi
- 11 spektral bant (termal dahil)
- Tarımsal kullanım: Uzun dönemli arazi kullanımı değişimi, iklim değişikliği etkisi izleme
2.3 MODIS (NASA)
- Günlük görüntüleme kapasitesi
- 250m-1km çözünürlük (düşük ancak zamansal süreklilik yüksek)
- 36 spektral bant
- Tarımsal kullanım: Kuraklık izleme, fenoloji analizi, bölgesel verim tahmini
2.4 Platform Karşılaştırması
| Platform | Mekansal Çöz. | Zamansal Çöz. | Spektral Bant | Maliyet | En Uygun Kullanım |
|---|---|---|---|---|---|
| Sentinel-2 | 10-60m | 5 gün | 13 | Ücretsiz | Tarla ölçeğinde mahsul izleme |
| Sentinel-1 | 5-20m | 6 gün | SAR (C-bant) | Ücretsiz | Bulutlu bölgelerde toprak nemi |
| Landsat-8/9 | 30m | 16 gün | 11 | Ücretsiz | Uzun dönemli değişim analizi |
| MODIS | 250-1000m | Günlük | 36 | Ücretsiz | Bölgesel/ulusal ölçek izleme |
| Planet | 3-5m | Günlük | 4-8 | Ücretli | Hassas tarım, yüksek çözünürlük |
| WorldView | 0.3-1.2m | İstek üzere | 8-16 | Pahalı | Tek bitki ölçeğinde analiz |
3. Spektral Bantlar ve Bitki İndeksleri
3.1 Elektromanyetik Spektrum ve Bitkiler
Bitkiler, elektromanyetik spektrumla benzersiz bir şekilde etkileşir: - Görünür ışık (400-700nm): Klorofil, mavi ve kırmızı ışığı emer; yeşili yansıtır - Kırmızı kenar (700-750nm): Sağlıklı bitkilerde keskin bir yansıım artışı (red edge) - Yakın kızılötesi (NIR, 750-1300nm): Yaprak hücre yapısı nedeniyle güçlü yansıım - Kısa dalga kızılötesi (SWIR, 1300-2500nm): Su içeriği bilgisi
Bu etkileşim, farklı bitki indekslerinin hesaplanmasına olanak tanır.
3.2 Temel Bitki İndeksleri
NDVI (Normalized Difference Vegetation Index):
$$\text{NDVI} = \frac{\text{NIR} - \text{Red}}{\text{NIR} + \text{Red}}$$
- Değer aralığı: -1 ile +1 arası
- Sağlıklı bitki: 0.6-0.9
- Cılız/stresli bitki: 0.2-0.4
- Su/toprak: < 0.1
- Kullanım: Bitki canlılığını izleme, verim tahmini, mahsul sınıflandırma
- Sınırlaması: Yoğun kanopi altında doygunluk (saturation) problemi
EVI (Enhanced Vegetation Index):
$$\text{EVI} = 2.5 \cdot \frac{\text{NIR} - \text{Red}}{\text{NIR} + 6 \cdot \text{Red} - 7.5 \cdot \text{Blue} + 1}$$
- NDVI'nın doygunluk sorununu azaltır
- Atmosferik ve toprak arka plan etkilerini minimize eder
- Yoğun bitki örtüsünde daha hassas
LAI (Leaf Area Index):
- Birim toprak alanı başına yaprak yüzeyi alanı (m^2/m^2)
- Sentinel-2'den biyofiziksel parametre olarak türetilir
- Tipik değerler: Mısır 3-6, Buğday 2-5, Orman 4-8
- Verim tahmini ve evapotranspirasyon modellerinde kritik girdi
SAVI (Soil-Adjusted Vegetation Index):
$$\text{SAVI} = \frac{(\text{NIR} - \text{Red})(1 + L)}{\text{NIR} + \text{Red} + L}$$
- $L$: Toprak parlaklık düzeltme faktörü (tipik olarak 0.5)
- Seyrek bitki örtüsünde (erken mevsim veya kurak bölgeler) NDVI'dan daha güvenilir
NDRE (Normalized Difference Red Edge):
$$\text{NDRE} = \frac{\text{NIR} - \text{RedEdge}}{\text{NIR} + \text{RedEdge}}$$
- Kırmızı kenar bandını kullanarak, klorofil içeriğine daha hassas
- Bitki azot durumu ve erken stres tespitinde üstün
3.3 İndeks Karşılaştırması
| İndeks | Avantaj | Dezavantaj | En Uygun Senaryo |
|---|---|---|---|
| NDVI | Basit, yaygın | Yoğun kanopide doygunluk | Genel bitki izleme |
| EVI | Doygunluk direnci | Mavi bant gerekli | Yoğun bitki (orman, mısır) |
| SAVI | Toprak etkisi düzeltmesi | Yoğun örtüde gereksiz | Kurak bölgeler, erken mevsim |
| LAI | Fiziksel anlam taşır | Ters çözümle hesaplanır | Verim modelleri |
| NDRE | Azot/strese hassas | Kırmızı kenar bant gerekli | Gübre yönetimi |
4. Google Earth Engine (GEE)
4.1 Nedir ve Neden Önemli?
Google Earth Engine, petabaytlarca uydu verisini bulut üzerinde işlemek için geliştirilmiş bir platformdur. Veri indirmeye gerek kalmadan, sunucu tarafında analiz yapılır.
Avantajları: - 40+ yıllık Landsat arşivi, tüm Sentinel verileri, MODIS vb. anında erişim - JavaScript ve Python API'si - Paralel işlem: Binlerce görüntü saniyeler içinde işlenir - Ücretsiz (araştırma ve eğitim için)
4.2 Örnek Kod: NDVI Zaman Serisi Hesaplama
// Google Earth Engine JavaScript API
// Türkiye sınırları içindeki buğday tarlalarında NDVI zaman serisi
var turkiye = ee.FeatureCollection('FAO/GAUL/2015/level0')
.filter(ee.Filter.eq('ADM0_NAME', 'Turkey'));
var sentinel2 = ee.ImageCollection('COPERNICUS/S2_SR')
.filterBounds(turkiye)
.filterDate('2025-03-01', '2025-08-31')
.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20));
// NDVI hesaplama fonksiyonu
var addNDVI = function(image) {
var ndvi = image.normalizedDifference(['B8', 'B4']).rename('NDVI');
return image.addBands(ndvi);
};
var ndviCollection = sentinel2.map(addNDVI);
// Aylık ortalama NDVI
var monthlyNDVI = ee.List.sequence(3, 8).map(function(month) {
return ndviCollection
.filter(ee.Filter.calendarRange(month, month, 'month'))
.select('NDVI')
.mean()
.set('month', month);
});
// Görselleştirme
Map.addLayer(ee.ImageCollection(monthlyNDVI).first(),
{min: 0, max: 0.8, palette: ['red', 'yellow', 'green']},
'NDVI Mart 2025');
5. SAR Radar: Bulutlu Havada Bile İzleme
5.1 SAR Temel Prensipleri
Sentetik Açıklıklı Radar (SAR), kendi mikrodalga sinyalini gönderip geri yansıyanı kaydeden aktif bir algılama sistemidir. Optik sensörlerden temel farkı:
- Bulut geçirgenlik: Mikrodalga dalgalar bulutları, sis ve hafif yağmuru deler
- Gece görüntüleme: Güneş ışığına ihtiyaç duymaz
- Yüzey yapısı: Toprak pürüzsüzlüğü ve nem içeriğine duyarlıdır
5.2 SAR ile Toprak Nemi Tahmini
Ben Abbes ve arkadaşları (2024), uzaktan algılama tabanlı toprak nemi tahmini konusunda kapsamlı bir derleme çalışması yayınlamıştır. Makale, fiziksel modeller, empirik yaklaşımlar ve makine öğrenimi yöntemlerinin entegrasyonunu incelemektedir.
Önemli Bulgular:
Toprak nemi (SM), tarımda sulama yönetimi, kuraklık izleme ve gıda güvenliği için kritik bir parametredir. Çalışmada incelenen yöntemler ve performansları:
| Yöntem | Veri Kaynağı | Çalışma Alanı | Performans |
|---|---|---|---|
| Random Forest + IEM | Sentinel-1A + Sentinel-2A | Güney Tunus | RMSE = 1.09%, R = 0.86 |
| Teacher-Student (CNN+RF) | Sentinel-1A + Sentinel-2A | Güney Tunus | RMSE = 0.8%, R = 0.96 |
| LSTM | SMAP | Carman, Kanada | 0.04 m^3/m^3 |
| ANN + CNN | SMOS | Küresel | RMSE: ANN=0.140, CNN=0.099 m^3/m^3 |
Çalışmada şu önemli noktalar vurgulanmıştır: - Sentinel-1 SAR verisinin toprak nemi tahmininde büyük potansiyeli vardır; ancak bitki örtüsü ve yüzey pürüzsüzlüğü sinyali etkiler - Makine öğrenimi yöntemleri, geleneksel fiziksel modellere (WCM, IEM) kıyasla daha yüksek doğruluk sağlamıştır - Sentinel-1 ve Sentinel-2 verisinin birlikte kullanılması (optik + radar füzyonu), tek başına kullanımdan daha iyi sonuçlar vermiştir - LSTM gibi zamansal modeller, toprak nemi zaman serisi tahmini için özellikle etkilidir
5.3 SAR Polarizasyon ve Tarımsal Anlamı
| Polarizasyon | Tarımsal Bilgi |
|---|---|
| VV | Toprak yüzey nemi, kısa bitkiler |
| VH | Bitki biyokütlesi, uzun bitkiler (mısır) |
| VV/VH oranı | Bitki büyüme aşaması |
6. Araştırma Örneği: Sentinel-2 ile Mısır Verim Tahmini
Desloires ve arkadaşları (2023), Sentinel-2 uydu görüntülerini kullanarak tarla ölçeğinde mısır verim tahminini makine öğrenimi yöntemleriyle gerçekleştirmiştir. Çalışmada 2017-2021 yılları arasında Iowa ve doğu Nebraska'daki 1164 tarla kullanılmıştır.
Veri ve Yöntem: - Sentinel-2'den türetilen bitki indeksleri: GNDVI, NDRE, NDWI - Biyofiziksel parametreler: LAI, FAPAR, fCOVER - Termal zaman (GDD) ile fenolojik normalleştirme - 3 büyüme evresi: vejetatif (V), üreme (R), yaşlanma (S) - Makine öğrenimi modelleri: Random Forest, XGBoost, sinir ağı
Yenilikçi Yaklaşım: Çoğu çalışmanın aksine, bu araştırma yıl-dışı (out-of-year) tahmin senaryosunu test etmiştir. Yani model, 2017-2020 verisiyle eğitilip 2021'i tahmin etmiştir. Bu, gerçek dünyadaki kullanım senaryosunu yansıtır: henüz hasat edilmemiş bir sezon için tahmin.
Sonuçlar: - LAI zirve değeri, verim tahmininde en etkili özellik olarak belirlenmiştir - Termal zamanla normalleştirme, takvim zamanı kullanımına göre tahmin doğruluğunu belirgin şekilde artırmıştır - Temmuz, tahmin için en kritik aydır (mısır için çiçeklenme ve tozlaşma dönemi) - Yıl-dışı senaryoda sıcaklık rejimleri arasındaki farklılık, modelin genelleştirme kapasitesini sınırlayan ana faktör olarak saptanmıştır
7. Uydu Verisi + ML İş Akışı
Tipik bir uzaktan algılama + makine öğrenimi projesi şu adımlardan oluşur:
1. Veri Edinme
|-- Sentinel-2 optik görüntüleri (GEE üzerinden)
|-- Sentinel-1 SAR verileri
|-- Meteorolojik veriler (sıcaklık, yağış)
2. Ön-İşleme
|-- Bulut maskeleme (SCL bantını kullanarak)
|-- Atmosferik düzeltme (Sen2Cor ile BOA reflektans)
|-- Mekansal birleştirme (20m bantları 10m'ye yeniden örnekleme)
3. Özellik Çıkarma
|-- Bitki indeksleri (NDVI, EVI, NDRE, SAVI)
|-- Biyofiziksel parametreler (LAI, fCOVER)
|-- Zamansal istatistikler (maks, min, ortalama, eğim)
|-- Fenolojik metrikler (sezon başlangıcı, zirve, bitiş)
4. Model Eğitimi
|-- Yer gerçekliği (ground truth) ile eşleştirme
|-- Özellik seçimi
|-- Model seçimi (RF, XGBoost, CNN, LSTM)
|-- Çapraz doğrulama
5. Tahmin ve Haritalama
|-- Piksel veya tarla bazında tahmin
|-- Sonuç haritalarının üretimi
|-- Doğruluk değerlendirmesi
8. Özet ve Önemli Noktalar
-
Sentinel-2, 10m çözünürlük ve 5 gün tekrar süresiyle tarımsal uzaktan algılamada en yaygın kullanılan ücretsiz uydu verisidir. 13 spektral bantı, çeşitli bitki indeksleri ve biyofiziksel parametrelerin hesaplanmasına olanak tanır.
-
Bitki indeksleri (NDVI, EVI, LAI, SAVI, NDRE) farklı amaçlara hizmet eder: NDVI genel izleme için, EVI yoğun örtüde, SAVI kurak bölgelerde, NDRE ise azot stresinde üstündür.
-
SAR radar (Sentinel-1), bulutlu iklim bölgelerinde optik verilerin kullanılamadığı durumlarda tarımsal izlemeyi mümkün kılar. Toprak nemi tahmininde ML yöntemleriyle entegre edildiğinde RMSE < %1 seviyesinde doğruluk elde edilmiştir.
-
Google Earth Engine, petabaytlarca uydu verisini bulut üzerinde işlemeye olanak tanıyan platform olarak tarımsal araştırmaları demokratikleştirmiştir.
-
Sentinel-2 ile mısır verim tahmini çalışmasında, fenolojik normalleştirme ve LAI zirve değeri, yıl-dışı tahmin senaryosunda kritik başarı faktörleri olarak belirlenmiştir.
Kaynaklar
- Ben Abbes, A., et al. (2024). Advances in remote sensing based soil moisture retrieval. Artificial Intelligence Review, 57, 224.
- Desloires, J., et al. (2023). Out-of-year corn yield prediction at field-scale using Sentinel-2 satellite imagery and ML methods. Computers and Electronics in Agriculture, 209, 107807.
- Drusch, M., et al. (2012). Sentinel-2: ESA's optical high-resolution mission for GMES operational services. Remote Sensing of Environment, 120, 25-36.
- Gorelick, N., et al. (2017). Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sensing of Environment, 202, 18-27.
- Rouse, J.W., et al. (1974). Monitoring vegetation systems in the Great Plains with ERTS. NASA Special Publication, 351, 309.