Преобразование растра в кадр данных и извлечение нужных значений в R

Question

Я хочу сделать батиметрическую карту, используя ggplot2.Я преобразовал растровый файл во фрейм данных, закрыл и преобразовал его в фрейм данных:

library(raster)

##Import a raster file
Bathymetry_dat <- raster("w001000.adf")

##Clop the raster file
Bathy_clopped <- crop(Bathymetry_dat, extent(84.11236, 108.4594, -4.046979, 24.09534))

##Convert it to a dataframe    
datframe_bathy<-as.data.frame(Bathy_clopped, xy = TRUE)

Затем я проверил значения (т.е. глубину в м) в фрейме данных:

> summary(datframe_bathy)
       x                y          w001000_COUNT    
 Min.   : 84.11   Min.   :-4.046   Min.   :   9945  
 1st Qu.: 90.20   1st Qu.: 2.987   1st Qu.:  81618  
 Median : 96.28   Median :10.021   Median : 168447  
 Mean   : 96.28   Mean   :10.021   Mean   : 210212  
 3rd Qu.:102.37   3rd Qu.:17.054   3rd Qu.: 336718  
 Max.   :108.45   Max.   :24.087   Max.   :1205362  
                                   NA's   :3449125 

Значения глубины (м) должны быть отрицательными и не должны быть такими большими.Затем я проверил файл батиметрии, импортированный в R.

> Bathy_clopped
class       : RasterLayer 
dimensions  : 3377, 2922, 9867594  (nrow, ncol, ncell)
resolution  : 0.008333333, 0.008333333  (x, y)
extent      : 84.10833, 108.4583, -4.05, 24.09167  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : in memory
names       : w001000 
values      : -6138, -1  (min, max)
attributes  :
           ID  COUNT
 from: -11584      1
 to  :     -1 676804

Есть два атрибута: ID и COUNT.Я предполагаю, что ID - это глубина (м), и я не знаю, что такое COUNT.И странно, что COUNT колеблется от 1 до 676804, но Summary () показывает значения от 9945 до 1205362.

Поэтому мой вопрос заключается в том, как преобразовать растровый файл в кадр данных, содержащий значения, которые я хочу?

С помощью значений "COUNT" я смог создать батиметрическую карту, но значения в легенде не верны.

Заранее спасибо.

Answer 1

Существует таблица атрибутов с вашим растровым слоем.Мы можем использовать levels для доступа к нему.

library(raster)

# See the attribute table
head(levels(Bathy_clopped)[[1]])
#       ID COUNT
# 1 -11584     1
# 2 -10944     1
# 3 -10907     1
# 4 -10900     1
# 5 -10879     1
# 6 -10878     1

Затем мы можем манипулировать этой таблицей атрибутов, заменив COUNT глубиной, равной ID.

# Get the attribute table
RAT <- levels(Bathy_clopped)[[1]]

# Replace COUNT with ID
RAT$Depth <- RAT$ID
RAT$COUNT <- NULL

# Replace the attribute table
Bathy_clopped2 <- Bathy_clopped
levels(Bathy_clopped2)[[1]] <- RAT

# Create a single layer based on the new RAT
Bathy_clopped2 <- deratify(Bathy_clopped2)

# Create a data frame
datframe_bathy2 <-as.data.frame(Bathy_clopped2, xy = TRUE)

Теперь значения в datframe_bathy2 соответствуют ожидаемым.

# Summarize the data frame
summary(datframe_bathy2)
#       x                y              COUNT        
# Min.   : 84.11   Min.   :-4.046   Min.   :-6138    
# 1st Qu.: 90.20   1st Qu.: 2.987   1st Qu.:-3800    
# Median : 96.28   Median :10.021   Median :-2428    
# Mean   : 96.28   Mean   :10.021   Mean   :-2154    
# 3rd Qu.:102.37   3rd Qu.:17.054   3rd Qu.:  -65    
# Max.   :108.45   Max.   :24.087   Max.   :   -1    
#                                   NA's   :3449125