Advertisement
Guest User

Untitled

a guest
Jun 27th, 2019
84
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.19 KB | None | 0 0
  1. raster = readGDAL("filename.TIFF")
  2. summary(raster)
  3. Object of class SpatialGridDataFrame
  4. Coordinates:
  5. min max
  6. x -180 180
  7. y -90 90
  8. Is projected: FALSE
  9. proj4string :
  10. [+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0]
  11. Grid attributes:
  12. cellcentre.offset cellsize cells.dim
  13. x -179.5 1 360
  14. y -89.5 1 180
  15. Data attributes:
  16. band1
  17. Min. : 1.0
  18. 1st Qu.:255.0
  19. Median :255.0
  20. Mean :226.9
  21. 3rd Qu.:255.0
  22. Max. :255.0
  23. NA's :245
  24.  
  25. head(data[7:8]) #Xp = longitude, Yp = latitude
  26. Xp Yp
  27. 7 -50.41777 1.33138
  28. 8 -59.25075 -19.16780
  29. 9 -67.25043 -27.83412
  30. 10 -54.91759 -15.00130
  31. 12 -67.58375 -28.00078
  32. 13 -65.08385 -34.33386
  33.  
  34. sp <- SpatialPoints(data[7:8])
  35. summary(sp)
  36. Object of class SpatialPoints
  37. Coordinates:
  38. min max
  39. Xp -76.08341 -35.25171
  40. Yp -38.66702 1.33138
  41. Is projected: NA
  42. proj4string : [NA]
  43. Number of points: 184
  44.  
  45. data$ndvi<-extract(raster, sp)
  46. Error in .subset(x, j) : invalid subscript type 'S4'
  47.  
  48. simr <- raster(ncol=36, nrow=18)
  49. simr[] <- 1:ncell(simr)
  50. simndvi<-extract(raster, sp)
  51. head(simndvi)
  52. [1] 301 373 408 373 408 444
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement