Jak analizować ramkę danych XML do R.


103

Próbowałem przeanalizować ramkę danych XML do R, ten link bardzo mi pomógł:

jak utworzyć ramkę danych R z pliku xml

Ale nadal nie byłem w stanie rozwiązać mojego problemu:

Oto mój kod:

data <- xmlParse("http://forecast.weather.gov/MapClick.php?lat=29.803&lon=-82.411&FcstType=digitalDWML")
xmlToDataFrame(nodes=getNodeSet(data1,"//data"))[c("location","time-layout")]
step1 <- xmlToDataFrame(nodes=getNodeSet(data1,"//location/point"))[c("latitude","longitude")]
step2 <- xmlToDataFrame(nodes=getNodeSet(data1,"//time-layout/start-valid-time"))
step3 <- xmlToDataFrame(nodes=getNodeSet(data1,"//parameters/temperature"))[c("type="hourly"")]

Ramka danych, którą chcę mieć, wygląda następująco:

latitude  longitude   start-valid-time   hourly_temperature
29.803     -82.411  2013-06-19T15:00:00-04:00    91
29.803     -82.411  2013-06-19T16:00:00-04:00    90

Utknąłem na tym xmlToDataFrame(), Każda pomoc byłaby bardzo mile widziana, dzięki.

Odpowiedzi:


103

Dane w formacie XML rzadko są zorganizowane w sposób umożliwiający działanie xmlToDataFramefunkcji. Lepiej jest wyodrębnić wszystko z list, a następnie powiązać listy w ramce danych:

require(XML)
data <- xmlParse("http://forecast.weather.gov/MapClick.php?lat=29.803&lon=-82.411&FcstType=digitalDWML")

xml_data <- xmlToList(data)

W przypadku danych przykładowych uzyskanie lokalizacji i czasu rozpoczęcia jest dość proste:

location <- as.list(xml_data[["data"]][["location"]][["point"]])

start_time <- unlist(xml_data[["data"]][["time-layout"]][
    names(xml_data[["data"]][["time-layout"]]) == "start-valid-time"])

Dane dotyczące temperatury są nieco bardziej skomplikowane. Najpierw musisz dostać się do węzła zawierającego listy temperatur. Następnie musisz wyodrębnić obie listy, przejrzeć każdą z nich i wybrać tę, która ma „godzinę” jako jedną z wartości. Następnie musisz wybrać tylko tę listę, ale zachować tylko wartości z etykietą „wartość”:

temps <- xml_data[["data"]][["parameters"]]
temps <- temps[names(temps) == "temperature"]
temps <- temps[sapply(temps, function(x) any(unlist(x) == "hourly"))]
temps <- unlist(temps[[1]][sapply(temps, names) == "value"])

out <- data.frame(
  as.list(location),
  "start_valid_time" = start_time,
  "hourly_temperature" = temps)

head(out)
  latitude longitude          start_valid_time hourly_temperature
1    29.81    -82.42 2013-06-19T16:00:00-04:00                 91
2    29.81    -82.42 2013-06-19T17:00:00-04:00                 90
3    29.81    -82.42 2013-06-19T18:00:00-04:00                 89
4    29.81    -82.42 2013-06-19T19:00:00-04:00                 85
5    29.81    -82.42 2013-06-19T20:00:00-04:00                 83
6    29.81    -82.42 2013-06-19T21:00:00-04:00                 80

94

Użyj xpath bardziej bezpośrednio, aby uzyskać wydajność i przejrzystość.

time_path <- "//start-valid-time"
temp_path <- "//temperature[@type='hourly']/value"

df <- data.frame(
    latitude=data[["number(//point/@latitude)"]],
    longitude=data[["number(//point/@longitude)"]],
    start_valid_time=sapply(data[time_path], xmlValue),
    hourly_temperature=as.integer(sapply(data[temp_path], as, "integer"))

prowadzący do

> head(df, 2)
  latitude longitude          start_valid_time hourly_temperature
1    29.81    -82.42 2014-02-14T18:00:00-05:00                 60
2    29.81    -82.42 2014-02-14T19:00:00-05:00                 55

12
To naprawdę powinna być akceptowana odpowiedź. Jest bardziej zwięzły, a xpath ma znacznie lepszą wydajność niż iteracja po listach.
SchaunW

40

Oto częściowe rozwiązanie przy użyciu xml2. Podzielenie roztworu na mniejsze części ogólnie ułatwia upewnienie się, że wszystko jest w jednej linii:

library(xml2)
data <- read_xml("http://forecast.weather.gov/MapClick.php?lat=29.803&lon=-82.411&FcstType=digitalDWML")

# Point locations
point <- data %>% xml_find_all("//point")
point %>% xml_attr("latitude") %>% as.numeric()
point %>% xml_attr("longitude") %>% as.numeric()

# Start time
data %>% 
  xml_find_all("//start-valid-time") %>% 
  xml_text()

# Temperature
data %>% 
  xml_find_all("//temperature[@type='hourly']/value") %>% 
  xml_text() %>% 
  as.integer()

8
Pomocna odpowiedź. Jeśli ktoś inny się na to natknie, oto link do samouczka Hadleya na temat korzystania z xml2: blog.rstudio.com/2015/04/21/xml2
Richard Erickson

9

Możesz wypróbować poniższy kod:

# Load the packages required to read XML files.
library("XML")
library("methods")

# Convert the input xml file to a data frame.
xmldataframe <- xmlToDataFrame("input.xml")
print(xmldataframe)
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.