grdfftgmt - Bulutta Çevrimiçi

Program:

ADI

grdfft - Dalga sayısı (veya frekans) alanındaki ızgaralarda matematiksel işlemler yapın

SİNOPSİS

grdfft ingrid [ ingrid2 ] dış dosya [ azimut ] [ seviye ] [ [ölçek|g] ] [r|x|y][w[k]] ] [
[r|x|y]params ] [ [ölçek|g] ] [
[f|q|s|nx/ny][+a|d|h|l][+e|n|m][+tgenişlik][+w[sonek]][+z[p]] ] [ ölçek ] [ [seviye] ] [ -fg
]

Not: Seçenek bayrağı ve ilişkili bağımsız değişkenler arasında boşluk bırakılamaz.

TANIM

grdfft 2-B ileri Hızlı Fourier Dönüşümünü alacak ve bir veya daha fazlasını gerçekleştirecek
uzaya geri dönmeden önce frekans alanındaki matematiksel işlemler
ihtisas. Yeni değerleri bir çıktıya yazmadan önce verileri ölçeklendirmek için bir seçenek sunulur
dosya. Izgaranın yatay boyutlarının metre cinsinden olduğu varsayılmıştır. Coğrafi
ızgaralar belirtilerek kullanılabilir. -fg dereceyi metreye ölçekleyen seçenek. eğer varsa
boyutları km cinsinden olan ızgaraları kullanarak bunu metre olarak değiştirebilirsiniz. grededit veya ölçeklendir
ile çıktı matematik.

GEREKLİ ARGÜMANLAR

