dimfiltergmt - Trực tuyến trên đám mây

CHƯƠNG TRÌNH:

TÊN

dimfilter - Lọc theo hướng các tệp có lưới 2-D trong miền không gian (hoặc thời gian)

SYNOPSIS

máy lọc mờ input_file.nc khoảng cách_cờ [chế độ] đầu ra_file.nc
[ vòng cổ ] [ tăng ] [ khu ] [] [[cấp]] [ -f]

Lưu ý: Không được phép có khoảng trống giữa cờ tùy chọn và các đối số liên quan.

MÔ TẢ

máy lọc mờ sẽ lọc một .nc tệp trong miền không gian (hoặc thời gian) bằng cách chia
lọc vòng tròn thành n_sector, áp dụng một trong các tích chập chính đã chọn hoặc
bộ lọc không tích chập cho từng khu vực và chọn kết quả cuối cùng theo
bộ lọc thứ cấp đã chọn. Nó tính toán khoảng cách bằng cách sử dụng hình học Descartes hoặc hình cầu.
Đầu ra .nc tùy chọn tệp có thể được tạo như một tiểu vùng của đầu vào và / hoặc với
người mới -Increment. Bằng cách này, người ta có thể có "không gian thừa" trong dữ liệu đầu vào để có
sẽ không có hiệu ứng cạnh cho lưới đầu ra. Nếu bộ lọc là thông thấp, thì đầu ra
có thể được lấy mẫu ít thường xuyên hơn so với đầu vào. -Q là dành cho chế độ phân tích lỗi và chỉ
yêu cầu tổng số cột trong tệp đầu vào, có chứa
độ sâu. Cuối cùng, người ta nên biết rằng máy lọc mờ sẽ không tạo ra một đầu ra suôn sẻ như các
bộ lọc không gian làm được vì nó trả về giá trị trung bình tối thiểu trong số N phương tiện của N các ngành. Các
đầu ra có thể thô trừ khi dữ liệu đầu vào không bị nhiễu. Do đó, một bộ lọc bổ sung
(ví dụ: Gaussian qua bộ lọc lưới) của dữ liệu được lọc DiM thường được khuyến nghị.

YÊU CẦU TRANH LUẬN

input_file.nc
Lưới dữ liệu được lọc.

-Dkhoảng cách_cờ
Khoảng cách cờ cho biết lưới (x, y) liên quan như thế nào đến bộ lọc chiều rộng, như sau:

cờ = 0: lưới (x, y) cùng đơn vị với chiều rộng, Khoảng cách Descartes. cờ = 1: lưới
(x, y) tính bằng độ, chiều rộng tính bằng km, khoảng cách Descartes. cờ = 2: lưới (x, y)
tính bằng độ, chiều rộng tính bằng km, dx được chia tỷ lệ theo cos (giữa y), khoảng cách Descartes.

Các tùy chọn trên là nhanh nhất vì chúng chỉ cho phép tính toán ma trận trọng lượng
Một lần. Ba tùy chọn tiếp theo chậm hơn vì chúng tính toán lại trọng lượng cho mỗi
vĩ độ.

cờ = 3: lưới (x, y) tính bằng độ, chiều rộng tính bằng km, dx được chia tỷ lệ theo cosine (y), Descartes
tính toán khoảng cách.

cờ = 4: lưới (x, y) tính bằng độ, chiều rộng tính bằng km, tính khoảng cách hình cầu.

-F[chế độ]
Đặt loại bộ lọc chính. Chọn trong số các bộ lọc tích chập và không tích chập.
Nối mã bộ lọc theo sau là đường kính đầy đủ chiều rộng. Tích chập có sẵn
bộ lọc là:

(b) Boxcar: Tất cả các trọng lượng đều bằng nhau.

(c) Cosine Arch: Trọng lượng tuân theo một đường cong hình cung cosine.

(g) Gaussian: Trọng số được cho bởi hàm Gauss.

Bộ lọc không chập là:

(m) Median: Trả về giá trị trung bình.

(p) Xác suất khả năng tối đa (công cụ ước tính chế độ): Trả về giá trị phương thức. Nếu hơn
nhiều hơn một chế độ được tìm thấy, chúng tôi trả về giá trị trung bình của chúng. Nối - hoặc + vào bộ lọc
width nếu bạn muốn trả về giá trị nhỏ nhất hoặc lớn nhất trong các giá trị phương thức.

-N
Đặt loại bộ lọc thứ cấp và số lượng các cung dây nơ. n_sector phải
là số nguyên và lớn hơn 0. Khi n_sector được đặt thành 1, bộ lọc phụ là
không hiệu quả. Các bộ lọc phụ có sẵn là:

(l) Thấp hơn: Trả về giá trị nhỏ nhất của tất cả các giá trị đã lọc.

(u) Trên: Trả về giá trị tối đa của tất cả các giá trị đã lọc.

(a) Trung bình: Trả về giá trị trung bình của tất cả các giá trị đã lọc.

(m) Trung vị: Trả về giá trị trung bình của tất cả các giá trị đã lọc.

