POINT CLOUD REGİSTRASYON VE ALIGNMENT HATALARINI AZALTMAK

Topkapı Sarayı Üçüncü Avlu kapısının üst söve seviyesinde, üç farklı taramadan gelen taş yüzeyler birbirine 2,4 cm sapıyordu. Saha ekibi 41 istasyon almıştı; ofiste registrasyon raporu açıldığında ilk on iki istasyon mean error 1,8 mm ile temiz oturmuş, on üçüncüden sonra hata yığılarak gidiyordu. Sebep tek bir karelajdı: avlu ortasındaki dört hedef küresinden ikisi gün ortasında güneşten ısınıp tripodlarında milimetrik kaymıştı, sonraki tüm istasyonlar bu kayan referansa kilitlenip yamulttu. Vakanın değerlendirmesi şu olmuştu: hatayı yapan algoritma değil, hatayı görmeyen iş akışıydı.

Lazer tarama projesinin %70'i tarama sonrası registrasyon ve hizalama aşamasında geçer; saha günü tüm projenin sadece üçte biridir. Yapı belgeleme süreçleri için sektörel rehber niteliğindeki USIBD doküman standartlarında tarama doğruluğu LOA (Level of Accuracy) sınıflarına bağlanır ve bu sınıf, sözleşme imzalanırken tolerans bütçesini baştan çiviler. LOA30 talep edilen bir röleve işinde 5 mm'lik bir registrasyon sapması teslim edilemez sonuç demektir; saha disiplini bu eşiği önceden bilerek kurulmalıdır. Bu yazıda registrasyon hatasının nereden geldiğini, ICP ailesi çözücülerin parametrelerini, target-based ve cloud-to-cloud yöntemlerin gerçek farkını ve TKGM koordinat sistemine bağlanmanın disiplinini ele alıyoruz.

REGİSTRASYON NE TÜR BİR PROBLEMDİR

Point cloud registrasyon problemi, iki ya da daha fazla noktasal koordinat setini ortak bir referans çerçevesinde hizalamayı arayan bir optimizasyon işidir. Her tarayıcı istasyonu kendi yerel koordinat sisteminde veri üretir; tarayıcının sıfırı her istasyonda farklıdır. İşin matematiksel özü, bir tarama bulutunu rijid bir dönüşümle (üç eksenli dönme + üç eksenli öteleme, toplam altı serbestlik derecesi) diğerinin üstüne oturtmaktır. Çözücü, korelasyonu kuran nokta çiftleri arasındaki kare uzaklığın toplamını minimize edecek dönüş matrisi R ve öteleme vektörü T'yi arar.

Bu altı parametre kulağa az gelir; gerçekte sayıların duyarlılığı kritiktir. Bir istasyon 100 metre öteden tarıyorsa, dönme matrisindeki 0,01 derecelik bir hata uç noktada 1,7 cm sapma üretir. İstasyonlar zincirleme bağlandığında bu sapma katlanır; "hata yayılımı" denen olgu, beşinci istasyondan sonra registrasyon raporunda dikkat çekici şekilde tırmanmaya başlar. Hatayı ondan azaltmanın yolu tek bir istasyona kilitlenmek değil, ağı dağıtık tutmak ve her segmente jeodezik bir bağlantı verebilmektir.

HATA KAYNAKLARI HANGİLERİDİR?

Saha registrasyon raporu okumayı bilmeden yöntem tartışmak verimsizdir. Beş ana hata kaynağı vardır ve her birinin sapma raporunda farklı bir belirtisi bulunur.

Bu beş kaynak içinden ilk üçü saha disipliniyle çözülür; sonradan yazılım içinde "düzeltmek" mümkün değildir. Yetersiz overlap'i hiçbir çözücü icat etmez; veri yoksa veri yoktur. Dördüncü kaynak, hata yayılımı, kapalı çevrim oluşturmakla ve jeodezik nokta bağlamakla kontrol altına alınır. Beşinci kaynak ofiste segmentasyonla çözülür ama kaybedilen tarama süresi geri gelmez.

TARGET-BASED VE CLOUD-TO-CLOUD HANGİSİ NE ZAMAN?

Registrasyon iki ana stratejiyle yapılır ve profesyonel iş akışı çoğu zaman ikisini birleştirir. Karar saha planlama aşamasında verilir, ofiste verilmez.

