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

CHƯƠNG TRÌNH:

TÊN

grdgradient - Tính toán đạo hàm định hướng hoặc gradient từ lưới

SYNOPSIS

độ dốc in_grdfile out_grdfile [ phương vị[/azim2]] [[a][c][o][n]] [
[s | p]azim / elev[/môi trường xung quanh/đồn/đặc biệt/tỏa sáng]] [ cờ ] [[e][t][amp] [/sigma[/bù đắp]]
] [ khu ] [ độ dốc ] [[cấp]] [ -fg ] [ -n]

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Ả

độ dốc có thể được sử dụng để tính đạo hàm có hướng theo một hướng nhất định (-A),
hoặc hướng (-S) [và độ lớn (-D)] của gradient vectơ của dữ liệu.

Giá trị ước tính trong hàng / cột đầu tiên / cuối cùng của đầu ra phụ thuộc vào các điều kiện biên (xem
-L).

YÊU CẦU TRANH LUẬN

in_grdfile
Tệp lưới 2-D để tính toán đạo hàm có hướng. (Xem CÁC ĐỊNH DẠNG TẬP TIN GRID
dưới đây).

-Gout_grdfile
Tên của tệp lưới đầu ra cho đạo hàm có hướng. (Xem CÁC ĐỊNH DẠNG TẬP TIN GRID
dưới đây).

CHỌN TRANH LUẬN

-Aphương vị[/azim2]
Hướng phương vị cho một đạo hàm có hướng; phương vị là góc theo x, y
máy bay đo theo độ dương theo chiều kim đồng hồ từ phía bắc (hướng + y) về phía
hướng đông (hướng + x). Phủ định của đạo hàm có hướng,
- [dz / dx * sin (phương vị) + dz / dy * cos (phương vị)], được tìm thấy; sự phủ định mang lại giá trị dương
khi dốc của z (x, y) xuống dốc trong phương vị hướng, ý nghĩa chính xác cho
che ánh sáng của hình ảnh (xem hình ảnh lưới và lưới xem) bởi một nguồn sáng
trên mặt phẳng x, y chiếu từ phương vị phương hướng. Tùy chọn, cung cấp hai
phương vị, -Aphương vị/azim2, trong trường hợp đó các gradient ở mỗi hướng này là
được tính toán và cái lớn hơn có độ lớn được giữ lại; điều này hữu ích cho
chiếu sáng dữ liệu với hai hướng của cấu trúc thẳng hàng, ví dụ: -A0/270
chiếu sáng từ phía bắc (trên cùng) và phía tây (bên trái).

-D [a] [c] [o] [n]
Tìm hướng của gradient dương (lên-dốc) của dữ liệu. Thay vào đó
tìm khía cạnh (hướng xuống dốc), sử dụng -Dạ. Theo mặc định, chỉ đường là
được đo theo chiều kim đồng hồ từ phía bắc, như phương vị in -A bên trên. Nối c sử dụng thông thường
Góc Descartes được đo ngược chiều kim đồng hồ so với chiều dương x (đông).
Nối o để báo cáo định hướng (0-180) chứ không phải hướng (0-360). Nối n
để thêm 90 độ cho tất cả các góc (ví dụ: để tạo ra các vết hằn cục bộ của bề mặt).

-E [s | p]azim / elev[/môi trường xung quanh/đồn/đặc biệt/tỏa sáng]
Tính toán độ tỏa sáng Lambertian thích hợp để sử dụng với hình ảnh lưới và lưới xem. Các
Lambertian Reflection giả định một bề mặt lý tưởng phản chiếu tất cả ánh sáng
đánh vào nó và bề mặt có vẻ sáng như nhau từ mọi hướng nhìn. phương vị
và độ cao là phương vị và độ cao của vectơ ánh sáng. Tùy chọn, cung cấp môi trường xung quanh
đồn đặc biệt tỏa sáng là các tham số kiểm soát các thuộc tính phản xạ
của bề mặt. Giá trị mặc định là: 0.55/0.6/0.4/10 Để lại một số giá trị
chưa được chạm, chỉ định = làm giá trị mới. Ví dụ -E60/30/=/0.5 thiết lập phương vị độ cao
và đồn đến 60, 30 và 0.5 và để lại các thông số phản xạ khác
không bị đụng chạm. Nối s để sử dụng một thuật toán Lambertian đơn giản hơn. Lưu ý rằng với biểu mẫu này
bạn chỉ phải cung cấp các tham số góc phương vị và độ cao. Nối p sử dụng
Xấp xỉ tuyến tính từng mảnh Peucker (thuật toán đơn giản hơn nhưng nhanh hơn; trong trường hợp này
các phương vị và độ cao được quấn cứng đến 315 và 45 độ. Điều này có nghĩa là ngay cả khi bạn
cung cấp các giá trị khác mà chúng sẽ bị bỏ qua.)

