Hello, I have following problem: I wish to extract z coordinates from a polyline shapefile. Using readOGR() i get following:
> data <- readOGR(path,'tesikov') OGR data source with driver: ESRI Shapefile Source: "C:/Users/Janoska/Documents/GPS", layer: "tesikov" with 1 features and 1 fields Feature type: wkbLineString with 3 dimensions Warning message: In readOGR("C:/Users/Janoska/Documents/GPS", "tesikov") : Z-dimension discarded Clearly, Z dimension is stored in the data, but readOGR does not read it. Hence, I followed this post: https://stat.ethz.ch/pipermail/r-sig-geo/2007-March/001847.html data <- coordinates(readShapeLines('path/tesikov.shp')) > str(data) List of 1 $ :List of 1 ..$ : num [1:335, 1:2] -540376 -540369 -540361 -540353 -540344 ... Resulting coordinates have only two dimensions. I am wondering, whether Z dimension can be extracted only from points? Any suggestion would be appreciated. Thanks, Zbynek -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/problem-extracting-Z-value-from-line-shapefile-tp7580600.html Sent from the R-sig-geo mailing list archive at Nabble.com. _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo