Ошибка чтения TIF файл «Географические данные CRS, не совместимой данных»

Question

Я пытаюсь прочитать в файле TIF следующим образом:

d = readGDAL('F182013.v4c.stable_lights.avg_vis.tif') 
Error in validityMethod(as(object, superClass)) : 
    Geographical CRS given to non-conformant data: -180.004166667 

Файл является night light data из американского Национального управления океанических и атмосферных исследований, так он должен быть действительным файлом. Возможно, проблема связана с численной проблемой с плавающей точкой.

Благодарен за любой совет, как обходным путем.

Answer 1

readGDAL вызывает функцию проверки ошибок validityMethod, которая проверяет, что долгота находится в пределах +/- 180 градусов для WGS84. Нефатальные ошибки и сообщения могут быть подавлены:

d = readGDAL('F182013.v4c.stable_lights.avg_vis.tif', silent = TRUE) 

так, что файл может, по крайней мере быть прочитана как SpatialGridDataFrame. Соответствующий географический диапазон объекта все еще находится за пределами +/- 180 град, но ошибка может быть в пределах допуска, который вы готовы принять.

bbox(d) 
     min  max 
x -180.00417 180.00417 
y -65.00417 75.00417 

Answer 2

Обходной путь:

require(tiff) 
r = raster(readTIFF(filepath, as.is = T)) 
extent(r) = extent(-180,180,-65,75) 
crs(r) = CRS("+proj=longlat +datum=WGS84")