Блог

Главная - Вопросы по программированию - Загрузка нескольких растров в R приводит к некоторым с атрибутами

Загрузка нескольких растров в R приводит к некоторым с атрибутами

#r #raster #r-raster

#r #растр #r-растр

Вопрос:

Я пытаюсь загрузить несколько растров (всего 87) в R для анализа. Растры загружаются непосредственно в стек с использованием этого кода:

 filelist <- list.files(path="mypath", pattern="tif$", full.names=TRUE)
stk1 <- raster::stack(filelist)

Данные в основном загружаются правильно, за исключением 9 проблемных полос, которые, похоже, имеют ту же проблему. Эти плохие полосы загружаются с «атрибутами». Смотрите пример ниже

goodband

 > goodband
class       : RasterLayer 
dimensions  : 1768, 3674, 6495632  (nrow, ncol, ncell)
resolution  : 40, 40  (x, y)
extent      : 135800.4, 282760.4, 22297.69, 93017.69  (xmin, xmax, ymin, ymax)
coord. ref. :  proj=tmerc  lat_0=49  lon_0=-2  k=0.9996012717  x_0=400000  y_0=-100000  datum=OSGB36  units=m  no_defs  ellps=airy  towgs84=446.448,-125.157,542.060,0.1502,0.2470,0.8421,-20.4894 
data source : goodband.tif 
names       : Good 
values      : 0.5, 36  (min, max)

плохая полоса

 > badband
class       : RasterLayer 
dimensions  : 1768, 3674, 6495632  (nrow, ncol, ncell)
resolution  : 40, 40  (x, y)
extent      : 135800.4, 282760.4, 22297.69, 93017.69  (xmin, xmax, ymin, ymax)
coord. ref. :  proj=tmerc  lat_0=49  lon_0=-2  k=0.9996012717  x_0=400000  y_0=-100000  datum=OSGB36  units=m  no_defs  ellps=airy  towgs84=446.448,-125.157,542.060,0.1502,0.2470,0.8421,-20.4894 
data source : badband.tif 
names       : Bad 
values      : 33, 6824  (min, max)
attributes  :
        ID OID Value Count
 from:   0   0    33   433
 to  : 529 529  6824  1307

Как вы можете видеть, полоса пропусков имеет дополнительные «атрибуты» внизу. При попытке проанализировать плохие диапазоны, они просто заполнены нулевыми значениями и не работают как хорошие диапазоны.

Не могли бы вы посоветовать, как лучше всего решить эту проблему?

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

1. Вы пробовали открывать «плохие полосы», например, с помощью QGIS, и посмотреть, все ли с ними в порядке? Если нет, то проблема не в rasterstack, а во входных данных (что мне кажется вероятным)

2. Да, я открыл в ArcGIS (где они были созданы), и они отлично работают. Я также воссоздал их и повторно расположил растры с той же проблемой

3. Возможно, странная вещь заключается в том, что ваш «плохой» диапазон значений имеет совершенно другой диапазон значений, чем «хороший» диапазон. Это правильно?

4. Данные являются геохимическими данными. Диапазоны не так уж и отличаются, они по-прежнему на два порядка от минимального до максимального

5. Да, но один равен от 0,5 до 36, а другой от 33 до 6824. Если они являются одним и тем же параметром, это кажется странным. Но опять же, вы, очевидно, знаете это лучше меня!