Yapay Açıklıklı Radar Tekniği

Hafta 10:
Yapay açıklık radar (Synthetic Aperture Radar (SAR)) Giriş
SAR ile DEM Oluşturulurken İzlenecek Temel
Adımlar
 Uygulama için kullanılacak Sentinel-1A ve Sentinel 1-B Single Look Complex IW
görüntüleri, Avrupa Uzay Ajansı’nın Copernicus Sentinel Data Hub adlı internet veri
tabanında ücretsiz olarak alınmıştır.
 Sentinel-1 verilerinin işlenmesi, Avrupa Uzay Ajansı’nın (ESA) geliştirdiği, bünyesinde
Sentinel-1 Toolbox, Sentinel-2 Toolbox ve Sentinel-3 Toolbox adlı toolbox’ların kombine
edildiği, SNAP (Sentinel Application Platform) adlı yazılım kullanılmıştır. Yazılım,
Avrupa Uzay Ajansı’nın web sitesinden ücretsiz olarak edinilebilir.
 S1A_IW_SLC__1SDV_20161128T035801_20161128T035828_014137_016D38_F010
S1B_IW_SLC__1SDV_20161122T035718_20161122T035745_003066_00536D_50A5
S1A’nın yörünge numarası 14137, S1B’nin yörünge numarası 3066’dır.
SAR ile DEM Oluşturulurken İzlenecek Temel
Adımlar
 1. InSAR Stack Overview
 2-Coregistrasyon(S1 TOPS Coregistration )
 3. İnterferogram Oluşturma ve Coherence Yorumlama
 4. Debursting
 5.Goldstein Faz Filtreleme
 6.Multilooking
 7.Phase Unwrapping (Faz Açma)
 8. Fazdan Yüksekliğe Geçiş (Phase to Elevation)
 9. Range-Doppler Terrain Correction (Arazi Düzeltmesi)
1. InSAR Stack Overview

Şekil 1. Verilerin
açılması ve harita
üzerindeki yeri
Dik baz mesafesi olan Bperp (perpendicular
baseline) değeri
 Veri setleri işlenmeye başlanmadan önce,
uyduların veri setlerinin alındığı andaki
konumları arasındaki dik baz mesafesi olan
Bperp (perpendicular baseline) değerinin
bilinmesi gerekir. Bu değeri öğrenmek için,
SNAP yazılımında Radar>Interferometric
sekmesi altında bulunan “InSAR Stack
Overview” fonksiyonu kullanılır.
 Elde edilen görüntü S1A ve S1B görüntü
çifti, 128.78 metre dik baz mesafesine
(Bperp) sahiptir, ki bu mesafe verimli bir
uygulama için yeterli bir mesafedir.
2-Coregistrasyon(S1 TOPS Coregistration )

 Coregistrasyon, zemin üzerindeki her bir noktanın hem master hem de slave
görüntüsündeki aynı (range, azimuth) piksel içinde bulunmasını sağlar.
 Elde edilen iki adet kompleks veri seti, bünyesinde 3 adet şerit (swath) ve her bir şerit için
VV ve VH polarizasyonlarını içerir. Şeritler, içinde birbirinden ayrılmış yatay şeritler
(burst) bulundurur. Her burst, birbirlerinden çizgiler ile ayrılmıştır
 Yapılacak interferometrik uygulamada, iki veya daha fazla veri setindeki görüntüler tek bir
veri içinde birleştirilmek üzere “stack” haline getirilmelidir.
 Uygulama için her iki veri setindeki aynı şerit ve polarizasyonlar seçilmelidir. Aksi
takdirde yazılım hata verecektir.
TOPSAR-Split işlemi: IW1 şeridindeki VH polarizasyonlu görüntü seçilmiştir.
Yörünge dosyasının seçimi

 Apply Orbit File ve Apply

Orbit File(2) sekmelerinde
yörünge verisi olarak
“Sentinel Precise”
seçilmelidir. Bu yardımcı
yörünge verisi, SAR verisinin
alımının yapıldığı andaki
uydu konumu hakkında bilgi
içerir. Sentinel-1 uydularının
yörünge verileri, yazılım
tarafından otomatik olarak
indirilir.
Back-Geocoding

 Back-Geocoding sekmesinde,
