Stat 133 Homework 10

In [2]:

    

library(DataComputing)
library(XML)


    

In [3]:

    

data = read.csv(file="earthquakes.csv")%>%
  filter(Magnitude>=4)
long = data$Longitude
lat = data$Latitude
head(data)


    

    
Out[3]:





DateTime
Latitude
Longitude
Depth
Magnitude
MagType
NbStations
Gap
Distance
RMS
Source
EventID

	
1
1968/03/21 21:54:59.94
37.0378
-121.7407
11.86
4.3
ML
29
133
6
0.1
NCSN
1001511
	
2
1968/05/30 08:03:01.06
38.1722
-123.2057
5
4.2
ML
38
324
91
0.3
NCSN
1001643
	
3
1969/09/25 13:21:15.06
35.8658
-120.841
6.15
4.09
Md
8
223
34
0.05
NCSN
1003106
	
4
1969/09/30 03:33:05.84
36.8468
-121.5635
5.49
4
Md
13
73
5
0.09
NCSN
1003117
	
5
1969/10/02 04:56:45.30
38.4978
-122.664
0.22
5.6
ML
38
104
52
0.22
NCSN
1003129
	
6
1969/10/02 06:19:56.39
38.45
-122.7535
5.14
5.7
ML
53
139
58
0.22
NCSN
1003132

In [4]:

    

Doc = newXMLDoc()
Root = newXMLNode("kml",namespaceDefinitions = "http://www.opengis.net/kml/2.2", doc = Doc)
Docmt = newXMLNode("Document", parent = Root)
Name = newXMLNode("Name", "Earthquakes", parent = Docmt)
Description = newXMLNode("Description", "4+ Earthquakes, 1966-present", parent = Docmt)


    

In [5]:

    

n = nrow(data)
Dtime = as.character(data$DateTime)
Dtimefix = gsub("/", "-", Dtime)
time = gsub(" ","T",Dtimefix)

for (i in 1:n)
{
  M = newXMLNode("Placemark", parent = Docmt)
  T = newXMLNode("Point", parent = M)
  cood = c(long[i],",",lat[i])
  newXMLNode("coordinates", cood, parent = T)
  S = newXMLNode("TimeStamp",parent=M)
  newXMLNode("when", time[i],"+08:00", parent=S)
}


saveXML(Doc, "earthquakes6.kml")


    

    
Out[5]:




'earthquakes6.kml'