# Rasterizing feature data

**URL:** <https://discourse.julialang.org/t/rasterizing-feature-data/59594>\
**Category:** Geo\
**Created:** [April 19, 2021, 1:53pm UTC](https://discourse.julialang.org/t/rasterizing-feature-data/59594 "2021-04-19T13:53:07Z")\
**Posts on this page:** 1\
**Showing post:** 12

<div class="post-metadata">

**Author:** ![mihalybaci](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/mihalybaci/32/13528_2.png) [@mihalybaci](https://discourse.julialang.org/u/mihalybaci)\
**Post date:** [April 20, 2021, 4:32pm UTC](https://discourse.julialang.org/t/rasterizing-feature-data/59594/12 "2021-04-20T16:32:48Z")

</div>

Thanks for the replies!

@joa-quim - This seems pretty easy, though I don’t have GMT installed because I couldn’t get past versioning issues with GMT and GDAL (mentioned in [this thread](https://discourse.julialang.org/t/gmt-jl-errors-when-trying-to-plot-problem-with-png/47860)).

@visr - This method seems easy enough as well, I just don’t have enough experience with GDAL (thought I could fumble around a bunch). I tried to save out the resulting raster with

```julia
AG.gdalrasterize(AG.Dataset(dataset.ptr), ["-of","MEM","-tr","0.05","0.05"]) do ds_raster
    # e.g. save the raster
    println(ds_raster)
    AG.write(ds_raster, "test.tif")
end

```

but no file gets saved. Is this the right way to go about it? The result of `println` tells me the resulting file is a 1x1 dataset, even when I use my real WKT, which seems odd.

---

_[View the full topic](https://discourse.julialang.org/t/rasterizing-feature-data/59594)._