def get_geotransform(plot_center_xy, args):
    """
    Get geotransform from plot center and plot dims.
    Structure: top_left_x, pix_width_in_meters, _, top_left_y, pix_height_in_meters (neg for north up picture)
    Geotransform reference : https://gdal.org/user/raster_data_model.html
    """

    return [
        plot_center_xy[0] - args.diam_meters // 2,  # xmin
        args.diam_meters / args.diam_pix,
        0,
        plot_center_xy[1] + args.diam_meters // 2,  # ymax
        0,
        -args.diam_meters / args.diam_pix,
        # negative b/c in geographic raster coordinates (0,0) is at top left
    ]