Làm thế nào tôi có thể đo diện tích từ tọa độ địa lý?
Câu trả lời:
PostGIS 1.5 đã giới thiệu một loại GEOGRAPHY mới . Các GEOGRAPHYloại cho phép tọa unprojected trên một hình cầu phải được lưu trữ trong một bảng PostGIS, và một số chức năng phân tích phải thực hiện khi họ.
ST_Area truy vấn có thể được thực hiện theo đa giác loại GEOGRAPHY để tính diện tích của chúng theo mét vuông.
Truy vấn sau đây đưa ra diện tích của tất cả các đa giác bằng cách sử dụng hình cầu (hiện tại chỉ WGS-84hỗ trợ hình cầu), giả sử chúng được lưu trữ bằng cách sử dụng GEOGRAPHYloại:
SELECT ST_Area(the_geom) FROM table_of_polygons;
Thuật toán được sử dụng để tính diện tích trên một hình cầu có thể được lấy từ mã nguồn .
Đây là một liên kết đến một số mã sẽ mang lại diện tích của một đa giác đơn giản (ban đầu từ Diễn đàn Gió thế giới): http://forum.worldwindcentral.com/showthread.php?t=20724 . Điều này giải quyết vấn đề trên một hình cầu, đại khái dựa trên mối quan hệ:
S = diện tích đa giác; theta là tổng các góc bên trong tính bằng radian; n là số đỉnh; r là bán kính của mặt cầu.
Xem thêm (nguồn hình ảnh công thức): http://www.geom.uiuc.edu/docs/reference/CRC-form Formula / node59.html
Tôi sẽ rất vui mừng khi thấy các liên kết và / hoặc mã cho khu vực đa giác trên một hình cầu bắt buộc.
Đây là nguồn cho phép tính đơn giản mà chúng tôi thực hiện trong OpenLayers. Phương pháp này xuất phát từ "Một số thuật toán cho đa giác trên một quả cầu" (Robert. G. Chamberlain và William H. Duquette, NASA JPL Publication 07-03). Mã được liên kết ở trên là để xác định diện tích của vòng tuyến tính (có tọa độ địa lý). Các khu vực cho Đa giác và MultiPolygons được tóm tắt từ các vòng.
var area = 0.0;
var len = ring.components && ring.components.length;
if (len > 2) {
var p1, p2;
for (var i=0; i<len-1; i++) {
p1 = ring.components[i];
p2 = ring.components[i+1];
area += OpenLayers.Util.rad(p2.x - p1.x) *
(2 + Math.sin(OpenLayers.Util.rad(p1.y)) +
Math.sin(OpenLayers.Util.rad(p2.y)));
}
area = area * 6378137.0 * 6378137.0 / 2.0;
}
Các thành phần vòng là hai mảng phần tử của các chuỗi x, y (lon, lat) trong đoạn mã trên. Phương thức OpenLayers.Util.rad chỉ chuyển đổi độ thành radian (deg * PI / 180).
Bạn sẽ cần chuyển đổi tọa độ địa lý của mình thành hình chiếu có hệ thống tọa độ cho phép bạn sử dụng toán học Cartesian để tính diện tích.
Tôi tin rằng UTM là phép chiếu tiêu chuẩn được chấp nhận, vì rất đơn giản để chọn một vùng dựa trên vĩ độ và kinh độ của bạn, và độ méo là tối thiểu, thậm chí trên các vùng. Vì vậy, nếu bạn có đa giác kích thước của Texas, bạn có thể sử dụng UTM Vùng 14 N và nó vẫn sẽ khá chính xác.
Nếu đa giác của bạn ở cực Bắc hoặc Nam, thì bạn nên sử dụng UPS thay vào đó, vì các phép chiếu UTM kém chính xác hơn các cực và bạn sẽ nhanh chóng đi qua chúng khi các ranh giới trở nên nhỏ hơn (vì chúng đi theo các đường kinh độ)
Khi các điểm của bạn nằm trong hệ tọa độ thân thiện với Cartesian, bạn có thể coi chúng như các đa giác trên lưới và tính diện tích.
Một lớp PolygonArea đã được thêm vào GeographicLib vào năm 2011-2007. Điều này sẽ tính diện tích hình elip thực sự của một đa giác có các cạnh là trắc địa. Không giống như PostGIS, phương thức này không đòi hỏi tích hợp số. Để biết thông tin (và liên kết đến bài báo có công thức xuất phát), xem
http://geographiclib.sf.net/html/ classGeographicLib_1_1PolygonAreaT.html
(Liên kết được cố định để phản ánh tổng quát hóa của PolygonArea đến một lớp mẫu.)