I have downloaded a GML-File from this Website.
Then I use the following code to process this file in R
:
library(rgdal)
test <- readOGR(dsn='WLV_GZP.gml')
I am printing the map like this:
test3 <- subset(test, qualitativeValue=='Rote Zone Lawine') #For speeding up the Test-print
plot(test3)
The result is turned 90 degrees to the right and flipped (East-West). How can I get this map to work? Any help is greatly appreciated.
edit: The result is the same when using the sf
library:
test4 <- sf::st_read('WLV_GZP.gml')
test4 <- subset(test4, qualitativeValue=='Rote Zone Lawine')
plot(test4)
In order to inspect the bbox I transform the map:
roads <- st_transform(test4, crs = 4326)
This yields a bbox:
Dimension: XY
Bounding box: xmin: -23.4819 ymin: 58.59924 xmax: -18.90591 ymax: 62.56157
Geod开发者_StackOverflow中文版etic CRS: WGS 84
Since the map displays Austria this is of course way off...
PS: I would like to post a PNG-File but I received an error when uploading the file.
精彩评论