0

私は8つの環境変数のraster.stackと米国南部のエコリージョンのシェープファイルを持っています

files<-list.files(path='E:/Ecoregions Models/Border Bioclim/',pattern='asc', full.names=TRUE)
predictors<-stack(files)

##subset the main shape file into 12 individual shapefiles for the regions of interest 
ER_11.1<-Level.2.ecoregs[Level.2.ecoregs$NA_L2CODE=="11.1",]...

関数を使用してextract、すべてのエコリージョンから 5,000 のランダム ポイントを抽出し、それらを別のオブジェクトとして保存しようとしています。

個々のエコリージョン オブジェクトをラスター ファイルに投影するのに問題があります。すでに投影を変更しましたが、まだ問題があります。

>bio_1
class       : RasterLayer 
dimensions  : 324, 444, 143856  (nrow, ncol, ncell)
resolution  : 0.08333333, 0.08333333  (x, y)
extent      : -125, -88, 18.5, 45.5  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=48 +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
data source : E:\Ecoregions Models\Border Bioclim\ModelVariables\bio_1.asc 
names       : bio_1 

> ER_11.1
class       : SpatialPolygonsDataFrame 
nfeatures   : 22 
extent      : -1989641, -1450736, -1748693, -208349.6  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=48 +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
nvariables  : 8
names       : NA_L2CODE,                NA_L2NAME, NA_L1CODE,                NA_L1NAME,                       NA_L2KEY,                     NA_L1KEY,   Shape_Leng,   Shape_Area 
min values  :      11.1, MEDITERRANEAN CALIFORNIA,        11, MEDITERRANEAN CALIFORNIA, 11.1  MEDITERRANEAN CALIFORNIA, 11  MEDITERRANEAN CALIFORNIA,     331.1174, 1.483274e+08 
max values  :      11.1, MEDITERRANEAN CALIFORNIA,        11, MEDITERRANEAN CALIFORNIA, 11.1  MEDITERRANEAN CALIFORNIA, 11  MEDITERRANEAN CALIFORNIA, 5417939.1114, 8.472524e+04 
>

ここでの問題は、さまざまな程度にあるに違いありません。ラスター パッケージでこれを解決する方法について何か提案はありますか?

4

1 に答える 1

1

bio_1 の CRS は明らかに間違っています。

resolution  : 0.08333333, 0.08333333  (x, y)
extent      : -125, -88, 18.5, 45.5  (xmin, xmax, ymin, ymax)

つまり、米国の大部分の経度/緯度座標ですが、使用している

coord. ref. : +proj=utm +zone=48 +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0

おそらく、データを変換 (投影) したいときに投影名を変更したためでしょうか? できることは次のとおりです。

bio_1 の CRS を元の値に戻します (推測します)

projection(bio_1) <- "+proj=longlat +datum=WGS84"

ポリゴンを同じCRSに変換します

library(rgdal)
ER <- spTransform(ER_11.1, CRS(projection(bio_1)) )

次に、値を抽出します。

v <- extract(bio_1, ER)
于 2013-07-04T20:20:09.350 に答える