Bài viết này trình bày kết quả xây dựng quy trình bán tự động chiết tách thông tin nhiệt độ bề mặt từ tư liệu ảnh viễn thám hồng ngoại nhiệt Landsat 5, 7 và 8.
Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất Tập 58, Kỳ (2017) 27-33 27 Xây dựng quy trình bấ n tự độ ng chiế t tấ ch thông tin nhiệ t độ bề mậ t từ ẩ nh viễ n thấ m hồ ng ngôậ i nhiệ t Lândsât sử dụ ng phầ n mề m Erdâs Imâginê Trịnh Lê Hùng 1,*, Vũ Dânh Tuyên 2, Vương Trọng Kha Học viện Kỹ thuật Quân sự, Việt Nam Trường Đại học Tài nguyên Môi trường Hà Nội, Việt Nam Khoa Trắc địa - Bản đồ Quản lý đất đai, Trường Đại học Mỏ - Địa chất, Việt Nam THÔNG TIN BÀI BÁO TĨM TẮT Q trình: Nhận 15/3/2017 Chấp nhận 18/6/2017 Đăng online 31/8/2017 Nhiệt độ bề mặt những yếu tố quan trọng nghiên cứu khí hậu tác động của nó đến môi trường sống của người Nhiệt độ bề mặt có thể được chiết tách từ tư liệu ảnh viễn thám hồng ngoại nhiệt, nhiên việc xác định nhiệt độ bề mặt các phần mềm xử lý ảnh thường nhiều thời gian phải thực nhiều tính toán riêng lẻ Bài báo trình bày kết xây dựng quy trình bán tự động chiết tách thông tin nhiệt độ bề mặt từ tư liệu ảnh viễn thám hồng ngoại nhiệt Landsat 5, Quy trình có ưu điểm đơn giản, giúp tiết kiệm thời gian công sức xây dựng đồ phân bố nhiệt độ bề mặt Từ khóa: Nhiệ t đợ bề mậ t Viễ n thấ m Hồ ng ngôậ i nhiệ t Landsat © 2017 Trường Đại học Mỏ - Địâ chất Tất quyền bảô đảm Mở đầu Nhiệ t độ mậ t đấ t lầ mộ t thông số đầ u vầ ô quân trộ ng cấ c mô hình nghiên cứu khí hậ u, đấ nh giấ hậ n hấ n, phấ t hiệ n đẩ ô nhiệ t đô thị Nhiệ t độ cố thể được chiế t tấ ch từ tư liệ u ẩ nh viễ n thấ m hồ ng ngôậ i nhiệ t Lândsât, Astêr, MODIS (Jin êt âl., 2015) Sô với cấ c phương phấ p truyề n thố ng dựâ số liệ u củ â cấ c trậ m quân trấ c, phương phấ p viễ n thấ m nhiệ t với nhiề u ưu điể m vượt trộ i diệ n tích phủ trù m rộ ng, thời giân cậ p nhậ t ngấ n,… đẫ được sử dụ ng hiệ u quẩ thế giới cũ ng ở Việ t Nâm trông nghiên cứu sự phân bố _ *Tác giả liên hệ E-mail: trinhlehung125@gmail.com nhiệ t độ bề mậ t Dô được cung cấ p hôầ n tôầ n miễ n phí với chu kỳ cậ p nhậ t 16 ngầ y, ẩ nh hồ ng ngôậ i nhiệ t Lândsât 5, vầ lầ mộ t nguồ n tư liệ u chủ đậ ô trông nghiên cứu phân bố nhiệ t độ mậ t đấ t Chô đế n nây, thế giới đẫ cố nhiề u nghiên cứu sử dụ ng tư liệ u ẩ nh viễ n thấ m hồ ng ngôậ i nhiệ t Lândsât vầ trông đấ nh giấ sự phân bố nhiệ t độ bề mậ t vầ mố i quân hệ với lớp phủ (Alipour et al., 2004; Balling and Brazel, 1988; Fei Yuan et al., 2007; Garcia Cueto et al., 2007; Hyung Moo Kim et al., 2005; Mâltick êt âl., 2008; Sôrbinô êt âl., 2004; Kumâr êt âl., 2012) Mộ t số nghiên cứu gầ n (Jin êt âl., 2015; Chên êt âl., 2014) đẫ sử dụ ng tư liệ u ẩ nh hồ ng ngôậ i nhiệ t Lândsât trông đấ nh giấ sự phân bố nhiệ t độ bề mậ t ở mộ t số khu vực tậ i Trung Quố c 28 Trịnh Lê Hùng nnk/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 27-33 Ở Việ t Nâm, đẫ cố mộ t số nghiên cứu ứng dụ ng phương phấ p viễ n thấ m nhiệ t nghiên cứu phân bố nhiệ t độ bề mậ t cấ c đô thị lớn Tư liệ u ẩ nh hồ ng ngôậ i nhiệ t Lândsât vầ Astêr đẫ được sử dụ ng trông đấ nh giấ mố i quân hệ giữâ nhiệ t độ bề mậ t vầ lớp phủ khu vực thầ nh phố Hồ Chí Minh (Trầ n Thị Vân vầ nnk, 2009), khu vực thầ nh phố Hầ Nộ i (Trịnh Lê Hù ng, 2014) Ả nh hồ ng ngôậ i nhiệ t Lândsât cũ ng được sử dụ ng trông cấ c nghiên cứu đấ nh giấ độ ẩ m đấ t, phân vù ng nguy hậ n hấ n (Trịnh Lê Hù ng, 2014), phấ t hiệ n vầ giấ m sấ t chấ y mổ (Trịnh Lê Hù ng, 2016) Tư liệ u ẩ nh hồ ng ngôậ i nhiệ t Lândsât7 ETM+ cũ ng được sử dụ ng nhầ m xấ c định nhiệ t độ bề mậ t phụ c vụ xây dựng mô hình dự bấ ô nguy chấ y rừng ở tỉnh Đấ k Lấ k (Lưu Thế Anh vầ nnk, 2014) Cấ c nghiên cứu cũ ng chô thấ y, sự chênh lệ ch giữâ nhiệ t độ chiế t tấ ch từ tư liệ u ẩ nh hồ ng ngôậ i nhiệ t Lândsât với số đô thực tế đậ t khôẩ ng từ 0,5 đế n 2,00C (Trầ n Thị Vân vầ nnk, 2009; Zhâng êt âl., 2006; Kumâr êt al., 2012) Xấ c định nhiệ t độ bề mậ t từ tư liệ u ẩ nh hồ ng ngôậ i nhiệ t Lândsât lầ mộ t công việ c tương đố i phức tậ p dô quấ trình xử lý phẩ i trẩ i quâ nhiề u bước riêng lể Đầ u tiên, kênh hồ ng ngôậ i nhiệ t ẩ nh Lândsât (kênh ẩ nh Lândsât vầ 7, kênh 10 hôậ c 11 ẩ nh Lândsât 8) được sử dụ ng để xấ c định nhiệ t độ bức xậ (brightnêss têmpêrâturê) thêô mô hình dô nhầ sẩ n xuấ t ẩ nh cung cấ p (Lândsât dâtâ usêrs hândbôôk, 2016; Lândsât dâtâ usêrs hândbôôk vêrsiôn 2.0, 2016), sâu đố nhiệ t độ bề mậ t được xấ c định sở hiệ u chỉnh mố i quân hệ với lớp phủ thông quâ độ phấ t xậ (êmissivity) (Trầ n Thị Vân vầ nnk, 2009; Yuân ând Bâuêr, 2007; Jin êt âl., 2015) Để tính tôấ n nhiệ t độ bề mậ t, mộ t số nghiên cứu đẫ xây dựng cấ c môdulê tính nhiệ t độ bầ ng cấ c ngôn ngữ lậ p trình khấ c nhâu (Trịnh Lê Hù ng, 2016; Zhâng êt âl., 2006) Ngôầ i râ, quấ trình tính tôấ n nhiệ t độ bề mậ t cũ ng cố thể được thực hiệ n bầ ng cấ ch sử dụ ng cấ c phầ n mề m xử lý ẩ nh với sự hỗ trợ củ â cấ c công cụ Bând Mâth (ENVI), Môdêlêr (ERDAS Imâginê) Mậ c dù vậ y, nế u sử dụ ng cấ c thâô tấ c đơn lể chô từng bước trông quy trình tính nhiệ t độ sễ tố n kế m nhiề u thời giân vầ công sức trông quấ trình xử lý ẩ nh Bầ i bấ ô nầ y trình bầ y kế t quẩ xây dựng quy trình bấ n tự độ ng chiế t tấ ch thông tin nhiệ t độ bề mậ t từ tư liệ u ẩ nh Lândsât 5, vầ với sự hỗ trợ củ â công cụ Môdêlêr trông phầ n mề m ERDAS Imâginê, giú p tiế t kiệ m thời giân trông quấ trình tính tôấ n Phương pháp nghiên cứu Nhiệ t độ bề mậ t (lând surfâcê têmpêrâturê) cố thể được xấ c định thêô mô hình dô NASA cung cấ p thêô cấ c bước sâu: Bước 1: Chuyển giá trị số nguyên của ảnh sang giá trị bức xạ phổ Để tính nhiệ t độ bề mậ t, đầ u tiên, giấ trị số nguyên củ â kênh hồ ng ngôậ i nhiệ t được chuyể n sâng giấ trị thực củ â bức xậ điệ n từ (𝑊𝑚−2 µ𝑚−1 ) Với ẩ nh LANDSAT vầ 7, giấ trị bức xậ phổ được xấ c định sâu (Lândsât sciêncê dâtâ usêrs handbook, 2016): L Lmax Lmin ( DN DN ) Lmin (1) DN max DN Trông đố : L - giấ trị bức xậ phổ ; Lmax, Lmin - giấ trị bức xậ phổ ứng với DNmax vầ DNmin (giấ trị nầ y được lấ y từ filê mêtâdâtâ trông dữ liệ u ẩ nh Lândsât) (Lândsât sciêncê dâtâ usêrs handbook, 2016); DNmax - giấ trị số lớn nhấ t (=255); DNmin - giấ trị số nhổ nhấ t (=1) Bảng Giá trị Lmax, Lmin đối với ảnh hồng ngoại nhiệt LANDSAT TM, ETM+ Kênh Tư liệu Lmax Lmin LANDSAT /ETM+ High 6.1 12,65 3,2 gain LANDSAT /ETM +Low 6.2 17,04 0,0 gain LANDSAT TM 15,503 1,238 Đố i với ẩ nh Lândsât 8, giấ trị bức xậ phổ được xấ c định thêô công thức sâu (Landsat science data users handbook version 2.0, 2016): L M L Qcal AL (2) Trông đố : ML, AL - hệ số chuyể n đổ i, được lấ y trông filê mêtâdâtâ ẩ nh Lândsât (giấ trị RADIANCE_MIUL_BAND_x vầ RADIANCE_ADD_BAND_x) (Bẩ ng 2); Qcal - Giấ trị số củ â ẩ nh (digitâl numbêr) Trịnh Lê Hùng nnk/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 27-33 Bảng Hệ số 𝑀𝐿 , 𝐴𝐿 đối với ảnh Landsat STT Tư liệu Kênh ML AL Landsat TIRS 10 3,3420.10-4 0,10000 Landsat TIRS 11 3,3420.10-4 0,10000 Bước 2: Xác định nhiệt độ bức xạ Ở bước tiế p thêô, giấ trị bức xậ phổ được sử dụ ng để tính nhiệ t độ bức xậ , hây cồ n gộ i lầ nhiệ t độ độ chố i (brightnêss têmpêrâturê) thêô công thức sâu (Landsat science data users handbook, 2016): K2 TB (3) K1 ln( 1) L Trông đố : 𝑇𝐵 - Nhiệ t độ bức xậ (K); 𝐾1 - Hầ ng sớ chủ n đở i (W/(m2.sr.µm)); 𝐾2 - Hầ ng số chuyể n đổ i (0K) Giấ trị 𝐾1 , 𝐾2 được cung cấ p trông filê mêtâdâtâ ẩ nh Lândsât (Bẩ ng 3) (Landsat science data users handbook, 2016; Landsat science data users handbook version 2.0, 2016) Bảng Hệ số K1, K2 đối với tư liệu ảnh hồng ngoại nhiệt Landsat TM, ETM+, Landsat STT Tư liệu Kênh K1 (W/(m2.sr.µm)) K2 (0K) Landsat 607,76 1260,56 TM Landsat 666,09 1282,71 ETM+ 10 774,89 1321,08 Landsat 11 480,89 1201,14 Bước 3: Xác định độ phát xạ bề mặt Nhiệ t độ bề mậ t cố quân hệ mậ t thiế t với lớp phủ , dô vậ y để tính nhiệ t độ bề mậ t cầ n phẩ i hiệ u chỉnh ẩ nh hưởng củ â lớp phủ thông quâ độ phấ t xậ (êmissivity) Độ phấ t xậ cố thể xấ c định bầ ng cấ ch sử dụ ng nhiệ t kế bức xậ chuẩ n thông quâ kế t quẩ đô cấ c mẫ u đậ i diệ n Tuy nhiên phương phấ p nầ y chỉ phù hợp với điề u kiệ n phồ ng thí nghiệ m vầ rấ t tố n kế m Để xấ c định độ phấ t xậ bề mậ t từ tư liệ u viễ n thấ m cố thể sử dụ ng phương phấ p dựâ kế t quẩ phân lôậ i lớp phủ (Zâng êt âl., 2006) hôậ c phương phấ p sử dụ ng chỉ số thực vậ t NDVI (Vâlôr ând Câsêllês, 1996; Vân dê Griênd ând Owên, 29 1993) Trông nghiên cứu nầ y, độ phấ t xậ bề mậ t được xấ c định thêô phương phấ p dô Vâlôr vầ Câsêllês (1996) đề xuấ t dựâ chỉ số thực vậ t NDVI (Nôrmâlizêd Difêrêncê Vêgêtâtiôn Indêx) Phương phấ p nầ y cố nhiề u ưu điể m sô với phương phấ p dựâ kế t quẩ phân lôậ i lớp phủ dô cố thể xấ c định độ phấ t xậ chi tiế t tới từng pixêl cũ ng cố thể ấ p dụ ng hiệ u quẩ chô cấ c khu vực với nhiề u kiể u bề mậ t khấ c nhâu (Trầ n Thị Vân vầ nnk, 2009) Thêô phương phấ p nầ y, độ phấ t xậ bề mậ t được xấ c định sâu (Vâlôr ând Câsêllês, 1996): v Pv s (1 Pv ) (4) Trông đố : ε - Độ phấ t xậ bề mậ t 𝜀𝑣 - Độ phấ t xậ đố i với thực vậ t 𝜀𝑠 - Độ phấ t xậ đố i với đấ t trố ng 𝑃𝑣 - Tỉ lệ thực vậ t trông mộ t pixêl ẩ nh, được xấ c định thêô công thức sâu (Jin êt âl., 2015): Pv ( NDVI NDVI ) NDVI max NDVI (5) Như vậ y, để xấ c định độ phấ t xậ cầ n phẩ i biế t trước độ phấ t xậ chô đấ t vầ thực vậ t Đầ u tiên, cầ n lấ y mộ t số dữ liệ u mẫ u đố i với đấ t vầ thực vậ t ẩ nh chỉ số NDVI, sâu đố độ phấ t xậ củ â đấ t vầ thực vậ t được xấ c định thêô công thức dô Vân dê Griênd vầ Owên đề xuấ t (Vân dê Griênd ând Owên, 1993): (6) 1,0094 0,047ln( NDVI ) Bước 4: Xác định nhiệt độ bề mặt Ở bước cuố i cù ng, nhiệ t độ bề mậ t (LST - lând surfâcê têmpêrâturê) được xấ c định thêô công thức sâu (Alipôur êt âl., 2004; Chên êtâl., 2014; Jin et al., 2015): LST TB T B ln (7) Trông đố : 𝜆 - giấ trị bước số ng trung tâm kênh hồ ng ngôậ i nhiệ t; h.c , σ - hầ ng số Stêfân Bôltzmânn (5,67.10-8 (Wm-2.K-4)); h - Hầ ng số Plânk (6,626.10-34J.sec); c - vậ n tố c ấ nh sấ ng (2,998.108 m/s) 30 Trịnh Lê Hùng nnk/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 27-33 Kết quả và thảo luận Trông nghiên cứu nầ y tiế n hầ nh xây dựng cấ c quy trình bấ n tự độ ng chiế t tấ ch thông tin nhiệ t độ bề mậ t từ tư liệ u ẩ nh hồ ng ngôậ i nhiệ t Lândsât cấ c thế hệ khấ c nhâu, bâô gồ m ẩ nh nhiệ t Lândsât 5, vầ Đầ u vầ ô bâô gồ m ẩ nh Lândsât vầ độ phấ t xậ bề mậ t, đầ u râ lầ nhiệ t độ bề mậ t Sở dĩ đầ u vầ ô củ â quy trình bâô gồ m kế t quẩ xấ c định độ phấ t xậ bề mậ t lầ dô cấ c khu vực khấ c nhâu cố đậ c điể m vậ t lý củ â bề mậ t khấ c nhâu, dô vậ y để nâng câô độ chính xấ c cầ n phẩ i xấ c định độ phấ t xậ chô từng khu vực cụ thể (Trầ n Thị Vân vầ nnk, 2009) Việ c nầ y cố thể thực hiệ n rấ t đơn giẩ n thêô công thức (4) Tư liệ u ẩ nh vệ tinh Lândsât đầ u vầ ô đẫ được tiề n xử lý nhầ m lôậ i bổ cấ c sâi số về mậ t hình hộ c vầ sâi số về phổ trước đưâ vầ ô quy trình Ví dụ mô hình tự độ ng chiế t tấ ch thông tin nhiệ t độ bề mậ t từ tư liệ u ẩ nh vệ tinh Lândsât TM được trình bầ y Hình Tương tự đố i với ẩ nh Lândsât 5, trông nghiên cứu cũ ng xây dựng cấ c quy trình xấ c định nhiệ t độ bề mậ t từ ẩ nh hồ ng ngôậ i nhiệ t Lândsât vầ sở thây đổ i mộ t số hệ số trông quấ trình tính tôấ n để phù hợp với cấ c tư liệ u ẩ nh Để đấ nh giấ kế t quẩ xấ c định nhiệ t độ bề mậ t bầ ng cấ c mô hình đẫ được xây dựng, trông bầ i bấ ô nầ y tiế n hầ nh thử nghiệ m với lôậ i ẩ nh hồ ng ngôậ i nhiệ t khấ c nhâu: Lândsât TM, Lândsât ETM+ vầ Lândsât Hình Quy trình bán tự động chiết tách thông tin nhiệt độ bề mặt từ tư liệu ảnh Landsat (a) (b) Hình Ảnh Landsat TM khu vực Lâm Hà (Lâm Đồng) ngày 07 - 02 - 2011 kết xác định nhiệt độ bề mặt Trịnh Lê Hùng nnk/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 27-33 Thử nghiệm 1: ảnh Landsat TM Ả nh Lândsât TM sử dụ ng để tính nhiệ t độ bề mậ t được chụ p ngầ y 07 - 02 - 2011 khu vực huyệ n Lâm Hầ , tỉnh Lâm Đồ ng (Hình 2â) Kế t quẩ xấ c định nhiệ t độ bề mậ t được trình bầ y hình 2b, trông đố nhiệ t độ thấ p nhấ t đậ t 25,420C, nhiệ t độ câô nhấ t đậ t 41,370C Thử nghiệm 2: Ảnh Landsat ETM+ Trông thử nghiệ m nầ y sử dụ ng tư liệ u ẩ nh Landsat ETM+ chụ p ngầ y 08 - 11 - 2007 khu vực thầ nh phố Thấ i Nguyên, tỉnh Thấ i Nguyên (Hình 3â) Kế t quẩ xấ c đinh nhiệ t độ bề mậ t được thể hiệ n hình 3b, trông đố nhiệ t độ thấ p nhấ t vầ câô nhấ t đậ t tương ứng lầ 22,040C vầ 40,710C Thử nghiệm 3: Ảnh Landsat Trông thử nghiệ m nầ y, tư liệ u viễ n thấ m được sử dụ ng để tính nhiệ t độ bề mậ t lầ ẩ nh Lândsât khu vực huyệ n Thậ ch Hầ , tỉnh Hầ Tĩnh chụ p ngầ y 18 - 06 - 2013 (Hình 4â) Kế t quẩ xấ c định nhiệ t độ bề mậ t được trình bầ y hình 4b, trông đố nhiệ t độ thấ p nhấ t đậ t 21,620C, nhiệ t độ câô nhấ t đậ t 38,670C Phân tích kế t quẩ nhậ n được trông bâ thử nghiệ m cố thể nhậ n thấ y, nhiệ t độ bề mậ t cố quân hệ mậ t thiế t với lớp phủ Nhiệ t độ thường đậ t câô tậ i cấ c khu vực không cố hôậ c cố ít thực vậ t chê phủ (khu vực đô thị, đấ t trố ng, đấ t nông nghiệ p), trông tậ i cấ c khu vực cố thực vậ t phủ dầ y hôậ c mậ t nước, nhiệ t độ thấ p rấ t nhiề u Chênh lệ ch giữâ nhiệ t độ câô nhấ t vầ thấ p nhấ t trông cấ c thử nghiệ m đề u đậ t 160C Để sô sấ nh, đấ nh giấ tố c độ tính tôấ n, trông nghiên cứu cũ ng thực hiệ n tính nhiệ t độ đố i với cấ c khu vực nghiên cứu sử dụ ng chương trình LST dô nhố m tấ c giẩ xây dựng sở ngôn ngữ lậ p trình C++ (Trịnh Lê Hù ng, 2014) Phân tích kế t quẩ nhậ n được chô thấ y, nhiệ t độ tính bầ ng quy trình xây dựng phầ n mề m ERDAS Imâginê vầ chương trình LST lầ nhâu, mậ c dù vậ y tố c độ tính tôấ n sử dụ ng chương trình LST lâu (a) (b) Hình Ảnh Landsat ETM+ ngày 08 - 11 - 2007 khu vực thành phố Thái Nguyên kết xác định nhiệt độ bề mặt (a) 31 (b) Hình Ảnh Landsat khu vực Thạch Hà, Hà Tĩnh ngày 18 - 06 - 2013 kết xác định nhiệt độ bề mặt 32 Trịnh Lê Hùng nnk/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 27-33 nhiề u (khôẩ ng giây sô với giây) Điề u nầ y cố thể giẩ i thích dô sự khố khân trông xử lý dữ liệ u ẩ nh vệ tinh với kích thước lớn, trông đố sử dụ ng quy trình xây dựng từ phầ n mề m ERDAS Imâginê chô phế p tậ n dụ ng khẩ nâng xử lý ẩ nh từ phầ n mề m nầ y Kết luận Ả nh viễ n thấ m hồ ng ngôậ i nhiệ t Lândsât lầ mộ t nguồ n tư liệ u phông phú vầ quý giấ trông chiế t tấ ch thông tin nhiệ t độ bề mậ t Với độ phân giẩ i không giân trung bình (60 - 120 m), thời giân cậ p nhậ t ngấ n (16 ngầ y), đậ c biệ t được cung cấ p hôầ n tôầ n miễ n phí (tậ i wêbsitê www.glovis.usgs.gov), ẩ nh hồ ng ngôậ i nhiệ t Lândsât đẫ vầ đâng được sử dụ ng rộ ng rẫ i vầ mâng lậ i hiệ u quẩ quân trộ ng trông nghiên cứu diễ n biế n nhiệ t độ bề mậ t Phân tích kế t quẩ nhậ n được thử nghiệ m tính tôấ n nhiệ t độ bề mậ t từ tư liệ u ẩ nh hồ ng ngôậ i nhiệ t Lândsât 5, vầ chô thấ y, quy trình xây dựng trông nghiên cứu nầ y cố thể tính tôấ n nhiệ t độ mộ t cấ ch đơn giẩ n, tiế t kiệ m thời giân với độ chính xấ c đẩ m bẩ ô, đồ ng thời không đồ i hổ i người sử dụ ng phẩ i cố kỹ nâng lậ p trình Kế t quẩ nhậ n được trông nghiên cứu cố thể sử dụ ng phụ c vụ công tấ c giấ m sấ t vầ đấ nh giấ phân bố nhiệ t độ bề mậ t cũ ng cố thể lầ m đầ u vầ ô chô cấ c mô hình dự bấ ô chấ y rừng, phấ t hiệ n vầ giấ m sấ t chấ y ngầ m cũ ng đấ nh giấ hậ n hấ n Tài liệu tham khảo Lưu Thế Anh, Trầ n Anh Tuấ n, Hôầ ng Thị Huyề n Ngộ c, Lê Bấ Biên, 2014 Ứng dụ ng tư liệ u ẩ nh viễ n thấ m vầ công nghệ GIS thầ nh lậ p bẩ n đồ nguy chấ y rừng tỉnh Đấ k Lấ k, Tạp chí Các khoa học về Trái đất 36(3), 252 - 261 Trầ n Thị Vân, Hôầ ng Thấ i Lân, Lê Vân Trung, 2009 Phương phấ p viễ n thấ m nhiệ t trông nghiên cứu phân bố nhiệ t độ bề mậ t đô thị, Tạp chí Các khoa học về Trái đất 31(2), 168 - 177 Trịnh Lê Hù ng, 2016 Xấ c định nhiệ t độ bề mậ t từ tư liệ u viễ n thấ m nhiệ t phụ c vụ phấ t hiệ n vầ giấ m sấ t chấ y ngầ m, thực nghiệ m khu vực mổ thân Nông Sơn, tỉnh Quẩ ng Nâm, Tạp chí Công nghiệp Mỏ 2, 73 - 76 Trịnh Lê Hù ng, 2014 Nghiên cứu sự phân bố nhiệ t độ bề mậ t bầ ng dữ liệ u ẩ nh vệ tinh đâ phổ LANDSAT, Tạp chí Các khoa học về Trái đất 36(1), 82 - 89 Trịnh Lê Hù ng, 2014 Ứng dụ ng dữ liệ u viễ n thấ m hồ ng ngôậ i nhiệ t LANDSAT nghiên cứu độ ẩ m đấ t sở chỉ số khô hậ n nhiệ t độ thực vậ t, Tạp chí Các khoa học về Trái Đất, 36(3), 262 270 Alipour T., Sarajian M R., Esmaseily A., 2004 Land surface temperature estimation from thermal band of LANDSAT sensor, case study: Alashtar city, The international archives of the Photogrammetry, Remote sensing and spatial information sciences, Vol XXXVIII-4/C7 Amit Kumar Verma, Sanjeev Kumar, 2015 Mapping fire hazard in Rajaji national park, future perspective of wildlife habitat conservation by using remote sensing and GIS, International Journal of Technical Research and Application 3(5) 5, 81 - 88 Balling R C., Brazel S W., 1988 High - resolution surface temperature patterns in a complex urban Terrain, Photogrammetric engineering and Remote sensing 54(9), 1289 - 1293 Chen Du, Huazhong Ren, Qiming Qin, Jinjie Meng, Jing Li, 2014 Split-Window algorithm for estimating land surface temperature from Landsat TIRS data, Geoscience and Remote Sensing, DOI: 10.1109/IGARSS.2014.6947256 Fei Yuan, Marvin E Bauer, 2007 Comparison of impervious surface area and normalized difference vegetation index as indicators of surface urban heat island effects in LANDSAT imagery, Remote sensing of Environment 106,375 - 386 Garcia Cueto O R., Jauregui Ostos E., Toudert D., Tejeda Martinez A., 2007 Detection of the urban heat island in Mexicali and its relationship with land use, Atmosfera 20(2), 111 - 131 Hongyuan Huo, Xiaoguang Jiang, Xianfeng Song, Zhao-Liang Li, Zhupya Ni, Caixia Gao, 2014 Detection of coal fire dynamics and propagation direction from multi-temporal nighttime Landsat SWIR and TIR data A case study on the Rujigou coalfield, northwest China, Remote Sensing 6, 1234 - 1259 Trịnh Lê Hùng nnk/Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 58 (4), 27-33 Hyung Moo Kim, Beob Kyun Kim, Kang Soo You, 2005 A statistic correlation analysis algorithm between land surface temperature and vegetation index, International journal of information processing systems 1(1), 102 - 106 Jinqu Zhang, Yunpeng Wang, Yan Li, 2006 A C++ program for retrieving land surface temperature from the data of Landsat TM/ETM+ band 6, Computers and Geosciences, 32, 1796 - 1805 Maltick Javed, Yogesh Kant, D B Bharath, 2008 Estimation of land surface temperature over Delhi using LANDSAT-7 ETM+, Journal Ind Geophysics Union 12(3), 131 - 140 Meijun Jin, Junming Li, Caili Wang, Ruilan Shang, 2015 A practical split-window algorithm for retrieving land surface temperature from Landsat data and a case study of an urban area in China, Remote Sensing 7, 4371 - 4390 NASA, Landsat science data users handbook, 2016 https://landsat.gsfc.nasa.gov/wpcontent/uploads/2016/08/Landsat7_Handbo ok.pdf 33 NASA, Landsat (L8) data users handbook version 2.0, 2016 https://landsat.usgs.gov/sites/default/files/ /Landsat8DataUsersHandbook.pdf Sorbino J A., Jimenez Munoz J C., Leonardo P., 2004 Land surface temperature retrieval from Landsat TM5, Remote Sensing of Environment 90(4), 434 - 440 Sundara Kumar K., Udaya Bhaskar P., Padmakumari K., 2012 Estimation of land surface temperature to study urban heat island effect using LANDSAT ETM+ image, International journal of Engineering Science and technology 4(2), 771 - 778 Valor E., Caselles V., 1996 Mapping land surface emissivity from NDVI Application to European African and South American areas, Remote sensing of Environment 57, 167 - 184 Van de Griend A A., Owen M., 1993 On the relationship between thermal emissivity and the normalized difference vegetation index for natural surface, International journal of remote sensing 14, 1119 - 1131 ABSTRACT A procedure for semi-automated calculating of land surface temperature using landsat thermal infrared image based on erdas imagine program Hung Le Trinh 1,*, Tuyen Danh Vu 2, Kha Trong Vuong Military Technical Academy, Vietnam Hanoi University of Natural Resources and Environment, Vietnam Faculty of Geomatics and Land Administration, Hanoi University of Mining and Geology, Vietnam Land surface temperature is one of the most important factors in climatology studies and human environment interactions Land surface temperature can be extracted from thermal infrared remote sensing data, however temperature determination using image processing software is difficult, due to individual calculations This paper present of development procedure for semi-automated calculating of land surface temperature using Landsat 5, and thermal infrared data The procedure has the advantage of simplicity, saving time, effort in mapping land surface temperature distribution Keywords: Procedure, semi-automated, Land surface temperature, Remote sensing, Thermal infrared, Landsat ... tự động chiết tách thông tin nhiệt độ bề mặt từ tư liệu ảnh Landsat (a) (b) Hình Ảnh Landsat TM khu vực Lâm Hà (Lâm Đồng) ngày 07 - 02 - 2011 kết xác định nhiệt độ bề mặt Trịnh Lê... (Bẩ ng 3) (Landsat science data users handbook, 2016; Landsat science data users handbook version 2.0, 2016) Bảng Hệ số K1, K2 đối với tư liệu ảnh hồng ngoại nhiệt Landsat TM, ETM+, Landsat STT... Giá trị Lmax, Lmin đối với ảnh hồng ngoại nhiệt LANDSAT TM, ETM+ Kênh Tư liệu Lmax Lmin LANDSAT /ETM+ High 6.1 12,65 3,2 gain LANDSAT /ETM +Low 6.2 17,04 0,0 gain LANDSAT TM 15,503 1,238 Đố