Comment créer un fichier KML à l'aide de R

J'ai écrit un script R pour obtenir un certain point de la carte de données (valeurs de Latitude et Longitude). Je suis capable de tracer dans R et de les visualiser. Mais maintenant, je veux générer un fichier KML à partir de ces données et de les afficher à l'aide de Google Earth. Afin que je puisse les partager avec des collègues et qu'ils puissent le voir sur Google Earth.

Quelle est la meilleure méthode /paquet pour ce faire ?

OriginalL'auteur Vijay Barve | 2011-10-18

  1. 18

    Vérifier la writeOGR fonction dans le rgdal paquet. Voici un exemple simple:

    library("sp")
library("rgdal")
data(meuse)
coordinates(meuse) <- c("x", "y")
proj4string(meuse) <- CRS("+init=epsg:28992")
meuse_ll <- spTransform(meuse, CRS("+proj=longlat +datum=WGS84"))
writeOGR(meuse_ll["zinc"], "meuse.kml", layer="zinc", driver="KML")

    Les objets exportés sont SpatialPointsDataFrame, SpatialLinesDataFrame, ou SpatialPolygonsDataFrame objets, tel que défini dans la sp paquet.

    R> class(meuse)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"

    Pour l'écriture avec le KML pilote, notez que les géométries doivent être coordonnées géographiques avec le datum WGS84.

    OriginalL'auteur rcs

  2. 4

    Je pense vaut la peine de mentionner le plotKML.

    Toutefois, pour faciliter le partage entre collègues que j'ai trouvé intéressant le mapview paquet basé sur dépliant paquet. On peut enregistrer une carte en tant que document HTML avec diverses options pour une carte d'arrière-plan; pas besoin de Google Earth et le HTML de la carte sur votre navigateur.

    Quelques exemples:

    library(sp)
library(rgdal)
library(raster)
library(plotKML)
library(mapview)

# A SpatialPointsDataFrame example
data(meuse)
coordinates(meuse) <- ~x+y
proj4string(meuse) <- CRS("+init=epsg:28992") # CRS Amersfoort (Netherlands)
# make a KML file from SpatialPointsDataFrame object
# will get a warning like "Reprojecting to +proj=longlat +datum=WGS84 ..."
# as it is expected to work with geographical coordinates with datum=WGS84, 
# but seems that it takes care of the reprojecting. 
plotKML::kml(meuse,
             file.name    = "meuse_cadium.kml",
             points_names = meuse$cadmium,
             colour    = "#FF0000",
             alpha     = 0.6,
             size      = 1,
             shape     = "http://maps.google.com/mapfiles/kml/pal2/icon18.png")
# Or, an easy to make interactive map with mapView()
mapView(meuse)

# A RasterLayer example   
data(meuse.grid)
gridded(meuse.grid) <- ~x+y
proj4string(meuse.grid) <- CRS("+init=epsg:28992")
dist_rst <- raster(meuse.grid["dist"])
# make a KML file from RasterLayer object
plotKML::kml(dist_rst,
             file.name    = "dist_rst.kml",
             colour_scale = SAGA_pal[[1]])
# Or, easy to make interactive map with mapView() - display raster and add the points
mapView(dist_rst, legend=TRUE) + meuse
# However, note that for bigger raster datasets mapView() might reduce from resolution

    Plus d'exemples avec plotKML ici, avec un tutoriel ici. Pour mapview, une intro peut être trouvé ici.

    OriginalL'auteur Valentin

  3. 0

    Si vous êtes prêts à sortir de R, il existe un programme gratuit appelé DNRGarmin peut prendre un fichier séparé par des virgules comme une .txt et les convertir en .kml pour l'importer dans google earth.

    Vous pouvez le trouver ici:

    http://www.dnr.state.mn.us/mis/gis/tools/arcview/extensions/DNRGarmin/DNRGarmin.html

    donc dans R:

    my.geo.data <- all.my.data[ c("unique.id", "lats", "longs")]

write.csv( my.geo.data, file = "myGeoData.txt")

    ouvrir DNRGarmin,

    Fichier -> Load from -> File -> myGeoData.txt Ensuite,
    Fichier -> Enregistrer pour -> File -> myGeoData.kml

    @rcs de conseils re: WGS84 s'applique pour cette réponse aussi.

    Bonne chance

    OriginalL'auteur Chris

  4. 0

    Si vous/votre /vos collègues savent QGIS, c'est un très bon moyen pour afficher des données dans Google Earth. QGIS a la particularité de montrer à Google Earth sur le fond de carte, puis vous pouvez ouvrir vos données spatiales et il sera affiché sur la carte de base. Bien sûr, il exige que vos données soient correctement projetées rcs dit.

    Ici, vous devez exporter vos points comme une forme de fichiers à l'aide de la maptools paquet et Spatiale des Points de paquet:

    library(maptools)
library(sp)

## define projection
myProjection <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

## your points in format dataframe
coordinates.df <- as.data.frame(MyCoordinates) 

## the number of points you have as dataframe
number <- as.data.frame(NumberOfPoints)

## convert points to Spatial Points Dataframe
myPoints.spdf <- SpatialPointsDataFrame(coordinates.df, number, proj4string = CRS(myProjection))

## save shapefile
writeSpatialShape(myPoints.spdf, "MyPointsName")

    Vos points peuvent maintenant être ouverts dans QGIS et l'afficher dans Google Earth. Dans QGIS, vos données peuvent également être enregistrés en tant que fichier kmz si nécessaire.

    OriginalL'auteur dtanon