ingrid Çalıştırılacak 2 boyutlu ikili ızgara dosyası. (Aşağıdaki IZGARA DOSYA FORMATLARI'na bakın). İçin
çapraz spektral işlemler, ayrıca ikinci ızgara dosyasını verir inrd2.

-Gdış dosya
Çıktı ızgara dosyasının veya 1-D spektrum tablosunun adını belirtin (bkz. -E). (Görmek
IZGARA DOSYA FORMATLARI aşağıdadır).

İSTEĞE BAĞLI ARGÜMANLAR

-Aazimut
Yönlü türevi alın azimut derece CW cinsinden ölçülen yön
kuzeyden.

-Cseviye
Yukarı doğru (için seviye > 0) veya aşağı doğru (için seviye < 0) alana devam et seviye
metre.

-NS[ölçek|g]
Alanı farklılaştırın, yani d(field)/dz alın. Bu çarpmaya eşdeğerdir
frekans alanında kr ile (kr radyal dalga sayısıdır). Bir ölçek ekle
ile çarp (kr * ölçek) Bunun yerine. Alternatif olarak, ekle g olduğunu belirtmek için
veriler metre cinsinden jeoid yüksekliklerdir ve çıktı mGal cinsinden yerçekimi anomalileri olmalıdır.
[Varsayılan ölçek yok].

-E[r|x|y][w[k]]
Radyal yönde güç spektrumunu tahmin edin [r]. Yer x or y hemen sonra
-E bunun yerine spektrumu x veya y yönünde hesaplamak için. ızgara dosyası yok
yaratıldı. Bir ızgara verilirse, f (yani, frekans veya dalga numarası), güç[f],
ve güçte [f] 1 standart sapma ile ayarlanan dosyaya yazılır -G [stdout]. Eğer
f ve 8 nicelik yazıyoruz: Xpower[f], Ypower[f], tutarlı
güç[f], gürültü gücü[f], faz[f], kabul[f], kazanç[f], tutarlılık[f]. Her biri
miktarı kendi 1-std dev hata tahmini takip eder, dolayısıyla çıktı 17'dir
sütunlar geniş. Ekle w frekans yerine dalga boyu yazmak için. Eğer ızgaranız
coğrafi daha fazla ekleyebilirsiniz k dalga boylarını metre [Varsayılan] ile ölçeklendirmek için
km.

-F[r|x|y]params
Verileri filtreleyin. Yer x or y hemen sonra -F filtrelemek için x or y sadece yön;
varsayılan izotropiktir [r]. Kosinüs konik bant geçişi, Gaussian arasında seçim yapın
bant geçiren filtre veya Butterworth bant geçiren filtre.

Kosinüs-konik:
Dört dalga boyu belirtin lc/lp/hp/hc doğru birimlerde (bkz. -fg) tasarlamak için
bant geçiren filtre: daha büyük dalga boyları lc veya daha az hc kesilecek,
dalga boyları daha büyük lp ve daha az hp geçirilecek ve dalga boyları
arasında kosinüs konik olacaktır. Örneğin, -F1000000 / 250000 / 50000 / 10000 -fg
bant geçirecek, dalga boylarını > 1000 km ve < 10 km kesecek, geçecek
250 km ile 50 km arası dalga boyları. Yüksek geçiren veya alçak geçiren filtre yapmak için,
için tire (-) verin hp/hc or lc/lp. Örneğin, -döviz-/-/50/10 alçak geçiş yapacak x,
geçen dalga boyları > 50 ve dalga boylarını reddeden < 10. -fy1000/250/-/-
yüksek geçiş yapacak y, geçen dalga boyları < 250 ve dalga boylarını reddeden > 1000.

Gauss bant geçişi:
eklemek lo/hi, doğru birimlerde iki dalga boyu (bkz. -fg) tasarlamak için
bant geçiren filtre Verilen dalga boylarında Gauss filtre ağırlıkları
0.5 olsun. Yüksek geçiren veya alçak geçiren filtre yapmak için, kısa çizgi (-) koyun. hi
or lo sırasıyla dalga boyu Örneğin, -F-/30, verileri bir
30'da yarı ağırlıklı Gauss filtresi, -F400/-
verileri.

Butterworth bant geçişi:
eklemek lo/hi/sipariş, doğru birimlerde iki dalga boyu (bkz. -fg) Ve
bir bant geçiren filtre tasarlamak için filtre sırası (bir tamsayı). verilen anda
dalga boylarında Butterworth filtre ağırlıkları 0.5 olacaktır. Yüksek geçiş yapmak için
veya alçak geçiren filtre için bir tire (-) verin. hi or lo dalga boyu
sırasıyla. Örneğin, -F-/30/2, 2. sırayı kullanarak verileri alçak geçirecek
30'da yarı ağırlığa sahip Butterworth filtresi, -F400/-/2 yüksek geçiş yapacak
veri.

-BEN[ölçek|g]
Alanı entegre edin, yani integral_over_z'yi (alan * dz) hesaplayın. Bu
frekans alanında kr ile bölmeye eşdeğerdir (kr, radyal dalga sayısıdır).
Bölmek için bir ölçek ekleyin (kr * ölçek) Bunun yerine. Alternatif olarak, ekle g için
veri kümenizin mGal'deki yerçekimi anormallikleri olduğunu ve çıktının geoid olması gerektiğini belirtin
metre cinsinden yükseklikler. [Varsayılan ölçek yok].

-N[f|q|s|nx/ny][+a|[+d|h|l][+e|n|m][+tgenişlik][+w[sonek]][+z[p]]
FFT için uygun ızgara boyutlarını seçin veya sorgulayın ve isteğe bağlı olarak ayarlayın
parametreler. FFT boyutunu kontrol edin:
-NF FFT'yi verilerin gerçek boyutlarını kullanmaya zorlar.

-Nq daha uygun boyutlar hakkında bilgi alacak, bunları rapor edecek ve devam edecek.

-N'ler isteğe bağlı boyutların bir listesini sunacak ve ardından çıkacaktır.

-Nnx/ny dizi boyutunda FFT yapacak nx/ny (>= ızgara dosya boyutu olmalıdır). Varsayılan
FFT'nin hızını ve doğruluğunu optimize eden boyutları >= verileri seçer. FFT ise
boyutlar > ızgara dosyası boyutları, veriler uzatılır ve sıfıra doğru incelir.

Verilerin trendden arındırılmasını kontrol edin: Doğrusal bir trendi kaldırmak için değiştiriciler ekleyin:
+d: Trend verileri, yani en uygun doğrusal trendi kaldırın [Varsayılan].

+a: Yalnızca ortalama değeri kaldırın.

+h: Yalnızca orta değeri kaldırın, yani 0.5 * (maks + min).

+l: Verileri rahat bırakın.

Verilerin genişletilmesini ve daraltılmasını kontrol edin: Uzantının nasıl kontrol edileceğini kontrol etmek için değiştiricileri kullanın.
ve daraltma yapılacaktır:
+e kenar noktası simetrisi uygulayarak ızgarayı genişletir [Varsayılan],

+m kenar ayna simetrisi uygulayarak ızgarayı genişletir

+n veri uzantısını kapatır.

Daraltma, veri kenarından FFT ızgara kenarına [%100] gerçekleştirilir. Değiştirmek
yoluyla bu yüzde +tgenişlik. Ne zaman +n yürürlüktedir, konik uygulanır
hiçbir uzantı mevcut olmadığından veri marjları yerine [0%].

Geçici sonuçların yazılmasını kontrol edin: Ayrıntılı araştırma için aşağıdakileri yazabilirsiniz.
ileri FFT'ye geçirilen ara ızgara; bu büyük ihtimalle olmuştur
Eğilimi bozulmuş, tüm kenarlar boyunca nokta simetrisi ile uzatılmış ve konik. Ekle
+w[sonek] hangi çıktı dosya adlarının oluşturulacağı (yani, ingrid_prefix.ext)
[konik], nerede ext dosya uzantınızdır. Son olarak, karmaşık ızgarayı kaydedebilirsiniz.
ekleyerek ileri FFT tarafından üretilen +z. Varsayılan olarak gerçeği yazıyoruz ve
hayali bileşenler ingrid_gerçek.ext ve ingrid_resim.ext. Ekle p kaydetmek için
bunun yerine dosyalara büyüklük ve fazın kutupsal biçimi ingrid_mag.ext ve
ingrid_faz.ext.

-Sölçek
Her öğeyi şu şekilde çarpın: ölçek uzay alanında (frekans alanından sonra
operasyonlar). [Varsayılan 1.0'dır].

-V[seviye] (Daha ...)
Ayrıntı düzeyini seçin [c].

-fg Coğrafi ızgaralar (boylam, enlem boyutları) metreye dönüştürülecek
mevcut elipsoid parametreleri kullanılarak bir "Düz Dünya" yaklaşımı yoluyla.

-^ or sadece -
Komutun sözdizimi hakkında kısa bir mesaj yazdırın, ardından çıkar (NOT: Windows'ta
sadece kullan -).

-+ or sadece +
Açıklama da dahil olmak üzere kapsamlı bir kullanım (yardım) mesajı yazdırın.
modüle özgü seçenek (ancak GMT ortak seçenekleri değil), ardından çıkar.

-? or yok hayır argümanlar
Seçeneklerin açıklamasını da içeren eksiksiz bir kullanım (yardım) mesajı yazdırın, ardından
çıkışlar.

--versiyon
GMT sürümünü yazdırın ve çıkın.

--gösteri-veri dizini
GMT paylaşım dizininin tam yolunu yazdırın ve çıkın.

GRID DOSYA BİÇİMLER

Varsayılan olarak GMT, ızgarayı COARDS şikayeti netCDF'de tek duyarlıklı yüzer olarak yazar
dosya formatı. Bununla birlikte, GMT, yaygın olarak kullanılan diğer birçok ızgarada ızgara dosyaları üretebilir.
dosya biçimlerini destekler ve ayrıca ızgaraların "paketlenmesini" kolaylaştırır, kayan noktayı yazar
1 veya 2 baytlık tamsayılar olarak veri. Kesinlik, ölçek ve ofseti belirtmek için kullanıcı şunları yapmalıdır:
son eki ekle =id[/ölçek/ofset[/nan]], nerede id ızgaranın iki harfli tanımlayıcısıdır
tip ve hassasiyet ve ölçek ve ofset isteğe bağlı ölçek faktörü ve ofset
tüm ızgara değerlerine uygulanır ve nan eksik verileri belirtmek için kullanılan değerdir. Durumunda
iki karakter id olarak sağlanmamaktadır. =/ölçek a'dan daha id=nf varsayılır. Ne zaman
ızgaraları okurken, format genellikle otomatik olarak tanınır. Değilse, aynı son ek
ızgara dosya adlarını girmek için eklenebilir. Görmek grddönüştürme ve Bölüm ızgara-dosya formatı
Daha fazla bilgi için GMT Teknik Referans ve Yemek Kitabı.

Birden çok ızgara içeren bir netCDF dosyasını okurken, GMT varsayılan olarak şunları okuyacaktır:
o dosyada bulabilen ilk 2 boyutlu ızgara. GMT'yi başka bir okumaya ikna etmek için
ızgara dosyasındaki çok boyutlu değişken, ekleme ?Değişkenadı dosya adına, nerede
Değişkenadı değişkenin adıdır. Özel anlamdan kaçmanız gerekebileceğini unutmayın.
of ? Shell programınızda önüne ters eğik çizgi koyarak veya
tırnaklar veya çift tırnaklar arasında dosya adı ve sonek. NS ?Değişkenadı eki de kullanılabilir
çıkış ızgaralarının varsayılandan farklı bir değişken adı belirtmesi için: "z". Görmek
grddönüştürme ve GMT Teknik bölümünün CF için değiştiricileri ve ızgara dosya formatı
Daha fazla bilgi için Referans ve Yemek Kitabı, özellikle 3-,
4 veya 5 boyutlu ızgaralar.

GRID MESAFE BİRİMLERİ

Izgara yatay birim olarak metreye sahip değilse, +ubirim giriş dosyasına
Belirtilen birimden metreye dönüştürülecek isim. Izgaranız coğrafiyse, dönüştürün
sağlayarak metrelere olan mesafeleri -fg yerine.

İLGİLİ HUSUSLAR

netCDF COARDS ızgaraları otomatik olarak coğrafi olarak tanınacaktır. Diğer ızgaralar için
dereceleri metreye dönüştürmek istediğiniz coğrafi ızgaraları seçin, -fg. eğer veri
her iki kutba da yakınsanız, ızgara dosyasını bir dikdörtgen üzerine yansıtmayı düşünmelisiniz.
kullanarak koordinat sistemi grd projesi

ÖRNEKLER

Mag_0.nc dosyasındaki deniz seviyesindeki manyetik anormallikleri 800 m seviyesine kadar yukarı doğru devam ettirmek için
Deniz seviyesinden yukarıda:

gmt grdfft mag_0.nc -C800 -V -Gmag_800.nc

Coğrafi bir ızgarada m (geoid.nc) cinsinden jeoid yüksekliklerini serbest hava yerçekimine dönüştürmek için
mGal'deki anormallikler:

gmt grdfft jeoid.nc -Dg -V -Ggrav.nc

mGal'deki (faa.nc) yerçekimi anormalliklerini dikey sapmalara dönüştürmek için (in
mikro-radyanlar) 038 yönünde, jeoid elde etmek için önce yerçekimini entegre etmeliyiz, sonra
yönlü türevi alın ve son olarak radyanları mikro radyanlara ölçeklendirin:

gmt grdfft faa.nc -Ig -A38 -S1e6 -V -Gdefl_38.nc

Yerçekimi anomalilerinin ikinci dikey türevleri, yerçekiminin eğriliği ile ilgilidir.
alan. Bunları iki kez farklılaştırarak mGal/m^2 olarak hesaplayabiliriz:

gmt grdfft yerçekimi.nc -D -D -V -Ggrav_2nd_derivative.nc

Ortak kayıtlı batimetri ve yerçekimi ızgaraları için çapraz spektral tahminleri hesaplamak ve
sonucu km cinsinden dalga boylarının fonksiyonları olarak rapor edin, deneyin

gmt grdfft bathymetry.nc yerçekimi.grd -Ewk -fg -V > cross_spectra.txt

Eğilim düşürme, nokta simetri yansıması ve sivrilmeden sonra FFT öncesi ızgarayı incelemek için
ham spektrumunun gerçek ve hayali bileşenlerini kaydetmenin yanı sıra
topo.nc'deki verileri deneyin

gmt grdfft topo.nc -N+w+z -fg -V

Artık topo_taper.nc, topo_real.nc ve topo_imag.nc'deki verilerin çizimlerini yapabilirsiniz.

onworks.net hizmetlerini kullanarak grdfftgmt'yi çevrimiçi kullanın