Landsat8 DN to TOA Reflectance (trừ các giá trị)


8

Tôi đang học cách thực hiện hiệu chỉnh khí quyển bằng tay, do đó không sử dụng các phương pháp tự động có sẵn trong ENVI

siêu dữ liệu (băng 1) cho cảnh Landsat 8 của tôi như sau;

RADIANCE_MULT_BAND_1 = 1.2344E-02

RADIANCE_ADD_BAND_1 = -61.72016

REFLECTANCE_MULT_BAND_1 = 2.0000E-05

REFLECTANCE_ADD_BAND_1 = -0.100000

CN_ELEVATION = 48.12237303

Vì vậy, khi tôi sử dụng toán toán băng tần để tính TOA Radiance, tôi gõ vào đây:

(0,012344 * B1) -61,72016

Và nó mang lại một số liệu thống kê quang phổ trông có vẻ đúng Quang phổ

Khi tôi sử dụng bất kỳ công thức toán học ban nhạc nào sau đây;

Độ phản xạ TOA (không bao gồm điều chỉnh độ cao mặt trời)

(B1 * 0,00002) -0,1

Độ phản xạ TOA (bao gồm cả điều chỉnh độ cao mặt trời)

((B1 * 0,00002) -0,1) / 0,74456226676389733207607359928648

Hoặc tính toán trực tiếp từ DN thô đến TOA Reflectance (bao gồm cả điều chỉnh độ cao mặt trời)

(((

Tôi kết thúc với một cái gì đó như thế này:

Phản xạ TOA

Có vẻ không chính xác.

Tôi đã googled và các thuật ngữ tìm kiếm của tôi là rác rưởi hoặc không ai sử dụng toán học ban nhạc thủ công để làm điều này.

[BIÊN TẬP]

Câu trả lời của Mikkel dưới đây với phép tính toán băng tần được sửa đổi (bỏ qua các bước trung gian) mang lại các thống kê sau:

nhập mô tả hình ảnh ở đây

Mà trông có vẻ thực tế hơn.

Tự hỏi nếu tôi nên thêm 13.4305 vào mỗi băng tần, hoặc bỏ qua mức tối thiểu âm?

[KẾT THÚC]

[EDIT # 2]

Mikkel đã được phát hiện lại, khi đề xuất rằng các giá trị -13.4305 được tách biệt với các khu vực "không có dữ liệu" xung quanh hình ảnh;

nhập mô tả hình ảnh ở đây

như bạn có thể thấy, chỉ bị cô lập ở rìa phía bắc của cảnh.

[KẾT THÚC SỐ 2]

Câu trả lời:


7

Một lỗ hổng trong cách tiếp cận của bạn. Bạn không cần phải thông qua DN để rạng rỡ.
Bạn có thể đi thẳng đến DN để phản xạ.
Chỉ cần dính vào ((B1*0.00002)-0.1)/0.74457226676389733207607359928648.


Điều này làm việc khá tốt! tôi cho rằng hoặc không bao giờ được giải thích rằng phản xạ TOA radiance / TOA excl. độ cao mặt trời không cần thiết các bước trong trình tự điều chỉnh.
kyuss

theo bài đăng đã chỉnh sửa ở trên, có cần thêm 13.4305 vào mỗi băng tần để đưa chúng lên tối thiểu = 0 hay không?
kyuss

@kyuss, hãy thử kiểm tra xem những vị trí -0.134305 đó nằm ở đâu. Tôi nghi ngờ rằng đó là khu vực không có dữ liệu xung quanh rìa ảnh.
Mikkel Lydholm Rasmussen

thật tuyệt vời, tôi thậm chí không kiểm tra bất cứ thứ gì ngoại trừ biên giới phía nam nơi có những giá trị tiêu cực. cảm ơn tất cả sự giúp đỡ của bạn
kyuss

@kyuss bây giờ bước tiếp theo trong quy trình của bạn là "đơn giản" để thực hiện hiệu chỉnh khí quyển.
Mikkel Lydholm Rasmussen

2

Đây là một bài đọc tốt về vật lý của hiệu chỉnh khí quyển và có thể giúp giải quyết vấn đề ... Tôi là một người ấn nút, tôi cảm thấy nỗi đau của bạn, ít nhất bạn có thể chỉ cần nhấp vào nút QUAC hoặc FLAASH tuyệt đẹp đó khi bạn ở bên ngoài của học viện.

http://modis.gsfc.nasa.gov/data/atbd/atbd_mod08.pdf


cảm ơn Steven. tôi đã - như bạn đã suy luận - đã học được nhiều trong quá trình học nhưng tài liệu modis đó có thể giúp lấp đầy khoảng trống hiểu biết của tôi: chuyển giao bức xạ.
kyuss
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.