У меня возникают проблемы при работе с файлом DEM в R, который доступен в виде файла ADF сетки ESRI (например, папка с различными файлами ADF, такими как hdr.adf, w001001.adf et c.)
Я могу использовать растровый пакет для чтения и построения файла:
dem<-raster("w001001.adf")
plot(dem)
Получившийся график выглядит хорошо, и я могу делать такие вещи, как обрезка файла. Однако я хочу перепроектировать файл из:
+proj=laea +lat_0=-100 +lon_0=6370997 +x_0=45 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
В:
+proj=aea +lat_1=49.0628 +lat_2=50.4997 +lon_0=-113.5986 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +towgs84=0,0,0
на основе другого растра. Два растра имеют разные экстенты и разрешения в дополнение к проекциям, и я обычно использую spatial_sync_raster из пакета spatial.tools R для обработки этого типа преобразования. Но в результате получается пустой растр. Я не могу понять, связано ли это с тем, что я неправильно прочитал данные ESRI adf (то есть они все в памяти) или потому что пространственный синтаксис c слишком сложен. Я попытался просто обрезать слой dem до приблизительной степени в соответствии со вторым растром и использовать команду projectRaster из пакета raster (таким образом, устраняя проблему с разрешением и просто пробуя прямую проекцию), но снова получаю пустой растр. Так что я думаю, что проблема заключается в импорте файла ADF. Несколько дней поиска в Google не выявили решений только для R (я бы не стал вызывать ArcGIS).
Я ценю любую помощь!
Количество битов кода:
dem_p<-projectRaster(dem,crs=target_proj)
dem_p<-spatial_sync_raster(dem,climate,method="bilinear")