ỨNG DỤNG VIỄN THÁM VÀ GIS TRONG THEO DÕI VÀ TÍNH TOÁN BIẾN ĐỘNG ĐƯỜNG BỜ KHU VỰC PHAN THIẾT

Tạp chí Khoa học Cơng nghệ biển T11 (2011) Số Tr - 13 ỨNG DỤNG VIỄN THÁM VÀ GIS TRONG THEO DÕI VÀ TÍNH TỐN BIẾN ðỘNG ðƯỜNG BỜ KHU VỰC PHAN THIẾT PHẠM THỊ PHƯƠNG THẢO Viện Hải dương học HỒ ðÌNH DUẨN Viện Tài nguyên, Môi trường Phát triển bền vững TP Huế ðẶNG VĂN TỎ Trường ðại học Khoa học Tự nhiên Tp HCM Tóm tắt: Bằng kỹ thuật đại từ công nghệ viễn thám GIS, việc theo dõi tính tốn tốc độ thay đổi đường bờ thực cách nhanh chóng hiệu Với nguồn ảnh Landsat ña thời gian, tranh sơ biến ñộng ñường bờ khu vực Phan Thiết giai ñoạn 1973 - 2004 ñược vẽ lại Bên cạnh đó, với hỗ trợ đắc lực từ phần mở rộng DSAS (Digital Shoreline Analysis System) công cụ GIS, tốc ñộ thay ñổi ñường bờ khu vực Hàm Tiến tính tốn giai đoạn trước sau có cơng trình ven bờ I GIỚI THIỆU Ngày nay, với tích hợp khéo léo công nghệ viễn thám GIS, việc theo dõi tính tốn biến động đường bờ thực nhanh chóng hiệu Ngồi chuyến thực ñịa dụng cụ thiết bị truyền thống, nay, thực đo vẽ ñường bờ thiết bị GPS, ảnh kỹ thuật số, hệ thống máy camera video [4, 5, 6, 12] Bên cạnh đó, đường bờ cịn thu phạm vi rộng lớn từ ảnh máy bay ảnh vệ tinh Từ liệu đó, việc đánh giá phân tích diễn biến đường bờ thực trực tiếp đưa vào cơng cụ GIS để tính tốn tốc độ thay đổi đường bờ Tuy cơng nghệ viễn thám GIS Việt Nam phát triển muộn nhiều so với giới nước ta ñã ñạt ñược số kết nghiên cứu định Về ứng dụng viễn thám phân tích ñánh giá biến ñộng ñường bờ, số nghiên cứu ñã ñược thực số khu vực, tập trung chủ yếu vùng đồng sơng Hồng sông Cửu Long, số khu vực cửa sông Lộc An, Cửa ðại, Thuận An Cụ thể như, với liệu ảnh Landsat, SPOT ảnh máy bay, nhóm tác giả Phạm Viết Cường cộng ñã nghiên cứu thay ñổi vùng ven bờ cửa sơng Hải Phịng từ năm 1975 đến năm 1988 [15] Biến động đường bờ cửa sơng chính: cửa Trà Lý, cửa Ba Lạt, cửa Văn Úc cửa ðáy thuộc dải ven biển đồng sơng Hồng từ năm 1926 đến 1995 nhóm tác giả Nguyễn Tứ Dần Nguyễn Thế Tiệp thực liệu ảnh SPOT, Landsat, đồ địa hình, ảnh máy bay [8] Các biến ñổi lịch sử khu vực cửa sơng Lộc An từ ảnh viễn thám đa thời gian từ năm 1953 ñến năm 2002 ñược thực nhóm tác giả Phạm Bách Việt từ ảnh Landsat ảnh máy bay [9] Với mục đích quản lý vùng ven bờ, tác giả Nguyễn Hạnh Quyên sử dụng kỹ thuật viễn thám khoanh vùng ñất ngập nước khu vực vịnh Hạ Long, ñồng thời cho kết thay ñổi rừng ngập mặn thay ñổi bờ biển khoảng thời gian khác từ 1988 ñến 2002 ảnh Landsat TM ETM+ [7] Việc phát thay ñổi bờ biển vùng ven bờ cửa sông Cửu Long thời gian 1989 đến 2004 thực nhóm tác giả Trần Thị Vân Trịnh Thị Bình [14] Về ứng dụng GIS tính tốn định lượng tốc độ thay ñổi ñường bờ, số kết ñã thu thập ñược Với liệu ñường bờ ñược số hóa từ liệu khảo sát đo đạc từ năm 1905 ñến 1992 khu vực Hải Hậu, Nam ðịnh, tác giả ðặng Văn Tỏ Phạm Thị Phương Thảo trình bày số kết tính tốn tốc độ thay đổi đường bờ từ phương pháp thống kê có sẵn phần mở rộng DSAS cơng cụ GIS [3] Tương tự, đường bờ lịch sử khu vực Phan Thiết, Bình Thuận từ năm 1973 ñến 2002 ñược rút trích từ ảnh Landsat tốc ñộ thay ñổi ñường bờ khu vực ñược tính tốn nhanh chóng phần mở rộng DSAS nhóm tác giả Phạm Thị Phương Thảo, Hồ ðinh Duẩn ðặng Văn Tỏ thực [10] Với phát triển cực nhanh dịch vụ du lịch từ năm 1995, Phan Thiết ñang phải ñối mặt với tượng xói lở nghiêm trọng tác động cơng trình ven bờ [2] Do đó, Phan Thiết chọn làm khu vực nghiên cứu (hình 1) Bài báo trình bày hai vấn đề: phân tích diễn biến ñường bờ theo thời gian giai ñoạn 1973 - 2004 kết tính tốn tốc độ thay đổi đường bờ hàng năm DSAS khu vực Phan Thiết trước sau có cơng trình ven bờ Hình 1: Khu vực nghiên cứu II ðIỀU KIỆN TỰ NHIÊN KHU VỰC NGHIÊN CỨU ðặc điểm khí tượng Phan Thiết chịu ảnh hưởng hai chế độ gió mùa: gió mùa ðơng Bắc vào mùa đơng gió mùa Tây Nam vào mùa hè Gió mùa ðơng Bắc bắt ñầu thổi từ tháng 11 kéo dài tới tháng năm sau Vận tốc gió trung bình mùa vào khoảng - 10 m/s Gió mùa Tây Nam bắt ñầu từ tháng ñến tháng Vận tốc gió trung bình khoảng - m/s Tháng - giai ñoạn chuyển tiếp từ gió mùa mùa đơng sang gió mùa mùa hè tháng - 10 giai ñoạn chuyển tiếp ngược lại [1] ðặc ñiểm thủy văn Do chịu tác động hệ thống gió mùa ðơng Bắc Tây Nam, Phan Thiết có hai chế độ sóng gió chủ yếu Về mùa đơng, sóng có hướng chiếm ưu hướng ðơng Bắc, có chuyển ðơng, độ cao sóng trung bình dao động từ 1,2 - 2,0 m Trong trường hợp gió có hướng Tây Bắc chiếm ưu sóng có hướng Nam Tây Nam, sóng với hướng có độ cao trung bình 0,7 - 1,5 m Tuy nhiên, chế độ sóng vào thời kỳ thường khơng ổn định bị ảnh hưởng bão áp thấp nhiệt đới, sóng ngồi khơi ven bờ thường quan sát thấy khơng đồng [1, 11] Thủy triều vùng ven biển Bình Thuận phức tạp nằm khu vực chuyển tiếp chế ñộ nhật triều khơng phía Bắc (tiêu biểu Qui Nhơn) bán nhật triều khơng phía Nam (Vũng Tàu) [1] Dao ñộng triều khoảng 2m Chế ñộ thủy triều gây dòng chảy mạnh ven bờ lên đến 50 - 70 cm/s [1] III PHƯƠNG PHÁP NGHIÊN CỨU Viễn thám ðể ñạt hiệu việc rút trích đường mực nước cho khu vực bờ biển vịnh Phan Thiết, phương pháp kết hợp phân tích giá trị ngưỡng ảnh tỷ số ñược áp dụng cho ảnh Landsat Ưu ñiểm phương pháp kết hợp loại bỏ nhiễu vùng nghiên cứu có độ phủ thực vật cao nhiễu vùng sóng vỡ Cơng thức cụ thể sau: - ðối với Landsat MSS: (B3+B4)/B1 - ðối với Landsat TM, ETM+: (B5+B7)/B2 có kết hợp B7 Do khơng có số liệu đo đạc thực tế ñồ ñịa hình ñáy tỷ lệ cao khu vực vịnh Phan Thiết ñể phục vụ cho việc hiệu chỉnh triều, độ phân giải khơng gian khơng cao (30 m) ảnh Landsat nên việc hiệu chỉnh triều bỏ qua Vì thế, kết đường mực nước rút trích xem đường bờ GIS Sau hồn thành tập liệu đường bờ, việc tính tốn tốc độ thay đổi đường bờ thực DSAS, phần mở rộng ArcView Thieler cộng viết ngôn ngữ Avenue phần mềm ArcView [13] Cơng việc tính tốn phân tích ñường bờ ñược tiến hành sau: Xác ñịnh đường chuẩn (baseline) đường bờ tính tốn (shoreline) Tạo tuyến cắt ngang vng góc bờ (transect) Tính tốn tốc độ thay đổi đường bờ Trong DSAS có nhiều phương pháp tính tốn thống kê khác nhau, phương pháp tính tốc độ thơng qua điểm ñầu – ñiểm cuối (End-Point Rate), phương pháp tính tốc ñộ trung bình AOR (Average Of Rate), phương pháp hồi quy tuyến tính (Linear Regression), phương pháp gập gãy (Jack-Knife Rate) Tùy theo chất lượng số lượng ñường bờ ta chọn phương pháp thống kê phù hợp Dựa liệu thu thập ñược, phương pháp hồi quy tuyến tính chọn để phân tích kết IV DỮ LIỆU SỬ DỤNG Nguồn liệu ảnh Landsat thu thập từ trang web Cơ quan ðịa chất Hoa Kỳ (USGS) (http://edcsns17.cr.usgs.gov/EarthExplorer/) Các ảnh ñã ñược nắn chỉnh theo hệ qui chiếu WGS-84 UTM, áp dụng cho vùng 49 Danh sách ảnh khu vực vịnh Phan Thiết ñược thu thập ñược ghi bảng Bảng 1: Các ảnh Landsat khu vực vịnh Phan Thiết STT Loại ảnh ðộ phân giải (m) Ngày chụp Năm qui ước Landsat MSS 57 01/01/1973 1973 Landsat MSS 57 31/01/1976 1976 Landsat TM 28.5 30/12/1990 1991 Landsat ETM+ 30 13/11/1999 2000 Landsat ETM+ 30 01/12/2000 2001 Landsat ETM+ 30 05/01/2002 2002 Landsat ETM+ 30 24/01/2003 2003 Landsat TM 30 20/12/2004 2005 Các ảnh Landsat ñược chụp khoảng thời gian từ tháng 11 ñến cuối tháng năm sau, trùng với thời kỳ gió mùa ðơng Bắc Do đó, việc tính tốn phân tích thay ñổi ñường bờ ñược thời kỳ gió mùa ðơng Bắc biến đổi qua năm V KẾT QUẢ VÀ THẢO LUẬN Nhằm tiện việc phân tích, khu vực nghiên cứu chia thành khu vực nhỏ từ trái sang phải ñược ñánh dấu từ đến hình Các khu vực 1, khu vực tập trung dân cư, số nơi có bờ vách ñá Khu vực ñoạn bờ ðồi Dương, nằm hai cửa sông: Cái Cà Ty ðây đoạn bờ có dịch vụ du lịch phát triển Khu vực ñoạn bờ Hàm Tiến, Mũi Né nơi tập trung hầu hết khách sạn khu nghỉ dưỡng cao cấp phục vụ du lịch ngồi nước Hình 2: Khu vực nghiên cứu Do khu vực 1, có đường bờ tương đối ổn định xói lở trình bày hình 3, vậy, khu vực khu vực ñược khảo sát nghiên cứu chi tiết Khu vực hay khu vực ðồi Dương gần cửa sơng nên đường bờ có biến ñộng ñáng kể qua nhiều năm từ 1973 ñến 2005 hình Ngồi tác động q trình ñộng lực ven bờ, ñường bờ khu vực chịu ảnh hưởng trình vận chuyển phù sa đổ từ sơng Sự diện kè cảng cửa sông Cái cửa sông Cà Ty góp phần tác động đến q trình vận chuyển trầm tích Tất yếu tố làm cho q trình biến đổi đường bờ nơi phức tạp Ở cửa sông Cà Ty, vào năm 1973, có xuất doi cát phía bờ Nam Trong năm 1976, vật liệu phù sa không bị đẩy qua khu vực phía Nam mà tích tụ lại bờ bắc, hình thành doi cát phía Bắc cửa sông ðến năm 1991, khu vực cửa sơng bắt đầu có biến đổi phức tạp với hai doi cát hai bên cửa sông, ñặc biệt hình thành doi cát dài bờ Bắc chắn ngang cửa sông Vào năm 90 kỷ trước, kè cảng phía Nam cửa sơng xây dựng, góp phần ổn định đoạn bờ Vật liệu tích tụ bên phía Nam kè Năm 2000, kè phía Bắc cửa sơng Cà Ty xây dựng, góp phần ngăn chặn q trình hình thành doi cát phía Bắc cửa sơng Hiện tượng xói lở bắt đầu xuất phía bên trái kè phía Nam cửa sơng Cà Ty năm 2002 Ngun nhân có q trình đổ vật liệu lấn biển công ty tư nhân, từ bắt đầu xuất cơng trình đê biển để bảo vệ đoạn bờ Vì thế, từ sau năm 2002, đường bờ khu vực cửa sơng Cà Ty tương ñối ổn ñịnh Khu vực Khu vực Khu vực Hình 3: Kết đường bờ khu vực 1, Ở cửa sông Cái, vào năm 1976 nơi bắt đầu hình thành doi cát nhỏ phía Bắc cửa sơng Doi cát tiếp tục phát triển kéo dài phía Nam cho ñến năm 1991, 2000, 2001 Nguyên nhân dịng chảy từ sơng đổ kết hợp với trình ven bờ tạo nên lắng ñọng phù sa khu vực gần cửa sông Năm 2002, khu vực cửa sơng mở lại, làm cho phía Nam cửa sơng hình thành doi cát xoắn dọc bờ Vào ñầu năm 2003, kè cảng phía Nam cửa sơng Cái chưa xây xong, nơi cịn tồn vài doi cát nhỏ Năm 2005, doi cát dần sau có diện kè phía Nam Khu vực cửa sơng cố định nhằm ổn định dần q trình biến đổi 1973 1973-1976 1976-1991 1991-2000 2000-2001 2001-2002 2002-2003 2003-2005 Hình 4: Sự thay ñổi ñường bờ khu vực (ðồi Dương) từ năm 1973 ñến 2005 Khu vực hay khu vực Hàm Tiến có đường bờ biến động tương đối ổn ñịnh so với trường hợp khu vực ðồi Dương Do khu vực nằm khuất sau Mũi Né nên khơng bị ảnh hưởng lớn gió ðơng Bắc vào cuối năm Về sóng biển, khu vực chủ yếu có trường sóng khúc xạ từ trường sóng ðơng Bắc ngồi khơi truyền vào Nhằm chống xói lở, khu vực xây dựng cơng trình bảo vệ bờ biển từ cuối thập niên 90 kỷ trước Do vậy, phân tích biến động ñường bờ khu vực này, hai giai ñoạn xác định có mặt cơng trình trình bày: giai đoạn trước năm 2000 giai ñoạn sau năm 2000 Từ giai ñoạn 1973 - 2000, dòng vật chất di chuyển từ Bắc xuống Nam (từ phải qua trái hình 5), khu vực phía Bắc đường bờ có xu hướng lùi vào đất liền khu vực phía Nam ngược lại lấn biển Hình 5: Sự thay đổi đường bờ khu vực (Hàm Tiến) từ năm 1973 ñến 2000 Từ giai đoạn 2000 - 2005, cơng trình bảo vệ ñường bờ khu khách sạn nghỉ dưỡng xuất hiện, quy luật biến ñộng ñường bờ nơi bắt đầu thay đổi (hình 6) Cụ thể, khu vực gần mũi đá Ơng ðịa bắt đầu bị xói lở nghiêm trọng giai đoạn 2003 2005 Giai ñoạn 2003 - 2005 Hình 6: Sự thay ñổi ñường bờ khu vực (Hàm Tiến) từ năm 2000 ñến 2005 Một số kết phân tích cho thấy khu vực ðồi Dương có thay đổi ñường bờ gần cửa sông phức tạp theo thời gian Việc tính tốn định lượng tốc độ xói lở bồi tụ không hiệu phần doi cát ðồi Dương tính tốn phần đường bờ Từ hạn chế đó, việc tính tốn tốc độ thay đổi đường bờ thực cho khu vực hay khu vực Hàm Tiến, bỏ qua khu vực ðồi Dương Sau ñây số kết tính tốn định lượng thay đổi đường bờ từ DSAS cho khu vực Mũi Né (hình 7, 8) Hình trình bày đường bờ từ 1973 ñến 2005 transect sử dụng DSAS Khu vực tính tốn dài 11.200 m có tổng cộng 112 transect Các transect ñược xây dựng dọc ñường bờ khu vực Hàm Tiến ñược ñánh số theo thứ tự tăng dần từ trái qua phải Khoảng cách transect 100 m ñộ dài transect 300 m Khu vực tính tốn riêng lẻ hai giai ñoạn: giai ñoạn 1973 - 2000 giai đoạn 2000 2005 Hình 7: Xây dựng transect khu vực (Hàm Tiến) từ năm 1973 ñến 2005 Hình trình bày biểu đồ biểu thị tốc ñộ thay ñổi ñường bờ khu vực Hàm Tiến tương ứng với giai đoạn 1973 - 2000 phía bên trái giai đoạn 2000 - 2005 phía bên phải Trên biểu ñồ, giá trị dương biểu thị cho tốc ñộ bồi tụ hàng năm giá trị âm biểu thị cho tốc độ xói lở hàng năm Kết tính tốn cho thấy giai đoạn 1973 2000, ñoạn bờ từ transect ñến transect 36 có xu hướng bồi hàng năm với tốc độ trung bình +2 m Từ transect 36 đến transect 112 có tốc độ xói lở trung bình hàng năm gần -2 m Tốc độ thay đổi đường bờ trung bình cho khu vực Hàm Tiến ñạt giá trị 0m/năm giai đoạn từ 1973 đến 2000 Nhìn chung, đường bờ có xu hướng cân Trong giai đoạn 2000 - 2005, đường bờ khu vực Hàm Tiến có xu hướng xói lở cao với tốc độ trung bình gần -4 m/năm 10 3 0 Tốc ñộ thay ñổi ñường bờ (m/năm) Tốc ñộ thay ñổi ñường bờ (m/năm) -3 -6 -9 -12 -3 -6 -9 -12 -15 -15 11 21 31 41 51 61 71 81 91 101 111 Số transect 11 21 31 41 51 61 71 81 91 101 111 Số transect Giai ñoạn 1973 - 2000 Giai ñoạn 2000 - 2005 Hình 8: Tốc độ thay đổi đường bờ khu vực Hàm Tiến giai ñoạn 1973 - 2000 giai ñoạn 2000 - 2005 VI KẾT LUẬN Với trợ giúp đắc lực cơng nghệ viễn thám GIS, việc theo dõi tính tốn thay đổi ñường bờ ñược thực nhanh chóng dễ dàng Kết phân tích ảnh cho thấy hai khu vực có nhiều biến đổi ðồi Dương Hàm Tiến ðối với ðồi Dương, khu vực có nhiều biến động hai cửa sơng Tuy nhiên, thay ñổi dần ổn ñịnh sau kè cảng ñược xây dựng ðối với Hàm Tiến, quy luật thay ñổi ñường bờ bắt ñầu từ cuối thập niên 90 kỷ trước phát triển dịch vụ du lịch bắt ñầu nở rộ Việc xây dựng đê biển khơng đồng khách sạn khu nghỉ dưỡng ñã làm ñường bờ nơi ñây diễn biến phức tạp Từ DSAS, tính tốn định lượng cho thấy, tốc ñộ thay ñổi ñường bờ hàng năm giai đoạn 1973-2000 xấp xỉ m, cịn giai ñoạn 2000 - 2005 gần -4 m Việc rút trích tự động bán tự động thơng tin từ ảnh viễn thám cho phép ta tiết kiệm thời gian nguồn nhân lực nhiều so với việc ño ñạc thực ñịa truyền thống trường Tuy nhiên, số thơng tin cần phải hiệu chỉnh thêm điểm khống chế mặt đất, thơng tin dao động triều, địa hình đáy… để có ñược nguồn liệu ñường bờ hữu ích xác Bên cạnh đó, cần phải có nguồn ảnh có độ phân giải cao để thấy biến động nhỏ cục Từ làm sở cho việc tính tốn tốc độ thay đổi đường bờ mơ hình số trị cơng nghệ GIS Các kết ñạt ñược phần cho thấy khả áp dụng công nghệ viễn thám GIS ñể ñánh giá biến ñộng ñường bờ cho khu vực khác Việt Nam 11 TÀI LIỆU THAM KHẢO Bùi Hồng Long cộng sự, 2000 ðiều tra điều kiện tự nhiên, mơi trường, nguồn lợi vùng ven bờ vịnh Phan Thiết xây dựng ñịnh hướng phát triển bền vững kinh tế - xã hội ñịa phương Viện Hải dương học Nha Trang Dang Vang To, 2006 Beach Erosion in Doi Duong - Phan Thiet Tourist Resort and Its Proposed Measure, Vietnam-Japan Estuary Workshop, Vietnam: 151-156 Dang Van To and Pham Thi Phuong Thao, 2008 A Shoreline Analysis using DSAS in Nam Dinh Coastal Area, International Journal of Geoinformatics, 4(1): 37-42 ðặng Văn Tỏ Phạm Thị Phương Thả, 2008 Xác ñịnh ñường mực nước từ ảnh số, Tuyển tập Báo cáo Khoa học Hội nghị Khoa học “Biển ðông-2007”, Viện Hải dương học Nha Trang: 591-600 Holman, R A., Sallenger, A H., Lippmann, T C and Haines, J W., 1993 The Application of Video Image Processing to The Study of Nearshore Processes, Oceanography, 6(3), pp 78-85 Li, R., 1997 Mobile Mapping: An Emerging Technology for Spatial Data Acquisition, Photogrammetric Engineering and Remote Sensing, 63(9): 1165-1169 Nguyen Hanh Quyen, Tran Minh Y, Le Thi Thu Hien, 2004 Using Remote Sensing Techniques for Coastal Zone Management in HaLong Bay (Vietnam), Proceedings of International Symposium on Geoinformatics for Spatial Infrastructure Development in Earth and Allied Sciences: 151-156 Nguyễn Tứ Dần, Nguyễn Thế Tiệp, 2003 Xu biến ñộng cửa sơng dải ven biển đồng sơng Hồng, Tạp chí Khoa học Cơng nghệ biển, T3, Số 1: 25-35 Pham Bach Viet, Truong Ngoc Tuong, Nguyen Thanh Minh, 2004 Using TimeSeries Remotely Sensed Data to Trace Historical Changes of LocAn River Mouth Area, Proceedings of International Symposium on Geoinformatics for Spatial Infrastructure Development in Earth and Allied Sciences: 172-177 10 Phạm Thị Phương Thảo, Hồ ðình Duẩn ðặng Văn Tỏ, 2009 Biến động ñường bờ khu vực Mũi Né, Kỷ yếu Kỷ niệm 25 năm ngày Thành lập Viện Cơ học Tin học Ứng dụng, Tp HCM: 334-337 11 Phan Văn Hoặc, 1986 Các đặc trưng khí tượng-hải-thủy văn vùng biển Phan Thiết 12 12 Tanaka, H., 2006 Monitoring of Short-term Morphology Change at A River Mouth, Vietnam-Japan Estuary Workshop, Vietnam: 1-6 13 Thieler, E R., Martine, D., and Ergul, A., 2003 The Digital Shoreline Analysis System, Version Shoreline Change Measurement Software Extension for ArcView, USGS Open-File Report 03-076 14 Tran Thi Van, Trinh Thi Binh, 2008 Shoreline Change Detection to Serve Sustainable Management of Coastal Zone in Cuu Long Estuary, International Symposium on Geoinformatics for Spatial Infrastructure Development in Earth and Allied Sciences: 351-356 15 http://www.aars-acrs.org/acrs/proceeding/ACRS1989/Papers/OCN89-1.htm 16 Pham Viet Cuong, Nguyen Hong Chau, Tran Minh Hien, 1989 Application of Remote Sensing Imagery for Investigation in the Haiphong Estuarine and Coastal Zone, Proceeding of Asian Association on Remote Sensing APPLICATION OF REMOTE SENSING AND GIS FOR CALCULATING THE COASTLINE CHANGES IN PHAN THIET PHAM THI PHUONG THAO, HO DINH DUAN, DANG VAN TO Summary: By integrating remote sensing with GIS, the monitoring and calculation of shoreline changes are carried out quickly and effectively With multi-temporal Landsat imageries, preliminary view is outlined about shoreline change in Phan Thiet from 1973 to 2004 In addition, by the extension of ArcView, DSAS, shoreline change rates are calculated in Ham Tien area before and after some coastal constructions were built Ngày nhận bài: 18 - 10 - 2010 Người nhận xét: TS Nguyễn Bá Xuân 13