Target-based registrasyon sahaya fiziksel referans nesneleri (sphere, checkerboard, retro-reflective hedef) yerleştirir ve istasyonları bu hedeflere göre kilitler. Her tarayıcı en az üç ortak hedefi gördüğünde rijid dönüşüm tek çözüme oturur. Hedef yerleşim disiplini şudur: her istasyon görüş alanında minimum dört küre, küreler 3D dağılımda (hepsi aynı düzlemde değil, farklı yüksekliklerde), bitişik istasyonlar arasında en az üç ortak küre. Hedef küresi malzemesi kritik: gün boyu güneşte kalan plastik küre 3-5 mm genleşir, sabah ölçülen referans öğleden sonra kaymış sayılır. Profesyonel proje magnetik tabanlı, gölgede tutulan veya köpük çekirdekli düşük-genleşme küreler kullanır.

Cloud-to-cloud registrasyon fiziksel hedef olmadan, iki taramanın doğal geometrisi (köşe, kenar, yüzey normalleri) üzerinden eşleştirme arar. ICP (Iterative Closest Point) ailesi algoritmalar tam bu işi yapar. Hızlı ve esnek ama yetersiz overlap'te ve düzgün geometrik özellik bulunmayan ortamlarda (boş depo, geniş cephe) yanılır. İki tarama %40 üstü örtüşüyorsa cloud-to-cloud genelde milimetre altı yakınsar, %20 altına düştüğünde sapar.

Hibrit yaklaşım sahada her iki ekonomiyi korur. Target-based birincil çözüm olarak kullanılır, raporun ana sayısı (mean target error) buradan gelir. Cloud-to-cloud doğrulama olarak çalıştırılır; eğer target-based 2 mm sonuç verdiyse cloud-to-cloud da 2 mm civarı çıkmalı. İki yöntemin farkı 3 mm'yi geçiyorsa, hedef yerleşiminde veya küre tutarsızlığında problem var demektir. Restorasyon ve LOA40 (1 mm) işlerde target-based zorunludur; LOA20 (1,5 cm) altı keşif işlerinde cloud-to-cloud tek başına yeterli olabilir.

ICP ALGORİTMASI VE PARAMETRELERİ

Cloud-to-cloud registrasyonun kalbi olan ICP, üç adımı belirli bir yakınsama eşiği sağlanana kadar tekrar eder. Birinci adım korelasyon: source bulutundaki her nokta için target bulutundaki en yakın nokta bulunur. İkinci adım dönüşüm tahmini: nokta çiftlerinin kare uzaklığını minimize eden R ve T hesaplanır (genellikle SVD ile). Üçüncü adım uygulama: source buluta R ve T uygulanır ve yeni RMSE ölçülür. Algoritma RMSE değişimi belirli bir eşiğin altına düşene veya max iterasyona ulaşılana kadar döner.

Bu temel akışın çıktısını belirleyen dört parametre vardır.

• Max correspondence distance: İki nokta arasındaki maksimum eşleşme mesafesi. Çok büyük olursa yanlış komşulukla eşleşir (yanlış yerel minimuma sıkışır), çok küçük olursa hiç eşleşme bulamaz. Tipik başlangıç bina röleve için 5-10 cm, kaba hizalamada 20-50 cm.
• Max iteration: Maksimum iterasyon sayısı. Çoğu uygulama 30-100 arası varsayılan tutar; karmaşık geometride 500-2000'e yükseltilir. Her iterasyon hesap maliyetidir.
• Convergence threshold: RMSE değişimi bu değerin altına düşerse algoritma durur. Tipik 1e-6, hassas röleve için 1e-8.
• Initial transformation: Kaba ön hizalama. ICP yerel optimizasyondur; başlangıç pozisyonu hedeften 1 metre ve 5 derece uzaktaysa yakınsar, 10 metre ve 45 derece uzaktaysa yanlış minimuma kilitlenir.

İki ana ICP varyantı vardır. Point-to-point ICP klasiktir; nokta çiftleri arası Öklid uzaklığını minimize eder. Point-to-plane ICP ise hedef yüzey normallerini kullanır ve düz yüzeylerde (duvar, döşeme, çatı) çok daha hızlı yakınsar. Bina rölevesinde point-to-plane neredeyse her zaman tercih edilir; doğal ortam ve organik form taramasında point-to-point daha güvenli kalır.

OVERLAP NE KADAR OLMALI?

Registrasyon hatasının en sık temel sebebi sahada yetersiz overlap'tir ve ofiste hiçbir çözüm getirmez. Komşu istasyonlar arası örtüşme oranı şu eşiklerle değerlendirilir: %30 altı zayıf (cloud-to-cloud yakınsamaz), %30-40 minimum kabul, %40-60 sağlam zon, %60 üstü gereksiz iş. Profesyonel saha planı %40-50 hedefler.

