मैं इसके अंदर निहित आकार और बिंदु लेने की कोशिश कर रहा हूं और आकारफाइल, बिंदुओं और फिर एक साजिश के साथ समाप्त होता हूं अंततः अल्फा पारदर्शिता की कुछ मात्रा के साथ शीर्ष पर एक quadrat विश्लेषण overlayed।एक गर्मी की तरह शैली में ggplot2 में ggplot2 में मानचित्र के शीर्ष पर क्वाड्रेट गणना कैसे करें
मैंने कोशिश की और इस उदाहरण है कि काम करना चाहिए के साथ आया था:
library(ggplot2)
library(maps)
library(maptools)
library(spatstat)
library(plyr)
# generate polygons
usa <- map('usa', plot=F)
IDs <- sapply(strsplit(usa$names, ':'), function(x) x[1])
crs <- CRS('+proj=longlat +ellps=WGS84')
usa.sp <-map2SpatialPolygons(usa, IDs=usa$names, proj4string=crs)
usa.pts <- fortify(usa.sp, region="id")
base.plot <- ggplot(usa.pts, aes(x=long, y=lat)) +
geom_path(aes(group=id)) +
coord_equal()
# generate point process
lon <- runif(5000, min(usa.pts$long, na.rm=TRUE), max(usa.pts$long, na.rm=TRUE))
lat <- rnorm(5000, mean(usa.pts$lat, na.rm=TRUE), sd=sd(usa.pts$lat, na.rm=TRUE)/2)
points <- data.frame(lon, lat)
points.sp <- SpatialPoints(points, proj4string=crs)
points <- points[which(!is.na(over(points.sp, usa.sp))),]
# the first plot
base.plot +
geom_point(data=points, aes(x=lon, y=lat), color="red", alpha=I(0.5))
तो फिर क्या मेरी शेपफ़ाइल और डेटा मोटे तौर पर कैसा लगेगा हैं:
अगले चरणों में मैं चाहूँगा विभिन्न सेल आकारों के quadrat विश्लेषण करने के लिए लेते हैं और अंक के साथ साजिश पर परिणाम साजिश। मान लीजिए कि मैं सफलतापूर्वक एक क्वाड्रैक्टकाउंट ऑब्जेक्ट प्राप्त करता हूं (डेटा को अभी सही प्रारूप में घुमाता हूं), मैं इसे एक हीटमैप-जैसे ओवरले के रूप में प्लॉट करने के बारे में कैसे जाउंगा? अलग नक्शा अनुमानों को सक्षम करने के
library(mapproj)
base.plot +
stat_bin2d(data=points, aes(x=lon, y=lat), bins=100) +
coord_map()
मैं आपको यह भी coord_map
का उपयोग नहीं बल्कि coord_equal
से, सुझाव है mapproj
पैकेज स्थापित करने के बाद:
शायद एक संबंधित सवाल यह है कि quadratcount वर्ग की छवि भाग लेते हैं, यह किसी शेपफ़ाइल जहां प्रत्येक कोशिका उचित गिनती मूल्य के साथ एक बहुभुज है, और साजिश है कि एक geom_tile परत के रूप में करने के लिए है? –