Ứng dụng gis trong quản lý sản xuất lúa và cảnh báo sớm tình hình nhiễm rầy nâu trên lúa

Đại Học Quốc Gia Tp.Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA TRƯƠNG CHÍ QUANG ỨNG DỤNG GIS TRONG QUẢN LÝ SẢN XUẤT LÚA VÀ CẢNH BÁO SỚM TÌNH HÌNH NHIỄM RẦY NÂU TRÊN LÚA Chuyên ngành: Bản đồ, viễn thám hệ thông tin địa lý LUẬN VĂN THẠC SĨ TP.HỒ CHÍ MINH, tháng năm 2009 ii CƠNG TRÌNH ĐƯỢC HỒN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Cán hướng dẫn khoa học: TS TRẦN TRỌNG ĐỨC Cán chấm nhận xét 1: TS HỒ ĐÌNH DUẪN Cán chấm nhận xét 2: TS CAO VĂN PHỤNG Luận văn thạc sĩ bảo vệ HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA, ngày 05 tháng 09 năm 2009 iii ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HOÀ Xà HỘI CHỦ NGHĨA VIỆT NAM Độc Lập - Tự Do - Hạnh Phúc -oOo Tp HCM, ngày tháng năm 2009 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: TRƯƠNG CHÍ QUANG Giới tính : Nam 6/ Nữ Ngày, tháng, năm sinh : 03-12-1978 Nơi sinh : Cần Thơ Chuyên ngành : Bản đồ, viễn thám hệ thơng tin địa lý Khố (Năm trúng tuyển) : 2007 1- TÊN ĐỀ TÀI: ỨNG DỤNG GIS TRONG QUẢN LÝ CANH TÁC LÚA VÀ CẢNH BÁO SỚM TÌNH HÌNH RẦY NÂU TRÊN LÚA 2- NHIỆM VỤ LUẬN VĂN: Nghiên cứu tình hình canh tác lúa dịch hại lúa tỉnh Đồng Tháp Xây dựng ứng dụng độc lập không phụ thuộc vào hệ GIS thương mại để quản lý liệu tạo đồ quản lý tình hình xuống giống, cấu giống lúa, tình hình rầy nâu lúa xã sở sử dụng sở liệu không gian để tổ chức liệu sử dụng thư viện nguồn mở để lập trình truy xuất liệu hình học Ứng dụng cho phép tạo đồ cảnh báo tình hình rầy nâu đầu vụ dựa giống lúa 3- NGÀY GIAO NHIỆM VỤ: 02-02-2009 4- NGÀY HOÀN THÀNH NHIỆM VỤ 03-07-2009 5- HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: TS TRẦN TRỌNG ĐỨC Nội dung đề cương Luận văn thạc sĩ Hội Đồng Chuyên Ngành thông qua CÁN BỘ HƯỚNG DẪN (Họ tên chữ ký) CHỦ NHIỆM BỘ MÔN QUẢN LÝ CHUYÊN NGÀNH (Họ tên chữ ký) iv LỜI CẢM TẠ Qua thời gian học tập rèn luyện trường Đại Học Bách Khoa TP Hồ Chí Minh, với nỗ lực phấn đấu thân với nhiệt tình giảng dạy hướng dẫn quý Thầy Cô trường đại học Bách Khoa Tp Hồ Chí Minh giúp tơi hồn tất khóa học cách tốt đẹp Tôi xin chân thành cảm tạ đến TS Trần Trọng Đức, môn Địa Tin Học, khoa Kỹ Thuật Xây Dựng Thầy tận tâm hướng dẫn, truyền đạt kiến thức, kinh nghiệm quý báu suốt trình học tập thực luận văn tốt nghiệp Xin cám ơn đến toàn thể quý Thầy Cô môn Địa Tin Học, khoa Kỹ Thuật Xây Dựng, trường Đại Học Bách Khoa Xin cám ơn Ban chủ nhiệm môn, Thầy Cô bạn đồng nghiệp môn Khoa Học Đất & QLĐĐ, Khoa Nông Nghiệp Sinh Học Ứng Dụng, trường Đại Học Cần Thơ tạo điều kiện hỗ trợ cho tơi tham dự hồn thành khoá học Xin dâng lời biết ơn sâu sắc đến ba má, anh chị em gia đình hỗ trợ tinh thần vật chất để an tâm hồn thành q trình học tập thực hiên tốt đề tài Cuối xin gởi lời cám ơn đến người vợ yêu quý động viên khích lệ tơi từ tơi bước chân vào khố học Một lần tơi xin trân trọng bày tỏ lịng biết ơn với tất cả! TRƯƠNG CHÍ QUANG v TÓM TẮT LUẬN VĂN THẠC SĨ Dịch rầy nâu bộc phát diện rộng vùng Đồng Bằng Sông Cửu Long xác định chủ yếu người dân canh tác liên tục sử dụng giống lúa nhiễm rầy diện tích lớn Do vậy, giám sát diện tích canh tác giống lúa giai đoạn phát triển lúa tiền đề để triển khai biện pháp phòng trừ dịch hại rầy nâu hiệu Để làm điều cần thiết xây dựng hệ thống quản lý liệu canh tác lúa Luận văn nghiên cứu xây dựng ứng dụng GIS ngôn ngữ lập trình Visual Studio.NET 2008 thư viện mã nguồn mở SharpMap Dữ liệu khơng gian thuộc tính lưu sở liệu không gian SQL Server 2008 nên dễ dàng thực thi tốn phân tích liệu khơng gian Ứng dụng có chức cho phép quản lý tiến độ gieo sạ, loại giống lúa diện tích gieo trồng; quản lý tình hình dịch hại địa phương theo thời gian; xây dựng đồ thời vụ gieo trồng, phân bố giống, tỉ lệ giống sử dụng theo khuyến cáo, dịch hại phát sinh theo thời điểm, đồ cảnh báo nguy nhiễm dịch hại rầy nâu Hệ thống thực kiểm chứng sử dụng liệu canh tác tỉnh Đồng Tháp Qua tìm mối tương quan tuyến tính diện tích canh tác giống kháng rầy diện tích nhiễm rầy cấp huyện Áp dụng phương trình hồi quy tuyến tính xây dựng đồ cảnh báo diện tích nhiễm rầy vụ Đông Xuân 20082009, so sánh với liệu nhiễm rầy thật tỉnh vụ cho thấy đồ cảnh báo giải thích hầu hết diện tích nhiễm rầy thật tỉnh Nếu triển khai sử dụng hệ thống công cụ trợ giúp đắc lực cho nhà quản lý trồng trọt bảo vệ thực vật cấp công tác quản lý mùa vụ, hỗ trợ cơng tác phịng tránh dịch hại hiệu Từ khoá: GIS sản xuất lúa, cảnh báo sớm rầy nâu, SQL Server 2008, SharpMap vi MỤC LỤC LỜI CẢM TẠ iv TÓM TẮT LUẬN VĂN THẠC SĨ v MỤC LỤC vi DANH SÁCH HÌNH MINH HOẠ ix DANH SÁCH BẢNG xii CHƯƠNG TỔNG QUAN 1.1 ĐẶT VẤN ĐỀ 1.2 MỤC TIÊU NGHIÊN CỨU 1.3 NỘI DUNG NGHIÊN CỨU 1.4 PHẠM VI NGHIÊN CỨU .3 CHƯƠNG TÌNH HÌNH NGHIÊN CỨU 2.1 CÔNG TÁC QUẢN LÝ NHÀ NƯỚC VỀ TRỒNG TRỌT 2.1.1 Cơ quan quản lý .4 2.1.2 Công tác quản lý 2.2 TÌNH HÌNH RẦY NÂU GÂY HẠI TRÊN LÚA 2.2.1 Tình hình rầy nâu gây hại lúa nước ta 2.2.2 Nguyên nhân bộc phát rầy nâu ĐBSCL .7 2.3 CÁC LOẠI DỰ BÁO DỊCH HẠI CÂY TRỒNG 2.4 TÌNH HÌNH NGHIÊN CỨU NGỒI NƯỚC 2.5 TÌNH HÌNH NGHIÊN CỨU TRONG NƯỚC .9 CHƯƠNG GIỚI THIỆU CÔNG NGHỆ VÀ PHƯƠNG PHÁP NGHIÊN CỨU 11 3.1 GIỚI THIỆU VỀ HỆ THỐNG THÔNG TIN ĐỊA LÝ 11 3.1.1 Hệ thống thông tin địa lý (GIS) .11 3.1.2 Các hợp phần GIS 11 3.1.3 Ứng dụng GIS 13 3.2 GIỚI THIỆU VỀ CƠ SỞ DỮ LIỆU KHÔNG GIAN 13 3.2.1 Định nghĩa .14 3.2.2 Kiểu liệu không gian 14 vii 3.2.3 Các hệ quản trị sở liệu không gian 16 3.2.4 Lựa chọn hệ quản trị sở liệu 18 3.3 HIỂN THỊ DỮ LIỆU KHÔNG GIAN 20 3.4 THỐNG KÊ HỒI QUY TUYẾN TÍNH 21 3.4.1 Hồi quy tuyến tính 21 3.4.2 Hồi quy tuyến tính đơn 22 3.4.3 Hệ số xác định hệ số tương quan 22 3.5 Kết luận lựa chọn công nghệ phương pháp nghiên cứu 22 CHƯƠNG PHÂN TÍCH VÀ THIẾT KẾ ỨNG DỤNG .24 4.1 PHÂN TÍCH YÊU CẦU .24 4.1.1 Yêu cầu chung .24 4.1.2 Yêu cầu chi tiết 24 4.2 PHÂN TÍCH CHỨC NĂNG 25 4.2.1 Nhóm chức hệ thống 25 4.2.2 Nhóm chức nhập xuất liệu .25 4.2.3 Nhóm chức phân tích, thống kê liệu 26 4.2.4 Nhóm chức thành lập đồ 26 4.3 THIẾT KẾ CƠ SỞ DỮ LIỆU 26 4.3.1 Dữ liệu không gian 26 4.3.2 Mơ tả liệu thuộc tính .27 4.3.3 Mơ hình thực thể - liên kết 27 4.3.4 Xây dựng sở liệu 28 4.4 THIẾT KẾ ỨNG DỤNG .34 4.4.1 Mơ hình tổ chức quản lý 34 4.4.2 Thiết kế giao diện 35 4.4.3 Kết nối đến sở liệu .35 4.4.4 Quản lý đăng nhập hệ thống 37 4.4.5 Mơ hình luồng liệu 38 4.5 NHẬN XÉT 61 viii CHƯƠNG ỨNG DỤNG 62 5.1 GIỚI THIỆU KHU VỰC ỨNG DỤNG 62 5.2 Yêu cầu phần cứng phần mềm 62 5.3 Ứng dụng .63 5.3.1 Thiết lập kết nối, quản lý người dùng 63 5.3.2 Nhập liệu 66 5.3.3 Báo cáo tình hình canh tác 70 5.3.4 Phân tích tương quan giống gieo trồng mức độ nhiễm rầy .72 5.3.5 Bản đồ tình hình canh tác 75 5.3.6 Bản đồ cảnh báo sớm diện tích nhiễm rầy nâu 79 CHƯƠNG KẾT LUẬN KIẾN NGHỊ 84 6.1 KẾT LUẬN 84 6.2 KIẾN NGHỊ 84 TÀI LIỆU THAM KHẢO .86 ix DANH SÁCH HÌNH MINH HOẠ Hình 2.1 Mơ hình phân cấp quản lý ngành nơng nghiệp Hình 2.2 Tình hình nhiễm rầy nâu tỉnh ĐBSCL (Tổng hợp từ nguồn TT BVTV phía Nam) Hình 3.1 Các thành phần GIS (Nguồn J.E Harmon and S.J Anderson, 2003).12 Hình 3.2 Sơ đồ phân cấp đối tượng dạng Geometry .15 Hình 4.1 Sơ đồ chức ứng dụng 25 Hình 4.2 Mơ hình thực thể kết hợp 28 Hình 4.3 Sơ đồ quan hệ bảng sở liệu 29 Hình 4.4 Mơ hình kết nối trao đổi liệu cấp 34 Hình 4.5 Xây dựng giao diện 35 Hình 4.6 Cửa sổ thiết lập kết nối sở liệu .36 Hình 4.7 Lưu đồ mã hoá lưu chuỗi kết nối 36 Hình 4.8 Thiết kế khung đăng nhập hệ thống 37 Hình 4.9 Lưu đồ kiểm tra người dùng 37 Hình 4.10 Mơ hình luồng liệu ứng dụng .38 Hình 4.11 Cửa sổ nhập diện tích kế hoạch xuống giống 39 Hình 4.12 Thiết kế giao diện nhập liệu xuống giống thu hoạch .40 Hình 4.13 Lưu đồ tổng quát lưu liệu xuống giống thu hoạch 41 Hình 4.14 Giao diện thống kê diện tích xuống giống .42 Hình 4.15 Lưu đồ thống kê diện tích xuống giống 43 Hình 4.16 Thống kê diện tích giai đoạn sinh trưởng lúa 44 Hình 4.17 Lưu đồ thống kê giai đoạn sinh trưởng lúa 45 Hình 4.18 Giao diện xây dựng đồ xuống giống 46 Hình 4.19 Bản đồ mẫu diện tích xuống giống 46 Hình 4.20 Lưu đồ trình xây dựng đồ xuống giống 47 Hình 4.21 Lưu đồ thực thủ tục tính diện tích chia khoảng biểu diễn 48 x Hình 4.22 Cửa sổ xây dựng diện tích giống lúa 49 Hình 4.23 Minh hoạ đồ phân bố loại giống lúa 49 Hình 4.24 Lưu đồ thực xử lý liệu xây dựng đồ phân bố giống lúa 50 Hình 4.25 Store Procedure xử lý liệu đồ giống lúa 50 Hình 4.26 Cửa sổ nhập liệu tình hình nhiễm dịch hại .51 Hình 4.27 Lưu đồ nhập liệu dịch hại 52 Hình 4.28 Cửa sổ báo cáo tình hình dịch hại 53 Hình 4.29 Lưu đồ xây dựng loại báo cáo nhiễm dịch hại 54 Hình 4.30 Thiết kế cửa sổ tạo đồ diên tích nhiễm rầy nâu .55 Hình 4.31 Bản đồ mẫu phân bố diện tích nhiễm rầy nâu 55 Hình 4.32 Lưu đồ xử lý liệu đồ phân bố sâu bệnh 56 Hình 4.33 Cửa sổ xây dựng đồ phân bố giống lúa nhiễm rầy 56 Hình 4.34 Lưu đồ bước xây dựng đồ phân bố giống nhiễm rầy 57 Hình 4.35 Bản đồ mẫu khu vực bị nhiễm rầy nâu qua nhiều vụ .58 Hình 4.36 Store procedure xây dựng liệu tình hình nhiễm rầy nhiều vụ 58 Hình 4.37 Cửa sổ xây dựng đồ cảnh báo rầy .59 Hình 4.38 Lưu đồ thủ tục xây dựng đồ cảnh báo .60 Hình 4.39 Bản đồ mẫu khu vực nguy nhiễm rầy nâu .61 Hình 5.1 Thiết lập kết nối đến sở liệu 63 Hình 5.2 Cửa sổ quản lý người dùng 64 Hình 5.3 Đăng nhập chương trình 65 Hình 5.4 Giao diện ứng dụng 66 Hình 5.5 Bản đồ hành tỉnh Đồng Tháp 67 Hình 5.6 Đưa shapefile vào sở liệu .68 Hình 5.7 Nhập cấu giống vụ Đơng Xuân 2008-2009 69 Hình 5.8 Nhập liệu tình hình rầy nâu tỉnh 69 Hình 5.9 Báo cáo tình hình xuống giống 70 PHỤ LỤC MỘT SỐ THỦ TỤC XỬ LÝ DỮ LIỆU BẢN ĐỒ THỦ TỤC TẠO BẢN ĐỒ XUỒNG GIỐNG USE [dbPlant] GO SET ANSI_NULLS ON SET QUOTED_IDENTIFIER ON GO ============================================= Author: Create date: Description: ============================================= CREATE PROCEDURE [dbo].[BandoXuonggiong] Add the parameters for the stored procedure here @nam nvarchar(50), @tenvu nvarchar(50), @Khoang float out, @min float out AS BEGIN Kiem tra xe bang tmp_Dientich1 da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich1' And type = 'U') drop table tmp_Dientich1 Kiem tra xe bang tmp_Dientich2 da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich2' And type = 'U') drop table tmp_Dientich2 Kiem tra xe bang tmp_Dientich3 da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich3' And type = 'U') drop table tmp_Dientich3 Kiem tra xe bang tmp_Dientich4 da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich4' And type = 'U') drop table tmp_Dientich4 Kiem tra xe bang tmp_Dientich4 da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich5' And type = 'U') drop table tmp_Dientich5 Tinh dien tich xuong giong cua tung xa vu va nam duoc nhap select maxa ,sum(dtxuonggiongps) as tongdientich into #tmp_Dientichxg from xuongiong where nam = @nam and tenvu=@tenvu group by maxa declare @DtMax as float declare @DtMin as float set @Dtmax = (select MAX(tongdientich) from #tmp_Dientichxg ) set @Dtmin = (select min(tongdientich) from #tmp_Dientichxg ) select @Khoang=(@DtMax-@DtMin)/5 set @min = @DtMin if (@Khoang>0) Chia dientich khoang bang begin select xa.*,tongdientich into tmp_Dientich1 from xa,#tmp_Dientichxg where xa.maxa=#tmp_Dientichxg.maxa and #tmp_Dientichxg.tongdientich @Khoang and #tmp_Dientichxg.tongdientich =@Khoang*2 and #tmp_Dientichxg.tongdientich =@Khoang*3 and #tmp_Dientichxg.tongdientich =@Khoang*4 end else set @khoang=0 drop table #tmp_Dientichxg END GO THỦ TỤC TẠO BẢN ĐỒ PHÂN BỐ DIỆN TÍCH GIỐNG NHIỄM RẦY USE [dbPlant] GO SET ANSI_NULLS ON GO SET QUOTED_IDENTIFIER ON GO ============================================= Author: Truong Chi Quang Create date: 23/6/2009 Description: Tao dư lieu ban dien tich cac giong nhiem, chia lam khoang bieu dien> ============================================= CREATE PROCEDURE [dbo].[BandoGiongNhiemray] Add the parameters for the stored procedure here @nam nvarchar(50), @tenvu nvarchar(50) AS BEGIN Kiem tra xe bang tmp_Dientich1 da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich1' And type = 'U') drop table tmp_Dientich1 Kiem tra xe bang tmp_Dientich2 da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich2' And type = 'U') drop table tmp_Dientich2 Kiem tra xe bang tmp_Dientich3 da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich3' And type = 'U') drop table tmp_Dientich3 Kiem tra xe bang tmp_Dientich4 da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich4' And type = 'U') drop table tmp_Dientich4 Kiem tra xe bang tmp_Dientich4 da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich5' And type = 'U') drop table tmp_Dientich5 Tính tong dien tich cac giong lua nhiem ray xã select maxa ,sum(dtxuonggiongps) as [dtgiong] into #tmp_DientichGiong from xuongiong where nam = @nam and tenvu=@tenvu and tengiong in (select distinct tengiong from giongnhiemray) group by maxa Dien tich xuong giong tung xa select maxa ,sum(dtxuonggiongps) as [dtxuonggiong] into #tmp_Dientichxg from xuongiong where nam = @nam and tenvu=@tenvu group by maxa Tim % dien tich canh tac giong lua nhiem ray select #tmp_Dientichxg.maxa, dtgiong,dtxuonggiong, dtgiong*100/dtxuonggiong as phantramxg into #tmp_DtGiongnhiem from #tmp_DientichGiong,#tmp_Dientichxg where #tmp_DientichGiong.maxa=#tmp_Dientichxg.maxa set @Khoang=(select @@ROWCOUNT ) select xa.*,phantramxg into tmp_Dientich1 from xa,#tmp_DtGiongnhiem where xa.maxa=#tmp_DtGiongnhiem.maxa and #tmp_DtGiongnhiem.phantramxg < 20 select xa.*,phantramxg into tmp_Dientich2 from xa,#tmp_DtGiongnhiem where xa.maxa=#tmp_DtGiongnhiem.maxa and #tmp_DtGiongnhiem.phantramxg >=20 and #tmp_DtGiongnhiem.phantramxg < 40 select xa.*,phantramxg into tmp_Dientich3 from xa,#tmp_DtGiongnhiem where xa.maxa=#tmp_DtGiongnhiem.maxa and #tmp_DtGiongnhiem.phantramxg >=40 and #tmp_DtGiongnhiem.phantramxg < 60 select xa.*,phantramxg into tmp_Dientich4 from xa,#tmp_DtGiongnhiem where xa.maxa=#tmp_DtGiongnhiem.maxa and #tmp_DtGiongnhiem.phantramxg >=60 and #tmp_DtGiongnhiem.phantramxg < 80 select xa.*,phantramxg into tmp_Dientich5 from xa,#tmp_DtGiongnhiem where xa.maxa=#tmp_DtGiongnhiem.maxa and #tmp_DtGiongnhiem.phantramxg >= 80 Xoá bảng tạm drop table #tmp_DientichGiong drop table #tmp_Dientichxg drop table #tmp_DtGiongnhiem END GO THỦ TỤC TẠO BẢN ĐỒ CẢNH BÁO SỚM RẦY NÂU USE [dbPlant] GO /****** Object: StoredProcedure [dbo].[BandoCanhbao] Script Date: 06/30/2009 16:58:19 ******/ SET ANSI_NULLS ON GO SET QUOTED_IDENTIFIER ON GO ============================================= Author: Truong Chi QUANG Create date: 23/6/2009 Description: Ban canh bao ray ============================================= CREATE PROCEDURE [dbo].[BandoCanhbao] Add the parameters for the stored procedure here @nam nvarchar(50), @tenvu nvarchar(50), @Hangso float = 0, @HesoX float =1, @TrongsoGiong as float= 1, @TrongsoLichsu as float=1, @Khoang float=0 out, @min float=0 out AS BEGIN Kiem tra xe bang tmpDulieu da co chua neu co thi xoa di Bang chua du lieu cua ban If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmpDulieu' And type = 'U') drop table tmpDulieu Kiem tra xe bang tblDientichxa_tmp da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich1' And type = 'U') drop table tmp_Dientich1 Kiem tra xe bang tblDientichxa_tmp da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich2' And type = 'U') drop table tmp_Dientich2 Kiem tra xe bang tblDientichxa_tmp da co chua neu co thi xoa di If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich3' And type = 'U') drop table tmp_Dientich3 Tinh dien tich xuong giong casc giong lua nhiem ray cua tung xa If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich4' And type = 'U') drop table tmp_Dientich4 If Exists (Select * from dbplant.dbo.sysobjects Where Name = 'tmp_Dientich5' And type = 'U') drop table tmp_Dientich5 Din tich xuong giong nhiem ray tung huyen select mahuyen ,sum(dtxuonggiongps) as [dtgiong] into #tmp_DientichGiong from xuongiong,xa where xuongiong.maxa=xa.maxa and nam = @nam and tenvu=@tenvu and tengiong in (select distinct tengiong from giongnhiemray) group by mahuyen Din tich xuong giong tung huyen select mahuyen ,sum(dtxuonggiongps) as [dtxuonggiong] into #tmp_Dientichxg from xuongiong,xa where xuongiong.maxa=xa.maxa and nam = @nam and tenvu=@tenvu group by mahuyen Tim % dien tich canh tac giong lua nhiem ray select #tmp_Dientichxg.mahuyen, dtgiong,dtxuonggiong, (@Hangso+(100-dtgiong*100/dtxuonggiong)* @TrongsoGiong*@HesoX) as phantramxg into #tmp_DtGiongnhiem from #tmp_DientichGiong,#tmp_Dientichxg where #tmp_DientichGiong.mahuyen=#tmp_Dientichxg.mahuyen - Tinh dien tich nhiem ray lon nhat cua cac xa moi vu select xa.maxa,xa.mahuyen, tenvu,nam,max(dtnnhe+dtntb+dtnnang) as maxdientich into #tmp_DientichRayxa from nhiemsaubenh,xa where nhiemsaubenh.maxa=xa.maxa and masb ='29' and nam @nam group by xa.maxa,xa.mahuyen,tenvu,nam Tinh tong dien tich nhiem ray lon nhat moi huyen vu select mahuyen,tenvu,nam,SUM(maxdientich) as tongdientich into #tmp_DientichRay from #tmp_DientichRayxa group by mahuyen,tenvu,nam dien tich trung binh nhiem ray cac xa nhieu vu select mahuyen, avg(tongdientich) as dttrungbinh into #tmp_dtRayTrungbinh from #tmp_DientichRay group by mahuyen -% dien tich nhiem ray = dtnhiem ray/dt xuong giong select #tmp_Dientichxg.mahuyen, dttrungbinh,dtxuonggiong, dttrungbinh*100*@TrongsoLichsu/dtxuonggiong as phantramnhiem into #tmp_DientichNhiemray from #tmp_dtRayTrungbinh,#tmp_Dientichxg where #tmp_dtRayTrungbinh.mahuyen=#tmp_Dientichxg.mahuyen select #tmp_DtGiongnhiem.mahuyen, phantramnhiem,phantramxg, (phantramnhiem+phantramxg)/2 as dientichcanhbao into #tmp_DientichCbao from #tmp_DtGiongnhiem LEFT OUTER JOIN #tmp_DientichNhiemray ON #tmp_DtGiongnhiem.mahuyen=#tmp_DientichNhiemray.mahuyen Cap nhât cac cot dientichcanh bao = null gia trị dtgiong canh tac nhiem ray da tinh bang pt tuong quan update #tmp_DientichCbao set dientichcanhbao= phantramxg where dientichcanhbao is NULL set @Khoang=(select @@ROWCOUNT ) select xa.*,dientichcanhbao into tmp_Dientich1 from xa,#tmp_DientichCbao where xa.mahuyen=#tmp_DientichCbao.mahuyen and #tmp_DientichCbao.dientichcanhbao < 20 select xa.*,dientichcanhbao into tmp_Dientich2 from xa,#tmp_DientichCbao where xa.mahuyen=#tmp_DientichCbao.mahuyen and #tmp_DientichCbao.dientichcanhbao >=20 and #tmp_DientichCbao.dientichcanhbao < 40 select xa.*,dientichcanhbao into tmp_Dientich3 from xa,#tmp_DientichCbao where xa.mahuyen=#tmp_DientichCbao.mahuyen and #tmp_DientichCbao.dientichcanhbao >=40 and #tmp_DientichCbao.dientichcanhbao < 60 select xa.*,dientichcanhbao into tmp_Dientich4 from xa,#tmp_DientichCbao where xa.mahuyen=#tmp_DientichCbao.mahuyen and #tmp_DientichCbao.dientichcanhbao >=60 and #tmp_DientichCbao.dientichcanhbao < 80 select xa.*,dientichcanhbao into tmp_Dientich5 from xa,#tmp_DientichCbao where xa.mahuyen=#tmp_DientichCbao.mahuyen and #tmp_DientichCbao.dientichcanhbao >= 80 Xuất liệu bảng tam để chương trình gọi hiển thị select #tmp_DientichCbao.* into tmpDulieu from #tmp_DientichCbao Xố bảng tạm q trình tính tốn drop table #tmp_DtGiongnhiem drop table #tmp_DientichNhiemray drop table #tmp_DientichRay drop table #tmp_dtRayTrungbinh drop table #tmp_DientichGiong drop table #tmp_Dientichxg drop table #tmp_DientichCbao END GO PHỤ LỤC LỆNH LẬP TRÌNH HIỂN THỊ BẢN ĐỒ Private Sub BandoDichhai () ' Chuan bi du lieu xay duwng ban phan bo dich hại 'Goi store procedure Tinh tong dien tich nhiem dich hai cua cac xã ' xử lý chi khoảng diện tich bang de hien thi ban 'Store tra ve ket qua gis trij chia khoang neu gis trij tra ve =0 tuc kong ti thay du lieuj ' Store ao table ket qua: tmp_Dientich1,tmp_Dientich2,tmp_Dientich3,tmp_Dientich4,tmp_Dientich5 Dim StoreName As String Dim v_Nam As String Dim t_Ketqua As New DataTable Dim st_dotdtra As String = "" Dim v_khoang As Double = 0.0 ' Luu gia trij cuar khoang dien tich Dim v_min As Integer = ' Luu gia nho nhat cua dien tich tra ve v_Nam = Trim(txtNam.Text) If chkTatca.Checked = True Then v_soVu = -1 'Else 'v_soVu = End If StoreName = "BandoDichhai_DT" '' Goi store procedure tao ban xuonggiong 'If Not ParaIsNULL() Then Dim tmp As String = "" conn = New SqlConnection conn.ConnectionString = con_str Dim SqlDataAdapter1 As New SqlDataAdapter() Dim CmdString As New SqlCommand(StoreName, conn) CmdString.CommandType = CommandType.StoredProcedure CmdString.Parameters.Add(New SqlParameter("@khoang", System.Data.SqlDbType.Int, 20, "khoang")) CmdString.Parameters("@khoang").Direction = Data.ParameterDirection.Output CmdString.Parameters.Add(New SqlParameter("@min", System.Data.SqlDbType.Int, 20, "min")) CmdString.Parameters("@min").Direction = Data.ParameterDirection.Output conn.Open() CmdString.ExecuteNonQuery() v_khoang = CmdString.Parameters("@khoang").Value v_min = CmdString.Parameters("@min").Value conn.Close() If v_khoang > Then pBandoXuonggiong(v_min, v_khoang) Else MessageBox.Show("Không tìm thấy liệu để xây dựng đồ dịch hại.") End If 'End If End Sub Private Sub pBandoXuonggiong(ByVal v_min As Int16, ByVal v_khoang As Int16) 'Lop xa FrmMain.map.Layers.Clear() Dim layCountries = New SharpMap.Layers.VectorLayer("Lop xa") layCountries.DataSource = New SharpMap.Data.Providers.Katmai(FrmMain.ConnStr, "xa", "geom", "ID") layCountries.Style.Fill = New SolidBrush(Color.White) FrmMain.map.Layers.Add(layCountries) layCountries.Style.Outline = System.Drawing.Pens.DimGray layCountries.Style.EnableOutline = True Dim layDientich1 = New SharpMap.Layers.VectorLayer("Dientich1") layDientich1.DataSource = New SharpMap.Data.Providers.Katmai(FrmMain.ConnStr, "tmp_Dientich1", "geom", "ID") layDientich1.Style.Fill = New SolidBrush(Color.PaleGreen) FrmMain.map.Layers.Add(layDientich1) '//Set the polygons to have a black outline layDientich1.Style.Outline = System.Drawing.Pens.DimGray layDientich1.Style.EnableOutline = True Dim layDientich2 = New SharpMap.Layers.VectorLayer("Dientich2") layDientich2.DataSource = New SharpMap.Data.Providers.Katmai(FrmMain.ConnStr, "tmp_Dientich2", "geom", "ID") layDientich2.Style.Fill = New SolidBrush(Color.SpringGreen) FrmMain.map.Layers.Add(layDientich2) '//Set the polygons to have a black outline layDientich2.Style.Outline = System.Drawing.Pens.DimGray layDientich2.Style.EnableOutline = True Dim layDientich3 = New SharpMap.Layers.VectorLayer("Dientich3") layDientich3.DataSource = New SharpMap.Data.Providers.Katmai(FrmMain.ConnStr, "tmp_Dientich3", "geom", "ID") layDientich3.Style.Fill = New SolidBrush(Color.LimeGreen) FrmMain.map.Layers.Add(layDientich3) '//Set the polygons to have a black outline layDientich3.Style.Outline = System.Drawing.Pens.DimGray layDientich3.Style.EnableOutline = True Dim layDientich4 = New SharpMap.Layers.VectorLayer("Dientich4") layDientich4.DataSource = New SharpMap.Data.Providers.Katmai(FrmMain.ConnStr, "tmp_Dientich4", "geom", "ID") layDientich4.Style.Fill = New SolidBrush(Color.ForestGreen) FrmMain.map.Layers.Add(layDientich4) '//Set the polygons to have a black outline layDientich4.Style.Outline = System.Drawing.Pens.DimGray layDientich4.Style.EnableOutline = True Dim layDientich5 = New SharpMap.Layers.VectorLayer("Dientich5") layDientich5.DataSource = New SharpMap.Data.Providers.Katmai(FrmMain.ConnStr, "tmp_Dientich5", "geom", "ID") layDientich5.Style.Fill = New SolidBrush(Color.DarkGreen) FrmMain.map.Layers.Add(layDientich5) '//Set the polygons to have a black outline layDientich5.Style.Outline = System.Drawing.Pens.DimGray layDientich5.Style.EnableOutline = True ' Lay yer Ranh huyen Dim layRanhuyen = New SharpMap.Layers.VectorLayer("Ranh huyen") layRanhuyen.DataSource = New SharpMap.Data.Providers.Katmai(FrmMain.ConnStr, "ranhhuyen", "geom", "ID") layRanhuyen.Style.Line.Color = Color.Blue layRanhuyen.Style.Line.Width = 'layRanhuyen.Style.Line = Pens.AliceBlue FrmMain.map.Layers.Add(layRanhuyen) ' Lay yer lopduong Dim layChanel = New SharpMap.Layers.VectorLayer("Lop duong") layChanel.DataSource = New SharpMap.Data.Providers.Katmai(FrmMain.ConnStr, "lopduong", "geom", "ID") layChanel.Style.Line.Color = Color.Red FrmMain.map.Layers.Add(layChanel) ' Lop khung ban Dim layKhung = New SharpMap.Layers.VectorLayer("Lop khung") layKhung.DataSource = New SharpMap.Data.Providers.Katmai(FrmMain.ConnStr, "khung", "geom", "ID") layKhung.Style.Line.Color = Color.Black FrmMain.map.Layers.Add(layKhung) FrmMain.map.BackColor = Color.White 'LightBlue capnhatChudan() FrmMain.lblChudan.Text = "Diện tích nhiễm (ha)" FrmMain.lbl1.Text = v_min.ToString & "-" & (v_min + v_khoang).ToString FrmMain.lbl2.Text = (v_min + v_khoang).ToString & "-" & (v_min + * v_khoang).ToString FrmMain.lbl3.Text = (v_min + * v_khoang).ToString & "-" & (v_min + * v_khoang).ToString FrmMain.lbl4.Text = (v_min + * v_khoang).ToString & "-" & (v_min + * v_khoang).ToString FrmMain.lbl5.Text = (v_min + * v_khoang).ToString & "-" & (v_min + * v_khoang).ToString FrmMain.MapImage.Refresh() End Sub PHỤ LỤC BẢNG TỔNG HỢP DỮ LIỆU RẦY NÂU VỤ TỈNH ĐỒNG THÁP VỤ ĐÔNG XUÂN 2008-2009 TT Huyện thị 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 Tân Hồng OM 1490 OM 2718 OM 3536 OMCS 2000 VNĐ 95-20 3.245,5 Tài nguyên 128 3.245,6 IR 64 OM 2517 OM 2518 OM 2717 OM 2514 JASMINE 10.819,0 VĐ 20 OM 4085 OM 4088 MTL 250 HĐ OM 4498 OM 4668 OM 4495 OM 5930 OM 6073 OM 6162 OM 4218 OM 4900 Hồng Ngự 597,3 20,0 249,1 742,4 2.200,1 1.200,6 200,1 230,1 3.611,9 1.778,9 296,2 300,2 30,0 56,0 130,1 323,2 0,0 1.750,9 678,4 Tam Nơng 511,3 Th Bình Cao Lãnh 2.210,1 Th Lấp Vò Mười 178,6 626,5 198,5 465,2 0,0 1.342,0 111,0 2.818,0 255,0 2.525,7 284,0 128,9 19,0 1.578,6 472,0 10.144,9 75,0 2.219,9 785,7 389,0 118,6 2.525,7 581,0 321,0 758,8 1.784,9 4.735,8 5.647,0 662,0 167,0 40,0 550,8 5,0 250,0 98,5 947,1 253,0 1.470,5 2.506,8 0,0 437,0 112,0 0,0 3.805,2 8.437,9 5.657,4 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 526,0 248,8 511,1 248,2 503,0 1.070,0 288,0 Lai Vung 385,0 Ch.Thành TX SaĐéc 49,2 643,8 179,5 310,0 TPCL 175,6 32,4 109,0 1.000,0 8.086,0 6.434,0 143,0 757,0 2.841,5 25,0 249,0 150,0 326,0 138,0 769,0 615,0 72,0 652,5 148,0 91,2 94,9 430,3 47,4 0,0 873,8 170,0 30,5 38,2 26 27 28 29 30 II III OM 2490 OM 6561 OM 5199 AS 996 OM 4655 Giống khác Nếp IR 50404 BUTYL (IR13240-108) 30,0 230,0 0,5 8,0 947,1 2.525,7 1.495,4 4,0 597,1 9.787,4 9.578,9 21263,0 30196,0 20966,0 31572,0 34808,0 100,0 2.163,7 178,1 175,1 2.163,7 6.484,4 21637,5 Tổng DT Xuống giống 947,1 1.357,7 0,0 2.456,4 1.656,8 1.319,0 900,0 40,0 217,0 1.697,3 4,0 6.283,2 10.162,8 14081,0 14066,0 654,1 173,3 284,9 8.586,5 1.368,0 2.008,5 346,6 1824,0 3347,0 13551,0 PHỤ LỤC TỔNG HỢP DỮ LIỆU NHIỄM RẦY NÂU ĐỢT 24/2/2009 VỤ ĐÔNG XUÂN TỈNH ĐỒNG THÁP MAXA TENXA TENHUYEN DTNNHE DTTB DT MD NHE MDTB MDNANG 3000-4000 11608 An B×nh Cao L·nh 0 700-1000 11616 Mü Long Cao L·nh 195 15 1500-3000 11609 Mü Thä Cao L·nh 470 1000-2000 11614 Bình Hng Trung 11607 Nhị Mỹ Cao LÃnh 100 0 700-1000 Cao L·nh 0 1000-2000 11606 Ba Sao Cao L·nh 993 1500-2000 11605 Ph¬ng Trμ Cao LÃnh 300 0 500-1000 11603 Phơng Thịnh Cao L·nh 225 0 1000-1500 11613 T©n Héi Trung Cao LÃnh 450 0 750-1500 11615 Bình Hng Tây Cao L·nh 200 0 1000-1500 11617 Mü HiÖp Cao L·nh 159 15 750-1500 11602 G¸o Giång Cao L·nh 500 0 500-1000 11601 Phong Mü Cao L·nh 75 0 750-1500 1500-2000 111002 Tân Bình Châu Thnh 200 30 750-1500 3000-5000 111010 Phó Hùu Ch©u Thμnh 33 750-1500 3000-5000 111009 An Khánh Châu Thnh 684 23 750-1500 3000-5000 111008 Hoμ T©n Ch©u Thμnh 0 10 111001 T©n Phó Trung Ch©u Thμnh 674 126 750-1500 3000-5000 1500-3000 11308 Thêng L¹c Hång Ngù 30 0 750-1500 11302 Thêng Thíi HËu A 11303 Thêng Thíi HËu B 11309 Thêng Thíi TiỊn 11310 Thêng Phíc Hång Ngù 30 0 750-1500 Hång Ngù 35 0 750-1500 Hång Ngù 25 10 50 750-1500 1500-3000 3000-5000 Hång Ngù 5 750-1500 1500-3000 3000-4000 11301 Thêng Phíc Hång Ngù 10 20 20 7500-1500 1500-3000 3000-5000 11306 An B×nh B Hång Ngù 45 10 750-1500 1500-3000 11314 An B×nh A Hång Ngù 70 15 750-1500 1500-3000 11305 Bình Thạnh Hồng Ngự 0 10 750-1500 1500-3000 11311 Long Kh¸nh A Hång Ngù 0 750-1500 11312 Long Kh¸nh B Hång Ngù 0 750-1500 1500-3000 11304 T©n Héi Hång Ngù 0 750-1500 1500-3000 11908 Long Th¾ng Lai Vung 50 1500-2000 11912 Phong Hoμ Lai Vung 120 1,5 750-1000 11910 T©n Hoμ Lai Vung 0 800-1300 11909 VÜnh Thíi Lai Vung 125 0 700-1000 8000-10000 5000-8000 11907 Hoμ Long Lai Vung 70 0 800-1000 11906 Hoμ Thμnh Lai Vung 20 2000-3000 11903 Long Hëu Lai Vung 30 2000-3000 11808 Héi An Đông Lấp Vò 40 30 0,5 1000-1500 1500-2000 11807 Vĩnh Thạnh Lấp Vò 25 11 1500-2500 11804 Long Hng A LÊp Vß 225 275 11109 P An Hoμ TX Sa Đéc 10 2500 11103 Xà Tân Quy Tây 11108 Tân Phú Đông TX Sa Đéc 150 140 2500 TX Sa §Ðc 0 30 11102 Tân Quy Đông TX Sa Đéc 88 11101 Tân Khánh Đông 11406 Phú Hiệp TX Sa Đéc 19 Tam N«ng 43 0 750-1500 11401 An Hoμ Tam N«ng 0 750-1500 11405 Phó Thμnh B Tam Nông 0 750-1500 11407 Phú Đức Tam Nông 34 0 750-1500 11408 Tân Công Sính Tam Nông 65 0 750-1500 11409 Ho Bình Tam Nông 32 0 750-1500 11410 Phó Cêng Tam N«ng 28 750-1500 11411 TT Trμm Chim Tam N«ng 11 0 750-1500 1000-1200 3000-5000 1500-2500 5000 3000-10000 750-1000 5000 2000-3000 1500-3000 11412 Phó Thä Tam N«ng 110 0 750-1500 11404 Phó Thμnh A Tam N«ng 15 0 750-1500 11201 Bình Phú Tân Hồng 260 2000-3000 3500-6000 11202 Tân Hộ Cơ Tân Hồng 190 2000-3000 3500-6000 11203 Thông Bình Tân Hồng 100 220 2000-3000 3500-6000 11205 T©n Thμnh B T©n Hång 260 40 2000-3000 3500-6000 11206 TT Sa Rμi T©n Hồng 190 2000-3000 3500-6000 11207 Tân Công Chí T©n Hång 450 50 2000-3000 3500-6000 11208 An Phíc T©n Hång 240 50 2000-3000 3500-6000 11209 T©n Phíc T©n Hång 310 20 2000-3000 3500-6000 30 20 750-1500 1500-3000 2000-3000 1500-3000 11507 TT Thanh B×nh Thanh Bình 11502 Tân Quới Thanh Bình 45 11503 Tân Huề Thanh Bình 10 750-1500 11504 Tân Long Thanh Bình 15 0 750-1500 11505 Tân Bình Thanh Bình 70 0 750-1000 11501 Tân Ho Thanh Bình 45 0 750-1500 11513 Tân Phú Thanh Bình 650 250 180 750-1500 1500-3000 11508 B×nh Thμnh Thanh B×nh 70 0 750-1500 1500-3000 11509 B×nh TÊn Thanh B×nh 160 0 750-1500 3000-7000 11510 An Phong Thanh B×nh 250 500 750-1500 1500-3000 11511 Phú Lợi Thanh Bình 1120 650 750-1500 1500-3000 11512 Tân Mỹ Thanh Bình 900 1200 300 750-1500 1500-3000 11506 Tân Thạnh Thanh Bình 150 0 750-1000 11702 Hng Thạnh Tháp Mời 150 500 750-1500 11703 Trờng Xuân Tháp Mời 0 >1500-2000 11705 Mü Hoμ Th¸p Mêi 0 750-1500 >1500-3000 11708 Tân Kiều Tháp Mời 300 180 750-1500 >1500-3000 11709 TT Mü An Th¸p Mêi 200 0 750-1500 11701 Thạnh Lợi Tháp Mời 1000 500 750-1500 >750-3000 11710 Mü An Th¸p Mêi 200 80 750-1500 >1500-3000 11001 P.11 TP Cao L·nh 12 20 20 750-1500 2000-3000 5000-8000 11002 Mü Ng·i TP Cao L·nh 60 85 18 750-1000 2000-3000 4000.-6000 11003 Mü T©n TP Cao L·nh 65 207 40 750-1000 2000-3000 4000-5000 11010 P.6 TP Cao L·nh 11004 Mü Trμ TP Cao L·nh 95 750-1000 4000-6000 4000-6000 >3000-5000 2000-3000 5000-8000 ... địa lý Khố (Năm trúng tuyển) : 2007 1- TÊN ĐỀ TÀI: ỨNG DỤNG GIS TRONG QUẢN LÝ CANH TÁC LÚA VÀ CẢNH BÁO SỚM TÌNH HÌNH RẦY NÂU TRÊN LÚA 2- NHIỆM VỤ LUẬN VĂN: Nghiên cứu tình hình canh tác lúa dịch... sánh nhờ vào khả biểu diễn phân bố, tìm lân cận khơng gian, phân tích liệu khơng gian… Từ vấn đề trên, luận văn ? ?Ứng dụng GIS quản lý sản xuất lúa cảnh báo sớm tình hình nhiễm rầy nâu lúa? ?? thực... phân bố giống lúa, đồ diện tích canh tác giống lúa nhiễm rầy, đồ phân bố diện tích nhiễm rầy nâu o Kết xuất đồ cảnh báo sớm tình hình rầy nâu đầu vụ dựa diện tích giống lúa nhiễm rầy nâu mối tương