2017-06-01 2 views
1

Ich benutze R, um in eine Formdatei zu lesen, um mit den sp-Paket-Polygonen von Ölfeldern (Länge Breite mit WGS84) und ihren jeweiligen Gebieten zu analysieren. Leider kenne ich die Einheit der Flächenleistung nicht. Z.B. die Flächenleistung beträgt im Durchschnitt 0,85 Einheiten (max. 4,34), was wahrscheinlich nicht in Quadratkilometern liegt, da dies für Ölfelder viel zu klein wäre.R sp: Einheit der Fläche von Polygon

Kennt jemand die Einheit der Flächenausgabe von Polygonen im SP-Paket? Danke vielmals!

Antwort

2

Um eine korrekte Flächenberechnung für ein Polygon in lat-lon-Koordinaten zu erhalten, wäre es besser, sie zuvor mit "spTransform" in eine metrische Projektion mit gleichem Flächeninhalt zu konvertieren. Alternativ können Sie Paket geosphere verwenden, die

tun können „Sphärische Trigonometrie für geographische Anwendungen. Das heißt, Entfernungen berechnen und lated Maßnahmen re- für Winkel (Länge/Breite) locations“

Zum Beispiel diese:

require(geosphere) 
areaPolygon(mypoly) 

(mit mypoly ein spatialPolygons Objekt ist) werden Sie in Quadratkilometern die Umgebung geben.

HTH.

+0

Vielen Dank, LoBu! – schimo

1

Wenn sf verwenden, können Sie tatsächlich die resultierenden Einheiten erhalten:

> library(sf) 
Linking to GEOS 3.5.1, GDAL 2.1.3, proj.4 4.9.2, lwgeom 2.3.2 r15302 
> demo(nc, ask = FALSE, echo = FALSE) 
Reading layer `nc.gpkg' from data source `/home/edzer/R/x86_64-pc-linux-gnu-library/3.4/sf/gpkg/nc.gpkg' using driver `GPKG' 
converted into: MULTIPOLYGON 
Simple feature collection with 100 features and 14 fields 
Attribute-geometry relationship: 0 constant, 8 aggregate, 6 identity 
geometry type: MULTIPOLYGON 
dimension:  XY 
bbox:   xmin: -84.32385 ymin: 33.88199 xmax: -75.45698 ymax: 36.58965 
epsg (SRID): 4267 
proj4string: +proj=longlat +datum=NAD27 +no_defs 
> st_area(nc[1:2,]) 
Units: m^2 
[1] 1137388604 611077263 
> units::set_units(st_area(nc[1:2,]), km^2) 
Units: km^2 
[1] 1137.3886 611.0773 

Für longlat Daten, wie hier verwendet es geosphere hinter den Kulissen.

Verwandte Themen