开发者

Map is flipped and turned when reading GML data with r

开发者 https://www.devze.com 2022-12-07 19:54 出处:网络
I have downloaded a GML-File from this Website. Then I use the following code to process this file in R:

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.

0

精彩评论

暂无评论...
验证码 换一张
取 消