ฉันมีไฟล์ในการฉายภาพ British National Grid:
Geometry: 3D Polygon
Feature Count: 5378
Extent: (9247.520209, 14785.170099) - (638149.173223, 1217788.569952)
Layer SRS WKT:
PROJCS["British_National_Grid",
GEOGCS["GCS_airy",
DATUM["OSGB_1936",
SPHEROID["Airy_1830",6377563.396,299.3249646]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",49],
PARAMETER["central_meridian",-2],
PARAMETER["scale_factor",0.9996012717],
PARAMETER["false_easting",400000],
PARAMETER["false_northing",-100000],
UNIT["Meter",1]]
cat: Integer (9.0)
ฉันสามารถใช้ GDAL / OGR เพื่อให้ได้พื้นที่รวมของรูปหลายเหลี่ยมทั้งหมดใน shapefile เป็นเฮกตาร์หรือไม่
ฉันสงสัยว่ามันเป็นไปได้ด้วย-sql
หรือเปล่า:
ogrinfo -sql "SELECT SUM(ST_Area(geom::geography)) FROM mytable" myshapefile.shp
ERROR 1: Undefined function 'ST_Area' used.
แต่พยายามที่ฉันได้รับ
ฉันเดาว่าฉันสามารถนำเข้า Shapefile ลงใน QGIS เพิ่มแอททริบิวของพื้นที่ลงในแต่ละรูปหลายเหลี่ยมแล้วรวมมัน แต่ฉันจะใช้เครื่องมือบรรทัดคำสั่งแทนถ้าเป็นไปได้
SUM_OGR_GEOM_AREA (Real) = 4459037129.50955
นี้จะช่วยให้ฉัน มันอยู่ในเฮกตาร์หรือหน่วยอื่น ๆ ไหม? และมันสำคัญว่าการฉายซอร์สโค้ดของฉันคืออะไร?