मेरे पास ब्रिटिश राष्ट्रीय ग्रिड प्रक्षेपण में एक आकृति है:
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 का उपयोग सभी पॉलीगोन के कुल क्षेत्रफल को शेपफाइल में, हेक्टेयर में प्राप्त करने के लिए कर सकता हूँ?
मैं सोच रहा था कि क्या यह संभव है -sql
, कुछ इस तरह से:
ogrinfo -sql "SELECT SUM(ST_Area(geom::geography)) FROM mytable" myshapefile.shp
लेकिन कोशिश है कि मुझे मिल जाए ERROR 1: Undefined function 'ST_Area' used.
।
मुझे लगता है कि मैं क्यूजीआईएस में शेपफाइल को आयात कर सकता हूं, प्रत्येक बहुभुज के लिए एक क्षेत्र विशेषता जोड़ सकता हूं, और फिर इसे योग कर सकता हूं, लेकिन यदि संभव हो तो मैं कमांड लाइन टूल का उपयोग करूंगा।
SUM_OGR_GEOM_AREA (Real) = 4459037129.50955
। क्या यह हेक्टेयर या किसी अन्य इकाई में है? और क्या इससे कोई फर्क पड़ता है कि मेरे स्रोत शेपफाइल में क्या प्रक्षेपण है?