(p) Chế độ: Trả về chế độ của tất cả các giá trị đã lọc.

-Gđầu ra_file.nc
đầu ra_file.nc là đầu ra của bộ lọc.

CHỌN TRANH LUẬN

-I x_inc [và tùy chọn y_inc] là phần tăng đầu ra. Nối m để chỉ phút,
or c để chỉ giây. Nếu mới x_inc, y_inc KHÔNG phải là bội số nguyên của
những cái cũ (trong dữ liệu đầu vào), quá trình lọc sẽ chậm hơn đáng kể. [Mặc định: Giống nhau
làm đầu vào.]

-R tây, đông, namvà phía bắc xác định Vùng của các điểm đầu ra. [Vỡ nợ:
Giống như đầu vào.]

-T Chuyển đổi đăng ký nút cho lưới đầu ra để trở nên ngược lại với
lưới đầu vào [Mặc định cho đăng ký giống như lưới đầu vào].

-Qvòng cổ vòng cổ là tổng số cột trong tệp bảng văn bản đầu vào. Đối với chế độ này,
nó mong đợi để đọc độ sâu bao gồm một số cột. Mỗi cột đại diện cho một
lưới lọc với chiều rộng bộ lọc, có thể lấy bằng grd2xyz -Z. Kết quả
sẽ là trung vị, MAD và trung bình. Vì vậy, cột có phương tiện trung gian được sử dụng để tạo
thành phần khu vực và cột có MAD được sử dụng để xử lý lỗi
phân tích.

-V [cấp] (hơn ...)
Chọn mức độ chi tiết [c].

-f [i | o]thông tin về colinfo (hơn ...)
Chỉ định kiểu dữ liệu của cột đầu vào và / hoặc cột đầu ra.

-^ or chỉ -
In một thông báo ngắn về cú pháp của lệnh, sau đó thoát ra (LƯU Ý: trên Windows
chỉ sử dụng -).

-+ or chỉ +
In thông báo sử dụng rộng rãi (trợ giúp), bao gồm giải thích về bất kỳ
tùy chọn dành riêng cho mô-đun (nhưng không phải tùy chọn chung theo GMT), sau đó thoát.

-? or Không đối số
In thông báo sử dụng hoàn chỉnh (trợ giúp), bao gồm giải thích về các tùy chọn, sau đó
lối thoát hiểm.

--phiên bản
In phiên bản GMT và thoát.

--show-datadir
In đường dẫn đầy đủ đến thư mục chia sẻ GMT và thoát.

GRID FILE M FORU ĐƠN

Theo mặc định, GMT ghi ra lưới khi độ chính xác duy nhất nổi trong netCDF khiếu nại COARDS
định dạng tệp. Tuy nhiên, GMT có thể tạo các tệp lưới trong nhiều lưới thông dụng khác
định dạng tệp và cũng tạo điều kiện cho cái gọi là "đóng gói" lưới, viết ra dấu phẩy động
dữ liệu dưới dạng số nguyên 1 hoặc 2 byte. Để chỉ định độ chính xác, tỷ lệ và độ lệch, người dùng nên
thêm hậu tố =id[/tỉ lệ/bù đắp[/nan]], ở đâu id là một mã định danh gồm hai chữ cái của lưới
loại và độ chính xác, và tỉ lệ và bù đắp là hệ số tỷ lệ tùy chọn và bù đắp
được áp dụng cho tất cả các giá trị lưới và nan là giá trị được sử dụng để chỉ ra dữ liệu bị thiếu. Trong trường hợp
hai nhân vật id không được cung cấp, như trong =/tỉ lệ hơn một id=nf được giả định. Khi nào
lưới đọc, định dạng thường được tự động nhận dạng. Nếu không, cùng một hậu tố
có thể được thêm vào tên tệp lưới đầu vào. Nhìn thấy chuyển đổi lưới và phần lưới-tệp-định dạng của
Tài liệu tham khảo Kỹ thuật và Sách dạy nấu ăn GMT để biết thêm thông tin.

Khi đọc tệp netCDF có chứa nhiều lưới, theo mặc định, GMT sẽ đọc
lưới 2 chiều đầu tiên có thể tìm thấy trong tệp đó. Để thu hút GMT đọc một
biến đa chiều trong tệp lưới, nối thêm ?tên sơn dầu đến tên tệp, ở đâu
tên sơn dầu là tên của biến. Lưu ý rằng bạn có thể cần phải thoát khỏi ý nghĩa đặc biệt
of ? trong chương trình shell của bạn bằng cách đặt dấu gạch chéo ngược trước nó hoặc bằng cách đặt
tên tệp và hậu tố giữa dấu ngoặc kép hoặc dấu ngoặc kép. Các ?tên sơn dầu hậu tố cũng có thể được sử dụng
cho lưới đầu ra để chỉ định một tên biến khác với tên mặc định: "z". Nhìn thấy
chuyển đổi lưới và các công cụ sửa đổi Phần-cho-CF và định dạng tệp lưới của GMT Technical
Tham khảo và Sách dạy nấu ăn để biết thêm thông tin, đặc biệt là về cách đọc các mối nối của 3,
Lưới 4 - hoặc 5 chiều.

