Создание спецлиба с использованием гиперспектрального изображения EO-1 Hyperion

Question

Я пытаюсь создать объект Speclib в R из образа EO-1 Hyperion, и мне нужна помощь в создании кода. Я создал растровый объект для своего изображения, но не могу понять, как работает код функции speclib для создания объекта speclib. Я предполагаю, что должны быть проблемы с моим объектом растрового кирпича, не соответствующим надлежащей структуре, но не уверен, как проверить / решить эту проблему. Кроме того, при проверке 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 полосы), а затем провожу анализ, какие полосы являются лучшими по отношению к переменной отклика.