2017-01-31 4 views
0

Ich habe das Raster (DEM)Wie überlagert man den Cernatin-Punktwert im Levelplot mit R?

class  : RasterLayer 
dimensions : 164, 96, 15744 (nrow, ncol, ncell) 
resolution : 0.0002777778, 0.0002777778 (x, y) 
extent  : 83.47915, 83.50582, 28.80114, 28.84669 (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : in memory 
names  : layer 
values  : -4.737531, 0.554273 (min, max) 

und ich möchte der Punkt

p = -3.9837500 -2.6327150 -1.4857878 -1.0200000 -0.7716000 

in der folgenden lat-long

lat=28.80046, 28.80381, 28.81314, 28.81314, 28.81144 
lon=83.49621, 83.49524, 83.49450, 83.49450, 83.49201 
+0

welche Art des Suchens 'sp.points' hast du getan? Sicherlich gibt es Beispiele für das Zeichnen von Punkten auf jeder grafischen Ausgabe, die aktuelle Pakete produzieren. Y. –

+0

Ich habe mit 'Plot' Befehl plotten und es funktioniert. Aber ich muss mit 'levelplot', dem Paket 'rasterVis' plotten. –

Antwort

0

eine spatialpointsdataframe Make-Werte zu überlagern;

p = -3.9837500 -2.6327150 -1.4857878 -1.0200000 -0.7716000 
lat=28.80046, 28.80381, 28.81314, 28.81314, 28.81144 
lon=83.49621, 83.49524, 83.49450, 83.49450, 83.49201 

df <- data.frame(lat = lat, lon= lon, z = p) 
coordinates(df) <- ~lon+lat 
crs(df) <- projection(yourRas) 

Und rasterVis Grundstück fügen mit

levelplot(yourRas) + 
    layer(sp.points(df, size=3)) 
Verwandte Themen