Tích hợp viễn thám và hệ thống thông tin địa lý đánh giá biến động đất đai giai đoạn 2005 2014 huyện thanh trì thành phố hà nội

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP & PTNT HỌC VIỆN NÔNG NGHIỆP VIỆT NAM NGUYỄN TUẤN ANH TÍCH HỢP VIỄN THÁM VÀ HỆ THỐNG THƠNG TIN ĐỊA LÝ ĐÁNH GIÁ BIẾN ĐỘNG ĐẤT ĐAI GIAI ĐOẠN 20052014 HUYỆN THANH TRÌ - THÀNH PHỐ HÀ NỘI CHUYÊN NGÀNH: QUẢN LÝ ĐẤT ĐAI MÃ SỐ: 60.85.01.03 NGƯỜI HƯỚNG DẪN KHOA HỌC: TS LÊ THỊ GIANG HÀ NỘI, NĂM 2015 LỜI CAM ĐOAN Tôi xin cam đoan kết nghiên cứu Các số liệu, kết nêu luận văn trung thực chưa cơng bố cơng trình Tôi xin cam đoan rằng, giúp đỡ trình thực luận văn cám ơn, thơng tin trích dẫn rõ nguồn gốc Hà Nội, ngày 29 tháng 09 năm 2015 Tác giả luận văn Nguyễn Tuấn Anh ii LỜI CẢM ƠN Trong q trình học tập hồn thành luận văn tốt nghiệp, nhận nhiều động viên, giúp đỡ thầy cơ, bạn bè gia đình Tôi muốn bày tỏ cảm ơn sâu sắc tới tất người Tơi xin bày tỏ lịng biết ơn với giáo, Tiến sĩ Lê Thị Giang, giảng viên Bộ môn Hệ thống thông tin đất đai, Khoa Quản lý đất đai,Học viện Nông nghiệp Việt Nam người định hướng cho lựa chọn đề tài, đưa nhận xét quý giá trực tiếp hướng dẫn suốt trình nghiên cứu hồn thiện khóa luận Xın cảm ơn thầy, cô giáo Khoa Quản lý đất đai, Học viện Nông nghiệp Việt Nam dạy bảo tận tình cho tơi suốt thời gian học tập trường Xin cảm ơn anh, chị phòng Tài nguyên – Mơi trường huyện Thanh Trì, Thành phố Hà Nội tạo điều kiện cho tơi q trình thực đề tài Xin cảm ơn gia đình bạn bè tạo điều kiện, giúp đỡ, động viên suốt q trình học tập hồn thành luận văn tốt nghiệp Hà Nội, ngày tháng năm 2015 Học viên thực Nguyễn Tuấn Anh iii MỤC LỤC Lời cam đoan ii Lời cảm ơn iii Mục lục iv Danh mục bảng vi Danh mục hình vii MỞ ĐẦU 1.1 Tính cấp thiết đề tài: .1 1.2 Mục đích đề tài .3 1.3 Yêu cầu đề tài .3 CHƯƠNG TỔNG QUAN TÀI LIỆU 1.1 Biến động đất đai 1.1.1 Khái quát chung đồ biến động sử dụng đất 1.1.2 Các phương pháp nghiên cứu biến động sử dụng đất 1.2 Viễn thám .8 1.2.1 Khái niệm chung công nghệ viễn thám 1.2.2 Lịch sử hình thành phát triển công nghệ viễn thám 10 1.2.3 Các loại vệ tinh viễn thám 13 1.2.4 Phương pháp xử lý thông tin viễn thám 20 1.3 Hệ thống thông tin địa lý (GIS) 23 1.3.1 Khái niệm GIS 23 1.3.2 Dữ liệu hệ thống thông tin địa lý 26 1.3.3 Phân tích xử lý số liệu GIS 26 1.4 Tình hình nghiên cứu ngồi nước 28 1.4.1 Ứng dụng nghiên cứu công nghệ viễn thám kết hợp với GIS giới 28 1.4.2 Ứng dụng nghiên cứu công nghệ viễn thám kết hợp với GIS Việt Nam 29 CHƯƠNG NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU 32 2.1 Đối tượng phạm vi nghiên cứu 32 iv 2.1.1 Đối tượng nghiên cứu 32 2.1.2 Phạm vi nghiên cứu 32 2.2 Nội dung nghiên cứu 32 2.3 Phương pháp nghiên cứu 33 2.3.1 Phương pháp điều tra, thu thập số liệu thứ cấp 33 2.3.2 Phương pháp thu thập số liệu sơ cấp 33 2.3.3 Phương pháp xử lý số liệu 33 CHƯƠNG KẾT QUẢ VÀ THẢO LUẬN 35 3.1 Điều kiện tự nhiên - kinh tế xã hội 35 3.1.1 Điều kiện tự nhiên 35 3.1.2 Điều kiện kinh tế - xã hội 41 3.1.3 Thực trạng công tác quản lý đất đai 45 3.1.4 Hiện trạng sử dụng đất địa bàn huyện Thanh Trì 46 3.2 Xây dựng đồ trạng sử dụng đất 48 3.2.1 Thu thập tư liệu 49 3.2.2 Tăng cường chất lượng ảnh 49 3.2.3 Xây dựng tệp mẫu vùng nghiên cứu 50 3.2.4 Phân loại ảnh 52 3.2.5 Lập đồ sử dụng đất 56 3.3 Xây dựng đồ đánh giá biến động đất đai 63 3.3.1 Xây dựng đồ biến động đất đai 63 3.3.2 Đánh giá biến động đất đai 66 KẾT LUẬN VÀ KIẾN NGHỊ 69 TÀI LIỆU THAM KHẢO 71 v DANH MỤC BẢNG Số bảng Tên bảng Trang 1.1 Tóm tắt phát triển viễn thám qua kiện 11 1.2 Các hệ thống vệ tinh Landsat 13 1.3 Các thông số kỹ thuật loại cảm 14 1.4 Các hệ vệ tinh SPOT 16 1.5 Các thông số kỹ thuật cảm vệ tinh SPOT 16 1.6 Đặc điểm ảnh vệ tinh VNREDSat-1 18 1.7 Bảng so sánh hai phương pháp giải đoán ảnh viễn thám 23 3.1 Hiện trạng sử dụng đất năm 2014 huyện Thanh Trì 47 3.2 Bảng xây dựng tệp mẫu ảnh năm 2005 50 3.3 Các loại hình sử dụng đất huyện Thanh Trì 50 3.4 Bảng xây dựng tệp mẫu ảnh năm 2010 51 3.5 Bảng xây dựng tệp mẫu ảnh năm 2014 51 3.6 Mẫu khoá giải đoán ảnh vệ tinh 52 3.7 Ma trận sai số phân loại ảnh năm 2014 56 3.8 Thống kê diện tích đất giải đốn năm 2005 57 3.9 Thống kê diện tích đất giải đốn năm 2010 59 3.10 Thống kê diện tích đất giải đoán năm 2014 61 3.11 Gán loại đất thành mã số riêng 63 3.12 Diện tích biến động loại đất giai đoạn 2005 - 2010 66 3.13 So sánh thay đổi diện tích loại đất từ 2005 – 2010 66 3.14 Diện tích biến động loại đất giai đoạn 2010 - 2014 67 3.15 So sánh thay đổi diện tích loại đất từ 2010 – 2014 68 vi DANH MỤC HÌNH Số hình Tên hình Trang 1.1 Thành lập đồ biến động phương pháp so sánh sau phân loại 1.2 Nghiên cứu biến động phương pháp phân loại ảnh đa thời gian 1.3 Nguyên lý thu nhận hình ảnh viễn thám 1.4 Ảnh vệ tinh Landsat miền đông Kazakhstan, độ phân giải 15m 15 1.5 Ảnh vệ tinh SPOT thành phố New York, Mỹ ngày 11/9/2001 17 1.6 Ảnh vệ tinh QuickBird, độ phân giải 1m 17 1.7 Ảnh chụp từ vệ tinh VNREDSat-1 khu vực sông Hồng với độ phân giải 2.5 m 19 1.8 Quá trình thu nhận, xử lý ứng dụng vệ tinh VNREDSat-1 19 1.9 Sơ đồ phần mềm hệ thống GIS 25 1.10 Nguyên lý chồng lắp đồ 27 1.11 Trạm thu ảnh vệ tinh &Trung tâm Quản lý Dữ liệu quốc gia 29 1.12 Nghiên cứu sóng thần 30 3.1 Biểu đồ cấu trạng loại đất năm 2014 48 3.2 Trình tự giải đốn ảnh viễn thám cơng nghệ số 48 3.3 Tư liệu ảnh viễn thám huyện Thanh Trì sau cắt ranh giới 49 3.4 Lấy mẫu cho loại hình sử dụng đất 50 3.5 Ảnh 2005 sau phân loại 54 3.6 Ảnh 2005 sau lọc nhiễu 54 3.7 Ảnh đánh giá độ xác kết phân loại 55 3.8 Biểu đồ biến động đất đai huyện Thanh Trì giai đoạn 2005 - 2010 67 3.9 Biểu đồ biến động đất đai huyện Thanh Trì giai đoạn 2010 - 2014 68 vii MỞ ĐẦU 1.1 Tính cấp thiết đề tài: Lịch sử nhân loại chứng minh tảng cho sống hoạt động sản xuất người bắt nguồn từ đất đai Đất đai tư liệu đặc biệt khơng thay được, thành phần quan trọng môi trường sống, địa bàn phân bố khu dân cư, xây dựng cơng trình văn hóa, xã hội, an ninh, quốc phịng Đất đai có tính chất đặc trưng khiến không giống tư liệu sản xuất khác Đất đai nguồn lực quan trọng hàng đầu cho phát triển kinh tế, kỹ thuật, sở hạ tầng xã hội ln cố định diện tích, vị trí khơng gian vơ hạn thời gian sử dụng quốc gia.Có thể nói đất đai vấn đề xuyến suốt thời đại va phải quản lý chặt chẽ đất lập phương án quy hoạch sử dụng đất nhằm nâng cao hiệu sản xuất xã hội điều tất yếu Để quản lý chặt chẽ quỹ đất ngày hạn chế việc nắm số liệu biến động đất đai thời kỳ thông qua số liệu thống kê, kiểm kê đất đai việc xây dựng đồ biến động sử dụng đất cần thiết Bản đồ biến động sử dụng đất thành lập từ đồ trạng sử dụng đất có song đồ trạng sử dụng đất thành lập năm lần vào năm kiểm kê đất đai nên khó thành lập đồ biến động sử dụng đất cho giai đoạn bất kỳ, thêm vào chất lượng đồ trạng sử dụng đất thành lập kỳ kiểm kê đất đai chưa cao Viễn thám (Remote sensing) xác định phương pháp nghiên cứu đối tượng, tượng thiết bị đặt cách đối tượng khoảng cách mà khơng địi hỏi tiếp xúc trực tiếp với đối tượng Thuật ngữ viễn thám sử dụng Mỹ vào năm 1960 bao hàm lĩnh vực đo ảnh, giải đoán ảnh, địa chất ảnh Ở Việt Nam, viễn thám quan tâm kể từ năm 1980 nước ta tham gia tổ chức vũ trụ quốc tế InterCosmos Đặc điểm quan trọng tư liệu viễn thám có ưu điểm giàu thông tin, chu kỳ thu nhận thơng tin ngắn, xử lý diện rộng phương pháp viễn thám ưu điểm hẳn phương pháp cổ điển khác nghiên cứu diễn biến trình xảy đất thời gian: thay đổi sử dụng đất, diễn biến rừng, trình xói mịn đất… Hệ thống thơng tin địa lý (GIS) tổ chức tổng thể hợp phần: thiết bị, phần mềm, liệu địa lý người điều hành, thiết kế hoạt động cách hiệu nhằm thu nhận, lưu trữ, điều khiển, phân tích, hiển thị toàn dạng liệu địa lý phục vụ giải lớp rộng lớn tốn ứng dụng liên quan đến vị trí địa lý bề mặt trái đất GIS sử dụng công cụ kỹ thuật để giải vấn đề cấp thiết thời đại tình trạng triệt phá rừng, xuống cấp môi trường, vấn đề thị hố, dự báo biến động khí hậu GIS cung cấp cơng cụ mạnh để xây dựng, tổ chức, xử lý quản lý liệu cung cấp thông tin trợ giúp cho chuyên gia GIS nhà quản lý việc định đắn, giải pháp hữu hiệu cho vấn đề Sự kết hợp viễn thám kỹ xử lý số liệu hệ thống thông tin địa lý (GIS) hệ thống định vị tồn cầu (GPS) giúp có cơng cụ hồn chỉnh để tìm hiểu thay đổi sử dụng đất, số lượng vị trí phân bố Đó phương pháp theo dõi thay đổi sử dụng đất ưu việt hẳn so với phương pháp truyền thống khơng gian, thời gian kinh phí, ta theo dõi diễn biến tự nhiên tác động người hàng chục năm trở lại Ngoài việc mang lại hiệu cao, xác phương pháp kết hợp ảnh viễn thám cơng nghệ GIS cịn dễ cập nhật Với ưu điểm đó, việc sử dụng tư liệu ảnh viễn thám để thành lập đồ thay đổi sử dụng đất trở thành phương pháp hiệu nghiên cứu thay đổi tài ngun đất, góp phần tích cực công tác quản lý đất đai, giúp Nhà Nước đưa sách tầm cỡ vĩ mơ hạn chế biến động đất đai Huyện Thanh Trì huyện nằm vị trí cửa ngõ phía Nam Thủ đơ, tuyến giao thông huyết mạch Trong năm gần với chế kinh tế thị trường nói chung phát triển huyện Thanh Trì nói riêng có bước phát triển mạnh mẽ, dẫn đến nhu cầu sử dụng đất cho mục đích khác khơng ngừng thay đổi, nhiên vấn đề đặt đất đai có hạn Xuất phát từ nhu cầu thực tế tính cấp bách việc xác định biến động đất đai địa bàn huyện Thanh Trì, đồng ý Khoa Quản lý đất đai, Khoa Sau đại học - trường Đại học Nông Nghiệp Hà Nội, tiến hành nghiên cứu đề tài: “Tích hợp Viễn thám Hệ thống thông tin địa lý đánh giá biến động đất đai giai đoạn 2005-2014 huyện Thanh Trì – Thành phố Hà Nội” 1.2 Mục đích đề tài - Xây dựng đồ biến động đất đai đánh giá biến động đất đai qua hai giai đoạn 2005 – 2010, 2010 - 2014 địa bàn huyện Thanh Trì – Thành phố Hà Nội công nghệ viễn thám kết hợp với GIS 1.3 Yêu cầu đề tài - Bản đồ trạng sử dụng đất năm 2005, 2010, 2014 kết hợp với phương pháp điều tra thực địa - Bản đồ xây dựng nằm hệ tọa độ VN-2000, đảm bảo độ xác theo quy phạm thành lập đồ Bộ Tài nguyên Môi trường KẾT LUẬN VÀ KIẾN NGHỊ Kết luận Thanh Trì huyện nằm ven nội thành thành phố Hà Nội, có diện tích tự nhiên 6.292,71 ha, cấu kinh tế huyện có chuyển dịch tích cực theo hướng tăng tỷ trọng Công nghiệp Thương mại - Dịch vụ, giảm dần tỷ trọng Nơng nghiệp Huyện Thanh trì có nguồn nhân lực dồi ,trình độ dân trí cao ,có q trình thị hóa nhanh nên kinh tế phát triển mạnh toàn diện tương lai khu vực thị thành phố Hà Nội Tích hợp tư liệu ảnh viễn thám hệ thống thông tin địa lý GIS xây dựng tệp liệu mẫu gồm loại hình sử dụng đất bao gồm: sông, đất mặt nước chuyên dùng, đất nông nghiệp, đất cơng trình đất giao thơng Bằng việc sử dụng kết hợp viễn thám GIS giải đoán ảnh viễn thám năm 2005, 2010 2014 Trên sở thành lập sơ đồ sử dụng đất huyện năm 2005, 2010 2014 đạt độ xác cao Bằng phương pháp số kỹ thuật GIS xây dựng đồ biến động đất giai đoạn 2005-2010 giai đoạn 2010 - 2014 huyện Thanh Trì qua đánh giá biến động đất đai hai giai đoạn 2005 – 2010 2010 – 2014 Trong giai đoạn 2005 – 2010, diện tích đất xây dựng tăng 339.55ha có 78.94 đất mặt nước, 339.11ha đất nơng chuyển sang mục đích xây dựng Nguyên nhân biến động trình thị hố nhanh chóng địa bàn phường nên nhiều cơng trình xây dựng để phục vụ nhu cầu nhân dân; diện tích đất nơng nghiệp giảm 236.07ha 339.11ha chuyển sang đất xây dựng, 253.14ha đất chuyển sang mặt nước chuyên dùng, 4.03 sang đất giao thông 8.11 sang đất sông suối Nguyên nhân biến động mở đường, giãn dân nên nhiều diện tích đất hàng năm bị chuyển mục đích Như vậy, thấy giai đoạn 2005 – 2010 diện tích đất nơng nghiệp có xu giảm dần, diện tích đất xây dựng tăng đáng kể, diện tích đất mặt nước chuyên dùng tăng dần Sự biến động tăng giảm loại đất chưa nhiều 69 Trong giai đoạn 2010 – 2014, diện tích đất xây dựng tăng 1,886.73ha có 394.46 đất mặt nước, 1180,35ha đất nơng chuyển sang mục đích xây dựng Nguyên nhân biến động q trình thị hố nhanh chóng địa bàn phường nên nhiều cơng trình xây dựng để phục vụ nhu cầu nhân dân, diện tích đất nơng nghiệp giảm 1,514.91ha 1180,35ha chuyển sang đất xây dựng, 394.46 đất chuyển sang mặt nước chuyên dùng, 11.17 sang đất giao thơng 5.29 sang đất sơng suối Như vậy, thấy giai đoạn 2010 - 2014 diện tích đất nơng nghiệp có xu giảm, diện tích đất xây dựng tăng đáng kể, diện tích đất mặt nước chuyên dùng tăng dần Sự biến động tăng giảm loại đất chưa nhiều đáng kể đặc biệt đất nông nghiệp đất xây dựng Kiến nghị - Phương pháp kết hợp công nghệ viễn thám với GIS nghiên cứu biến động đất đai có độ tin cậy cao Tuy nhiên để xác định xác biến động nên chọn ảnh thu nhận mùa năm, điều khó thực nước thuộc vùng khí hậu nhiệt đới Việt Nam, mây sương mù phổ biến ảnh cần cố độ phân giải cao không việc giải đốn ảnh gặp khó khăn gây ảnh hưởng tới kết giải đoán nhiên ảnh viễn thám có độ phân giải cao có giá thành cao - Trong giai đoạn 2005 – 2014 tốc độ đô thị hố huyện Thanh Trì diễn mạnh nhanh dẫn tới biến động đất đai rõ nét hai loại đất cơng trình đất nơng nghiệp Chính quyền địa phương cần có sách cụ thể hỗ trợ việc làm, vốn để người dân ổn định sản xuất nâng cao đời sống 70 TÀI LIỆU THAM KHẢO Phần tiếng việt Phạm Văn Cự, Ứng dụng hệ thông tin địa lý viễn thám quản lý môi trường tài nguyên Việt Nam thực trạng, thuận lợi thách thức, Hà Nội 2008 Lê Thị Giang, Nguyễn Thị Châu Long (2010), Ứng dụng Hệ thống thông tin địa lý (GIS) quản lý thông tin bất động sản thị trấn Trâu Quỳ - Gia Lâm, Hà Nội, Tạp chí Khoa học Phát triển tập 8, số 1: 130 – 139 Lê Thị Giang, Đào Châu Thu (2002) Tìm hiểu thay đổi sử dụng đất huyện yên Châu tỉnh Sơn La công nghệ GIS viễn thám, Tạp chí Khoa học đất số 17 Bảo Huy, GIS viễn thám quản lý tài nguyên rừng môi trường, NXB Tổng hợp thành phố Hồ Chí Minh Nguyễn Thị Thu Hiền, Phạm Vọng Thành, Nguyễn Khắc Thời (2014) Đánh giá biến động sử dụng đất/lớp phủ huyên Tiên Yên tỉnh Quảng Nình giai đoạn 2000 – 2010, Tạp chí Khoa học Phát triển, tập (1), 43 – 51 Nguyễn Thị Thu Hiền, Nguyễn Khắc Thời (2011) Sử dụng ảnh vệ tinh công nghệ GIS thành lập đồ biến động đất đô thị thành phố Buôn Ma Thuột giai đoạn 2000 – 2009, Tạp chí Khoa học phát triển, số Đặng Thị Mỹ Lan (2008), Tài liệu hướng dẫn sử dụng ArcGIS bản, Viện nghiên cứu Địa Phịng Tài ngun Mơi trường huyện Thanh Trì, Báo cáo thống kê đất đai năm từ 2010 đến 2014 10 Trần Thị Băng Tâm (2006) Giáo trình Hệ thống thơng tin địa lý, trường Đại học Nông nghiệp Hà Nội Nguyễn Khắc Thời, Trần Quốc Vinh, Lê Thị Giang, Nguyễn Thị Thu Hiền (2008) Ứng dụng kỹ thuật Viễn thám công nghệ GIS để xác định biến động đất đai tiến trình thị hóa huyện Gia Lâm quận Long Biên thành phố Hà Nội, Tạp chí khoa học Đất, năm 2008 11 Nguyễn Khắc Thời cộng (2010) Giáo trình viễn thám, trường Đại học Nơng nghiệp Hà Nội 12 Phạm Vọng Thành (2013) Giáo trình Kỹ thuật viễn thám, trường Đại học Nông nghiệp Hà Nội 13 Trung tâm Viễn thám Quốc gia - Ứng dụng công nghệ viễn thám giám sát tài nguyên thiên nhiên môi trường Việt Nam 14 Trung tâm Viễn thám Quốc gia - Ứng dụng công nghệ viễn thám giám sát tài nguyên thiên nhiên môi trường Việt Nam 15 Trung tâm Viễn thám Quốc gia - Ứng dụng công nghệ viễn thám để quản lý dải ven biển 16 Trung tâm Viễn thám Quốc gia - Giới thiệu công nghệ viễn thám hệ thống thông tin địa lý ứng dụng quản lý tổng hợp đới bờ Việt Nam 17 Trung tâm Viễn thám Quốc gia - Các cơng trình thành lập chỉnh đồ địa hình 71 18 Uỷ ban nhân dân huyện Thanh Trì, Quy hoạch tổng thể phát triển kinh tế - xã hội huyện Thanh Trì đến năm 2020 19 Uỷ ban nhân dân huyện Thanh Trì, Báo cáo sử dụng đất huyện Thanh Trì năm 2010 đến 2014 20 Uỷ ban nhân dân huyện Thanh Trì, Báo cáo tình hình kinh tế - xã hội huyện Thanh Trì năm 2010 đến 2014 21 Trần Quốc Vinh (2004) Ứng dụng kỹ thuật Viễn thám GIS tìm hiểu thay đổi sử dụng đất nông – lâm nghiệp huyện Con Cuông – tỉnh Nghệ An, Tạp chí Khoa học đất số 21 22 Trần Quốc Vinh, Lê Thị Giang (2006) Sử dụng phương pháp giải đoán ảnh mắt thường phương pháp giải đoán số việc xây dựng đồ sử dụng đất huyện Thanh Trì, thành phố Hà Nội, Tạp chí Khoa học đất số 25 23 Trần Quốc Vinh, Đặng Hùng Võ, Đào Châu Thu (2011) Ứng dụng viễn thám hệ thống thông tin địa lý đánh giá xói mịn đất đồi gị huyện Tam Nơng tỉnh Phú Thọ, Tạp chí Khoa học Phát triển 2011, Tập 9, số 5, 823 – 833 Ấn phẩm điện tử 24 Viện hàn lâm khoa học công nghệ việt nam, Vệ tinh VNREDSat – 1: sau tháng hoạt động quỹ đạo, từ http://www.vast.ac.vn/tin-tuc-su-kien/tin-khoahoc/trong-nuoc/1668-ve-tinh-vnredsat-1-sau-3-thang-hoat-dong-tren-quy-dao 25 Phạm Khánh Chi, Giới thiệu chung công nghệ viễn thám, từ http://www.geoviet.vn/goc-ky-thuat/vn/400/464/293/0/gioi-thieu-chung-vecong-nghe-vien-tham.aspx 72 PHỤ LỤC Số liệu thực địa GPS cầm tay phục vụ đánh giá độ xác giải đốn ảnh năm 2014 STT Tọa độ X Tọa độ Y Tham chiếu 591281.225 2314914.838 591196.558 2315264.089 591256.875 2314968.76 590883.15 2315544.23 590661.561 2315702.981 626421.985 2225781.271 590664.868 2225873.875 590152.237 2227421.691 589983.565 2227785.493 10 589361.793 2229223.892 11 589275.803 2229680.300 12 589173.277 2229508.320 13 589282.418 2229374.343 14 589206.35 2229691.844 15 588564.734 2230710.492 16 589226.194 2225106.372 17 589292.34 2223849.599 18 589398.173 2218718.102 19 588773.094 2218288.153 20 588554.812 2220186.542 21 587794.133 2219604.458 22 588144.707 2218618.883 23 585703.921 2219366.332 ii 24 585667.541 2220517.493 25 585660.926 2220777.843 26 585313.66 2221406.494 27 585445.852 2230324.823 28 584873.789 2230239.355 29 584758.034 2231168.822 30 584338.007 2231051.304 31 584142.876 2232066.240 32 584646.115 2230912.418 33 584663.048 2231467.962 34 584794.282 2232141.025 35 584379.414 2234460.138 36 583733.83 2235524.791 37 583672.446 2234207.282 38 583462.896 2234819.458 39 583390.929 2233941.119 40 583632.229 2231931.585 41 583754.996 2234020.986 42 585233.094 2234327.056 43 585259.552 2236828.992 44 585259.552 2236629.370 45 584773.38 2236616.062 46 583956.08 2236948.766 47 584006.88 2238066.653 48 583812.146 2237880.338 49 583727.48 2238732.061 50 583725.363 2238918.376 iii ii iii iv v vi vii viii ix