Блог

Главная - Вопросы по программированию - Как преобразовать код растеризации() в gdal_rasterize() в R для растеризации точечного слоя источника данных OSM на основе количества в каждой ячейке?

Как преобразовать код растеризации() в gdal_rasterize() в R для растеризации точечного слоя источника данных OSM на основе количества в каждой ячейке?

#r #raster #gdal #r-raster #rasterizing

Вопрос:

Существует .osm источник данных, который я хочу растеризировать на основе количества точек в каждой ячейке. Я могу сделать это, сначала прочитав .osm файл в sf объект, а затем используя raster::rasterize() функцию. Однако, как отмечается в Геокомпьютере в R, это происходит медленно.

Похоже gdalUtils::gdal_rasterize() , что вместо этого можно использовать гораздо более быструю скорость, я не мог понять, как переписать подобную операцию для требуемого синтаксиса gdalUtils::gdal_rasterize() .

 raster::rasterize(x = input_sf, y = raster_template, field = 1, fun = "count")

До сих пор я пытался выполнить подобную команду, которая создавала только пустой растр без каких-либо данных.

 gdalUtils::gdal_rasterize(src_datasource = "xxx.osm",
                      l = "points",
                      dst_filename = "xxx.tif",
                      a = "osm_id",
                      te = c(xmin, ymin, xmax, ymax),
                      tr = c(50, 50),
                      verbose = T)

Должен ли я использовать этот sql аргумент? Я прочитал документацию OGR SQL здесь, но не понял ни одного полезного совета. Как я могу выполнить ту же операцию подсчета или любую другую пространственную операцию с помощью gdalUtils::gdal_rasterize() утилиты?

Метки: