dimfiltergmt - ออนไลน์ในคลาวด์

โครงการ:

ชื่อ

dimfilter - การกรองทิศทางของไฟล์กริด 2 มิติในโดเมนอวกาศ (หรือเวลา)

เรื่องย่อ

ตัวกรองแสง input_file.nc Distance_flag [โหมด] output_file.nc
[ ปลอกคอ ] [ เพิ่มขึ้น ] [ ภูมิภาค -ระดับ- -f]

หมายเหตุ ไม่อนุญาตให้มีช่องว่างระหว่างแฟล็กตัวเลือกและอาร์กิวเมนต์ที่เกี่ยวข้อง

DESCRIPTION

ตัวกรองแสง จะกรอง a .nc ไฟล์ในโดเมนช่องว่าง (หรือเวลา) โดยหารที่กำหนด
กรองวงกลมเข้า n_ภาค, ใช้หนึ่งใน convolution หลักที่เลือกหรือ
ตัวกรองที่ไม่บิดเบี้ยวในแต่ละภาคส่วนและเลือกผลลัพธ์สุดท้ายตาม
ตัวกรองรองที่เลือก มันคำนวณระยะทางโดยใช้เรขาคณิตคาร์ทีเซียนหรือทรงกลม
ผลลัพธ์ .nc file สามารถสร้างเป็น subregion ของ input และ/หรือ a
ใหม่ -Iเพิ่มขึ้น ด้วยวิธีนี้อาจมี "ช่องว่างเพิ่มเติม" ในข้อมูลที่ป้อนเข้าเพื่อให้มี
จะไม่มีเอฟเฟกต์ขอบสำหรับกริดเอาต์พุต หากตัวกรองเป็นแบบ low-pass แสดงว่าเอาต์พุต
อาจสุ่มตัวอย่างน้อยกว่าอินพุต -Q ใช้สำหรับโหมดการวิเคราะห์ข้อผิดพลาดและเท่านั้น
ต้องการจำนวนคอลัมน์ทั้งหมดในไฟล์อินพุตซึ่งมีตัวกรอง
ความลึก สุดท้ายนี้ควรจะรู้ว่า ตัวกรองแสง จะไม่ให้ผลผลิตที่ราบรื่นเหมือนอื่นๆ
ตัวกรองเชิงพื้นที่ทำเพราะมันส่งกลับค่ามัธยฐานขั้นต่ำจาก N ค่ามัธยฐานของ N ภาค NS
เอาต์พุตอาจหยาบเว้นแต่ข้อมูลอินพุตจะปราศจากสัญญาณรบกวน ดังนั้นการกรองเพิ่มเติม
(เช่น Gaussian via ตัวกรองกริด) ของข้อมูลที่กรอง DiM โดยทั่วไปแนะนำ

ที่จำเป็น อาร์กิวเมนต์

input_file.nc
ตารางข้อมูลที่จะกรอง

-DDistance_flag
ระยะทาง ธง บอกว่ากริด (x,y) เกี่ยวข้องกับตัวกรองอย่างไร ความกว้างดังนี้

ธง = 0: กริด (x,y) หน่วยเดียวกับ ความกว้าง, ระยะทางคาร์ทีเซียน. ธง = 1: ตาราง
(x,y) เป็นองศา ความกว้าง เป็นกิโลเมตร ระยะทางคาร์ทีเซียน ธง = 2: ตาราง (x,y)
เป็นองศา ความกว้าง ในหน่วยกม. dx ถูกขยายโดย cos (กลาง y) ระยะทางคาร์ทีเซียน

ตัวเลือกด้านบนนั้นเร็วที่สุดเพราะอนุญาตให้คำนวณเมทริกซ์น้ำหนักเท่านั้น
ครั้งหนึ่ง. สามตัวเลือกถัดไปจะช้ากว่าเพราะจะคำนวณน้ำหนักใหม่สำหรับแต่ละตัวเลือก
ละติจูด.