ĐỊA LÝ VÀ THỜI GIAN COORDINATE

Khi loại lưới đầu ra là netCDF, tọa độ sẽ được gắn nhãn "kinh độ",
"vĩ độ" hoặc "thời gian" dựa trên các thuộc tính của dữ liệu đầu vào hoặc lưới (nếu có) hoặc trên
-f or -R tùy chọn. Ví dụ, cả hai -f0x -f1t và -R90w / 90e / 0t / 3t sẽ dẫn đến một
lưới kinh độ / thời gian. Khi tọa độ x, y hoặc z là thời gian, nó sẽ được lưu trữ trong lưới
là thời gian tương đối kể từ kỷ nguyên được chỉ định bởi TIME_UNIT và TIME_EPOCH trong gmt.conf hồ sơ
hoặc trên dòng lệnh. Ngoài ra, đơn vị thuộc tính của biến thời gian sẽ cho biết
cả đơn vị này và kỷ nguyên.

VÍ DỤ

Giả sử rằng north_pacific_dbdb5.nc là một tệp có độ sâu 5 phút từ 140E đến 260E và
0N đến 50N và bạn muốn tìm các giá trị trung bình trong bán kính 300km (600km đầy
chiều rộng) của các điểm đầu ra, mà bạn chọn là từ 150E đến 250E và 10N đến 40N, và
bạn muốn các giá trị đầu ra mỗi 0.5 độ. Để ngăn các phương tiện không bị sai lệch bởi
mặt phẳng dốc, bạn muốn chia vòng tròn bộ lọc thành 6 cung và chọn
giá trị thấp nhất trong số 6 giá trị trung bình. Sử dụng tính toán khoảng cách hình cầu, bạn cần:

gmt dimfilter North_pacific_dbdb5.nc -Gfiltered_pacific.nc -Fm600 -D4 \
-Nl6 -R150 / 250/10/40 -I0.5 -V

Giả sử rằng cape_verde.nc là một tệp có độ sâu 0.5 phút từ 32W đến 15W và 8N đến
25N và bạn muốn xóa các đối tượng địa lý có quy mô độ dài nhỏ để xác định độ phồng trong
vùng mở rộng từ 27.5W đến 20.5W và 12.5N đến 19.5N và bạn muốn giá trị đầu ra mỗi
2 phút. Bằng cách sử dụng tính toán khoảng cách theo trường phái Cartesian, bạn cần:

gmt dimfilter 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

Giả sử rằng bạn đã tìm thấy một loạt các chiều rộng bộ lọc cho một khu vực nhất định và bạn đã lọc
dữ liệu đo độ sâu đã cho bằng cách sử dụng phạm vi chiều rộng bộ lọc (ví dụ: f100.nc f110.nc f120.nc
f130.nc) và bạn muốn xác định xu hướng khu vực bằng cách sử dụng phạm vi chiều rộng bộ lọc và
bạn muốn có được ước tính độ lệch tuyệt đối trung bình (MAD) tại mỗi điểm dữ liệu. Sau đó, bạn
sẽ cần phải làm:

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
dán f100.d f110.d f120.d f130.d> độ sâu.d
gmt độ sâu bộ lọc mờ. d -Q4> output.z

GIỚI HẠN

Khi làm việc với lưới địa lý (vĩ độ, vĩ độ), cả ba bộ lọc tích chập (boxcar,
cosine Arch và gaussian) sẽ chuẩn hóa đúng các trọng số của bộ lọc cho biến thể trong
kích thước hộp lưới với vĩ độ và xác định chính xác những nút nào cần thiết cho
tích chập khi bộ lọc "vòng tròn" vượt qua ranh giới tuần hoàn (0-360) hoặc chứa
cực địa lý. Tuy nhiên, các bộ lọc không gian, chẳng hạn như bộ lọc chế độ và trung vị, không sử dụng
trọng số và do đó chỉ nên được sử dụng trên lưới Descartes (hoặc ở vĩ độ rất thấp) mà thôi.
Nếu bạn muốn áp dụng các bộ lọc không gian như vậy, bạn nên chiếu dữ liệu của mình lên một vùng bằng
chiếu và chạy bộ lọc mờ trên lưới Descartes kết quả.

SCRIPT BẢN MẪU

Dim.template.sh là một tập lệnh khung sườn có thể được sử dụng để thiết lập một DiM hoàn chỉnh
phân tích, bao gồm cả phân tích MAD.

THAM KHẢO

Kim, S.-S. và Wessel, P. (2008), Lọc trung vị theo hướng cho phần dư theo vùng
Tách Bathymetry, Địa hóa học. Địa vật lý. hệ thống địa chất., 9, Q03005,
doi: 10.1029 / 2007GC001850.

Sử dụng dimfiltergmt trực tuyến bằng các dịch vụ onworks.net