Remove Points fall within the same pixel

This practical is to solve GPS points that are very close to each other in the same imagery pixel.

Load library

Read raster and shapefile point

point<-shapefile("F:\\Research\\Research_Cooperation\\ILRI\\UTM_Shapefile.shp")

point$Sum_Cx<-as.numeric(point$Sum_Cx)

path<-"F:\\Research\\Research_Cooperation\\ILRI\\output.tif"

myraster<-stack(path)

Extract the positional number of points in raster

mynumber<-cellFromXY(myraster, point)

Merging mynumber to point and convert to sf geodataframe

point$Raster_ID<-mynumber

library(sf)

Mypoint<-st_as_sf(point)

Mypoint$Raster_ID<-as.factor(Mypoint$Raster_ID)

Summary and merge

# Calculate the total of Sum_Cx grouped by Raster_ID
Mytotal<- Mypoint %>% group_by(Raster_ID)%>% summarise(Sum_total=sum(Sum_Cx))

# Remove duplicated Raster_ID 

Mypoint_Uni_Raster_ID<-Mypoint[!duplicated(Mypoint$Raster_ID),]

# Join two datasets together

Mypoint_final<-st_join(Mytotal,Mypoint_Uni_Raster_ID,by="Raster_ID")

head(Mypoint_final)

Simple feature collection with 6 features and 17 fields
geometry type:  POINT
dimension:      XY
bbox:           xmin: 566512.2 ymin: 2337469 xmax: 567087.4 ymax: 2337871
CRS:            +proj=utm +zone=48 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0