Đang tải... (xem toàn văn)
Mục tiêu và nhiệm vụ của luận vănMục tiêu:Ứng dụng tư liệu Radar và tư liệu Quang học để đánh giá nhanh mức độ ảnh hưởng của lũ lụt lên lớp phủ mặt đất.Nhiệm vụ: Tổng quan về phương pháp tiền xử lý dữ liệu Radarsat Xử lý dữ liệu Radarsat để chiết tách thông tin vùng ngập lũ Xử lý dữ liệu Landsat để chiết tách thông tin về lớp phủ bề mặt Tích hợp GIS để phân tích đánh giá mức độ ảnh hưởng lên từng đối tượng cụ thể thuộc lớp phủ bề mặtTính mới của luận vănỨng dụng công nghệ Viễn thám và GIS để đưa ra báo cáo nhanh về hiện trạng và mức độ ảnh hưởng của lũ lên các đối tượng chính thuộc lớp phủ bề mặt.
mục lục Danh mục từ viết tắt danh mục hình vẽ bảng biểu .4 phần mở đầu Chơng sở phơng pháp luận ứng dụng ảnh viễn thám nghiên cứu ngập lụt .11 1.1 Đặc điểm ảnh Quang học .11 1.2 Các đặc tính ảnh Radar .11 1.2.1 Cơ chế tán xạ Radar 12 1.2.2 Phân cực 13 1.2.3 Đặc tính phản xạ nớc ảnh Radar 14 1.3 Khả thông tin ảnh Quang học lớp phủ mặt đất 15 1.4 Khả thông tin ảnh Radar chiết tách vùng ngập lụt 17 1.5 Lợi ích kết hợp ảnh Quang học ảnh Radar 17 Chơng xử lý ảnh rADAR 19 2.1 ảnh vệ tinh RADARSAT-1 19 2.2 Định chuẩn ảnh RADARSAT 20 2.3 Xử lý nhiễu ảnh Radar .24 2.3.1 Mô hình nhiễu 25 2.3.2 Nhóm phơng pháp lọc cục (Adaptive Filter) 26 2.3.3 Lọc đa thời gian (Multitemporal Filter) 29 2.4 Nắn chỉnh hình học ảnh Radar 32 2.4.1 Những biến dạng hình học ảnh Radar 32 2.4.2 Phơng pháp nắn chỉnh hình học .35 2.5 Chiết tách vùng ngập từ ảnh Radar 38 Chơng nghiên cứu quy trình tích hợp ảnh radar ảnh quang học để đánh giá ảnh hởng tình trạng ngập lụt đến lớp phủ bề mặt 40 3.1 Nghiên cứu Quy trình tích hợp thông tin từ ảnh Radar ảnh Quang học .40 3.1.1 Sơ đồ quy trình 40 3.1.2 Mô tả quy trình 42 3.2 Nghiên cứu phơng pháp đánh giá ảnh hởng ngập lụt tới đối tợng lớp phủ bề mặt 45 3.2.1 Hệ thống phân loại 45 3.2.2 Các tiêu đánh giá mức độ ngập 46 3.2.3 Lập đồ mức độ ảnh hởng 46 3.3 Thực nghiệm đánh giá ảnh hởng tình trạng ngập lụt đến lớp phủ bề mặt 46 3.3.1 T liệu sử dụng 46 3.3.2 Chiết tách thông tin lớp phủ mặt đất từ ảnh Quang học 47 3.3.3 Xử lý ảnh RADAR chiết tách thông tin vùng ngập lụt 48 3.3.4 Chồng ghép thông tin Quang học RADAR 51 3.4 Phân tích đánh giá kết thực nghiệm 55 Kết luận .61 Kiến nghị 62 Tài liệu tham khảo 63 Phụ lục 65 Danh mục từ viết tắt CDPF CESBIO CNES dB DEM ENVI EROS ERS IRS Cơ quan Xử lý liệu Canada Trung tâm nghiên cứu Sinh trực thuộc Trung tâm Nghiên cứu Không gian Pháp Trung tâm nghiên cứu Không gian Pháp Deci-Ben Mô hình số độ cao Phần mềm thơng mại xử lý ảnh viễn thám Vệ tinh quan sát trái đất Israel Vệ tinh viễn thám Châu Âu Vệ tinh viễn thám ấn Độ LANDSAT V tinh ti nguyên ca M MIKE-11 Mụ hỡnh tớnh toỏn thy vn, thy lc ca an Mch NOAA RADAR RADARSAT SAR Vệ tinh nghiên cứu đại dơng khí tợng Mỹ Thiết bị dò tìm xác định khoảng cách sóng Raio Vệ tinh quan sát trái đất Canada Radar độ mở tổng hợp danh mục hình vẽ bảng biểu A Hình vẽ Hình 1-1 Dải tần số hoạt động Radar 13 Hình 1-2 Các bề mặt tán xạ Radar 13 Hình 1-3 Các kiểu tán xạ bề mặt khác 14 Hình 1-4 Các kiểu tán xạ môi trờng điện môi khác 14 Hình 1-5 Các kiểu phân cực viễn thám Radar 14 Hình 1-6 Cơ chế tán xạ bề mặt đối tợng ảnh Radar 15 Hình 1-7 Sự khác biệt tán xạ nớc phân cực HH VV 17 Hình 2-1 Các chế độ chụp ảnh vệ tinh RADARSAT-1 21 Hình 2-2 Sự tạo thành nhiễu ảnh Radar 28 Hình 2-3 Mối tơng quan độ lệch chuẩn giá trị trung bình cục 29 Hình 2-4 Hiện tợng co ngắn ảnh Radar 36 Hình 2-5 Hiện tợng biến dạng co ngắn ảnh Radar 37 Hình 2-6 Các tợng biến dạng ảnh Radar địa hình gây 37 Hình 2-7 Hiện tợng co ngắn phía trớc 38 Hình 2-8 Hiện tợng chồng đè ảnh Radar 39 Hình 2-9 Hiện tợng bóng ảnh Radar 39 Hình 2-10 ảnh hởng chênh cao địa hình tới vị trí điểm thực địa 43 Hình 3-1 ảnh Radarsat trớc lọc 55 Hình 3-2 ảnh Radarsat sau sử dụng phép lọc đa thời gian 56 Hình 3-3 ảnh Radarsat sau đợc lọc phơng pháp lọc Lee 56 Hình 3-4 ảnh Radarsat sau đợc lọc phơng pháp lọc Frost 57 Hình 3-5 Hiện trạng ngập thời điểm lũ ngày 06/11/1999 58 Hình 3-6 Hiện trạng lũ rút ngày 15/11/1999 58 Hình 3-7 Kết phân loại lớp phủ bề mặt từ ảnh Landsat ngày 02/06/2001 59 Hình 3-8 Mức độ ngập thời điểm lũ ngày 06/11/1999 60 Hình 3-9 Mức độ rút lũ thời điểm 15/11/1999 60 Hình 3-10 Bản đồ mức độ thiệt hại ngập thời điểm lũ 06/11/1999 66 B Bảng biểu Bảng 2-1 Các sản phẩm ảnh Radarsat-1 22 Bảng 2-2 Nguyên lý phân bố Radiometric Data Record 23 Bảng 3-1 Mức độ ảnh hởng đối tợng lũ gây 65 Biểu đồ 3-1 Diện tích đối tợng thời điểm lũ 61 Biểu đồ 3-2 Số liệu diện tích mặt nớc 62 Biểu đồ 3-3 Diện tích đối tợng bị ngập ngày 06/11/1999 63 Biểu đồ 3-4 Biểu đồ thể mức độ ảnh hởng theo diện tích đối tợng xảy lũ 63 phần mở đầu Tính cấp thiết luận văn Lũ, Lụt thiên tai có nguồn gốc khí hậu Thủy, hoả, đạo, tặc bốn loại tai ơng mà từ ngàn đời xa đợc ông cha ta coi tai họa nguy hiểm đe dọa ngời, thủy tai đợc xếp hạng hàng đầu Trong thập kỷ gần đây, lũ lụt ngày nhiều gây thiệt hại nặng nề không lãnh thổ Việt Nam nhiều nớc giới nh Tây Âu (1998, 2000), Bangladesh (2001), Trung Quốc (1998), Lịch sử nớc ta không ghi lại đầy đủ số liệu trận lụt lớn kỷ trớc, nhng biết nhiều trận lụt dẫn đến nạn đói chết hai triệu ngời năm 1945 trận lũ lịch sử năm 1971 gây vỡ đê làm ngập lụt nghiêm trọng nhiều tỉnh đồng Bắc Bộ Tại miền Trung, lũ lụt nghiêm trọng diện rộng xảy năm 1964 tỉnh từ Quảng Bình vào Phú Yên; lũ lịch sử làm vỡ đê gây ngập lụt lớn sông Cả, La năm 1978, năm 1998, 1999, lũ lụt lịch sử liên tiếp sông thuộc tỉnh từ Quảng Bình đến Khánh Hòa gây thiệt hại nghiêm trọng ngời tài sản Và liên tiếp tháng 10, 11 năm 2007 vừa qua nhân dân miền Trung lại phải liên tiếp hứng chịu trận lụt với quy mô lớn hẳn trận lụt lịch sử năm 1999 gây thiệt hại nặng nề ngời Ước tính khoảng 200 ngời chết, hàng trăm ngời bị tích, thiệt hại tài sản khoảng 9509 tỉ đồng (theo http://www.dartmouth.edu/floods/archives/2007sum) Đã có nhiều dự án, đề tài nghiên cứu khoa học nớc lũ lụt biện pháp phòng tránh Có thể kể số đề tài dự án đợc thực nớc nh: Đề tài 1: " Điều tra nghiên cứu cảnh báo ngập lụt phục vụ phòng tránh thiên tai lu vực sông Miền Trung" Viện Khí tợng Thủy văn thực năm 1999 -2002 Phơng pháp thực hiện: Sử dụng số liệu thực đo điều tra bổ sung sau sử dụng mô hình DEM để xây dựng đồ ngập Đề tài 2: " Nghiên cứu sở khoa học cho giải pháp tổng thể dự báo phòng tránh ngập lụt tỉnh Miền Trung"" Viện Địa lý, Trung tâm Khoa học tự nhiên Công nghệ Quốc gia thực năm 2000 -2004 Phơng pháp thực hiện: Sử dụng số liệu thực đo điều tra bổ sung Đề tài 3: " Nghiên cứu xây dựng tập đồ ngập lụt tỉnh Thừa Thiên Huế " Viện Địa lý, Trung tâm Khoa học tự nhiên Công nghệ Quốc gia thực năm 1999-2001 Phơng pháp thực hiện: Sử dụng số liệu thực đo điều tra bổ sung Đề tài 4: " Xây dựng đồ phân vùng ngập lụt phơng án cảnh báo nguy ngập lụt hạ lu sông Hơng, sông Bồ tỉnh Thừa Thiên Huế" Đài Khí tợng Thủy văn Khu vực Trung Trung Bộ thực năm 1999-2001 Phơng pháp thực hiện: Sử dụng số liệu thực đo điều tra bổ sung có kết hợp với mô hình số độ cao để lập đồ ngập Và số nghiên cứu nớc giới nh: Bangladesh xây dựng thành công hệ thống giám sát cảnh báo ngập lụt sở sử dụng mô hình thuỷ văn thuỷ lực MIKE-11 (của Đan Mạch) dới trợ giúp UNDP/WMO kết hợp với sử dụng t liệu viễn thám GMS, NOAA-12 NOAA-14 Hệ thống giám sát cảnh báo ngập lụt đợc áp dụng cho vùng lãnh thổ rộng 82000 km2, đoạn dài 7270 km sông, 195 nhánh, sử dụng 30 trạm giám sát Trung quốc xây dựng hệ thống giám sát cảnh báo ngập lụt sở sử dụng t liệu viễn thám FY-II, OLR, GPCP, ERS-II, SSM/I ấn Độ bắt đầu xây dựng hệ thống giám sát cảnh báo ngập lụt từ năm 1959 cho lu vực sông Hằng Hiện ấn độ có 145 trung tâm dự báo, 500 trạm khí tợng, 350 trạm thuỷ văn phục vụ cho vùng lu vực rộng 240000 km2, sử dụng khả thông tin t liệu ảnh vệ tinh IRS, TM Landsat-5, ERS, Radarsat Một số nớc thuộc Châu Phi sử dụng mô hình thuỷ văn FEWS NET kết hợp với hệ thống thông tin địa lý GIS để xây dựng hệ thống giám sát cảnh báo ngập lụt cho 5600 vùng hạ lu với trợ giúp xây dựng tổ chức USGS/EROS Nhìn chung, cách thức tiếp cận thực đề tài nớc chủ yếu dựa sở phơng pháp truyền thống sử dụng số liệu thực đo điều tra thực địa bổ sung kết hợp với mô hình số độ cao để chiết tách vết lũ Các kết thu đợc hầu hết đồ ngập lụt, đồ trạng lũ chu kỳ khác nhau, cha có số liệu chi tiết vùng ngập đánh giá nhanh ảnh hởng thiệt hại mà lũ lụt gây Còn đề tài nớc ứng dụng công nghệ viễn thám để nghiên cứu nhiều, song chủ yếu tập trung vào xây dựng mô hình thủy văn, thủy lực để giám sát cảnh báo ngập lụt mà cha đa số liệu cụ thể đánh giá nhanh tình trạng ngập lụt Nhằm tận dụng khả u việt công nghệ viễn thám, đặc biệt viễn thám Radar mang lại, luận văn đề cập cụ thể phân tích lũ lụt dựa vào t liệu SAR t liệu quang học để đa phơng pháp đánh giá nhanh mức độ ảnh hởng thiệt hại lũ lụt nhằm bổ sung phần cha đợc rõ ràng đề tài nớc nêu đề tài nghiên cứu nớc cha đề cập tới Mục tiêu nhiệm vụ luận văn Mục tiêu: ứng dụng t liệu Radar t liệu Quang học để đánh giá nhanh mức độ ảnh hởng lũ lụt lên lớp phủ mặt đất Nhiệm vụ: - Tổng quan phơng pháp tiền xử lý liệu Radarsat - Xử lý liệu Radarsat để chiết tách thông tin vùng ngập lũ - Xử lý liệu Landsat để chiết tách thông tin lớp phủ bề mặt - Tích hợp GIS để phân tích đánh giá mức độ ảnh hởng lên đối tợng cụ thể thuộc lớp phủ bề mặt Tính luận văn ứng dụng công nghệ Viễn thám GIS để đa báo cáo nhanh trạng mức độ ảnh hởng lũ lên đối tợng thuộc lớp phủ bề mặt Vùng nghiên cứu: Căn vào mức độ thiệt hại lũ lụt xảy miền Trung ác liệt nguy hiểm Chủ yếu dạng lũ quét, lũ lên nhanh, rút nhanh, tốc độ lớn, chảy mạnh, ngập dài ngày đồng nhỏ hẹp ven biển (khi thời gian lũ trùng với lúc thuỷ triều lên) mà nơi lại trung tâm dân c có kinh tế phát triển miền Trung, nên hậu mà lũ lụt gây cho tỉnh miền Trung không nhỏ ngời Chính nguyên nhân nêu tác giả định chọn Thừa Thiên Huế khu vực nghiên cứu mà cụ thể Thành phố Huế Đặc điểm điều kiện tự nhiên Vùng nghiên cứu có tọa độ từ vĩ tuyến 16 020 B đến 17015 B từ kinh tuyến 107005 Đ đến 108015 Đ, bao gồm hai lu vực sông Sông Hơng sông lớn tỉnh Thừa Thiên Huế, nằm từ vĩ tuyến 16 000 B đến 16045 B từ kinh tuyến 107000 Đ đến 108015 Đ, phía tây dãy Trờng Sơn, phía Bắc núi Bạch Mã, phía Nam nằm tiếp giáp thành phố Đà Nẵng phía đông giáp biển Đông Sự phân bố dòng chảy bất qui tắc phù hợp với phân bố ma rào Mùa lũ kéo dài từ tháng 10 đến tháng 12 lợng nớc chiếm từ 50 đến 80% tổng số lợng nớc hàng năm, có dao động mạnh năm Lợng nớc vào mùa lũ năm có nhiều trận lũ gấp ba lần so với năm có trận lũ trạm thuỷ văn Thờng Nhật, lu lợng dòng chảy dao động từ 18.5 l/skm2 (mùa khô) tới 228 l/skm2 (mùa lũ) Vùng bờ biển tỉnh Thừa Thiên Huế bị chi phối hệ thống phá Phá phá Tam Giang Phá chạy dọc bờ biển Huế, có diện tích 22000 Phía đông, phá bị ngăn tách với biển cồn cát với hai khe, Thuận An Từ Hiền phía tây phá cánh đồng lúa cửa sông Đây vùng có thắng cảnh tài nguyên sinh học Dân c c trú khai thác tài nguyên sinh học phá làm nông nghiệp Vì độ dốc nhỏ lợng nớc dâng cao từ sông ma, vùng có khả lũ lớn Các dải đá ngầm ven biển với độ cao lên tới 20m gây khó khăn cho nớc chảy tự trờng hợp lũ nghiêm trọng, làm vùng giống nh hồ lớn Đặc điểm kinh tế xã hội Thừa Thiên Huế nằm trục giao thông quan trọng xuyên Bắc-Nam, trục hành lang Đông-Tây nối Thái Lan - Lào - Việt Nam theo đờng Thừa Thiên Huế vào vị trí trung độ nớc, nằm thành phố Hà Nội thành phố Hồ Chí Minh hai trung tâm lớn hai vùng kinh tế phát triển nớc ta, nơi giao thoa điều kiện tự nhiên - kinh tế - xã hội hai miền Nam - Bắc Thừa Thiên Huế trung tâm văn hoá, du lịch, trung tâm giáo dục đào tạo, y tế lớn nớc cực phát triển kinh tế quan trọng vùng kinh tế trọng điểm miền Trung Bờ biển tỉnh dài 120 km, có cảng Thuận An vịnh Chân Mây với độ sâu 18 - 20m đủ điều kiện xây dựng cảng nớc sâu với công suất lớn, có cảng hàng không Phú Bài nằm đờng quốc lộ 1A đờng sắt xuyên Việt chạy dọc theo tỉnh, có 81 km biên giới với Lào Vị trí địa lý nh tạo điều kiện thuận lợi cho Thừa Thiên Huế phát triển sản xuất hàng hoá mở rộng giao lu kinh tế - xã hội với tỉnh nớc quốc tế Tổng hợp thông tin bão, lũ lụt Thừa Thiên Huế năm 1999 Theo Trung tâm Khí tợng Thủy văn quốc gia thống kê thiên tai lũ lụt miền Trung nói chung Thừa Thiên Huế nói riêng ác liệt nguy hiểm nhất, chủ yếu dạng lũ quét, lũ lên nhanh rút nhanh Còn nhớ trận lũ lụt lịch sử xảy năm 1999 miền Trung ma nhiều, tập trung khoảng thời gian ngắn, gây lũ lụt đột ngột vùng núi thuộc tỉnh Quảng Ngãi lũ lớn lu vực sông từ Thừa Thiên - Huế đến Khánh Hòa Đặc biệt, có trận lũ lớn sông tỉnh Quảng Ngãi Mức nớc lũ hầu hết sông đạt mức đỉnh mức báo động III Nhất tỉnh Quảng Ngãi, xảy lũ lịch sử sông Trà Khúc Trên sông Bình Định, Quảng Nam Thừa Thiên - Huế, có lũ lớn mức báo 10 động III Lũ gây ngập lụt diện rộng Chi tiết lũ ngập lụt tỉnh Thừa Thiên Huế nh sau: Huế, ma lớn có nơi 2000 mm vòng ngày từ đến tháng 11 Lũ lên nhanh, với cờng suất có tới 1m nớc trong giờ, biên độ lũ lớn 4-5 m, chí 8-10 m thợng nguồn sông Hơng, sông Bồ, gây lũ quét số nơi Hơn 90% khu dân c, kể vùng gò đồi bị ngập chìm nớc lũ tới 1-4 m -9 ngày Thành phố Huế bị ngập nặng nề, hấu hết nhà cửa, đờng phố chìm biển nớc mênh mông, giao thông liên lạc bị ngừng trệ Đặc biệt, nớc lũ sông Hơng phá vỡ phá Tam Giang, cửa biển Thuận An bị đổi dòng, tạo thành cửa mới, chia cắt cô lập vùng biển, ảnh hởng nghiêm trọng đến môi trờng sinh thái vùng phá Tam Giang vùng cầu Hai Tại hai cửa Hoà Duân T Hiền rộng tới 500-600 m, lũ trôi hộ dân với toàn nhà cửa, tài sản Cụ thể 40 xã huyện Hơng Trà, Quảng Điền, Phú Vang, Hơng Thuỷ Phú Lộc Thành phố Huế bị ngập lụt nghiêm trọng: 54000 nhà bị nhấn chìm nớc, vài địa phơng Phú Lộc chịu ảnh hởng lũ lên đột ngột với 76 nhà ngập khoảng 1.5m (trong có nhà bị trôi) buộc phải di chuyển Cấu trúc luận văn Cấu trúc luận văn bao gồm: Phần mở đầu, Nội dung luận văn, Phần kết luận, Kiến nghị, Tài liệu tham khảo Các phụ lục Trong đó, nội dung luận văn gồm chơng: Chơng 1: Cơ sở phơng pháp luận ứng dụng ảnh viễn thám nghiên cứu ngập lụt Chơng 2: Xử lý ảnh RADAR Chơng 3: Nghiên cứu quy trình tích hợp ảnh Radar ảnh Quang học để đánh giá ảnh hởng tình trạng ngập lụt tới trạng lớp phủ 55 3.4 Phân tích đánh giá kết thực nghiệm Phần phân tích đánh giá kết sử dụng phần thực nghiệm luận văn tổng hợp phơng pháp đánh giá dựa vào tiêu chí nh diện tích thời gian ngập cộng với việc phân tích số liệu tính toán thống kê thu đợc kết hợp với phơng pháp tích hợp thông tin công cụ GIS để đánh giá ảnh hởng lũ lụt tới lớp phủ bề mặt Trong phần phân tích đánh giá này, t liệu sử dụng gồm thời điểm: thời điểm 02/6/2001 cung cấp thông tin lớp phủ bề mặt thành phố Huế, hạn chế mặt t liệu nên đợc dùng để cung cấp số liệu nh thời điểm lũ, thời điểm 06/11/1999 thời điểm xảy lũ thời điểm 15/11/1999 thời điểm rút lũ Biểu đồ 3-1: Diện tích đối tợng thời điểm lũ Dới biểu đồ số liệu mô tả diện tích mặt nớc (không tính phần biển) thời điểm lũ 02/6/2001 (chính tổng diện tích đối tợng thủy hệ đất nuôi trồng thủy sản), thời điểm xảy lũ 06/11/1999 thời điểm lũ rút 15/11/1999 56 Biểu đồ 3-2: Số liệu diện tích mặt nớc Căn vào thời điểm diện tích ngập, sử dụng công cụ GIS để tính toán diện tích bị ảnh hởng lũ đến đối tợng Thuật toán sử dụng để xác định trạng ngập đợc mô tả nh sau: Thực chất trình xử lý xác định trạng vùng ngập toán toán tử logic (hợp giao lớp đối tợng) Đầu vào vùng ngập lớp phủ bề mặt, trình xử lý tính toán lớp đối tợng có vị trí mặt không gian Khi diện tích vùng ngập đợc chồng lên đợc chia cắt theo lớp đối tợng lớp phủ bề mặt Kết đầu lớp chứa chứa đựng giá trị thuộc tính lớp đối tợng đầu vào Kết tính toán trạng vùng ngập đợc thể biểu đồ 3-3 dới đây: 57 Biểu đồ 3-3: Diện tích đối tợng bị ngập ngày 06/11/1999 Xét mức độ ảnh hởng rừng đối tợng khó bị ảnh hởng nằm địa hình núi cao dốc nên có ma lũ lợng nớc bị thấm phần thoát cách nhanh chóng Do đó, đối tợng rừng không đợc đa vào để đánh giá Dới biểu đồ mô tình trạng đối tợng lũ xảy lũ Biểu đồ 3-4: Biểu đồ thể mức độ ảnh hởng theo diện tích đối tợng xảy lũ 58 Từ số liệu diện tích trạng ngập đối tợng đợc tính toán cộng với số liệu độ cao ngập (đợc tham khảo từ báo cáo trận lụt lịch sử Huế năm 1999), sau dựa vào tiêu chí đánh giá theo tầm quan trọng đối tợng, diện tích thời gian ngập để lập bảng đánh giá nhanh mức độ ảnh h ởng lũ đến bề mặt lớp phủ Theo biểu đồ 3-4 thấy lũ xảy ảnh hởng nặng nề nhất, gây hậu trực tiếp mặt kinh tế vùng nông nghiệp kết hợp nuôi trồng thủy sản Với 16400 bị ngập chìm (chiếm 93% tổng số đất nông - thủy sản), với vùng trũng thủy sản coi nh bị trắng bị trôi, vùng đợc đắp cao bị ảnh hởng nhiều nớc tràn vào gây nhạt hóa nguồn nớc ảnh hởng trực tiếp đến loại thủy sản đợc nuôi trồng Nếu tính sơ mang lại vài trăm triệu đồng tổn thất lớn mặt kinh tế mà ngời dân nói riêng Thành phố Huế nói chung phải chịu Do mức độ ảnh hởng đến đất nông - thủy sản liệt vào mức (mức bị ảnh hởng nghiêm trọng) Mặc dù đất dân c bị ảnh hởng chiếm 20.5% tổng diện tích đất dân c nhng thực tế vùng bị ngập nhiều lại vùng nội thành, nơi có mật độ dân c đông, ngời lại có tầm quan trọng lớn ngời dân bị ảnh hởng trực tiếp đến tính mạng thiệt hại mang tổn thất lớn Dựa vào số liệu thời điểm lũ rút hẳn ngày 15/11/1999 (tức sau 10 ngày) tính toán diện tích vùng ngập rút ngày gần 1200 ha, độ cao trung bình Nội đô thành phố Huế khoảng 1.5m nh trung bình ngày nớc rút đợc 15 cm Với số ngày bị ngập lâu nh mà tốc độ rút lại chậm nên ngời dân có nguy nhà cao, bên cạnh giao thông ngập chìm biển nớc làm cắt đờng vận chuyển cứu hộ lơng thực cho ngời dân nên mức độ thiệt hại tính mạng khó tránh khỏi Nh vậy, mức độ thiệt hại dân c đợc liệt mức ảnh hởng nghiêm trọng Với hệ thống giao thông bị ngập chìm lâu nh (có thể tính nhanh tốc độ rút quốc lộ chậm nhanh ngày muộn ngày), thơng mại bị tê liệt tạm thời, quan, xí nghiệp, nhà xởng bị đình trệ gây ảnh hởng nặng nề tới kinh tế Vì mức độ thiệt hại đến hệ thống đờng giao thông mà chủ yếu đờng đợc liệt vào mức ảnh hởng nghiêm trọng Đất nông nghiệp với diện tích bị ảnh hởng 7000 (chiếm 22.5% tổng diện tích đất nông nghiệp) mức độ thiệt hại theo diện tích cha phải lớn nhng thời gian ngập lâu nh gây trắng lúa, màu mà làm ảnh 59 hởng tới chất lợng đất để canh tác sau lũ rút Do mức độ thiệt đất nông nghiệp xếp vào mức (mức bị ảnh hởng trung bình) Còn loại đất khác chủ yếu cát bãi bồi, diện tích bị ảnh hởng 2980 (chiếm 39.4% tổng số loại đất khác), xét mức độ quan trọng loại đất cha có tầm quan trọng mức độ ảnh đối tuợng đa vào mức (mức bị ảnh hởng nhẹ) mức (mức không bị ảnh hởng) Qua phân cấp mức độ ảnh hởng nh sau: - Mức 0: Không bị ảnh hởng - Mức 1: Mức độ ảnh hởng nhẹ - Mức 2: Mức độ ảnh hởng trung bình - Mức 3: Mức độ ảnh hởng nghiêm trọng Dới kết đánh giá mức độ ảnh hởng lũ đến lớp phủ bề mặt: Tên đối tợng Diện tích ngập (Ha) Chiều dài ngập (Km) Đất nông - thủy sản Đất canh tác Dân c Đất khác Giao thông 16,412.22 7,275.42 477.63 2,980.35 Độ sâu ngập (m) Thời gian ngập (ngày) 670.83 2m 2m 1.5 m 1.5 m 2m Mức độ ảnh hởng ngập theo thời gian 10 10 10 10 10 Mức Mức Mức Mức Mức Bảng 3-1: Mức độ ảnh hởng đối tợng lũ gây Dựa vào bảng tiêu chí minh họa mức độ ảnh hởng đối tợng xảy lũ Chi tiết đợc minh họa Bản đồ mức độ ảnh hởngtrong phần phụ lục: 60 Hình 3-10: Bản đồ mức độ ảnh hởng ngập thời điểm lũ 06/11/1999 Cho đến cha có tiêu chí cụ thể để đánh giá mức độ ảnh hởng lũ lụt Những thiệt hại mang tính chất xã hội thờng khó đánh giá, nhng lại trầm trọng nhất, lâu dài dới tác động lũ, lụt Trong số hậu trớc mắt, thấy lũ, lụt làm nhiều nhà cửa, nơi c trú nhân dân, công trình công cộng, dân sinh bị phá huỷ h hại làm nẩy sinh hậu xã hội khó lờng hết đợc, chảng hạn nh việc giải chỗ ở, sản xuất nông nghiệp, ổn định đời sống, chống dịch bệnh, cải tạo môi trờng nớc, đất, Lũ, lụt phá hoại sở vật chất kỹ thuật, làm cạn kiệt nguồn lơng thực, thực phẩm, phá hoại mùa màng, đất đai, làm cho không canh tác đợc, gây nạn đói, ổn định đời sống văn hóa, xã hội, môi trờng, diện rộng gây hậu lâu dài, mà phải nhiều năm sau khắc phục đợc 61 Kết luận Trên sở nghiên cứu thử nghiệm luận văn, rút vài kết luận nh sau: Lũ lụt tợng thiên tai thờng xuyên xảy Việt Nam gây hậu nghiêm trọng nên cần có công cụ giám sát để dự báo nh đánh giá nhanh thiệt hại nhằm khắc phục hậu Để đạt đợc mục tiêu nói luận văn đề xuất sử dụng ảnh viễn thám, cụ thể kết hợp ảnh Radar Quang học ảnh Radar t liệu Việt Nam có kỹ thuật xử lý phức tạp, nhiên kết nghiên cứu luận văn đề xuất đợc quy trình xử lý ảnh Radar để chiết tách thông tin vùng ngập Quy trình nêu phần xử lý ảnh Radar quy trình cần thiết ứng dụng ảnh Radar nghiên cứu Trong định chuẩn ảnh khâu phức tạp nhng quan trọng cho việc phân loại có sở đối tợng có phản hồi tơng tự gần Tham số dùng để định chuẩn cho loại vệ tinh khác nhau, thông thờng quan quản lý vệ tinh cung cấp mua ảnh Nh ảnh Radarsat quan CDPF (Canadian Data Processing Facility) cung cấp Loại bỏ nhiễu ảnh Radar công đoạn đặc biệt quan trọng, định phần lớn cho thành công kết nghiên cứu thu đợc Nếu sử dụng ảnh đa thời gian cần thiết phải sử dụng phép lọc đa thời gian trớc sử dụng nhóm lọc cục Nếu sử dụng thời điểm ảnh nên sử dụng phép lọc Frost trớc tiến hành xử lý số T liệu W2 (Wide Beam Positions) ảnh Radarsat góc tới xấp xỉ 30 cộng với phân cực HH (phân cực đặc trng cho đặc tính phản hồi nớc ảnh Radar) t liệu thích hợp cho việc nghiên cứu lũ Luận văn xây dựng đợc hệ phân loại lớp phủ để sử dụng ảnh Quang học với kỹ thuật xử lý thông tin đa phổ để chiết tách đối tợng lớp phủ vùng nghiên cứu Thực nghiệm quy trình kết hợp t liệu Radar Quang học cho kết Bản đồ mức độ ảnh hởng lũ lụt lớp phủ mặt đất Thừa Thiên Huế Mặc dù đợc thử nghiệm vùng, song quy trình ứng dụng để đánh giá ảnh hởng trân lũ khác lớp phủ bề mặt Việt Nam 62 Kiến nghị Để giám sát nhanh nghiên cứu lũ cần sử dụng loạt ảnh Radar đa thời gian Cần kết hợp dự báo lũ với lập trình thu ảnh Trạm thu ảnh vệ tinh thuộc dự án Giám sát Tài nguyên thiên Môi trờng Bộ Tài nguyên Môi trờng để thu đợc ảnh thời điểm "nhạy cảm" lũ Cần kết hợp với mô hình thủy văn, thủy lực nhằm tính cao trình ngập, đa mức độ ngập cụ thể cho vùng, đồng thời dự báo vùng có khả bị ảnh hởng để đa biện pháp phòng chống cứu hộ kịp thời Cần nghiên cứu xây dựng chi tiết tiêu đánh giá mức độ thiệt hại cho đối tợng thuộc lớp phủ bề mặt để phục vụ nhanh chóng công tác đánh giá thiệt hại lũ lụt gây Khi đánh giá nhanh mức độ thiệt hại ảnh hởng lũ lụt cần phải có đầy đủ sở liệu GIS DEM với độ xác cao 63 Tài liệu tham khảo Tiếng Việt Phạm Văn Cự Ferdinand Bonn (2006) Giáo trình Viễn thám Radar Nguyễn Xuân Lâm (2003), "Nghiên cứu kinh nghiệm Thái Lan ứng dụng công nghệ Viễn thám phục vụ công tác quản lý tài nguyên thiên nhiên môi trờng Việt Nam, trớc hết tài nguyên đất nớc", Đề tài nghiên cứu khoa học theo Nghị định th hợp tác Việt Nam-Thái Lan Lê Bắc Huỳnh nnk (2000) Nghiên cứu thiên tai lũ, Dự án UNDP VIE/97/002- Diasaster Management Unit, Chơng trình phát triển Liên hiệp quốc, Việt Nam Trung tâm ứng dụng Công nghệ Vũ trụ (2001), Xây dựng đồ trạng sử dụng đất từ ảnh vệ tinh Radar phục vụ nghiên cứu đánh giá biến động trạng sử dụng đất đồng sông Cửu Long, Báo cáo đề tài, Viện Vật lý, Trung tâm Khoa học tự nhiên Công nghệ Quốc gia Chu Hải Tùng (2006), "Nghiên cứu khả ứng dụng kết hợp ảnh vệ tinh Radar Quang học để thành lập số lớp thông tin lớp phủ mặt đất" (2006), Đề tài khoa học cấp Bộ, Trung tâm Viễn thám, Bộ Tài nguyên Môi trờng Nguyễn Đình Dơng (2005), Nghiên cứu mô hình đánh giá phục giám sát tài nguyên môi trờng Hệ thông tin địa lý, Báo cáo chuyên đề, Viện Địa lý, Viện Khoa học Công nghệ Việt Nam Phòng Viễn thám Hệ thông tin địa lý (2000), Bớc đầu tiếp cận phơng pháp xử lý ảnh Radar qua việc theo dõi sinh trởng lúa ô mẫu vùng Hóc Môn - Củ Chi, Báo cáo nghiên cứu khoa học, Phân Viện Địa lý TP Hồ Chí Minh, Trung tâm Khoa học tự nhiên Công nghệ Quốc gia Trần Tuấn Đạt (2007), "Nghiên cứu ứng dụng công nghệ Viễn thám Hệ thống thông tin địa lý quản lý tổng hợp lu vực sông" (2007), Đề tài khoa học cấp sở, Trung tâm Viễn thám, Bộ Tài nguyên Môi trờng Trang thông tin điện tử UBND Tỉnh Thừa Thiên Huế, http://www.thuathienhue.gov.vn 10 Trang thông tin điện tử Quản lý đê điều phòng chống lụt bão Việt Nam, http://www.ccfsc.org.vn Tiếng Anh 64 Canadian Space Agency, "Extraction of Beta Nought and Sigma Nought", RADARSAT CDPF Products (2000) RADARSAT International, "Using your Radarsat Product", Radarsat Illuminated Chris Oliver., Shaun Quegan (1997), "Principles of SAR Image Formation", "Fundamental Properties of SAR Images", Understanding Synthetic Aperture Radar Images, SCITECH Publishing, Inc, pp 11-42, 75-99 Alexandre Bouvet, Instructions to use the multi-channel SAR filter, CESBIO, CNES, France Hervé Yesou, "Flood Mapping and Monitoring", Flood monitoring training course, ASAR Training course in Remote Sensing Center, MONRE, Hanoi, Vietnam (2008) Hoang Minh Hien, Tran Nhu Trung, Wim Looijen and Kees Hulsbergen, Flood Vulnerability Analysis and Mapping in Vietnam Robbert Misdorp, Hua Chien Thang, Nguyen Xuan Lam , Using Remote Sensing Data for Coastal TT-Hue Province, Vietnam, Providing information for Intergrated Coastal Zone Management Priero Villegas (2004), Flood Modelling in Perfume River Basin, Hue Province, Vietnam, Master of Science Thesis, International Institute for Geoinformation Science and Earth Observation, The Netherlands, pp 37-39, 41-47 United Nations Development Programme (UNDP), Immediate report on flood and rain in the central provinces from Thua Thien-Hue to Khanh Hoa -06 Dec 1999 10 Mekong River Commission Website, Monitoring floods with Radarsat, Examples of Earth Observation Data Products and Applications, http://www.mrcmekong.org 11 Dartmouth Flood Observatory Website, http://dartmouth.edu/~floods 65 Phụ lục 66 Bảng Radiometric Data Record ảnh Radarsat thời điểm 06/11/1999 0.000000E+00 0.000000E+00 0.000000E+00 _2 8450E+08 1.9626400E+08 1.9654339E+08 1.9682330E+08 1.9710360E+08 1.9738430E+08 1.9766499E+08 1.9794610E+08 1.9822760E+08 1.9850950E+08 1.9879190E+08 1.9907350E+08 1.9932061E+08 1.9956870E+08 1.9981630E+08 2.0006461E+08 2.0031290E+08 2.0056150E+08 2.0081040E+08 2.0105960E+08 2.0130880E+08 2.0155859E+08 2.0180840E+08 2.0205850E+08 2.0230960E+08 2.0259630E+08 2.0288301E+08 2.0317010E+08 2.0345730E+08 2.0374520E+08 2.0403320E+08 2.0432160E+08 2.0461000E+08 2.0489920E+08 2.0518840E+08 2.0547810E+08 2.0576810E+08 2.0605830E+08 2.0631320E+08 2.0656781E+08 2.0682240E+08 2.0707760E+08 2.0733280E+08 2.0758819E+08 2.0784370E+08 2.0809981E+08 2.0835619E+08 2.0861261E+08 2.0886930E+08 2.0912630E+08 2.0938320E+08 2.0967370E+08 2.0996819E+08 2.1026310E+08 2.1055810E+08 2.1085421E+08 2.1114960E+08 2.1144621E+08 2.1174250E+08 2.1203910E+08 2.1233661E+08 2.1263410E+08 2.1293200E+08 2.1323040E+08 2.1350061E+08 2.1376221E+08 2.1402370E+08 2.1428590E+08 2.1454810E+08 2.1481021E+08 2.1507299E+08 2.1533610E+08 2.1559920E+08 2.1586259E+08 2.1612630E+08 2.1639040E+08 2.1665430E+08 2.1691859E+08 2.1718360E+08 2.1744861E+08 2.1771350E+08 2.1797910E+08 2.1824461E+08 2.1851040E+08 2.1877661E+08 2.1904310E+08 2.1930960E+08 2.1957670E+08 2.1984381E+08 2.2011080E+08 2.2039101E+08 2.2069710E+08 2.2100360E+08 2.2131050E+08 2.2161779E+08 2.2192520E+08 2.2223341E+08 2.2254170E+08 2.2285000E+08 2.2315910E+08 2.2346819E+08 2.2377770E+08 2.2408770E+08 2.2439699E+08 2.2466840E+08 2.2494080E+08 2.2521290E+08 2.2548520E+08 2.2575830E+08 2.2603139E+08 2.2630430E+08 2.2657760E+08 2.2685160E+08 2.2712560E+08 2.2739981E+08 2.2767440E+08 2.2794899E+08 2.2825110E+08 2.2856610E+08 2.2888110E+08 2.2919661E+08 2.2951210E+08 2.2982840E+08 2.3014470E+08 2.3046150E+08 2.3077830E+08 2.3109600E+08 2.3141370E+08 2.3173181E+08 2.3205040E+08 2.3235710E+08 2.3263610E+08 2.3291581E+08 2.3319539E+08 2.3347539E+08 2.3375570E+08 2.3403630E+08 2.3431690E+08 2.3459779E+08 2.3487901E+08 2.3516061E+08 2.3544250E+08 2.3572440E+08 2.3600650E+08 2.3628950E+08 2.3657190E+08 2.3685510E+08 2.3713821E+08 2.3742210E+08 2.3770590E+08 2.3799000E+08 2.3827440E+08 2.3855880E+08 2.3884360E+08 2.3912861E+08 2.3941400E+08 2.3969981E+08 2.3998550E+08 2.4027190E+08 2.4055821E+08 2.4084490E+08 2.4113200E+08 2.4141890E+08 2.4170621E+08 2.4199379E+08 2.4228221E+08 2.4257010E+08 2.4285840E+08 2.4314739E+08 2.4343630E+08 2.4372560E+08 2.4404819E+08 2.4437981E+08 2.4471150E+08 2.4504310E+08 2.4537570E+08 2.4570781E+08 2.4604090E+08 2.4637430E+08 2.4670781E+08 2.4704179E+08 2.4737619E+08 2.4771101E+08 2.4804630E+08 2.4837290E+08 2.4866670E+08 2.4896090E+08 2.4925501E+08 2.4954979E+08 2.4984461E+08 2.5013970E+08 2.5043520E+08 2.5073101E+08 2.5102680E+08 2.5132280E+08 2.5161970E+08 2.5191600E+08 2.5221320E+08 2.5251021E+08 2.5280760E+08 2.5310570E+08 2.5340339E+08 2.5370179E+08 2.5400019E+08 2.5429930E+08 2.5459830E+08 2.5489770E+08 2.5519699E+08 2.5549710E+08 2.5579710E+08 2.5609750E+08 2.5639821E+08 2.5669880E+08 2.5700019E+08 2.5730200E+08 2.5760360E+08 2.5790520E+08 2.5820760E+08 2.5851030E+08 2.5881290E+08 2.5911640E+08 2.5941970E+08 2.5972339E+08 2.6002701E+08 2.6033101E+08 2.6063850E+08 2.6098670E+08 2.6133490E+08 2.6168400E+08 2.6203320E+08 2.6238280E+08 2.6273290E+08 2.6308301E+08 2.6343400E+08 2.6378499E+08 2.6413650E+08 2.6448800E+08 2.6484040E+08 2.6519240E+08 2.6553510E+08 2.6584400E+08 2.6615341E+08 2.6646310E+08 2.6677320E+08 2.6708310E+08 2.6739341E+08 2.6770410E+08 2.6801510E+08 2.6832610E+08 2.6863779E+08 2.6895002E+08 2.6926150E+08 2.6957430E+08 2.6988659E+08 2.7019920E+08 2.7051219E+08 2.7082560E+08 2.7113971E+08 2.7145331E+08 2.7176781E+08 67 2.7208211E+08 2.7365830E+08 2.7524160E+08 2.7683162E+08 2.7842842E+08 2.8003200E+08 2.8164230E+08 2.8325939E+08 2.8488371E+08 2.8651440E+08 2.8815238E+08 2.8979718E+08 2.9144880E+08 2.9320659E+08 2.9511139E+08 2.9702541E+08 2.9881219E+08 3.0050202E+08 3.0219920E+08 3.0390330E+08 3.0561440E+08 3.0733290E+08 3.0905840E+08 3.1079091E+08 3.1253030E+08 3.1427670E+08 3.1603059E+08 3.1779101E+08 3.1955891E+08 3.2133440E+08 3.2311642E+08 3.2490541E+08 3.2670202E+08 3.2850621E+08 3.3031690E+08 3.3213469E+08 3.3396010E+08 3.3579261E+08 3.3763219E+08 3.3947949E+08 3.4133402E+08 3.4319549E+08 3.4506410E+08 3.4693990E+08 3.4882339E+08 3.5060442E+08 2.7239680E+08 2.7397478E+08 2.7555901E+08 2.7715040E+08 2.7874851E+08 2.8035341E+08 2.8196509E+08 2.8358349E+08 2.8520922E+08 2.8684179E+08 2.8848109E+08 2.9012720E+08 2.9178019E+08 2.9358701E+08 2.9549318E+08 2.9740970E+08 2.9914931E+08 3.0084109E+08 3.0253958E+08 3.0424518E+08 3.0595770E+08 3.0767760E+08 3.0940390E+08 3.1113830E+08 3.1287910E+08 3.1462691E+08 3.1638170E+08 3.1814400E+08 3.1991331E+08 3.2168970E+08 3.2347360E+08 3.2526461E+08 3.2706259E+08 3.2886771E+08 3.3067981E+08 3.3249958E+08 3.3432640E+08 3.3616038E+08 3.3800141E+08 3.3984950E+08 3.4170531E+08 3.4356829E+08 3.4543830E+08 3.4731610E+08 3.4920109E+08 -25.0200005 2.7271190E+08 2.7302730E+08 2.7429069E+08 2.7460749E+08 2.7587670E+08 2.7619482E+08 2.7746941E+08 2.7778890E+08 2.7906899E+08 2.7938979E+08 2.8067520E+08 2.8099741E+08 2.8228819E+08 2.8261171E+08 2.8390851E+08 2.8423341E+08 2.8553510E+08 2.8586131E+08 2.8716899E+08 2.8749610E+08 2.8880970E+08 2.8913811E+08 2.9045722E+08 2.9078749E+08 2.9211149E+08 2.9244730E+08 2.9396739E+08 2.9434819E+08 2.9587549E+08 2.9625882E+08 2.9779398E+08 2.9813789E+08 2.9948691E+08 2.9982490E+08 3.0118000E+08 3.0151939E+08 3.0288000E+08 3.0322070E+08 3.0458691E+08 3.0492899E+08 3.0630131E+08 3.0664480E+08 3.0802211E+08 3.0836749E+08 3.0975040E+08 3.1009722E+08 3.1148560E+08 3.1183322E+08 3.1322781E+08 3.1357741E+08 3.1497750E+08 3.1532800E+08 3.1673430E+08 3.1708611E+08 3.1849741E+08 3.1885069E+08 3.2026819E+08 3.2062278E+08 3.2204589E+08 3.2240250E+08 3.2383130E+08 3.2418870E+08 3.2562362E+08 3.2598310E+08 3.2742301E+08 3.2778390E+08 3.2922950E+08 3.2959181E+08 3.3104310E+08 3.3140669E+08 3.3286429E+08 3.3322931E+08 3.3469251E+08 3.3505901E+08 3.3652790E+08 3.3689581E+08 3.3837030E+08 3.3873958E+08 3.4022042E+08 3.4059120E+08 3.4207770E+08 3.4244989E+08 3.4394202E+08 3.4431571E+08 3.4581360E+08 3.4618861E+08 3.4769280E+08 3.4806918E+08 3.4957859E+08 3.4995520E+08 0.0000000E+00 2.7334259E+08 2.7492461E+08 2.7651331E+08 2.7810829E+08 2.7971091E+08 2.8131939E+08 2.8293562E+08 2.8455811E+08 2.8618739E+08 2.8782400E+08 2.8946749E+08 2.9111821E+08 2.9282669E+08 2.9472950E+08 2.9664160E+08 2.9847430E+08 3.0016320E+08 3.0185910E+08 3.0356179E+08 3.0527210E+08 3.0698870E+08 3.0871280E+08 3.1044381E+08 3.1218189E+08 3.1392678E+08 3.1567939E+08 3.1743830E+08 3.1920490E+08 3.2097840E+08 3.2275901E+08 3.2454710E+08 3.2634230E+08 3.2814522E+08 3.2995440E+08 3.3177082E+08 3.3359482E+08 3.3542589E+08 3.3726410E+08 3.3910941E+08 3.4096240E+08 3.4282250E+08 3.4468970E+08 3.4656400E+08 3.4844611E+08 3.5028000E+08 Bảng Radiometric Data Record ảnh Radarsat thời điểm 15/11/1999 0.000000E+00 0.000000E+00 0.000000E+00 _2 7310E+08 1.0532530E+08 1.0547770E+08 1.0563010E+08 1.0578290E+08 1.0593580E+08 1.0608910E+08 1.0623430E+08 1.0637070E+08 1.0650720E+08 1.0664400E+08 1.0678090E+08 1.0691790E+08 1.0705510E+08 1.0719230E+08 1.0732980E+08 1.0746760E+08 1.0760530E+08 1.0774320E+08 1.0788130E+08 1.0801950E+08 1.0815770E+08 1.0829630E+08 1.0843510E+08 1.0857390E+08 68 1.0871280E+08 1.0941020E+08 1.1011130E+08 1.1081610E+08 1.1152490E+08 1.1223750E+08 1.1295380E+08 1.1367420E+08 1.1439840E+08 1.1512640E+08 1.1585840E+08 1.1659420E+08 1.1733390E+08 1.1807790E+08 1.1882560E+08 1.1957740E+08 1.2033320E+08 1.2109300E+08 1.2185680E+08 1.2262490E+08 1.2339680E+08 1.2417290E+08 1.2495330E+08 1.2573760E+08 1.2652580E+08 1.2731840E+08 1.2811510E+08 1.2891610E+08 1.2972130E+08 1.3053050E+08 1.3134410E+08 1.3216180E+08 1.3298390E+08 1.3381020E+08 1.3464070E+08 1.3547560E+08 1.3631480E+08 1.3715840E+08 1.3800600E+08 1.3885840E+08 1.3971501E+08 1.4049870E+08 1.4125520E+08 1.4201490E+08 1.4285610E+08 1.4373290E+08 1.4461400E+08 1.4549981E+08 1.4639000E+08 1.4728450E+08 1.4818360E+08 1.4908730E+08 1.4999539E+08 1.5090830E+08 1.5182570E+08 1.0885220E+08 1.0955000E+08 1.1025200E+08 1.1095760E+08 1.1166720E+08 1.1238040E+08 1.1309750E+08 1.1381860E+08 1.1454360E+08 1.1527250E+08 1.1600520E+08 1.1674180E+08 1.1748240E+08 1.1822720E+08 1.1897580E+08 1.1972830E+08 1.2048490E+08 1.2124570E+08 1.2201020E+08 1.2277910E+08 1.2355170E+08 1.2432880E+08 1.2510960E+08 1.2589480E+08 1.2668410E+08 1.2747730E+08 1.2827500E+08 1.2907660E+08 1.2988270E+08 1.3069290E+08 1.3150730E+08 1.3232580E+08 1.3314870E+08 1.3397580E+08 1.3480730E+08 1.3564301E+08 1.3648301E+08 1.3732739E+08 1.3817630E+08 1.3902920E+08 1.3988670E+08 1.4064979E+08 1.4140690E+08 1.4216739E+08 1.4303110E+08 1.4390870E+08 1.4479101E+08 1.4567730E+08 1.4656861E+08 1.4746400E+08 1.4836410E+08 1.4926840E+08 1.5017760E+08 1.5109139E+08 1.5200960E+08 1.0899130E+08 1.0969020E+08 1.1039270E+08 1.1109920E+08 1.1180950E+08 1.1252360E+08 1.1324150E+08 1.1396330E+08 1.1468890E+08 1.1541860E+08 1.1615220E+08 1.1688960E+08 1.1763110E+08 1.1837640E+08 1.1912570E+08 1.1987920E+08 1.2063650E+08 1.2139810E+08 1.2216350E+08 1.2293340E+08 1.2370700E+08 1.2448460E+08 1.2526650E+08 1.2605230E+08 1.2684240E+08 1.2763660E+08 1.2843490E+08 1.2923770E+08 1.3004440E+08 1.3085540E+08 1.3167050E+08 1.3249000E+08 1.3331400E+08 1.3414180E+08 1.3497410E+08 1.3581090E+08 1.3665160E+08 1.3749680E+08 1.3834630E+08 1.3920050E+08 1.4004650E+08 1.4080090E+08 1.4155880E+08 1.4233230E+08 1.4320640E+08 1.4408480E+08 1.4496781E+08 1.4585530E+08 1.4674710E+08 1.4764360E+08 1.4854450E+08 1.4944990E+08 1.5036000E+08 1.5127459E+08 1.5219370E+08 1.0913090E+08 1.0983040E+08 1.1053360E+08 1.1124090E+08 1.1195200E+08 1.1266680E+08 1.1338560E+08 1.1410810E+08 1.1483450E+08 1.1556510E+08 1.1629950E+08 1.1703750E+08 1.1777980E+08 1.1852600E+08 1.1927610E+08 1.2003050E+08 1.2078860E+08 1.2155080E+08 1.2231730E+08 1.2308770E+08 1.2386200E+08 1.2464060E+08 1.2542320E+08 1.2621010E+08 1.2700080E+08 1.2779590E+08 1.2859520E+08 1.2939880E+08 1.3020620E+08 1.3101820E+08 1.3183420E+08 1.3265450E+08 1.3347920E+08 1.3430810E+08 1.3514120E+08 1.3597850E+08 1.3682019E+08 1.3766650E+08 1.3851699E+08 1.3937181E+08 1.4019720E+08 1.4095210E+08 1.4171059E+08 1.4250670E+08 1.4338160E+08 1.4426110E+08 1.4514490E+08 1.4603350E+08 1.4692610E+08 1.4782350E+08 1.4872520E+08 1.4963150E+08 1.5054240E+08 1.5145810E+08 1.5237830E+08 1.0927040E+08 1.0997080E+08 1.1067490E+08 1.1138270E+08 1.1209470E+08 1.1281020E+08 1.1352970E+08 1.1425320E+08 1.1498040E+08 1.1571150E+08 1.1644690E+08 1.1718580E+08 1.1792870E+08 1.1867560E+08 1.1942680E+08 1.2018180E+08 1.2094080E+08 1.2170380E+08 1.2247100E+08 1.2324220E+08 1.2401750E+08 1.2479680E+08 1.2558040E+08 1.2636780E+08 1.2715950E+08 1.2795540E+08 1.2875540E+08 1.2955980E+08 1.3036820E+08 1.3118100E+08 1.3199800E+08 1.3281930E+08 1.3364460E+08 1.3447430E+08 1.3530840E+08 1.3614650E+08 1.3698930E+08 1.3783610E+08 1.3868770E+08 1.3954330E+08 1.4034770E+08 1.4110370E+08 1.4186280E+08 1.4268130E+08 1.4355720E+08 1.4443730E+08 1.4532230E+08 1.4621160E+08 1.4710541E+08 1.4800350E+08 1.4890610E+08 1.4981350E+08 1.5072530E+08 1.5164179E+08 1.5256280E+08 69 1.5274760E+08 1.5360330E+08 1.5441770E+08 1.5523530E+08 1.5612190E+08 1.5706490E+08 1.5801280E+08 1.5896530E+08 1.5992240E+08 1.6088410E+08 1.6179910E+08 1.6264830E+08 1.6350120E+08 1.6440030E+08 1.6538381E+08 1.6637190E+08 1.6736510E+08 1.6836301E+08 1.6936570E+08 1.7035341E+08 1.7123910E+08 1.7212819E+08 1.7302680E+08 1.7405181E+08 1.7508150E+08 1.7611650E+08 1.7715630E+08 1.7820099E+08 1.7925059E+08 1.8020130E+08 1.8112770E+08 1.8205800E+08 1.8307821E+08 1.8415070E+08 1.8522850E+08 1.8624490E+08 1.8719581E+08 1.8815090E+08 1.5293261E+08 1.5376610E+08 1.5458080E+08 1.5539901E+08 1.5631021E+08 1.5725410E+08 1.5820280E+08 1.5915640E+08 1.6011440E+08 1.6107699E+08 1.6196850E+08 1.6281870E+08 1.6367210E+08 1.6459661E+08 1.6558090E+08 1.6657021E+08 1.6756430E+08 1.6856310E+08 1.6956659E+08 1.7053030E+08 1.7141670E+08 1.7230640E+08 1.7323120E+08 1.7425730E+08 1.7528800E+08 1.7632390E+08 1.7736490E+08 1.7841059E+08 1.7946141E+08 1.8038630E+08 1.8131360E+08 1.8224450E+08 1.8329219E+08 1.8436590E+08 1.8544470E+08 1.8643480E+08 1.8738670E+08 -23.5000000 1.5311650E+08 1.5327880E+08 1.5392880E+08 1.5409139E+08 1.5474410E+08 1.5490781E+08 1.5556320E+08 1.5574610E+08 1.5649870E+08 1.5668739E+08 1.5744370E+08 1.5763299E+08 1.5839299E+08 1.5858370E+08 1.5934750E+08 1.5953910E+08 1.6030659E+08 1.6049880E+08 1.6127040E+08 1.6146000E+08 1.6213821E+08 1.6230810E+08 1.6298899E+08 1.6315939E+08 1.6384330E+08 1.6401440E+08 1.6479299E+08 1.6498970E+08 1.6577850E+08 1.6597610E+08 1.6676870E+08 1.6696720E+08 1.6776360E+08 1.6796330E+08 1.6876330E+08 1.6896390E+08 1.6976810E+08 1.6996950E+08 1.7070720E+08 1.7088450E+08 1.7159440E+08 1.7177210E+08 1.7248480E+08 1.7266330E+08 1.7343610E+08 1.7364099E+08 1.7446290E+08 1.7466899E+08 1.7549480E+08 1.7570179E+08 1.7653160E+08 1.7673981E+08 1.7757350E+08 1.7778230E+08 1.7862040E+08 1.7883040E+08 1.7964720E+08 1.7983190E+08 1.8057139E+08 1.8075680E+08 1.8149939E+08 1.8168541E+08 1.8243699E+08 1.8265059E+08 1.8350680E+08 1.8372120E+08 1.8458150E+08 1.8479661E+08 1.8566080E+08 1.8586541E+08 1.8662480E+08 1.8681510E+08 1.8757730E+08 1.8776819E+08 0.0000000E+00 1.5344099E+08 1.5425450E+08 1.5507120E+08 1.5593390E+08 1.5687600E+08 1.5782280E+08 1.5877440E+08 1.5973059E+08 1.6069150E+08 1.6162950E+08 1.6247810E+08 1.6333030E+08 1.6420410E+08 1.6518659E+08 1.6617390E+08 1.6716590E+08 1.6816290E+08 1.6916470E+08 1.7017120E+08 1.7106170E+08 1.7195021E+08 1.7284210E+08 1.7384610E+08 1.7487530E+08 1.7590899E+08 1.7694760E+08 1.7799170E+08 1.7904040E+08 1.8001650E+08 1.8094230E+08 1.8187160E+08 1.8286440E+08 1.8393570E+08 1.8501261E+08 1.8605520E+08 1.8700520E+08 1.8795960E+08