ธง = 3: ตาราง (x,y) เป็นองศา ความกว้าง ในหน่วยกม. dx ที่ปรับขนาดโดยโคไซน์ (y) คาร์ทีเซียน
การคำนวณระยะทาง

ธง = 4: ตาราง (x,y) เป็นองศา ความกว้าง เป็นกม. การคำนวณระยะทางทรงกลม

-F[โหมด]
ตั้งค่าประเภทตัวกรองหลัก เลือกระหว่างตัวกรองแบบบิดและไม่บิดเบี้ยว
ต่อท้ายรหัสตัวกรองตามด้วยเส้นผ่านศูนย์กลางเต็ม ความกว้าง. การบิดที่มีอยู่
ตัวกรองคือ:

(b) Boxcar: น้ำหนักทั้งหมดเท่ากัน

(c) โคไซน์อาร์ค: ตุ้มน้ำหนักเป็นไปตามเส้นโค้งโคไซน์อาร์ค

(g) เกาส์เซียน: น้ำหนักถูกกำหนดโดยฟังก์ชันเกาส์เซียน

ตัวกรองแบบไม่บิดคือ:

(m) ค่ามัธยฐาน: ส่งกลับค่ามัธยฐาน

(p) ความน่าจะเป็นที่เป็นไปได้สูงสุด (ตัวประมาณโหมด): ส่งกลับค่าโมดอล ถ้ามากกว่านี้
พบมากกว่าหนึ่งโหมดที่เราส่งคืนค่าเฉลี่ย ต่อท้าย - หรือ + กับตัวกรอง
width หากคุณต้องการส่งคืนค่ากิริยาที่เล็กที่สุดหรือใหญ่ที่สุด

-N
ตั้งค่าประเภทตัวกรองรองและจำนวนส่วนหูกระต่าย n_ภาค ต้อง
เป็นจำนวนเต็มและมากกว่า 0 เมื่อ n_ภาค ถูกตั้งค่าเป็น 1 ตัวกรองรองคือ
ไม่มีประสิทธิภาพ ตัวกรองรองที่มีอยู่คือ:

(l) ต่ำกว่า: ส่งคืนค่าต่ำสุดของค่าที่กรองทั้งหมด

(u) บน: ส่งกลับค่าสูงสุดของค่าที่กรองทั้งหมด

(a) ค่าเฉลี่ย: ส่งกลับค่าเฉลี่ยของค่าที่กรองทั้งหมด

(m) ค่ามัธยฐาน: ส่งกลับค่ามัธยฐานของค่าที่กรองทั้งหมด

(p) โหมด: กลับโหมดของค่าที่กรองทั้งหมด

-Goutput_file.nc
output_file.nc คือเอาต์พุตของตัวกรอง

ตัวเลือก อาร์กิวเมนต์

-I x_inc [และเป็นทางเลือก y_inc] คือการเพิ่มเอาต์พุต ผนวก m เพื่อระบุนาที
or c เพื่อระบุวินาที ถ้าใหม่ x_inc, y_inc ไม่ใช่ผลคูณของจำนวนเต็มของ
เก่า (ในข้อมูลที่ป้อน) การกรองจะช้ากว่ามาก [ค่าเริ่มต้น: Same
เป็นข้อมูลเข้า]

-R ทิศตะวันตก, ทางทิศตะวันออก, ใต้และ ทางทิศเหนือ กำหนดขอบเขตของจุดส่งออก [ค่าเริ่มต้น:
เหมือนกับอินพุต.]

-T สลับการลงทะเบียนโหนดสำหรับกริดเอาต์พุตเพื่อให้ตรงข้ามกับ
ตารางอินพุต [ค่าเริ่มต้นให้การลงทะเบียนเหมือนกับกริดอินพุต]