istenilen türdeki dijital yükseklik
modeli (DEM) ve enterpolasyon
metodu seçilebilir. Yaygın olarak
kullanılan ve SNAP yazılımı
tarafından otomatik olarak indirilen
SRTM (Shuttle Radar Topography
Mission) yükseklik modellerinin
dünya genelindeki çözünürlükleri 1
arc-second için yaklaşık olarak
30 metre, 3 arc-seconds için
yaklaşık olarak 90 metredir.
Çakıştırılan Görüntünün Saklanması

 Write sekmesinde, bir

çıktı klasörü tanımlanır
ve çıktı verisi için
istenilen herhangi bir
isim verilebilir.
Ardından “Run”
butonuna tıklanarak
coregistrasyon işlemi
gerçekleştirilir.

Coregister verisinin kaydedilmesi

Coregister edilip stack haline getirilmiş veri
 Yörünge düzeltmesi yapılmış ve
stack haline getirilmiş veri,
aşağıdaki şekilde yazılımda açılır.
3. İnterferogram Oluşturma ve Coherence Yorumlama

 İnterferogram master görüntüsünün, slave

görüntüsünün karmaşık eşleniği ile çapraz
çarpılmasıyla oluşturulur. Burada, iki
görüntünün genlikleri çarpılırken, faz iki
görüntü arasındaki faz farkını temsil eder.
4. Debursting

 İnterferogram görüntüsündeki
burst’lerin tek bir görüntü oluşturmak
üzere birleştirilmesi için debursting
işlemi uygulanmalıdır.
Deburst işlemi uygulanmış
interferometrik faz
görüntüsü
 5.Goldstein Faz Filtreleme
İnterferometrik faz,
 Zamansal dekorelasyon
 Geometrik dekorelasyon
 Hacim yayılımı
 Görüntü işlemede oluşan hata
 gibi etkilerden kaynaklanan gürültü sebebi ile bozulmaya uğrayabilir. [12]
 İki görüntü arasındaki tutarlılıkta azaldığında, interference desenlerinde kayıplar olur.
Fazın doğru bir şekilde unwrap edilebilmesi için, faz üzerinde filtre uygulanarak
sinyal/gürültü oranının artırılması gerekir. Bu sayede interferogramda mevcut olan şerit
(fringe) desenleri daha belirgin hale getirilir.
 Goldstein faz filtreleme ile, interferometrik şerit desenleri (fringes) daha keskin hale gelir.
Bu keskinleştirme, spektrumda şeritlerden dolayı kaynaklanan tepe noktasının filtrelenmesi
ile tepe noktasına daha yüksek bir rölatif ağırlık verilmesi sonucu gerçekleştirilir.
Goldstein faz
filtresi
uygulanmış
interferogram
6.Multilooking

 Bütün radar görüntüleri, belirli bir derecede, adına “radar speckle”

denilen etkiyle birlikte görünür. Speckle etkisi, görüntüde “tuz ve
karabiber” dokusu şeklinde gözükür. Bunun nedeni, her bir
çözünürlük hücresinde ortaya çıkan, çoklu saçılma ile dönüş yapan
sinyallerin rastgele yapıcı ve yıkıcı interference’larıdır. Örnek
olarak, büyük ve çim kaplı araziler gibi homojen yapıya sahip
bölgeler speckle etkisi olmadığında genellikle açık tondaki piksel
değerlerine sahip olurlar (A). Fakat gerçekte, her bir çözünürlük
hücresindeki çim yapraklarından yansıyan sinyaller, görüntüdeki
bazı piksellerin ortalama tona göre daha açık veya daha koyu
renklerde oluşmasına sebep olur (B).
  Speckle etkisi, multilooking denilen işlem ile giderilebilir.
Multilooking, tek bakışa sahip (single look) radar ışınının
(A) daha dar yapıdaki alt ışınlara bölünmesi prensibine
dayanır. Her alt ışın ile, sensör tarafından alımı yapılan
bölgede bağımsız bir bakış (look) oluşur. Bu bakışlar
içinde kalan bölgelerde speckle etkisine yine rastlanacaktır.
Fakat oluşturulan her bakış içindeki speckle’ların toplanıp
ortalamasının alınması ile, oluşturulacak olan görüntüdeki
speckle etkisi azaltılmış olur. Multilooking genellikle
verinin uydu tarafından algılandığı sırada yapılır.

Multilooking
Multilooking uygulanmış interferometrik faz görüntüsü
7.Phase Unwrapping (Faz Açma)

 İnterferogram verisinde, interferometrik faz belirsiz durumdadır ve yalnızca bir döngü

(2π) içerisinde bilinir. İnterferometrik faz bilgisinin topoğrafik yükseklik bilgisiyle
ilişkilendirmek için, faz görüntüsünde unwrapping (açma) işlemi gerçekleştirilmelidir.

