Load a world raster with ArchGDAL and crop australia

loadraster.jl

using ArchGDAL

function readtiff(file)
    img = ArchGDAL.registerdrivers() do
        ArchGDAL.read(file) do dataset
            ArchGDAL.read(dataset)
        end
    end
    img ./= maximum(img)
    img .= max.(0.0, img)
    img = img[:,:,1]'
end

cropaust(x) = x[950:1350, 3100:3600] # Australia

world_suitability = readtiff("your/path/limited_growth_2017_06.tif") 
suit = cropaust(world_suitability)