i.pansharpengrass - Trực tuyến trên đám mây

CHƯƠNG TRÌNH:

TÊN

i.pansharpen - Các thuật toán tổng hợp hình ảnh để làm sắc nét đa mặt với độ phân giải cao
kênh panchromatic

TỪ KHÓA

hình ảnh, kết hợp, làm sắc nét, Brovey, IHS, HIS, PCA

SYNOPSIS

i.pansharpen
i.pansharpen --Cứu giúp
i.pansharpen [-sl] đỏ=tên màu xanh lá=tên màu xanh da trời=tên pan=tên đầu ra=tên cơ sở phương pháp=chuỗi
[-ghi đè] [-giúp đỡ] [-dài dòng] [-yên tĩnh] [-ui]

Cờ:
-s
Xử lý nối tiếp thay vì xử lý song song

-l
Cân bằng lại kênh màu xanh lam cho LANDSAT

- ghi đè
Cho phép các tệp đầu ra ghi đè lên các tệp hiện có

--Cứu giúp
In tóm tắt sử dụng

--dài dòng
Đầu ra mô-đun dài dòng

--Yên lặng
Đầu ra mô-đun yên tĩnh

--ui
Buộc khởi chạy hộp thoại GUI

Tham số:
đỏ=tên [yêu cầu]
Tên của bản đồ raster được sử dụng cho

màu xanh lá=tên [yêu cầu]
Tên của bản đồ raster được sử dụng cho

màu xanh da trời=tên [yêu cầu]
Tên của bản đồ raster được sử dụng cho

pan=tên [yêu cầu]
Tên của bản đồ raster được sử dụng cho kênh panchromatic có độ phân giải cao

đầu ra=tên cơ sở [yêu cầu]
Tên cho (các) bản đồ raster basename đầu ra

phương pháp=chuỗi [yêu cầu]
Phương pháp mài chảo
Tùy chọn: người lớn, tôi, pca
Mặc định: ihs

MÔ TẢ

i.pansharpen sử dụng dải panchromatic có độ phân giải cao từ hình ảnh đa kính để
làm sắc nét 3 dải phân giải thấp hơn. Sau đó, 3 băng tần có độ phân giải thấp hơn có thể được kết hợp thành
hình ảnh màu RGB ở độ phân giải cao hơn (chi tiết hơn) so với mức có thể bằng cách sử dụng
3 dải ban đầu. Ví dụ, Landsat ETM có dải quang phổ độ phân giải thấp 1 (xanh lam), 2
(xanh lục), 3 (đỏ), 4 (gần IR), 5 (giữa IR) và 7 (giữa IR) ở độ phân giải 30m và cao
độ phân giải dải thơm 8 ở độ phân giải 15m. Độ sắc nét của chảo cho phép các dải tần 3-2-1 (hoặc
các kết hợp khác của dải phân giải 30m như 4-3-2 hoặc 5-4-2) để được kết hợp thành 15m
độ phân giải hình ảnh màu.
i.pansharpen cung cấp lựa chọn ba thuật toán 'làm sắc nét bằng chảo' khác nhau: IHS, Brovey,
và PCA.
Trong cáp IHS pan mài, 3 dải ban đầu có độ phân giải thấp hơn, được chọn là đỏ, xanh lá cây và
các kênh màu xanh lam để tạo hình ảnh tổng hợp RGB, được chuyển thành IHS (cường độ,
màu sắc và độ bão hòa) không gian màu. Dải panchromatic sau đó được thay thế cho
kênh cường độ (I), kết hợp với các kênh màu gốc (H) và bão hòa (S), và
được chuyển đổi trở lại không gian màu RGB ở độ phân giải cao hơn của dải màu panchromatic. Các
thuật toán cho điều này có thể được biểu diễn dưới dạng: RGB -> IHS -> [pan] HS -> RGB.
Với một Brovey pan mài, mỗi dải trong số 3 dải phân giải thấp hơn và dải panchromatic
được kết hợp bằng cách sử dụng thuật toán sau để tính toán 3 băng tần mới ở mức cao hơn
độ phân giải (ví dụ cho băng tần 1):
ban nhạc 1
ban nhạc mới1 = ----------------------- * panband
dải1 + dải2 + dải3
In PCA pan mài, một phân tích thành phần chính được thực hiện trên 3 thấp hơn ban đầu
các dải phân giải để tạo ra 3 hình ảnh thành phần chính (PC1, PC2 và PC3) và
các eigenvectors liên kết (EV), như vậy:
ban nhạc1 ban nhạc2 ban nhạc3
PC1: EV1-1 EV1-2 EV1-3
PC2: EV2-1 EV2-2 EV2-3
PC3: EV3-1 EV3-2 EV3-3
và
PC1 = EV1-1 * band1 + EV1-2 * band2 + EV1-3 * band3 - trung bình (dải 1,2,3)
Một PCA nghịch đảo sau đó được thực hiện, thay thế dải panchromatic cho PC1. Để làm điều này,
ma trận eigenvectors được đảo ngược (trong trường hợp này là chuyển vị), hình ảnh PC được
được nhân với các eigenvector với dải panchromatic được thay thế cho PC1, và giá trị trung bình của
mỗi dải được thêm vào mỗi dải hình ảnh đã biến đổi bằng cách sử dụng thuật toán sau (ví dụ
cho ban nhạc 1):
band1 '= pan * EV1-1 + PC2 * EV2-1 + PC3 * EV3-1 + trung bình (band1)
Việc phân bổ các kênh phụ thuộc vào vệ tinh. Ví dụ về hình ảnh vệ tinh
với các dải quang phổ có độ phân giải cao và các dải quang phổ có độ phân giải thấp hơn bao gồm
Landsat 7 ETM, QuickBird và SPOT.

