使用EO-1 Hyperion高光谱图像创建speclib

Question

我正在尝试从EO-1 Hyperion图像在R中创建一个speclib对象，并且需要帮助来创建代码。我已经创建了图像的rasterbrick对象，但是似乎无法弄清楚用于创建speclib对象的speclib函数代码。我假设我的rasterbrick对象肯定存在问题，不符合正确的结构，但不确定如何检查/解决此问题。此外，在检查str（图像）时，光栅砖对象缺少某些光谱信息，包括波长，值等。不知道为什么会这样。

Cau_Jun18 <- brick('Cau_Jun18_164Bands_FLAASH_SurRef_0_1_WGS84_UTM19N_24GCPs.dat') # Load image raster

str(Cau_Jun18)
Formal class 'RasterBrick' [package "raster"] with 12 slots
  ..@ file    :Formal class '.RasterFile' [package "raster"] with 13 slots
  .. .. ..@ name        : chr "E:\\1_R_PROJECTS\\5_EO-1_Hyperion_Analysis_in_R\\Cau_Jun18_164Bands_FLAASH_SurRef_0_1_WGS84_UTM19N_24GCPs.dat"
  .. .. ..@ datanotation: chr "FLT4S"
  .. .. ..@ byteorder   : chr "little"
  .. .. ..@ nodatavalue : num -Inf
  .. .. ..@ NAchanged   : logi FALSE
  .. .. ..@ nbands      : int 164
  .. .. ..@ bandorder   : chr "BIL"
  .. .. ..@ offset      : int 0
  .. .. ..@ toptobottom : logi TRUE
  .. .. ..@ blockrows   : int [1:164] 1 1 1 1 1 1 1 1 1 1 ...
  .. .. ..@ blockcols   : int [1:164] 979 979 979 979 979 979 979 979 979 979 ...
  .. .. ..@ driver      : chr "gdal"
  .. .. ..@ open        : logi FALSE
  ..@ data    :Formal class '.MultipleRasterData' [package "raster"] with 14 slots
  .. .. ..@ values    : logi[0 , 0 ] 
  .. .. ..@ offset    : num 0
  .. .. ..@ gain      : num 1
  .. .. ..@ inmemory  : logi FALSE
  .. .. ..@ fromdisk  : logi TRUE
  .. .. ..@ nlayers   : int 164
  .. .. ..@ dropped   : NULL
  .. .. ..@ isfactor  : logi FALSE
  .. .. ..@ attributes: list()
  .. .. ..@ haveminmax: logi FALSE
  .. .. ..@ min       : num [1:164] Inf Inf Inf Inf Inf ...
  .. .. ..@ max       : num [1:164] -Inf -Inf -Inf -Inf -Inf ...
  .. .. ..@ unit      : chr ""
  .. .. ..@ names     : chr [1:164] "B8..426.820010.Nanometers." "B9..436.989990.Nanometers." "B10..447.170010.Nanometers." "B11..457.340000.Nanometers." ...
  ..@ legend  :Formal class '.RasterLegend' [package "raster"] with 5 slots
  .. .. ..@ type      : chr(0) 
  .. .. ..@ values    : logi(0) 
  .. .. ..@ color     : logi(0) 
  .. .. ..@ names     : logi(0) 
  .. .. ..@ colortable: logi(0) 
  ..@ title   : chr(0) 
  ..@ extent  :Formal class 'Extent' [package "raster"] with 4 slots
  .. .. ..@ xmin: num 614669
  .. .. ..@ xmax: num 644039
  .. .. ..@ ymin: num 5298630
  .. .. ..@ ymax: num 5397150
  ..@ rotated : logi FALSE
  ..@ rotation:Formal class '.Rotation' [package "raster"] with 2 slots
  .. .. ..@ geotrans: num(0) 
  .. .. ..@ transfun:function ()  
  ..@ ncols   : int 979
  ..@ nrows   : int 3284
  ..@ crs     :Formal class 'CRS' [package "sp"] with 1 slot
  .. .. ..@ projargs: chr "+proj=utm +zone=19 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"
  ..@ history : list()
  ..@ z       : list()

我正在尝试使用所有频段组合（具有164个频段）创建各种植被指数，然后对与响应变量相关的频段进行最佳分析。