Я хочу працювати з растровими даними в R, але у мене є проблеми з експортом мого растрового набору даних з ArcMap. За допомогою "експорту даних" / GRID я отримую файл із розширенням * .aux. Що мені потрібно - * .grd.
Хтось може мені допомогти?
Я хочу працювати з растровими даними в R, але у мене є проблеми з експортом мого растрового набору даних з ArcMap. За допомогою "експорту даних" / GRID я отримую файл із розширенням * .aux. Що мені потрібно - * .grd.
Хтось може мені допомогти?
Відповіді:
Відповідь не експортуйте, бо в цьому немає потреби. Ви можете безпосередньо читати дані сітки Arc * за допомогою растрового пакету (коли бібліотека rgdal доступна). А якщо ви експортуєте, то краще уникати ascii (використовуйте, наприклад, GTiff).
library(raster)
library(rgdal)
#For an arc/grid (albem_s1 is the folder!):
r <- raster("G:/USGS/DEM/7_5min/VA/albem_s1")
#For a geotiff
r <- raster("G:/USGS/DEM/7_5min/VA/albem_s1.tif")
Більш зручне рішення для читання надає @RobertH в інших місцях цієї теми . Ви можете прочитати набір даних Arc * grid безпосередньо, використовуючи той самий підхід, описаний тут: просто посилайтеся на папку, в якій знаходяться файли даних сітки в raster
команді.
У Arc * використовуйте Raster to ASCII
інструмент (знайдений у Conversion Tools|From Raster
), щоб експортувати сітку у формат ASCII.
У R
, завантажте raster
і rgdal
пакети, а потім прочитайте растр за допомогою raster
функції, як у
r <- raster("G:/USGS/DEM/7_5min/VA/albem_s1.txt")
plot(r)
Щоб створити .grd версію, використовуйте writeRaster
, як в
writeRaster(r, "G:/USGS/DEM/7_5min/VA/albem_s1.grd")
Відповідь консолі опише, що знаходиться у .grd частині виводу (дані знаходяться у файлі .gri), як у
class : RasterLayer
filename : G:/USGS/DEM/7_5min/VA/albem_s1.grd
nrow : 1415
ncol : 1133
ncell : 1603195
min value : 70
max value : 960
projection : +proj=utm +zone=17 +ellps=GRS80 +datum=NAD83 +units=m +no_defs +towgs84=0,0,0
xmin : 686280
xmax : 720270
ymin : 4179990
ymax : 4222440
xres : 30
yres : 30
*.asc
? Я думаю, що обидва працюють, але, можливо, це може бути проблемою, якщо ви спробуєте відкрити *.txt
файл за допомогою ArcGIS.
.txt
а не .asc
. Я думав, що це може бути контр-інтуїтивно зрозумілим.
Не видається, що у файлу ESRI GRID є один файл із розширенням .grd. Натомість, здається, це формат з декількома папками / мультифайлами, який ArcMap сприймає як один файл (подібний до того, як працюють shapefiles та ArcInfo). У будь-якому випадку файл .aux - це файл-помічник у форматі GRID, який зберігає інформацію, яка не може бути збережена в самому растрі. Ось кілька посилань, які говорять про це:
Про допоміжні файли http://webhelp.esri.com/arcgiSDEsktop/9.3/index.cfm?TopicName= About_auxiliary_files
Коли створюється .aux файл? http://support.esri.com/en/knowledgebase/techarticles/detail/29821
конвертувати ваші растрові дані в ascii і зберегти в будь-якому вашому каталозі. то n R
library(rgdal)
libray(sp)
data<-readGDAL("C:/filepath/file")
це прочитає ваш сітковий файл ascii, ви також можете перевірити свої дані за допомогою команд varioud, таких як
summary(data)
попереджувально збережіть ваші растрові дані як TIFF, потім, використовуючи той самий пакет (rgdal), він зчитує файли tiff, які ви будете читати як
data<-readGDAL("C:/filepath/file.tiff")