-Lcờ Tình trạng ranh giới cờ có lẽ x or y or xy cho biết dữ liệu là định kỳ trong phạm vi
x hoặc y hoặc cả hai, hoặc cờ có lẽ g cho biết điều kiện địa lý (x và y là
lon và lat). [Mặc định sử dụng các điều kiện "tự nhiên" (đạo hàm riêng cấp hai chuẩn
cạnh bằng XNUMX).]

-Mạng lưới][amp] [/sigma[/bù đắp]]
Bình thường hóa. [Mặc định: không chuẩn hóa.] Các gradient thực tế g được bù đắp và
được chia tỷ lệ để tạo ra các gradient chuẩn hóa gn với cường độ đầu ra tối đa là amp.
If amp không được cung cấp, mặc định amp = 1. Nếu bù đắp không được đưa ra, nó được đặt thành
Trung bình của g. -N sản lượng gn = amp * (g - bù đắp) / max (abs (g - bù đắp)). -E
chuẩn hóa bằng cách sử dụng phân phối Laplace tích lũy mang lại gn = amp * (1.0 -
exp (câu hỏi(2) * (g - bù đắp)/ sigma)) ở đâu sigma được ước tính bằng cách sử dụng định mức L1 của
(g - bù đắp) nếu nó không được đưa ra. -Không chuẩn hóa bằng cách sử dụng Cauchy tích lũy
năng suất phân phối gn = (2 * amp / PI) * atan ((g - bù đắp)/ sigma) Ở đâu sigma
được ước tính bằng cách sử dụng định mức L2 của (g - bù đắp) nếu nó không được đưa ra.

-NS[đơn vị]xmin/xmax/ymin/ymax[NS] (hơn ...)
Chỉ định khu vực quan tâm. Sử dụng -R tùy chọn sẽ chọn một phần phụ của
in_grdfile lưới điện. Nếu tiểu mục này vượt quá ranh giới của lưới, chỉ
vùng chung sẽ được trích xuất.

-Sđộ dốc
Tên của tệp lưới đầu ra với độ lớn vô hướng của vectơ gradient. Đòi hỏi -D
nhưng làm cho -G không bắt buộc.

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

-fg Lưới địa lý (kích thước kinh độ, vĩ độ) sẽ được chuyển đổi thành mét
thông qua phép gần đúng "Trái đất phẳng" sử dụng các thông số ellipsoid hiện tại.

-n [b | c | l | n] [+ a] [+ bBC] [+ c] [+ tngưỡng] (hơn ...)
Chọn chế độ nội suy cho lưới.

-^ 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 KHOẢNG CÁCH CÁC ĐƠN VỊ

Nếu lưới không có mét làm đơn vị nằm ngang, hãy nối thêm +uđơn vị đến tệp đầu vào
tên để chuyển đổi từ đơn vị xác định sang mét. Nếu lưới của bạn là địa lý, hãy chuyển đổi
khoảng cách đến mét bằng cách cung cấp -fg thay thế.

GỢI Ý

Nếu bạn không biết những gì -N các tùy chọn sử dụng để tạo tệp cường độ cho hình ảnh lưới or
lưới xem, một lần thử đầu tiên tốt là -E0.6.

Thông thường 255 sắc thái là quá đủ cho các mục đích trực quan. Bạn có thể tiết kiệm 75% đĩa
khoảng trắng bằng cách thêm = nb / a vào tên tệp đầu ra out_grdfile.

Nếu bạn muốn tạo một số bản đồ được chiếu sáng về các tiểu vùng của một tập dữ liệu lớn, và bạn
cần các hiệu ứng chiếu sáng nhất quán trên tất cả các bản đồ, hãy sử dụng -N lựa chọn và
cung cấp cùng một giá trị của sigma và bù đắp đến độ dốc cho mỗi bản đồ. Một dự đoán tốt là
bù đắp = 0 và sigma tìm thấy bởi thông tin lưới -L2 or -L1 được áp dụng cho grd gradient không chuẩn hóa.

Nếu bạn chỉ cần x- hoặc y- dẫn xuất của lưới, sử dụng môn toán.

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.

VÍ DỤ

Để tạo tệp để làm sáng dữ liệu trong geoid.nc bằng cách sử dụng gradient chuẩn hóa trong
phạm vi [-0.6,0.6] bắt chước các nguồn sáng theo hướng bắc và tây:

gmt grdgradient geoid.nc -A0 / 270 -Ggradients.nc = nb / a -Ne0.6 -V

Để tìm hướng phương vị của vải đáy biển trong tệp topo.nc:

gmt grdgradient topo.nc -Dno -Gazimuths.nc -V

THAM KHẢO

Horn, BKP, Hill-Shading and the Reflectance Map, Proceedings of the IEEE, Vol. 69, Không.
1, tháng 1981 năm 14, trang 47-XNUMX. (http://people.csail.mit.edu/bkph/papers/Hill-Shading.pdf)

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