Блог

Главная - Вопросы по программированию - Присвоение CRS RaterBrick в R

Присвоение CRS RaterBrick в R

#r #raster #coordinate-systems #proj

#r #растр #системы координат #проект

Вопрос:

Я попытался преобразовать это df (c.data) в RasterBrick и установить значение crs EPSG:3035 ; Ламберт, равное площади Европы с кодом. Код создает растр, но не устанавливает CRS

 c.data
# A tibble: 361,478 x 6
         x       y   pop women mean_age hh_size
     <int>   <int> <int> <int>    <int>   <int>
 1 4334500 2684500    NA    NA       NA      NA
 2 4335500 2684500    NA    NA       NA      NA
 3 4336500 2684500    NA    NA       NA      NA
 4 4337500 2684500    NA    NA       NA      NA
 5 4338500 2684500    NA    NA       NA      NA
 6 4339500 2684500    NA    NA       NA      NA
 7 4334500 2685500    NA    NA       NA      NA
 8 4335500 2685500    NA    NA       NA      NA
 9 4336500 2685500    NA    NA       NA      NA
10 4337500 2685500    NA    NA       NA      NA
# ... with 361,468 more rows
  
 c.data_ras <- rasterFromXYZ(c.data, crs = st_crs(3035)$proj4string)
  
     c.data_ras
class      : RasterBrick 
dimensions : 868, 642, 557256, 4  (nrow, ncol, ncell, nlayers)
resolution : 1000, 1000  (x, y)
extent     : 4031000, 4673000, 2684000, 3552000  (xmin, xmax, ymin, ymax)
crs        : NA 
source     : memory
names      : pop, women, mean_age, hh_size 
min values :   1,     1,        1,       1 
max values :   6,     5,        5,       5

Вместо этого я получаю это предупреждение

Предупреждающее сообщение:

 In .newCRS(value) :
 proj=laea  lat_0=52  lon_0=10  x_0=4321000  y_0=3210000  ellps=GRS80  units=m  no_defs is not a valid PROJ.4 CRS string

Пожалуйста, как мне это сделать?

Спасибо

Спасибо @lime и @sequoia Я попробовал оба предложения, но я все еще получаю те же ошибки

 c.data_ras <- rasterFromXYZ(c.data, crs = st_crs(" init=epsg:3035")$proj4string)

Warning message:
    In .newCRS(value) :
       proj=laea  lat_0=52  lon_0=10  x_0=4321000  y_0=3210000  ellps=GRS80  units=m  no_defs is not a valid PROJ.4 CRS string

и

 c.data_ras <- rasterFromXYZ(c.data, crs = st_crs("EPSG:3035")$proj4string)
Warning message:
In .newCRS(value) :
   proj=laea  lat_0=52  lon_0=10  x_0=4321000  y_0=3210000  ellps=GRS80  units=m  no_defs is not a valid PROJ.4 CRS string

Комментарии:

1. заменить st_crs(3035) на st_crs("EPSG:3035")

2. … или с , crs = " init=epsg:3035" помощью . Я попробовал с образцами данных со страницы rasterFromXYZ справки, это сработало

3. упс, я хотел сказать st_crs("ESRI:3035") , что это должно сработать, и это намного быстрее, у меня была проблема с этим некоторое время назад, и это исправило ее.

4. К сожалению, ни один из них не сработал для меня.

5. попробуйте crs( c.data_ras) <- " proj=laea lat_0=52 lon_0=10 x_0=4321000 y_0=3210000 ellps=GRS80 units=m no_defs"