使用ppp函数绘图错误

Question

我在策划德克萨斯州休斯顿炼油厂的位置。我用以下代码定义了学习区域。

 P4S.latlon <- CRS("+proj=longlat +datum=WGS84 ")
 county.lines <-readShapePoly("tl_2017_us_county/tl_2017_us_county.shp", verbose = T, proj4string = P4S.latlon)

 harris <- subset(county.lines, county.lines$GEOID == 48201)
 brazoria <- subset(county.lines, county.lines$GEOID == 48039)
 galveston <- subset(county.lines, county.lines$GEOID == 48167)
 chambers <- subset(county.lines, county.lines$GEOID == 48071)
 liberty <- subset(county.lines, county.lines$GEOID == 48291)
 all.counties <- rbind(harris, brazoria, galveston, chambers, liberty)

 ShapePoly <- as(all.counties, "SpatialPolygons")
 study.area <- as(ShapePoly, "owin")

绘制研究区域的图是按计划进行的。

为了绘制炼油厂数据，我使用了以下代码

refinery <- data.frame(refinery)
attach(refinery)
refinery <- ppp(coords.x1, coords.x2, window = study.area)

炼油厂的数据最初是一个shapefile。我将其转换为数据框架，以便使用ppp()。

出现的问题是，当我阴谋炼油厂没有点，只有形状文件。

plot(refinery, pch = 20, col = "firebrick1")

​

​

我正在为一个类项目编写这个脚本，并且通常在空间项目中使用ArcGIS。我不知道是什么导致了这一点，因为我已经对ppp函数进行了广泛的研究，并且已经用不同的数据集复制了这段代码。我想把炼油厂划成对抗县形状的点。

Answer 1

我不能肯定地说，因为您没有提供一个可复制的示例(没有下载数据的链接，没有加载所使用的R包等)，但我认为这仅仅是由于着色:参数col控制研究区域的颜色(窗口)。要绘制给定颜色的点，请使用cols，例如cols = "firebrick1"。

注意：--你是直接在lon中工作的，所以你的点实际上是在一个球面(椭球)上，而在R中很少能直接对这些坐标进行分析(甚至绘图也略有错误，因为1 deg长与1 deg lat是不同的)。如果您继续分析，您确实需要将您的坐标投影到平面空间(例如，使用sf::st_transform()或sp::spTransform())。

Answer 2

在绘制类ppp的对象时，将执行函数plot.ppp。plot.ppp的帮助文件解释了点的颜色由参数cols控制。