How can I rasterize the polygons of a Shapefile in Julia?

Specific Domains Geo
17 
import GDAL
using ArchGDAL; const AG = ArchGDAL

raster = AG.registerdrivers() do
    AG.read("/vsizip/gadm36_ARG_shp.zip/gadm36_ARG_1.shp") do source_ds
        dest_ds = AG.Dataset(GDAL.rasterize(
            "gadm36_ARG_1.mem",
            Ptr{GDAL.GDALDatasetH}(C_NULL),
            source_ds.ptr,
            GDAL.rasterizeoptionsnew([
                "-of", "MEM",
                # "-l", "gadm36_ARG_1",   # The layer(s) used for input features.
                "-ts", "250", "250",      # target size: "-ts xsize ysize". Set it to whatever you want.
                "-ot", "Byte"
            ], C_NULL), C_NULL))
        result = AG.read(dest_ds)
        AG.destroy(dest_ds)
        result
    end
end