Unwrapping işlemi uygulanmış fazın

uygulanmamış faz ile kıyası
 Faz açma işlemi, komşu pikseller arasındaki faz farkını birbirlerine entegre ederek bu
belirsizliği giderir. İnterferogramdaki iki nokta arasındaki faz değişimi, gerçek yükseklik
değişiminin elde edilmesini sağlayan bir ölçüm sağlar. Uygulamadan en iyi sonuçları elde
etmek için, unwrapping işleminden önce multilooking ve faz filtremele işlemlerinin
gerçekleştirilmiş olması önerilir. Unwrapping işlemi uygulanan faz görüntüsünün kalitesi
ve güvenilirliği çoğunlukla girdi verisinin oluşturulmasında kullanılan master ve slave
görüntüleri arasındaki tutarlılığa (coherency) bağlıdır. Yalnızca tutarlılığın yüksek olduğu
bölgelerde güvenilir sonuçlar beklenmelidir.
 Unwrapping işlemi SNAP yazılımı
dahilinde değil, Linux tabanlı SNAPHU
adında bir algortima kullanılarak
gerçekleştirilir. SNAPHU, Stanford
Üniversitesi’nde Curtis Chen ve Howard
Zebker tarafından geliştirilmiş olan, faz
açma için kullanılan istatistiksel maliyetli
bir ağ akış algoritmasıdır.

Unwrapping işlemi
uygulanmış veri içinde yer
alan bantlar
Unwrapping işlemi uygulanmış faz görüntüsü
8. Fazdan Yüksekliğe Geçiş (Phase to Elevation)

 Unwrapping işlemi uygulanmış interfereometrik faz görüntüsünün, dijital yükseklik

modeline (DEM) dönüştürülmesi adımıdır.

İnterferometrik faz
görüntüsünden elde
edilen DEM
 9. Range-Doppler Terrain Correction (Arazi Düzeltmesi)

 Arazi Düzeltmesi, geometrik bozulmaları aynı bölgeye ait başka bir dijital yükseklik
modeli (DEM) kullanarak düzeltip bir harita projeksiyonlu bir veri üreterek görüntüyü
coğrafyayla kodlar. Coğrafi kodlama, eğik mesafe (slant range) veya zemin mesafesi
(ground range) geometrisindeki bir görüntüyü bir harita koordinat sistemine dönüştürür.
Arazi düzeltmesinde ekstra bir dijital yükseklik modeli kullanmasının sebebi
foreshortening, layover ve shadow (gölgelenme) gibi doğal geometrik etkileri gidermektir.
Range-Doppler arazi düzeltmesi
Range-Doppler arazi düzeltmesi uygulanmış DEM
uygulanmış verinin içindeki
elevation_VH bandı
  Arazi düzeltmesi
getirilmiş dijital yükseklik
modelindeli yükseklik
değerleri, Pixel Info
sekmesi açıkken görüntü
üzerinde farenin
gezdirilmesiyle
incelenebilir.

Dijital yükseklik modelindeki yükseklik değerleri

 Üretilen yükseklik modelinin gerçek topografya ile uyum
sağlayıp sağlamadığını görmek için, görüntü Google
Earth’te açılmak üzere kaydedilebilir. Bunun için,
görüntü üzerinde sağ tıklanarak Export View as Google
Earth KMZ seçilir ve görüntü KMZ uzantısında
kaydedilir. Ardından kaydedilen görüntü, Google Earth’te
açılıp incelenebilir.
 Üretilen yükseklik modelinin gerçek topografya ile uyum sağlayıp sağlamadığını görmek
için, görüntü Google Earth’te açılmak üzere kaydedilebilir. Bunun için, görüntü üzerinde
sağ tıklanarak Export View as Google Earth KMZ seçilir ve görüntü KMZ uzantısında
kaydedilir. Ardından kaydedilen görüntü, Google Earth’te açılıp incelenebilir.
Adımlar

 1. InSAR Stack Overview

 2-Coregistrasyon(S1 TOPS Coregistration )
 3. İnterferogram Oluşturma ve Coherence Yorumlama
 4. Debursting
 5.Goldstein Faz Filtreleme
 6.Multilooking
 7.Phase Unwrapping (Faz Açma)
 8. Fazdan Yüksekliğe Geçiş (Phase to Elevation)
 9. Range-Doppler Terrain Correction (Arazi Düzeltmesi)