Извлечение пикселей высочайшего качества MYD09GA - PullRequest
Новая
       23

Извлечение пикселей высочайшего качества MYD09GA

0 голосов
/ 30 марта 2020

На основании приведенного ниже описания слоя MYD09GA QC_500m мне нужно извлечь только данные самого высокого качества band1. Мой текущий код также приведен ниже. После применения маски, когда я пытаюсь сгенерировать диаграмму временных рядов для Band1, я получаю сообщение об ошибке: Ошибка при создании диаграммы: столбцы данных для оси # 0 не могут иметь тип string. Что-то не так с моей функцией маски? или мне нужно сделать какое-либо преобразование данных для генерации временных рядов? Может кто-нибудь помочь, пожалуйста. 

var mod = ee.ImageCollection("MODIS/006/MYD09GA"),
geometry1 = /* color: #d63000 */ee.Geometry.Point([10.644168108701706, 64.4188479841206]);

//Selecting Band 1 and QC layers
var dataset = mod.filterDate('2013-01-01', '2013-12-31').select('sur_refl_b01','QC_500m')
var single = dataset.first()

//create mask to extract only highest quality band 1 data 
var filter = function(image){ 
var mask = single.select('QC_500m').eq(0)
return ee.Image(image).updateMask(mask)
}
var maskapply = dataset.map(filter)
var chart = ui.Chart.image.seriesByRegion({
imageCollection: maskapply, 
regions: geometry1, 
reducer: ee.Reducer.mean(),
scale: 1000,
band: 'sur_refl_b01',
})
print(chart)

MYD09GA LAYERS info QC_500m: качество полосы отражения 500 м Бит 2-5: диапазон данных качества полосы 1, четыре бита 0000 (0): высшее качество 1000 (8): детектор мертвых сигналов; данные, интерполированные в L1B 1001 (9): солнечный зенит ≥ 86 градусов 1010 (10): солнечный зенит ≥ 85 и <86 градусов 1011 (11): пропущенный вход 1100 (12): внутренняя постоянная, используемая вместо климатологических данных как минимум для одна атмосфера c константа 1101 (13): коррекция за пределы пикселя ограничена до предельно допустимого значения 1110 (14): данные L1B неисправны 1111 (15): не обработано из-за глубокого океана или облаков </p>

Мой код это:

...