Selam jeofizik dostlari, akademik tarafindan bir destek aramam gerekiyor. Hacettepe yuksek lisans danismanligini yurutugum bir tezde, Cankiri-Corum havzasinda iki helikopter manyetik survey verisi var, MTA arsivinden 1992 ve 2018 tarihli iki ayri ucus, ayni bolgede ortakli. Iki dataseti tek konturu uzerinde gostermek istiyoruz, mineralizasyon takibi acisindan.

Geosoft Oasis montaj 2024 ile basliyorum ama veri yapisi tutmuyor: 1992 verisi MGA AGSO formatinda eski binary, 2018 verisi GDB modern. Once MAGMAP olu reduktif filtre uygulayip iki veriyi ayni datum'a cekmem gerek (RTP — reduction to pole), sonra montaj. Ucus yuksekligi farki 1992'de 80 m, 2018'de 60 m, upward continuation ile yukseklige normalize edicem.

Sorum: (1) Oasis'in MAGMAP yeni surumunde RTP filtresi inclination/declination otomatik mi cekiyor IGRF-13'ten, yoksa elle mi giriliyor? Cankiri 2018 icin I=58.2°, D=4.8° elle hesaplamistim ama MAGMAP'in ic referansi farkli olabilir. (2) Iki survey arasinda baseline drift duzeltmesi nasil yapiyorsunuz, ozellikle 26 yil arayla?

Ezgi hocam, Geosoft Oasis 7 yildir kullaniyorum, arsiv data ile son donemli data karsilastirmasi sik isim. MAGMAP RTP filtresinde IGRF-13 entegrasyonu var ama otomatik degil — MAGMAP > Filter Setup penceresinde "Get IGRF" butonu var, latitude longitude ve tarih giriyorsun, otomatik I ve D dolduruyor. Senin elle hesapladigin I=58.2 D=4.8 buyuk ihtimal MAGMAP'in cikaracagiyla aynidir, %1 sapma normal kabul edilir.

Baseline drift 26 yil arayla buyuk konu. 1992'de cesium magnetometre, 2018'de Geometrics G-823A muhtemelen, sensitivity ve baseline kalibrasyonu farkli. Yapilacak: ortak alandaki istasyon noktalarini bul, 1992 ve 2018 arasinda overlap zone medyan farkini hesapla, bu fark sabit drift olarak 1992 datasina ekle. Yaklasik 5-15 nT mertebesinde olur tipik.

Ezgi hocam, Cankiri-Corum havzasinda biz de 2019'da yer manyetik calismasi yapmistik (helikopter degil, yer-bazli proton magnetometre), MTA'nin 2018 helikopter verisiyle karsilastirma denedik. Bizim deneyim: helikopter 60 m vs yer-bazli ~1.5 m yukseklik farki cok buyuk, upward continuation ile normalize edince derin kaynaklar guzel uyusuyor ama yuzlek kucuk anomaliler kayboluyor. Senin durumda iki helikopter survey'in ucus yuksekligi yakin (80 vs 60), upward continuation ile rahat normalize edilir.

Diger oneri: kontur uzerinde yapilacak ana islem second vertical derivative (2VD) veya analytic signal, ham TMI verisinden daha yorumlanabilir cikti verir, edge detection icin. Geosoft'un MAGMAP'inde her ikisi de tek tikla.

jeofizik_tolga yazdi:

ortak alandaki istasyon noktalarini bul, 1992 ve 2018 arasinda overlap zone medyan farkini hesapla

Tolga, overlap zone median farkı yontemi pratik ve net, tesekkurler. Iki survey ortak alani yaklasik 24 km2, GDB'de polygon ile maskeleyip iki verinin aritmetik median farkini cikartiyim. Bu sabit offset olarak 1992 datasina ekleyince baseline drift duzelmis olur.

Fatma hocam, 2VD ve analytic signal sinyal yorumu icin onerinizi alacagim. Bizim tezde mineralizasyon zonlarini tespit etmek hedef, edge detection idealdir.

Ezgi hocam, MTA arsiv mudurlugu ile sorabiliyorsanız 1992 datasinin original survey raporunu da isteyin. O raporda baseline kalibrasyon noktalari ve magnetometre serial numarasi olur. Eger calibration noktalari hala saha olarak erisilebilir ise (cogu kez fiziksel marker yok artık ama bazi yerlerde var) ek dogrulama yapilabilir.

1992-2018 arasi sekuler degisim de hesaba katilmali: IGRF'e gore Cankiri'da bu surede TMI ortalama 1100-1300 nT artmis (sekuler trend). Iki dataseti karsilastirirken bu sekuler drift de cikartilmali, sadece geological signal kalsin.

