saya harus mengubahnya menjadi UTM untuk membuat buffer berfungsi.
wets<-readOGR(dsn=".",layer="shapefile")
r.raster <- raster()
extent(r.raster) <- extent(wets)
res(r.raster) <- 100
wets.r <- rasterize(wet,r.raster)
plot(wets.r)
wetsbuf<-buffer(wets.r,width=500)
Selama pembuatan buffer yang merupakan baris kode terakhir, ia memberikan peringatan ini:
Warning message:
In couldBeLonLat(x) :
raster has a longitude/latitude CRS, but coordinates do not match that
inilah informasinya
summary(wets.r)
layer
Min. 1
1st Qu. 1
Median 2
3rd Qu. 9
Max. 11
NA's 52629
summary(wets)
Object of class SpatialPolygonsDataFrame
Coordinates:
min max
x 683705 714088.8
y 4326266 4343768.0
Is projected: TRUE
proj4string :
[+proj=tmerc +lat_0=0 +lon_0=24 +k=0.9996 +x_0=500000 +y_0=0 +datum=GGRS87
+units=m +no_defs +ellps=GRS80 +towgs84=-199.87,74.79,246.62]
Data attributes:
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.0 2.5 5.0 5.0 7.5 10.0
wets.r
class : RasterLayer
dimensions : 175, 304, 53200 (nrow, ncol, ncell)
resolution : 100, 100 (x, y)
extent : 683705, 714105, 4326268, 4343768 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
data source : in memory
names : layer
values : 1, 11 (min, max)
attributes :
ID FID
from: 1 0
to : 11 10
Saya harus mengubah proyeksi agar memungkinkan untuk melakukan buffer.
Jawaban:
Ini adalah bagaimana Anda dapat memproyeksikan ulang raster di R menggunakan paket raster . Dalam contoh ini, input geotiff berada dalam sistem koordinat geografis NAD83 dan saya memproyeksikan kembali ke sistem koordinat yang diproyeksikan NAD 83 UTM 15. Referensi yang baik untuk proyeksi format Proj4, yang digunakan oleh RGDAL, dapat ditemukan di spatialreference.org .
sumber