Я пытаюсь исправить повернутый растр с помощью пакета R Raster.Ниже приведена реализация, в которой mnf_filename является растровым файлом, который я пытаюсь прочитать.
mnf <- brick(mnf_filename)
Выдает следующее предупреждение:
Этот файл имеет ротацию. Поддержка таких файлов ограничена, и результаты обработки данных могут быть неправильными.Действуйте с осторожностью и рассмотрите возможность использования функции «выпрямить»
Запуск функции выпрямления mnf_r <- rectify(mnf)
Выдает следующую ошибку:
Ошибка в.intersectExtent (x, y, validate = TRUE): объекты не пересекаются. Дополнительно: Предупреждающие сообщения: 1: In .local (x, ...): min y больше, чем max y 2: In .local (x, ...): min y больше, чем max y
И при выпрямлении всего одного слоя в стеке растров, например, mnf_r <- rectify(subset(mnf, 1) генерирует следующую ошибку:
Ошибка вvalidityMethod (объект): недопустимый экстент: ymin> = ymax
Я не могу сказать, что не так.Ниже приведены выходные данные экстента (mnf):
класс: экстент
xmin: 720770.2
xmax: 752353
ymin: 4328909
ymax: 4361487
Ниже выводится gdalinfo(mnf_filename):
[1] "Драйвер: GTiff / GeoTIFF"
[2] "Файлы: X: \gcwu \ AVIRIS \ f150429t01p00r07_refl_hpc18_v1_topo_American_MNF2.tif "[3]" Размер: 701, 1821 "
[4]" Система координат: "
[5]" PROJCS [\ "WGS 84 / UTM зона"
[6]" GEOGCS [\ "WGS 84 \", "
[7]" DATUM [\ "WGS_1984 \", "
[8]" СФЕРОИД [\ "WGS 84 \",6378137,298.257223563, "
[9]" AUTHORITY [\ "EPSG \", \ "7030 \"]], "
[10]" AUTHORITY [\ "EPSG \", \ "6326 \"]], "
[11]" PRIMEM [\ "Greenwich \", 0, "
[12]" AUTHORITY [\ "EPSG \", \ "8901 \"]], "
[13] "UNIT [\" Степень \ ", 0.0174532925199433,"
[14] "AUTHORITY [\" EPSG \ ", \" 9122 \ "]],"
[15] "AUTHORITY [\" EPSG \ "", \" 4326 \ "]],"
[16] "PROJECTION [\" Transverse_Mercator \ "],"
[17] "PARAMETER [\" latitude_of_origin \ ", 0],"
[18] "PARAMETER [\ "central_meridian \", - 123], "
[19]" PARAMETER [\ "scale_factor \", 0.9996], "
[20]" PARAMETER [\ "false_easting \", 500000],"
[21]" PARAMETER [\ "false_northing \", 0], "
[22]" UNIT [\ "meter \", 1, "
[23]" AUTHORITY [\ "EPSG\ ", \" 9001 \ "]],"
[24] "ОСЬ [\" Восток \ ", ВОСТОК],"
[25] "ОСЬ [\" Север \ ", СЕВЕР],"
[26] "AUTHORITY [\" EPSG \ ", \" 32610 \ "]]"
[27] "GeoTransform ="
[28] "720770.204, 13.16436670051244, -12.27597040459283"
[29] "4352881.37, -12.27597040459283, -13.16436670051244"
[30] "Метаданные:"
[31] "AREA_OR_POINT = Area"
[32] "TIFFTAG_XRESOLUTION = 1"
[33]"TIFFTAG_YRESOLUTION = 1"
[34] "Метаданные структуры изображения:"
[35] "INTERLEAVE = BAND"
[36] "Угловые координаты:"
[37] "Верхний левый (720770.204), 4352881.370) (120d26'23.54 \ W, 39d17'49.87 \ "N)"
[38] "Слева внизу (698415.662, 4328909.058) (120d42'21.35 \ W, 39d 5'12.32 \" N) "
[39] "Верхний правый" (729998,425, 4344275,915) (120 д.)20 '9,24 \ "W, 39d13' 2,37 \" N) "
[40]" Нижний правый (707643,883, 4320303.603) (120d36 '6,99 \ W, 39d 0'25,66 \ "N)"
[41] "Центр (714207.043, 4336592.486) (120d31'15.89" W, 39d 9 '7.89 \ "N)"