grdfiltergmt - Trực tuyến trên Đám mây

CHƯƠNG TRÌNH:

TÊN

grdfilter - Lọc lưới trong miền không gian (hoặc thời gian)

SYNOPSIS

bộ lọc lưới ingrid khoảng cách_cờ xchiều rộng[/chiều rộng2][bổ ngữ] khô khan [ tăng ] [ i|p|r ]
[ 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Ả

bộ lọc lưới sẽ lọc tệp lưới trong miền thời gian bằng cách sử dụng một trong các tích chập đã chọn
hoặc các bộ lọc đẳng hướng hoặc hình chữ nhật không tích chập và tính toán khoảng cách bằng cách sử dụng Descartes
hoặc Hình học hình cầu. Tệp lưới đầu ra có thể được tùy chọn tạo dưới dạng một vùng con
của đầu vào (thông qua -R) và / hoặc với mức tăng mới (thông qua -I) hoặc đăng ký (qua -T). Trong này
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ác cạnh sẽ không được sử dụng và
đầu ra có thể nằm trong một nửa chiều rộng của các cạnh đầu vào. Nếu bộ lọc thông qua thấp,
thì đầu ra có thể ít được lấy mẫu hơn đầu vào.

YÊU CẦU TRANH LUẬN

ingrid Tệp lưới của các điểm sẽ được lọc. (Xem CÁC ĐỊNH DẠNG TẬP TIN GRID bên dưới).

-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ờ = p: grid (px, py) với chiều rộng một số pixel lẻ; Khoảng cách Descartes.

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), Descartes
khoảng cách.

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.

cờ = 5: lưới (x, y) trong Mercator -Jm1 đơn vị img, chiều rộng tính bằng km, Khoảng cách hình cầu
phép tính.

-Fxchiều rộng[/chiều rộng2][bổ ngữ]
Đặt loại bộ lọc. Chọn trong số các bộ lọc tích chập và không tích chập. Sử dụng bất kỳ
mã bộ lọc x (liệt kê bên dưới) theo sau là đường kính đầy đủ chiều rộng. Điều này mang lại một
bộ lọc đẳng hướng; nối thêm /chiều rộng2 cho một bộ lọc hình chữ nhật (yêu cầu -Dp or -0).
Một số bộ lọc cho phép các đối số và công cụ sửa đổi tùy chọn.

Bộ lọc chuyển đổi (và mã của chúng) 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 đưa ra bởi hàm Gaussian, trong đó chiều rộng là 6 lần
Gaussian sigma truyền thống.

(f) Tùy chỉnh: Trọng số được cung cấp bởi các giá trị tính toán trước trong lưới trọng lượng bộ lọc
hồ sơ trọng lượng, mà phải có kích thước lẻ; cũng yêu cầu -0 và khoảng cách đầu ra
phải khớp với khoảng cách đầu vào hoặc là bội số nguyên.

(o) Toán tử: Trọng số được đưa ra bởi các giá trị được tính toán trước trong lưới trọng số của bộ lọc
hồ sơ trọng lượng, mà phải có kích thước lẻ; cũng yêu cầu -0 và khoảng cách đầu ra
phải khớp với khoảng cách đầu vào hoặc là bội số nguyên. Các trọng số được giả định là tổng thành
XNUMX để không có sự tích lũy tổng trọng lượng và quá trình chuẩn hóa sẽ được thực hiện.

Bộ lọc không chập (và mã của chúng) là:

(m) Median: Trả về giá trị trung bình. Để chọn một phần nối thêm lượng tử khác +qlượng tử in
phạm vi 0-1 [Mặc định là 0.5, tức là trung vị].

(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 +l or +u nếu bạn thích
muốn trả về giá trị thấp nhất hoặc trên cùng của các giá trị phương thức.

(h) Chế độ biểu đồ (một công cụ ước tính chế độ khác): Trả về giá trị phương thức làm trung tâm
của đỉnh vượt trội trong biểu đồ. Nối /băng thông để chỉ định binning
khoảng thời gian. Sử dụng công cụ sửa đổi +c để căn giữa các thùng trên bội số băng thông [Mặc định có
các cạnh thùng là bội số của băng thông]. Nếu 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 +l or +u nếu bạn muốn trả lại
dưới cùng hoặc trên cùng của các giá trị phương thức.

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

(L) Thấp hơn: Chỉ trả về giá trị tối thiểu của tất cả các giá trị dương.

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

(U) Trên: Chỉ trả về giá trị tối đa hoặc tất cả các giá trị âm.

Trong trường hợp L|U có thể không có dữ liệu nào vượt qua kiểm tra dấu hiệu ban đầu; trong
trường hợp đó bộ lọc sẽ trả về NaN.

-Gkhô khan
khô khan là tệp lưới đầu ra của bộ lọc. (Xem CÁC ĐỊNH DẠNG TẬP TIN GRID bên dưới).

CHỌN TRANH LUẬN

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

-Ni | p | r
Xác định cách các giá trị NaN trong lưới đầu vào ảnh hưởng đến đầu ra được lọc: Nối i đến
bỏ qua tất cả các NaN trong phép tính giá trị được lọc [Mặc định], r giống như i
ngoại trừ nếu nút đầu vào là NaN thì nút đầu ra sẽ được đặt thành NaN (chỉ
áp dụng nếu cả hai lưới được đồng đăng ký), và p điều này sẽ buộc giá trị được lọc
là NaN nếu tìm thấy bất kỳ nút lưới nào có giá trị NaN bên trong vòng tròn bộ lọc.

-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].

-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_etopo5.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
chiều rộng đầy đủ) 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 cứ 0.5 độ. Sử dụng tính toán khoảng cách hình cầu,
bạn cần:

gmt grdfilter North_pacific_etopo5.nc -Gfiltered_pacific.nc -Fm600 \
-D4 -R150 / 250/10/40 -I0.5 -V

Thay vào đó, nếu chúng ta muốn có kết quả vượt qua cao thì người ta có thể thực hiện kết quả vượt qua thấp tương ứng
bộ lọc sử dụng khoảng lưới thô vì grdfilter sẽ lấy mẫu lại kết quả giống nhau
độ phân giải dưới dạng lưới đầu vào để chúng tôi có thể tính toán phần dư, ví dụ:

gmt grdfilter North_pacific_etopo5.nc -Gresidual_pacific.nc -Fm-600 \
-D4 -R150 / 250/10/40 -I0.5 -V

Ở đây, lưới remod_pacific.nc sẽ có cùng độ phân giải 5 phút như ban đầu.

Để lọc tập dữ liệu trong ri gợn sóng.nc bằng cách sử dụng exp bộ lọc Gaussian dị hướng tùy chỉnh
(-0.5 * r ^ 2) có khoảng cách r từ tâm được cho bởi (2x ^ 2 + y ^ 2 -2xy) / 6, với chính
trục ở một góc 63 độ với phương ngang, hãy thử

gmt grdmath -R-10/10 / -10 / 10 -I1 X 2 POW 2 MUL Y 2 POW ADD XY MUL 2 MUL \
SUB 6 DIV NEG 2 DIV EXP DUP SUM DIV = gfilter.nc
gmt grdfilter ri gợn sóng.nc -Ffgfilter.nc -D0 -Gsmooth.nc -V

GIỚI HẠN

Để sử dụng -D5 tùy chọn lưới Mercator đầu vào phải được tạo bởi img2mercgrd bằng cách sử dụng -C
tùy chọn để gốc của các giá trị y là Xích đạo (tức là, x = y = 0 tương ứng với lon =
vĩ độ = 0).

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