Tƣơng tự nhƣ phần 1.2.3 đã chứng minh, ta có tỷ lệ giữa diện tích phẳng và diện tích nghiêng là: Snghiêng= cos 1 × Sphẳng (2.2)
Nhƣ vậy, để tính tốn diện tích của tam giác trên mặt phẳng nghiêng khi đã biết diện tích tam giác là hình chiếu trên mặt phẳng ngang và góc nghiêng, ta hồn tồn có thể tính dựa theo cơng thức nêu trên.
Trong trƣờng hợp tổng quát khi cần tính diện tích tam giác trên mặt phẳng nghiêng khơng có cạnh b nằm ngang (hình 2.4). ABC là tam giác trên mặt phẳng
nghiêng một góc α so với mặt phẳng ngang, AB’C’ là hình chiếu của tam giác ABC trên mặt phẳng nằm ngang. Kéo dài cạnh BC và B’C’ cắt nhau ở điểm D.
Hình 2.4: Tam giác nằm nghiêng trong trường hợp tổng quát
Ta có 𝑆𝐴𝐵𝐶 = 𝑆𝐴𝐷𝐶 − 𝑆𝐴𝐷𝐵 (*) 𝑆𝐴𝐵′𝐶′ = 𝑆𝐴𝐷𝐶′− 𝑆𝐴𝐷𝐵′ 𝑆𝐴𝐷𝐶 = cos 1 × 𝑆𝐴𝐷𝐶′
(tam giác nghiêng và tam giác nằm ngang có chung cạnh đáy, thế vào *)
𝑆𝐴𝐷𝐵 =
cos 1
× 𝑆𝐴𝐷𝐵′
(tam giác nghiêng và tam giác nằm ngang có chung cạnh đáy, thế vào *)
𝑆𝐴𝐵𝐶 𝑆𝐴 𝐵 ′ 𝐶′ = cos 1 Tức là 𝑆𝑛𝑔 ℎ𝑖ê𝑛𝑔 = cos 1 × 𝑆𝑝ℎẳ𝑛𝑔 (2.3)
Nhƣ vậy trong thực tế, diện tích thực của thửa đất (giả sử là mặt phẳng) nằm nghiêng ln lớn hơn diện tích đo đƣợc trên mặt phẳng bản đồ (thực chất là hình chiếu) với độ lớn 1/cosα lần, với α là góc dốc của thửa đất. Sự chênh lệch đó lớn hay nhỏ cịn phụ thuộc vào bề mặt của địa hình.
c.Tính diện tích đa giác nằm trên mặt phẳng nghiêng
Xét bài tốn tính diện tích của một đa giác phẳng nằm nghiêng. Diện tích đa giác sẽ bằng tổng diện tích các tam giác cấu thành nên đa giác đó.
𝑆đ𝑎 𝑔𝑖á𝑐 = 𝑆𝑡𝑎𝑚 𝑔𝑖á𝑐 = 𝑆𝑡𝑎𝑚 𝑔𝑖á𝑐𝑖 𝑛 𝑖=1 (2.4) Mặt khác theo cơng thức 2.3 ta có: 𝑆𝑡ℎự𝑐 = 1 cos 𝛼 𝑆𝑐ℎ𝑖ế𝑢
Với Sthực là diện tích thực của tam giác, Schiếu là diện tích của tam giác là hình
chiếu của tam giác cần tính trên mặt phẳng nằm ngang, α là góc nghiêng của tam giác. Áp dụng vào công thức 2.4, ta có: Sđa giác = 1 𝑐𝑜𝑠𝛼 𝑛 𝑖=1 𝑆𝑐ℎ𝑖ế𝑢𝑖 (2.5)
Với α là góc nghiêng của đa giác, n là số tam giác trong đa giác.
2.3.2. Quy trình tính tốn diện tích thửa đất bằng GIS
Với cách đo đạc thành lập bản đồ địa chính theo phƣơng pháp truyền thống, các thửa đất trên thực tế đƣợc đo vẽ dƣới dạng hình chiếu của nó trên mặt phẳng bản đồ. Vì vậy các kết quả này có độ chính xác cao đối với các khu vực đồng bằng địa hình bằng phẳng. Tuy nhiên, tại các khu vực đất đai có bề mặt địa hình dốc và gồ ghề thì sẽ dẫn đến sai số về diện tích và sai số này lớn hay nhỏ còn phụ thuộc vào độ chênh cao của địa hình. Ngày nay, với sự bùng nổ của công nghệ thông tin và khả năng ứng dụng rộng rãi của nó trong cơng tác quản lý đất đai thì việc tính tốn diện tích thực của các thửa đất đã đơn giản và thuận tiện hơn.
Quy trình tính tốn diện tích thực của thửa đất bằng GIS thể hiện nhƣ hình 2.5.
a. Xác định mục tiêu, nhiệm vụ
Khi tiến hành xây dựng bài tốn xác định diện tích thực cho một khu vực, cần xác định rõ mục đích của bài tốn, tiến hành điều tra tìm hiểu rõ khu vực từ đó
đƣa ra yêu cầu nhiệm vụ và xác định phạm vi thực hiện các nghiên cứu theo mục đích của bài tốn và các nhiệm vụ đề ra.
Hình 2.5: Quy trình tính tốn diện tích thực của thửa đất b. Điều tra thu thập thông tin và xử lý dữ liệu b. Điều tra thu thập thông tin và xử lý dữ liệu
Với đối tƣợng nghiên cứu là các thửa đất và dạng địa hình đối với từng thửa đất riêng biệt, theo yêu cầu tính tốn, dữ liệu đầu vào rất quan trọng là lớp điểm độ cao hoặc lớp đƣờng bình độ. Các dữ liệu này có thể thu thập trực tiếp từ bản đồ địa hình hoặc bản đồ hiện trạng sử dụng đất.
c. Tạo mơ hình số độ cao
Mơ hình số độ cao DEM - Digital Elevation Model đƣợc hiểu là một tập hợp Xác định mục tiêu, nhiệm vụ
Điều tra, thu thập thông tin và xử lý dữ liệu
Tạo mơ hình số độ cao (DEM)
Xây dựng lớp dữ liệu độ dốc
Xây dựng lớp dữ liệu ∑1/cosαi
Dữ liệu diện tích phẳng của thửa đất
Dữ liệu diện tích thực của thửa đất
Xây dựng lớp dữ liệu α
các điểm độ cao đƣợc lấy mẫu theo phƣơng pháp số để mơ tả bề mặt địa hình mà ở đó với mỗi điểm độ cao đƣợc lấy mẫu thì chỉ có duy nhất một tọa độ mặt phẳng tƣơng ứng. DEM thể hiện bề mặt địa hình trực quan dƣới dạng raster và cho phép tính tốn độ dốc của khu vực, một yếu tố ảnh hƣởng lớn đến diện tích của thửa đất. DEM có thể đƣợc tạo từ dữ liệu điểm độ cao hoặc đƣờng bình độ thơng qua một trong số các phƣơng pháp nội suy cơ bản nhƣ nội suy Spline, nội suy trị trung bình trọng số IDW hay nội suy Kriging.
- Phép nội suy Spline: là một hàm đa thức liên tục theo đoạn, nó có ƣu điểm là tính tốn nội suy nhanh chóng và có khả năng biểu diễn đƣợc các điểm đặc trƣng của khu vực nhỏ. Phƣơng pháp này dựa trên ba hay nhiều điểm đã biết.
- Phép nội suy trị trung bình trọng số: là phƣơng pháp xác định các giá trị chƣa biết bằng cách tính trị trung bình trọng số khoảng cách các giá trị trong vùng lân cận của mỗi pixel đã biết giá trị. Độ chính xác của phép nội suy này phụ thuộc vào các giá trị xung quanh nó và khơng phụ thuộc vào địa hình.
- Phép nội suy Kriging: là một quá trình xử lý chuyên sâu. Các tốc độ thực hiện phụ thuộc vào số điểm trong các tập dữ liệu đầu vào và kích thƣớc của cửa sổ tìm kiếm. Đây là phƣơng pháp nội suy dựa trên sự kết hợp giữa phƣơng pháp “phân tích bề mặt” và “trung bình trọng số”. “Phân tích bề mặt” tìm ra một phƣơng trình tốn học diễn tả xu hƣớng tổng quát của bề mặt nhƣng khơng tính đến các tính bất quy luật cục bộ. Nội suy cục bộ dùng để tính độ lệch từ xu hƣớng tồn cầu do sự không theo quy luật của khu vực. Phƣơng pháp “trung bình trọng số” dùng để tính sự biến thiên này. Trọng số đƣợc xác định bởi xu hƣớng của những độ lệch bề mặt toàn cầu và các điểm quan sát.
Chất lƣợng của phép nội suy phụ thuộc vào số lƣợng và sự phân bố của các điểm đã biết, vào độ chính xác của các giá trị dữ liệu ghi nhận và hàm toán học đƣợc chọn. Khi tiến hành nội suy ta sử dụng các phƣơng pháp khác nhau, và so sánh để xác định sự chênh lệch giữa các phƣơng pháp và sử dụng phƣơng pháp phù hợp nhất cho khu vực.
d. Xây dựng lớp dữ liệu độ dốc
đƣợc độ dốc tƣơng ứng đối với từng chi tiết trên lớp độ cao đó. Với cơng cụ Slope trong 3D Analyst, ta có thể dễ dàng thực hiện đƣợc yêu cầu này. Kết thúc q trình này ta có đƣợc các giá trị độ dốc trong từng pixel.
Nhìn chung các thơng số địa hình thƣờng đƣợc tính tốn thơng qua hai thơng số chính là độ dốc và hƣớng dốc. Đặc biệt trong công tác nghiên cứu, phân vùng cảnh quan, xác định sự phân hóa lãnh thổ theo các quy luật địa lý thì độ dốc là một căn cứ cơ bản và quan trọng. Các cấp độ dốc trên bản đồ độ dốc đƣợc phân chia và quy định màu sắc rõ ràng, trực quan giúp cho ngƣời sử dụng có thể dễ dàng thấy đƣợc sự biến thiên độ cao của địa hình.
e. Xây dựng lớp dữ liệu α, ∑1/cosαi
Sau khi thu đƣợc lớp dữ liệu độ dốc, các thông tin về độ dốc cần tiếp tục đƣợc thực hiện việc tính tốn giá trị góc nghiêng α của địa hình đƣợc biểu thị trong từng pixel riêng biệt, từ các giá trị góc α ta sẽ tính đƣợc giá trị 1/cosα tƣơng ứng.
Một trong các tính năng ƣu việt của GIS là giúp ngƣời sử dụng có thể nhanh chóng chuyển đổi định dạng của dữ liệu. Từ dữ liệu độ dốc hiển thị dƣới dạng độ, nhờ chức năng Raster Calculator trong cơng cụ 3D Analyst, ta có đƣợc dữ liệu dạng radian một cách dễ dàng. Cũng nhờ chức năng này đề tài tính tốn đƣợc giá trị
1/cosαi cho mỗi pixel.
f. Xây dựng lớp dữ liệu diện tích phẳng
Trên bản đồ số, mỗi thửa đất đƣợc biểu diễn bởi tập hợp các pixel. Vì vậy diện tích từng thửa đất sẽ đƣợc tính dựa vào kích thƣớc và số lƣợng của các pixcel. Thông qua công cụ Zonal Statistics as Table, ta sẽ thực hiện đƣợc việc thống kê các pixel và cho kết quả tính tổng các pixel đó với từng thửa đất. Tuy nhiên diện tích đƣợc tính tốn này sẽ bao gồm cả các pixel có phần nằm trong và ngồi thửa đất, dẫn đến các kết quả tính tốn cho độ chính xác chƣa cao.
g. Lớp diện tích thực của thửa đất
Ta thấy rằng, tỉ lệ giữa diện tích thực tế và diện tích đƣợc tính tốn trên bản đồ số của một thửa đất có dạng hình học trên mặt phẳng (ví dụ là hình chiếu của thửa đất nằm nghiêng một góc nào đó) ln bằng tỉ lệ:
𝐷𝑇𝑑𝑜𝑐𝑡ℎ𝑢𝑐 = 𝐷𝑇𝑑𝑜𝑐𝑝𝑖𝑥𝑒𝑙x𝐷𝑇𝑝ℎ𝑎𝑛𝑔 𝑡ℎ𝑢𝑐
𝐷𝑇𝑝ℎ𝑎𝑛𝑔𝑝𝑖𝑥𝑒𝑙 (2.6)
Áp dụng cơng thức này trong tính tốn ta có đƣợc diện tích thửa đất chỉ bao gồm các pixel nằm trọn trong các thửa đất đó. Sử dụng cơng cụ Field Calculator với phép tốn thích hợp để có đƣợc lớp diện tích thực của các thửa đất.
2.4. Tính tốn thửa đất có bề mặt mẫu
2.4.1. Tính tốn diện tích mặt cong trên hình cầu
Để kiểm chứng phƣơng pháp tính tốn diện tích thực của thửa đất theo lý thuyết, đề tài tiến hành thử nghiệm tính tốn trên mơ hình mơ phỏng của thửa đất có bề mặt địa hình khơng bằng phẳng là một thửa đất có bề mặt dạng hình cầu.
Xét mặt cầu cụt đƣợc giới hạn bởi hai góc α và β nhƣ hình 2.6 [18]:
Hình 2.6: Mặt cầu được giới hạn bởi hai góc α và β
Mặt cong ta cần tính diện tích đƣợc giới hạn dƣới bởi góc β và giới hạn trên bởi góc (α+β). Giả sử góc α đủ nhỏ để dây cung thuộc góc α là đƣờng thẳng, khi đó bài tốn tính diện tích mặt cong sẽ trở thành bài tốn tính diện tích xung quanh của hình nón cụt (hình 2.7).
Với r1 và r2 lần lƣợt là bán kính của hai đáy, l là độ dài đƣờng sinh.
Hình 2.7: Hình nón cụt tách ra từ hình cầu khi góc α đủ nhỏ [18]
Độ dài đƣờng sinh l sẽ đƣợc tính theo:
𝑙 = 𝑟1− 𝑟2 2 + ℎ2
Nếu góc α là đủ nhỏ, ta có thể tính đƣợc x theo định lý Sin trong tam giác:
𝑥 = 𝑅𝛼 = 2𝑅𝛼 2 = 2𝑅 sin 𝛼 2 Khi đó ta cũng có: 𝑟1 = 𝑅 cos 𝛽 và 𝑟2 = 𝑅𝑐𝑜𝑠 𝛽 + 𝛼
Nhƣ vậy, diện tích mặt cầu giới hạn bởi 2 góc β và (β+α) khi đó sẽ trở thành:
𝑆𝑥𝑞 = 𝜋𝑅 cos 𝛽 + cos 𝛽 + 𝛼 𝑅 cos 𝛽 − cos 𝛽 + 𝛼 2 + 𝑥2 (2.8)
x ≈ h do α rất nhỏ. Do đó:
𝑆𝑥𝑞 = 𝜋𝑅 cos 𝛽 + cos 𝛽 + 𝛼 𝑅 cos 𝛽 − cos 𝛽 + 𝛼 2 + 2𝑅 sin𝛼 2
2
𝑆𝑥𝑞 = 𝜋𝑅2 cos 𝛽 + cos 𝛽 + 𝛼 cos 𝛽 − cos 𝛽 + 𝛼 2 + 2 sin𝛼 2
2
(2.9) Áp dụng phần trên vào tính phần diện tích S của hình cầu này ta chia góc α làm k góc bằng nhau, k góc này xác định trên mặt cầu (k+1) đƣờng trịn có bán kính lần lƣợt R0, R1,…Rk nhƣ trên hình 2.8. Tiếp theo tính tổng diện tích xung quanh của
thì tổng diện tích này sẽ đạt tới diện tích của phần mặt cầu đang cần tìm.
Hình 2.8: Mặt cầu được giới hạn bởi hai góc α và β được chia thành k góc bằng nhau [18]
Diện tích xung quanh của hình nón cụt đƣợc biểu diễn theo cơng thức:
𝑆𝑥𝑞 =1 2𝑥2𝜋 𝑅1 + 𝑅2 = 𝜋𝑥 𝑅1 + 𝑅2 (2.10) Do góc 𝛼 2𝑘 rất nhỏ nên ta có đƣợc: 𝑥 = 2𝑅𝑠𝑖𝑛 𝛼 2𝑘 Và 𝑅𝑖 = 𝑅𝑐𝑜𝑠 𝛽 + 𝑖𝛼 𝑘 với 𝑖 = 0, 𝑘
Khi đó diện tích mặt cầu giới hạn giữa hai đƣờng trịn giới hạn bởi góc α và
β sẽ trở thành dạng khai triển chuỗi:
𝑆 = lim 𝑘→∞ 𝜋𝑥 𝑅𝑚 + 𝑅𝑚 −1 𝑘 𝑚 =1 Mà ta có đƣợc:
𝜋𝑥 𝑅𝑚 + 𝑅𝑚 −1 𝑘 𝑚 =1 = 2𝜋𝑅2𝑠𝑖𝑛 𝛼 2𝑘 𝑐𝑜𝑠𝛽 + 2 𝑐𝑜𝑠 𝛽 +𝑚𝑎 𝑘 + 𝑐𝑜𝑠 𝛼 + 𝛽 𝑘−1 𝑚 =1 Suy ra: 𝑆 = lim 𝑘→∞ 2𝜋𝑅2𝑠𝑖𝑛 𝛼 2𝑘 𝑐𝑜𝑠𝛽 + 2 𝑐𝑜𝑠 𝛽 +𝑚𝑎 𝑘 + 𝑐𝑜𝑠 𝛼 + 𝛽 𝑘−1 𝑚 =1 𝑆 = 2𝜋𝑅2𝑐𝑜𝑠 𝛽 +𝛼 2 𝑠𝑖𝑛 𝛼 2 = 2𝜋𝑅2 𝑠𝑖𝑛 𝛽 + 𝛼 − 𝑠𝑖𝑛𝛽 (2.11) 2.4.2. Kết quả thử nghiệm tính tốn
Đề tài thử nghiệm tính tốn diện tích có bề mặt dạng hình cầu thơng qua hai phƣơng pháp khác nhau là tính tốn bằng cơng thức tốn học đã đƣa ra ở trên và phƣơng pháp sử dụng phần mềm GIS theo quy trình đã đề xuất.
Với phần mềm GIS, ta thực hiện theo các bƣớc nhƣ sau:
1. Từ ArcGIS tạo một vịng trịn có bán kính R=500m với tâm là một điểm có tọa độ xác định (Xo=500000; Yo=2323000). Diện tích của hình trịn là 78.54 ha.
2. Từ ArcGIS sử dụng công cụ Create Random Point tạo một số lƣợng lớn các điểm nằm trong vịng trịn đó, cụ thể đề tài sử dụng 5000 điểm.
3. Với cơng cụ Add XY data, thêm thuộc tính tọa độ X,Y cho các điểm vừa tạo đƣợc.
4. Từ bảng thuộc tính của lớp điểm trên, tạo một trƣờng mới có tên DoCao
với định dạng Double, sử dụng cơng cụ Field Calculator ta có độ cao của 5000 điểm vừa tạo. Cơng thức tính là ℎ = 𝑅2 − 𝑥 − 𝑥𝑜 2 − 𝑦 − 𝑦𝑜 2 với 𝑥𝑜, 𝑦𝑜 là tâm của vòng trịn.
5. Kết thúc q trình này các điểm sẽ tạo thành một nửa hình cầu có bán kính R=500m với 5000 điểm độ cao nhƣ hình 2.9.
Hình 2.9: Mơ phỏng thửa đất dạng hình cầu
Áp dụng quy trình tính tốn theo hình 2.5 và kết quả thu đƣợc là diện tích hình cầu thể hiện ở hình 2.10.
Hình 2.10: Kết quả tính diện tích theo GIS
Cũng với thửa đất dạng hình cầu nhƣ trên, đề tài thực hiện tính tốn theo cơng thức (2.11) với góc α=90o và góc β=0o
nhƣ sau: S= 2𝜋𝑅2(sin(α+β)-sinβ) =2.π.5002 = 1570796 m2
Diện tích mặt cầu khi tính tốn bằng GIS cho kết quả S~142.1ha, tuy nhiên kết quả tính bằng cơng thức tốn học cho kết quả S~157 ha. Nhƣ vậy có thể thấy sự sai lệch diện tích giữa hai cách tính là khá lớn, mặc dù đã tốt hơn nhiều so với diện
tích hình chiếu của nó trên mặt phẳng (diện tích của hình trịn 78.54 ha nhƣ đã nói ở trên).
Để đánh giá sự thay đổi diện tích thực của thửa đất khi thay đổi số lƣợng điểm độ cao, đề tài thực hiện thay đổi số lƣợng điểm trong mặt cầu từ 5000 điểm đến 120.000 điểm. Sử dụng quy trình tính tốn diện tích thực bằng GIS sẽ có đƣợc các số liệu diện tích thay đổi theo số lƣợng các điểm độ cao và đƣợc biểu diễn bằng biểu đồ thống kê nhƣ trên hình 2.11.
Hình 2.11: Biểu đồ so sánh diện tích thực theo số lượng điểm tính tốn theo hai phương pháp
Nhìn vào biểu đồ trên hình 2.11 có thể nhận thấy: diện tích của mặt cầu biến thiên phụ thuộc vào số lƣợng điểm độ cao nằm trong mặt cầu đó. Với mặt cầu này khi tính tốn bằng cơng thức tốn học, ta có thể thấy diện tích tính đƣợc khơng phụ thuộc vào số lƣợng điểm nằm trên nó mà chỉ phụ thuộc vào bán kính R của đƣờng tròn. 130 135 140 145 150 155 160 5 10 20 30 40 50 60 70 80 90 100 110 120 DIỆN TÍCH THỰC TÍNH BẰNG CƠNG CỤ GIS
Số lượng điểm độ cao (đv: nghìn điểm)
Diện tích