r - स्थानिक अंक की ग्रिड कैसे बनाएं



grid geospatial (1)

आप sp::makegrid का उपयोग कर एक ग्रिड बना सकते हैं

library(sp)
library(rgdal)
library(raster)

# load some spatial data. Administrative Boundary
us <- getData('GADM', country = 'US', level = 1)
us$NAME_1
colorado <- us[us$NAME_1 == "Colorado",]


# check the CRS to know which map units are used
proj4string(colorado)
# "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

# Create a grid of points within the bbox of the SpatialPolygonsDataFrame 
# colorado with decimal degrees as map units
grid <- makegrid(colorado, cellsize = 0.1)

# grid is a data.frame. To change it to a spatial data set we have to
grid <- SpatialPoints(grid, proj4string = CRS(proj4string(colorado)))

plot(colorado)
plot(grid, pch = ".", add = T)

और ऑस्ट्रिया के साथ एक और उदाहरण (जीएडीएम कोड 'AUT' )।

ग्रिड केवल बहुभुज के भीतर इंगित करता है:

बहुभुज के भीतर केवल अंक निकालने के लिए, इस तरह के स्थान के आधार पर अंकों को कम करने के लिए `[` का उपयोग करें:

grid <- grid[colorado, ]
library(reshape2)
library(data.table)
library(dplyr)
library(magrittr)
library(ggplot2)
library(scales)
library(gstat)
library(DescTools)
library(sp)



#I want a colorado grid#
data("colorado.grid")

#making cordinates into spatial points dataframe#
coordinates(Gold_tracer_kri) <- ~ long_orig + lat_orig
#attempt at kriging but no grid#
lzn.kriged <- krige(Au ~ 1, Gold_tracer_kri, colorado.grid, model=lzn.fit)

lzn.kriged %>% as.data.frame %>%
ggplot(aes(long_orig=long_orig, lat_orig=lat_orig)) +     geom_tile(aes(fill=var1.pred)) + coord_equal() +
scale_fill_gradient(low = "yellow", high="red") +
scale_x_continuous(labels=comma) + scale_y_continuous(labels=comma) +
theme_bw()

स्थानिक डोमेन लोड करना

data("meuse.grid")

मैं आर में क्रिगिंग विधियों का उपयोग करने की कोशिश कर रहा हूं लेकिन मेरे डेटा के लिए ग्रिड नहीं ढूंढने के कारण मैं फंस गया हूं। मेरे डेटा को कोलोराडो के सभी अमेरिकी राज्यों को ऊपर ले जाता है और मैं अपने डेटा को अंतराल करने के लिए एक ग्रिड प्राप्त करने में सक्षम होना चाहता हूं। Iuse.grid के समान है जिसे मैं निम्नलिखित उदाहरण में उपयोग किया जाता है I

किसी भी सहायता की सराहना की जाएगी





geospatial