İstasyon planı CAD üzerinde önceden çizilir, sahada yapılmaz. Disiplinli plan üç soruyu cevaplar: hangi istasyondan hangi noktaları görebiliyorum, hangi istasyonlar arası overlap kaç metre, kapalı çevrim hangi istasyonlardan oluşuyor. Kapalı çevrim özellikle önemlidir: A-B-C-D-A şeklinde bir tarama döngüsü kurulduğunda, son istasyonda RMSE hata yayılımını anlık ölçer. Açık zincir registrasyon (A-B-C-D) son istasyonda hatayı saklar; kapalı çevrim hatayı görünür kılar ve dağıtmaya izin verir.

Tipik proje istasyon yoğunlukları şu seviyededir. Tek katlı küçük villa (150 m²): 12-18 istasyon. Çok katlı tarihi yapı, Ayasofya tip büyük strüktür: oda başı 1-3, büyük hacim için 8-12 istasyon, toplam 80-200. Bursa Demirtaş OSB ölçeğinde bir fabrika hattı: 60-120 istasyon. Belirleyici metrik m² değil, görüş engelleri ve detay yoğunluğudur. Her kolonun arkası, her duvar köşesi, her makine altı ayrı istasyon ister.

GNSS BAĞLANTISI VE TKGM KOORDİNATI

Türkiye'deki tüm jeodezik tabanlı işler ITRF96 referans çerçevesinde tanımlı TUTGA (Türkiye Ulusal Temel GPS Ağı) noktalarına ya da TUSAGA-Aktif istasyonlarına bağlanır. TKGM'nin tescilli koordinat sistemi bu çerçeveyi kullanır; restorasyon, kentsel dönüşüm, imar planı koordinasyonu işlerinde teslim edilen point cloud bu koordinat sistemine bağlı olmalıdır. Tarayıcı yerel sisteminde kalan bir tarama, sonraki BIM aşamasında imar planıyla çakıştırılamaz, kadastro pafta ile karşılaştırılamaz.

Bağlantı iki yoldan kurulur. Doğrudan georeferanslama yönteminde tarayıcı üzerine sabitlenmiş bir GNSS anteni vardır; tarayıcı her istasyonda kendi koordinatını TUSAGA-Aktif düzeltmesiyle doğrudan ölçer. Bu yöntem tek antenli sistemde sadece konum verir, yönelimi vermez; çift antenli kurulumda yönelim de ölçülür. Doğruluk genelde 2-5 cm, hassas RTK kurulumunda 1-2 cm civarındadır.

Dolaylı georeferanslama yönteminde sahada üç ya da daha fazla yer kontrol noktası seçilir; her noktada total station veya GNSS rover ile TKGM koordinatı kayıt edilir, aynı noktalar tarama bulutunda da hedef olarak ölçülür. Registrasyon önce tarayıcı yerel sisteminde tamamlanır, sonra "Helmert dönüşümü" ile bütün bulut TKGM çerçevesine taşınır. Doğruluk kontrol noktalarının kalitesiyle sınırlıdır; profesyonel projede 1-2 cm'ye iner. Restorasyon ve tarihi yapı belgelemede bu yöntem standarttır çünkü tek antenli doğrudan yöntemin doğruluğu yetmez.

RMSE NASIL OKUNUR?

Registrasyon yazılımları çıktı olarak bir sapma raporu üretir ve bu raporu okumayı bilmeden iş teslim edilmez. Üç sayı belirleyicidir.

• Mean error / overall RMSE: Tüm istasyonlar veya tüm hedef çiftleri arası ortalama sapma. Bina röleve için < 3 mm, restorasyon için < 1,5 mm hedef. Bu sayı projenin kalite imzasıdır.
• Max error: En kötü istasyon veya hedef çiftinin sapması. Mean'in 2 katından büyükse hata dağılımı homojen değil demektir; sorun belirli bir köşeye sıkışmış, bulup düzeltmek gerekir.
• Overlap percentage / link count: Komşu taramalar arası örtüşme oranı veya ortak nokta sayısı. %30 altı zayıf, %50 üstü sağlam. Düşük link sayısı ICP'nin yalnız yakınsadığı, sonuçların güvenilirliğinin azaldığı anlamına gelir.

USIBD LOA sınıflarının pratik karşılığı şudur: LOA10 (5 cm tolerans, kaba kayıt veya yaklaşık planlama), LOA20 (1,5 cm, MEP koordinasyonu için keşif), LOA30 (5 mm, BIM koordinasyon ve as-built modelleme), LOA40 (1 mm, restorasyon ve hassas dokümantasyon). Sözleşmede LOA30 yazıyor ama saha planı %25 overlap ve hedef küre yetersiz kurulmuşsa, raporda mean error 4 mm çıkar ve teslim reddedilir. LOA seçimi saha günü değil sözleşme günü yapılır.