Cengiz abi, sekuler trend uyarisi cok onemli, hic dusunmemistim. 26 yil ciddi sekuler degisim, IGRF-13 ile her iki tarih icin TMI prediction cikartip aradaki farki dataseti ayrik bicimde temizlemem gerekecek. Yoksa geological signal'a sekuler drift karistirmis olurum. Bu noktayi yarinki danisman gorusmemde gundeme alacagim.

MTA arsiv mudurlugune dilekce attim 3 hafta once 1992 raporu icin, henuz cevap yok. Bugun tekrar takip edecegim.

Ezgi, Geosoft alternatifi de masaya yatir: Pythonda fatiando-terra ve harmonica kutuphaneleri ucretsiz, RTP, upward continuation, derivatives hepsi orada. Akademik yayin reproducibility acisindan Python tabanli script paylasmak Geosoft GDB binary paylasmaktan cok daha temiz. Geosoft ile final harita cikar, ama metodoloji icin Python wrapper yazmayi dusun.

Ezgi hocam, ben GPR'ciyim ama manyetik tarafindan da ortakliga giriyorum bazen. Oasis'in son surumunde "Geophysics Workflow Builder" var, Python script'i grafiksel iceren, tekil ve coklu islemleri zincirleyebiliyorsun. Senin sekuler trend cikarma + RTP + median offset + 2VD zinciri tek workflow olur, batch isletilir. 1992 ve 2018 verisi icin paralel iki workflow tanimlarsan tezde reproducibility temiz cikar.

GPR_dilara yazdi:

son surumunde "Geophysics Workflow Builder" var, Python script'i grafiksel iceren

Dilara, Workflow Builder'i acmamistim, hemen denedim — gercekten tek tikla zincirlenebiliyor. Umur'un dedigi fatiando-terra ile sonra paralel test yapacagim, Python scripti supplementary materyal olarak yayinda paylasilacak.

Ezgi, son bir oneri: kontur cikti olusturduktan sonra mineral oneri etkisi acisindan SVD (analytic signal'in 3D versiyonu, source depth estimation) yap. Cankiri-Corum'da Cu-Mo mineralizasyonu derinlikte oldugu bolgelerde SVD anomalin merkezi derinlik tahmini verir. Salem (2007) makalesi standart referans.

Orhan, SVD ve Salem 2007 makalesi notunu aldim, depth estimation tezin son bolumune girer. Cu-Mo prospect'i bizim tez konusu ile ortakli, hassas bir argument oluyor.

1992'de helikopter manyetik survey yapip dosya saklamis MTA, 26 yil sonra ayni bolgeye yine ucup yine helikopter parasi yakmis. Devletin parayla helikopter ucurmasini sevdigi belli, yan yana iki dataset cikinca herkes "vay be" diyor.

Ezgi hocam, BIM tarafiyla ilgisi sinirli ama bir oneri: kontur cikti haritalarini IFC4 spatial data extension ile 3D BIM modele entegre etmek mumkun, ozellikle eger Cu-Mo prospect'i sondaj ve isletme planlamasinda kullanilacaksa. Geosoft'tan dxf export et, sonra BIMcollab veya Trimble Connect ile sondaj kuyularinin yanina overlay olarak yerlestir. Tez disinda industri-uygulama icin pratik.

Ezgi, Workflow Builder'i + fatiando-terra paralel testi karari mukemmel. Tez raporunda metodoloji bolumu zenginlesir, savunmada reviewer "neden Geosoft, neden acik kaynak" diye sorabilir, sen "ikisiyle de paralel kontrol ettim, sonuclar uyusuyor (RMS difference X)" diye cevaplarsan disiplin gosterir.

Bir teknik nokta: 1992 binary AGSO formatinin Geosoft'a import'unda DAT2GDB veya MAGSURVEY plugin gerekebilir, modern Oasis'te eski binary destek bazen kayip oluyor. Hata alirsan MTA'nin 1992 raporundan veriyi ASCII'ye export ettirmek alternatif yol.

Toparliyorum: (1) IGRF-13'ten 1992 ve 2018 icin sekuler trend hesaplanip ayri ayri cikartilacak, (2) RTP MAGMAP otomatik, latitude/longitude/tarih girilerek, (3) overlap median fark ile baseline offset duzelmesi, (4) Upward continuation ile 80m-60m ucus yuksekligi normalize, (5) 2VD + analytic signal + SVD edge detection ve depth estimation, (6) Workflow Builder + fatiando-terra paralel test, (7) MTA arsiv 1992 raporu takibi.

Cok teknik bir tartisma cikti, tez giris bolumune dipnot olarak forum metodoloji notu eklemek istiyorum, izniniz olursa atif yaparak. Tesekkurler herkese.