-Qปลอกคอ ปลอกคอ คือจำนวนคอลัมน์ทั้งหมดในไฟล์ตารางข้อความที่ป้อน สำหรับโหมดนี้
มันคาดว่าจะอ่านความลึกประกอบด้วยหลายคอลัมน์ แต่ละคอลัมน์แสดงถึงa
ตะแกรงกรองที่มีความกว้างของตัวกรองซึ่งสามารถหาได้โดย grd2xyz -Z. ผลลัพธ์
จะเป็นค่ามัธยฐาน MAD และค่าเฉลี่ย ดังนั้น คอลัมน์ที่มีค่ามัธยฐานจึงถูกใช้เพื่อสร้าง
ส่วนประกอบภูมิภาคและคอลัมน์ที่มี MADs ใช้เพื่อดำเนินการผิดพลาด
การวิเคราะห์

-วี[ระดับ] (มากกว่า ... )
เลือกระดับการใช้คำฟุ่มเฟือย [c]

-ฉ[ฉัน|o]คอลอินโฟ (มากกว่า ... )
ระบุชนิดข้อมูลของคอลัมน์อินพุตและ/หรือเอาต์พุต

-^ or เพียงแค่ -
พิมพ์ข้อความสั้น ๆ เกี่ยวกับไวยากรณ์ของคำสั่ง จากนั้นออก (หมายเหตุ: บน Windows
ใช้เพียง -).

-+ or เพียงแค่ +
พิมพ์ข้อความการใช้งาน (ช่วยเหลือ) อย่างกว้างขวาง รวมถึงคำอธิบายของทุก ๆ
ตัวเลือกเฉพาะโมดูล (แต่ไม่ใช่ตัวเลือกทั่วไปของ GMT) จากนั้นออก

-? or ไม่ ข้อโต้แย้ง
พิมพ์ข้อความการใช้งาน (วิธีใช้) ที่สมบูรณ์ รวมทั้งคำอธิบายตัวเลือก จากนั้น
ทางออก

--รุ่น
พิมพ์เวอร์ชัน GMT และออก

--show-datadir
พิมพ์เส้นทางแบบเต็มไปยังไดเรกทอรีแชร์ GMT และออก

GRID ไฟล์ รูปแบบ

โดยค่าเริ่มต้น GMT จะเขียนกริดออกเป็นความแม่นยำเดียวใน netCDF . ร้องเรียน COARDS
รูปแบบไฟล์. อย่างไรก็ตาม GMT สามารถสร้างไฟล์กริดในกริดอื่น ๆ ที่ใช้กันทั่วไปได้
รูปแบบไฟล์และยังอำนวยความสะดวกที่เรียกว่า "การบรรจุ" ของกริดการเขียนจุดลอยตัว
ข้อมูลเป็นจำนวนเต็ม 1 หรือ 2 ไบต์ ในการระบุความแม่นยำ มาตราส่วน และออฟเซ็ต ผู้ใช้ควร
เพิ่มคำต่อท้าย =id[/ขนาด/ชดเชย[/น่าน]], ที่ไหน id เป็นตัวระบุสองตัวอักษรของกริด
ชนิดและความแม่นยำ และ ขนาด และ ชดเชย เป็นสเกลแฟกเตอร์เสริมและออฟเซ็ตเป็น
นำไปใช้กับค่ากริดทั้งหมดและ น่าน คือค่าที่ใช้ระบุข้อมูลที่ขาดหายไป เผื่อ
ตัวละครทั้งสอง id ไม่ได้ระบุไว้เช่นใน =/ขนาด กว่า a id=nf จะถือว่า เมื่อไหร่
อ่านตาราง รูปแบบโดยทั่วไปจะรู้จักโดยอัตโนมัติ ถ้าไม่ใช่ คำต่อท้ายเดียวกัน
สามารถเพิ่มลงในชื่อไฟล์กริดอินพุต ดู แปลงกริด และรูปแบบไฟล์กริดมาตราของ
ข้อมูลอ้างอิงทางเทคนิค GMT และตำราอาหารสำหรับข้อมูลเพิ่มเติม