BİR TARİHİ YAPI RÖLEVESİ SENARYOSU

İstanbul Tarihi Yarımada'daki orta ölçekli bir Osmanlı kervansarayının röleve projesi şu disiplinle yürüdü. Sözleşmede LOA30 (5 mm) hedef belirlendi. Saha planlama CAD üzerinde önceden yapıldı; 76 istasyon, ortalama %45 overlap, dört kapalı çevrim (avlu, koridor, üst kat, çatı seviyesi) tasarlandı. Sahaya 28 düşük-genleşme küre yerleştirildi, her istasyon görüş alanında en az dört küre kalacak şekilde dağıtıldı. Avluya üç yer kontrol noktası kondu, koordinatları RTK GNSS ile TKGM ITRF96 sisteminde kayıt edildi.

İlk tarama günü sabah 07:00-12:00 arası yapıldı, ikinci gün öğleden sonra 14:00-19:00. Bu süre dağılımı bilinçliydi: gün ortası güneşi taş yüzeylerde termal genleşme yaratır, sabah ve geç öğlen tarama yapan ekip 3-5 mm'lik termal sapmayı bypass etti. Registrasyon iki aşamada işledi. Önce target-based birincil çözüm üretildi; mean error 1,9 mm, max error 3,4 mm çıktı. Sonra cloud-to-cloud point-to-plane ICP doğrulama çalıştırıldı; mean error 2,1 mm. İki sonuç arası 0,2 mm tutarlılık, kabul edilebilir bantta. Helmert dönüşümüyle bütün bulut TKGM koordinat sistemine taşındı; üç kontrol noktası üstündeki kalan sapma 1,4 cm civarı, RTK doğruluğunun teorik tavanına oturdu.

Modele aktarım için bulut RCP üretildi, BIM yazılımında Insert Point Cloud akışıyla bağlandı, View Range içinde Display by Elevation modunda renk kodlu kesit alındı. Mimari ekip duvar deformasyonlarını, üst kornişin 3,1 cm dışa eğikliğini, kemer açıklığındaki 1,8 cm asimetriyi anında okudu. Bu seviye somutluk LOA30'un teslim ettiği değerdir; LOA10 ile aynı yapı taranmış olsaydı bu detayların hiçbiri görünmezdi. Konuyu uçtan uca sıralı bir programla çalışmak isteyen okurlar için point cloud eğitimi sayfası modülleri tek tek listeler.

HATA YAYILIMI NASIL KONTROL EDİLİR?

Açık zincir registrasyonda hata istasyon başına yığılır; A-B-C-D-E-F-G şeklinde yedi istasyon zincirleme bağlandığında, yedinci istasyondaki konum hatası ilk istasyonun 5-7 katıdır. Kapalı çevrim bu hatayı görünür kılar ve dağıtmaya izin verir; A-B-C-D-A şeklinde döngü kurulduğunda, son A istasyonunun ilk A ile arasındaki sapma yayılımın doğrudan ölçüsüdür. Yazılım bu sapmayı "loop closure error" olarak raporlar ve çevrim üzerinde dengeli dağıtım yapar.

İkinci kontrol katmanı yer kontrol noktalarıdır. Geniş bir projede her 15-20 istasyonda bir GNSS veya total station ile koordinatlandırılmış kontrol noktası tarayıcı görüş alanına konur. Bu nokta registrasyon zincirini jeodezik gerçeğe bağlar; zincir uç noktası kayarsa kontrol noktasıyla çarpışır ve sapma anında görülür. Bursa Demirtaş OSB ölçeğinde bir fabrika hattı taramasında bu disiplin uygulanmazsa, hattın sonundaki son makineye gelen sapma 6-10 cm'ye çıkar; BIM koordinasyonu için kullanılamaz hale gelir.

Üçüncü katman saha kalibrasyon kontrolüdür. Tanıdık ve değişmez bir referans (çelik kolon kenarları arası mesafe, betonarme kiriş açıklığı) hem fiziksel olarak şerit metre ile ölçülür hem tarama bulutunda noktadan-noktaya ölçü alınır. İki değer 3 mm tolerans içinde uyuşuyorsa ağ sağlam; sapıyorsa hız öznitelikleri veya hedef yerleşiminde sorun var. Bu üç kontrolü düzenli yapan saha ekibinde 100 istasyonluk taramanın registrasyonu 2-3 gün sürer; uygulamayan ekipte aynı iş iki haftaya yayılır ve %20-30 yeniden tarama gerekir.