Konwersja danych punktowych na siatkową ramkę danych do analizy histogramu przy użyciu R?

Jestem bardzo nowy w korzystaniu z danych GIS i mam jedynie niewielkie doświadczenie z R. Czytałem o tym, jak analizować dane przestrzenne za pomocą książki PDF spatial-analyst.net, więc nie jestem całkowicie zagubiony, ale pomyślałem, że mógłbym opisać mój problem i ludzie mogą sugerować pomysły.

Mam zestaw danych z około 2000 pomiarami przy różnych współrzędnych szerokości / długości, chociaż prawdopodobnie podzielę ten zestaw danych, ponieważ dane były gromadzone przez 3 lata, a warunki zmieniały się z czasem. Nazwijmy mierzoną zmienną „IP”.

Chcę utworzyć mapę IP na całym badanym obszarze za pomocą Kriginga lub innej metody interpolacji na przykładowych danych. Następnie chcę utworzyć histogram mierzący wielkość terenu w różnych segmentach IP. Będę też musiał stworzyć histogram, który pokazuje liczbę próbek w każdym wiadrze (zauważ, że próbka może mieć wyższe lub niższe rzeczywiste IP niż to, co przewiduje kriging dla jego ziemi).

Śledzę, jak ładować dane do SpatialPointsDataFrame i przeprowadzam analizę kriging, w której mam problem z tym, jak przekonwertować te dane na siatkową ramkę danych, aby móc przeprowadzić analizę histogramu.

Wszelkie sugestie dotyczące konwertowania punktów na siatki?

Masz rację ... to całkiem proste! Pakiet „rastrowy” ma całkiem proste sposoby radzenia sobie z tworzeniem i manipulowaniem rastrami.

library(maptools)
library(raster)

# Load your point shapefile (with IP values in an IP field):
pts <- readShapePoints("pts.shp")

# Create a raster, give it the same extent as the points
# and define rows and columns:

rast <- raster()
extent(rast) <- extent(pts) # this might be unnecessary
ncol(rast) <- 20 # this is one way of assigning cell size / resolution
nrow(rast) <- 20

# And then ... rasterize it! This creates a grid version 
# of your points using the cells of rast, values from the IP field:
rast2 <- rasterize(pts, rast, pts$IP, fun=mean) 

Możesz przypisać rozmiar siatki i rozdzielczość na wiele sposobów - zapoznaj się z dokumentacją pakietu rastrowego.

Wartości komórek rastrowych z rasteryzacji można obliczyć za pomocą funkcji - „średnia” w powyższym przykładzie. Upewnij się, że umieściłeś to: w przeciwnym razie po prostu używa wartości IP od ostatniego punktu, w którym się znajduje!

Z pliku CSV:

pts <- read.csv("IP.csv")
coordinates(pts) <- ~lon+lat
rast <- raster(ncol = 10, nrow = 10)
extent(rast) <- extent(pts)
rasterize(pts, rast, pts$IP, fun = mean)

Pakiet plotKML ma funkcję o nazwie vect2rast. Ta funkcja zasadniczo rozszerza rasterizefunkcję dostępną w pakiecie rastrowym. Zaletą vect2rast; polega jednak na tym, że nie wymaga wprowadzania danych ze strony użytkownika, tj. automatycznie określa rozmiar komórki siatki i ramki granicznej na podstawie właściwości zestawu danych wejściowych. Rozmiar komórki siatki jest szacowany na podstawie gęstości / rozmiaru obiektów na mapie ( nndistfunkcja w pakiecie spatstat).

library(plotKML)
Rast2 <- vect2rast(pts)

# for large data sets use SAGA GIS:
Rast2 <- vect2rast(pts, method = "SAGA")