เมื่ออ่านไฟล์ netCDF ที่มีหลายกริด GMT จะอ่านตามค่าเริ่มต้น
ตาราง 2 มิติแรกที่หาได้ในไฟล์นั้น เพื่อเกลี้ยกล่อม GMT ให้อ่านอีก
ตัวแปรหลายมิติในไฟล์กริด ผนวก ?นามสกุล ไปที่ชื่อไฟล์โดยที่
นามสกุล เป็นชื่อของตัวแปร โปรดทราบว่าคุณอาจต้องหนีความหมายพิเศษ
of ? ในโปรแกรมเชลล์ของคุณโดยใส่แบ็กสแลชไว้ข้างหน้า หรือโดยการวาง
ชื่อไฟล์และคำต่อท้ายระหว่างเครื่องหมายคำพูดหรือเครื่องหมายคำพูดคู่ NS ?นามสกุล ต่อท้ายก็ใช้ได้นะ
สำหรับกริดเอาต์พุตเพื่อระบุชื่อตัวแปรที่แตกต่างจากค่าเริ่มต้น: "z" ดู
แปลงกริด และส่วนดัดแปลงสำหรับ CF และรูปแบบไฟล์กริดของ GMT Technical
อ้างอิงและตำราอาหารสำหรับข้อมูลเพิ่มเติม โดยเฉพาะอย่างยิ่งในการอ่าน splices ของ 3-,
กริด 4 หรือ 5 มิติ

ภูมิศาสตร์ AND เวลา พิกัด

เมื่อประเภทกริดเอาท์พุตเป็น netCDF พิกัดจะมีป้ายกำกับว่า "ลองจิจูด"
"ละติจูด" หรือ "เวลา" ตามคุณสมบัติของข้อมูลที่ป้อนหรือกริด (ถ้ามี) หรือตาม
-f or -R ตัวเลือก. ตัวอย่างเช่น ทั้งสอง -f0x -f1t และ -R90w/90e/0t/3t จะส่งผลให้ a
เส้นลองจิจูด/เวลา เมื่อพิกัด x, y หรือ z เป็นเวลา มันจะถูกเก็บไว้ในกริด
เป็นเวลาสัมพัทธ์ตั้งแต่ยุคตามที่ระบุโดย TIME_UNIT และ TIME_EPOCH ใน gmt.conf ไฟล์
หรือบนบรรทัดคำสั่ง นอกจากนี้ หน่วย คุณลักษณะของตัวแปรเวลาจะระบุ
ทั้งหน่วยนี้และยุค

ตัวอย่าง

สมมติว่า north_pacific_dbdb5.nc เป็นไฟล์การวัดความยาว 5 นาทีจาก 140E ถึง 260E และ
0N ถึง 50N และคุณต้องการหาค่ามัธยฐานของค่าภายในรัศมี 300 กม. (600 กม. เต็ม
ความกว้าง) ของจุดเอาท์พุต ซึ่งคุณสามารถเลือกได้ตั้งแต่ 150E ถึง 250E และ 10N ถึง 40N และ
คุณต้องการค่าเอาต์พุตทุกๆ 0.5 องศา เพื่อป้องกันไม่ให้ค่ามัธยฐานถูกลำเอียงโดย
ระนาบลาดเอียง คุณต้องการแบ่งวงกลมกรองออกเป็น 6 ส่วน แล้วเลือก
ค่าต่ำสุดระหว่าง 6 ค่ามัธยฐาน เมื่อใช้การคำนวณระยะทรงกลม คุณจะต้อง:

gmt หรี่ตัวกรอง north_pacific_dbdb5.nc -Gfiltered_pacific.nc -Fm600 -D4 \
-Nl6 -R150/250/10/40 -I0.5 -V

