Я импортирую и экспортирую цифровые модели рельефа АСТЕР в R, используя пакет raster. Эти данные имеют географическую систему координат c, так что экстент растра и его ячейки указаны в градусах.
Я считаю, что при импорте и экспорте одного и того же набора данных должно быть некоторое арифметическое значение c. Например, я импортирую игрушечную сцену АСТЕР ASTGTMV003_N48W122_dem.tif в R (я загрузил набор данных здесь ) и хочу узнать ее разрешение:
# example scene downloaded from EarthExplorer (https://earthexplorer.usgs.gov/)
r <- raster("ASTGTMV003_N48W122_dem.tif")
res(r)
И R сообщает мне, что разрешение (x, y): 0.0002777778 0.0002777778 QGIS, однако, предполагает, что тот же набор данных имеет гораздо больше цифр: 0.0002777777777777779944,-0.0002777777777777779944.
Итак, вот поворот: я теперь экспортирую этот набор данных как GeoTIFF и импортируйте его снова и оставьте все параметры без изменений. R скажет мне, что ничего не изменилось:
writeRaster(r, "test.tif") # export GTiff
r2 <- raster("test.tif") # import again
res(r) == res(r2) # which returns TRUE TRUE
Тем не менее, в ГИС я обнаружил, что мой новый test.tif имеет разрешение 0.0002777777778394882884,-0.0002777777778394882884, которое отличается от оригинальной сцены АСТЕР. В то время как R кажется прощающим в этом отношении, ГИС говорит мне, что экстенты и размеры ячеек НЕ одинаковы. Я воспроизвел этот пример с другими сценами АСТЕР и SRTM, и проблема остается. Тем не менее, цифры всегда отличаются от 13-го числа git и далее.
Я предполагаю, что проблема связана с неравным количеством строк и столбцов (3601, 3601). По крайней мере, используя данные ALOS3D 30 м с 3600 строками и столбцами, я не смог воспроизвести эту проблему.
Кто-нибудь сталкивался с подобной проблемой с экстентами и разрешениями при экспорте и импорте данных SRTM или ASTER?