GHI CHÚ

Mô-đun hiện chỉ hoạt động đối với hình ảnh 8-bit.
Lệnh tạm thời thay đổi vùng tính toán thành độ phân giải cao của
dải panchromatic trong quá trình tính toán làm sắc nét, sau đó khôi phục vùng trước đó
cài đặt. Tọa độ vùng hiện tại (và giá trị rỗng) được tôn trọng. Cao
độ phân giải hình ảnh panchromatic là biểu đồ phù hợp với dải mà nó được thay thế trước
thay thế (tức là, kênh cường độ để làm sắc nét IHS, dải tần có độ phân giải thấp được chọn
cho mỗi kênh màu với độ sắc nét Brovey và hình ảnh PC1 để làm sắc nét PCA).
Theo mặc định, lệnh sẽ cố gắng sử dụng xử lý song song, sử dụng tối đa 3 lõi
đồng thời. Cờ -s sẽ vô hiệu hóa xử lý song song, nhưng sử dụng
biểu thức r.mapcalc để giảm I / O đĩa.
Ba kênh đầu ra được làm sắc nét lia có thể được kết hợp với d.rgb or r.compozit. Màu sắc
có thể được tùy chọn tối ưu hóa với i.colors.enhance. Trong khi hình ảnh màu kết quả sẽ
ở độ phân giải cao hơn trong mọi trường hợp, 3 thuật toán làm sắc nét pan khác nhau về mặt
của phản ứng quang phổ.

VÍ DỤ

Độ sắc nét lia của hình ảnh Landsat từ Boulder, Colorado, Hoa Kỳ:
Tổ hợp # R, G, B ở cự ly 30m
g.region raster=p034r032_7dt20010924_z13_10 -p
d.rgb b=p034r032_7dt20010924_z13_10 g=lp034r032_7dt20010924_z13_20
r=p034r032_7dt20010924_z13_30
# i.pansharpen với thuật toán IHS
i.pansharpen red=p034r032_7dt20010924_z13_30 green=p034r032_7dt20010924_z13_20
blue=p034r032_7dt20010924_z13_10 pan=p034r032_7dp20010924_z13_80
output = ihs321 method = ihs
# hiển thị ở 15m
g.region raster = ihs321_blue -p
d.rgb b = ihs321_blue g = ihs321_green r = ihs321_red
Kết quả:

R, G, B tổng hợp of hạ cánh ở at 30m R, G, B tổng hợp of Brovey mài hình ảnh at 15m

R, G, B tổng hợp of IHS mài hình ảnh at 15m R, G, B tổng hợp of PCA mài hình ảnh at 15m "

Ví dụ: LANDSAT ETM + (Landsat 7), tập dữ liệu mẫu North Carolina:
# ban đầu ở 28 phút
g.region raster = lsat7_2002_10 -p
d.mon wx0.
d.rgb b=lsat7_2002_10 g=lsat7_2002_20 r=lsat7_2002_30
# i.pansharpen với thuật toán IHS
i.pansharpen red = lsat7_2002_30 @ PERMANENT \
green=lsat7_2002_20 blue=lsat7_2002_10 \
pan = lsat7_2002_80 method = ihs \
output = lsat7_2002_ihs
# hiển thị ở 14.25m
g.region raster = lsat7_2002_ihs_red -p
d.xóa
d.rgb r = lsat7_2002_ihs_red g = lsat7_2002_ihs_green b = lsat7_2002_ihs_blue
# so sánh trước / sau (hỗ trợ RGB trong "Nâng cao"):
g.gui.mapswipe
# tùy chọn cân bằng màu:
i.colors.enhance r = lsat7_2002_ihs_red g = lsat7_2002_ihs_green b = lsat7_2002_ihs_blue

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