สมมติว่า Cape_verde.nc เป็นไฟล์การวัดขนาด 0.5 นาที จาก 32W ถึง 15W และ 8N ถึง
25N และคุณต้องการลบคุณลักษณะขนาดเล็กเพื่อกำหนดการขยายใน
พื้นที่ขยายจาก 27.5W เป็น 20.5W และ 12.5N ถึง 19.5N และคุณต้องการค่าเอาต์พุตทุกๆ
2 นาที การใช้การคำนวณระยะทางคาร์ทีเซียนคุณต้อง:

gmt ตัวกรองหรี่ cape_verde.nc -Gt.nc -Fm220 -Nl8 -D2 -R-27.5/-20.5/12.5/19.5 -I2m -V
gmt grdfilter t.nc -Gcape_swell.nc -Fg50 -D2 -V

สมมติว่าคุณพบช่วงความกว้างของตัวกรองสำหรับพื้นที่ที่กำหนด และคุณกรอง
กำหนดข้อมูล Bathymetric โดยใช้ช่วงความกว้างของตัวกรอง (เช่น f100.nc f110.nc f120.nc
f130.nc) และคุณต้องการกำหนดแนวโน้มระดับภูมิภาคโดยใช้ช่วงความกว้างของตัวกรอง และ
คุณต้องการรับค่าประมาณการเบี่ยงเบนสัมบูรณ์มัธยฐาน (MAD) ที่จุดข้อมูลแต่ละจุด แล้วคุณ
จะต้องทำ:

gmt grd2xyz f100.nc -Z > f100.d
gmt grd2xyz f110.nc -Z > f110.d
gmt grd2xyz f120.nc -Z > f120.d
gmt grd2xyz f130.nc -Z > f130.d
วาง f100.d f110.d f120.d f130.d > ความลึก.d
gmt dimfilter depths.d -Q4 > output.z

ข้อ จำกัด

เมื่อทำงานกับกริดทางภูมิศาสตร์ (lat, lon) ตัวกรองการบิดทั้งสามตัว (boxcar,
โคไซน์อาร์คและเกาส์เซียน) จะทำให้น้ำหนักตัวกรองเป็นปกติสำหรับการเปลี่ยนแปลงใน
ขนาดกริดบ็อกซ์พร้อมละติจูด และกำหนดอย่างถูกต้องว่าโหนดใดที่จำเป็นสำหรับ
การบิดเบี้ยวเมื่อตัวกรอง "วงกลม" ข้ามขอบเขตเป็นระยะ (0-360) หรือมี
เสาทางภูมิศาสตร์ อย่างไรก็ตาม ตัวกรองเชิงพื้นที่ เช่น ตัวกรองค่ามัธยฐานและโหมด ห้ามใช้
น้ำหนักจึงควรใช้กับกริดคาร์ทีเซียนเท่านั้น (หรือที่ละติจูดต่ำมาก) เท่านั้น
หากคุณต้องการใช้ตัวกรองเชิงพื้นที่ดังกล่าว คุณควรฉายข้อมูลของคุณไปยังพื้นที่ที่เท่ากัน
ฉายภาพและรัน dimfilter บนกริดคาร์ทีเซียนที่เป็นผลลัพธ์

สคริปต์ เทมเพลต

dim.template.sh เป็นสคริปต์เชลล์โครงร่างที่สามารถใช้เพื่อตั้งค่า DiM . ที่สมบูรณ์
การวิเคราะห์รวมทั้งการวิเคราะห์ MAD

อ้างอิง

Kim, S.-S. และ Wessel, P. (2008) Directional Median Filtering for Regional-Residual
การแยก Bathymetry, จีโอเคม. ธรณีฟิสิกส์ จีโอซิสท์, 9,Q03005,
ดอย:10.1029/2007GC001850.

ใช้ dimfiltergmt ออนไลน์โดยใช้บริการ onworks.net