У меня есть прямоугольная angular сетка, представляющая большую «область исследования» со средней глубиной каждой ячейки сетки:
preddataA2<- data.frame(x=gridObjA2$Spatial[,1], y=gridObjA2$Spatial[,2], area=500*500, DEPTH=Depthgrid_A2$Data$Info)
класс для gridObjA1 = rAIStoolsGRID; class для Depthgrid_A1 = rAIStoolsGRID
Однако я хотел бы пересечь эту большую сетку с моей фактической областью съемки, которая составляет SpatialPolygonsDataFrame:
survey.areaA2 <- spTransform(survey.areaA2, CRSobj=CRS(mainproj))
В идеале я хотел бы иметь сетку 500 * 500 для моей конкретной области исследования, так как я буду использовать функцию predict для прогнозирования пространственного распределения из обобщенной аддитивной модели.
Я относительно новичок в R, поэтому я не совсем уверен, существует ли уже код для этого.
Спасибо