R'de Çokgen şekil dosyası kullanarak Rasterden Raster çıkarın

R için yeni ve raster paketini kullanıyorum. Varolan bir raster dosyasından çokgenler ayıklanırken bir sorunum var. Kullanırsam

extract(raster, poly_shape)

Rasterdeki fonksiyon her zaman verilerle bir liste oluşturur. Gerçekten istediğim, ArcGIS ile tekrar yükleyebileceğim başka bir raster dosyasını ayıklamak. Biraz daha okuduktan sonra kırpma fonksiyonunun gerçekten ihtiyacım olduğunu düşünüyorum. Ama bu işlevi kullanmaya çalıştığımda

crop(raster, poly_shape)

Bu hatayı alıyorum:

Error in .local(x, y, ...) : extents do not overlap
In addition: Warning message:
In intersect(extent(x), extent(y)) : Objects do not overlap

Raster ve poly_shape dosyaları her iki işlev için aynıdır. Burada neyin yanlış olabileceğini söyleyebilir misin? Kırpma işlevinin liste değil başka bir tarama oluşturması doğru mu?

EDIT : scope () işlevi benim için çalışmıyor. Ben hala aynı hatayı alıyorum. Ama 2 veri kümesinin çakıştığından eminim! İle

extract(raster, poly_shape)

Ondan doğru verileri alıyorum. Sadece bir liste gibi ve istediğim gibi bir tarama değil. Daha önce ArcGIS'deki veri kümelerini yükledim ve çok iyi uyuyorlar, bu yüzden projeksiyonu kontrol etmedim. Şimdi denedim

projection(raster) # "+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs"
projection(poly_shape) # "+proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs"

ve projeksiyonların uymadığını görebilirsiniz. Extract işlevi, dosyaları otomatik olarak doğru şekilde dönüştürebiliyor gibi görünüyor. Bunu biliyorum çünkü aşağıdakileri yaptım:

• ArcGIS'de de R'de çıkarılan poligonun tam kısmını kestim
• Ekstre edilen R poligonunun (liste) tüm değerlerinin toplamını hesapladım
• ArcGIS'de kestiğim tüm tarama hücrelerinin toplamını hesapladım

2 aynı sonucu var bu yüzden sonuç ayıklama işlevi doğru çalışması olması gerektiğini düşünüyorum. Şimdi sanırım 2 seçeneğim var:

1. Raster'ı tekrar çıkarılan listeden çıkarmanın bir yoluna ihtiyacım var veya
2. 2 veri kümesinin (raster + poly_shape) aynı önlemi kullanması gerekir ve kırpma işlevi çalışmalıdır

Burada ne yapmayı önerirsiniz?

Çıkartma işlevi tam olması gerektiği gibi davranıyor. Kırpma işlevini çokgenin kapsamını kullanmaya zorlayabilir ve ardından nesneyi çokgen alanını temsil eden tam raster döndürmek için maskeleyebilirsiniz. Hatayı almaya devam ederseniz, aslında verileriniz çakışmaz.

Lütfen R'nin "anında" projeksiyon yapmadığını unutmayın, bu yüzden projeksiyonlarınızı kontrol edin. "Kapsam ()" işlevini kullanarak uzantılarınızın çakışıp çakışmadığını kontrol edebilirsiniz.

Burada çokgen genişliğini kullanarak kırpma ve ardından "rasterleştirilmiş" çokgen kullanarak elde edilen raster maskeleme örneği.

# Add required packages
require(raster)
require(rgeos)
require(sp)

# Create some data using meuse 
data(meuse)
  coordinates(meuse) <- ~x+y
    proj4string(meuse) <- CRS("+init=epsg:28992")
data(meuse.grid)
  coordinates(meuse.grid) = ~x+y
    proj4string(meuse.grid) <- CRS("+init=epsg:28992")
      gridded(meuse.grid) = TRUE    
        r <- raster(meuse.grid) 
          r[] <- runif(ncell(r))

# Create a polygon
f <- gBuffer(meuse[10,], byid=FALSE, id=NULL, width=250, 
                         joinStyle="ROUND", quadsegs=10)   

# Plot full raster and polygon                       
plot(r)
  plot(f,add=T)

# Crop using extent, rasterize polygon and finally, create poly-raster
#          **** This is the code that you are after ****  
cr <- crop(r, extent(f), snap="out")                    
  fr <- rasterize(f, cr)   
    lr <- mask(x=cr, mask=fr)

# Plot results
plot(lr)
  plot(f,add=T)

Aslında aradığım şey işlevdi mask().

mask(raster, poly_shape)

hatasız çalışır ve aradığım şeyi döndürür.

Kapsamı gayet iyi çalışıyor ... Sanırım sizin dereceniz Xmin, Xmax, Ymin ve Ymax, rasterinizin X ve Y'sinden farklı - yani tam tersi