def bounds_for_naip(raster_dataset, rows, cols): """Clip the NAIP to 0 to cols, 0 to rows.""" left_x, right_x, top_y, bottom_y = \ NAIP_PIXEL_BUFFER, cols - NAIP_PIXEL_BUFFER, NAIP_PIXEL_BUFFER, rows - NAIP_PIXEL_BUFFER sw = pixel_to_lat_lon(raster_dataset, left_x, bottom_y) ne = pixel_to_lat_lon(raster_dataset, right_x, top_y) return {'sw': sw, 'ne': ne}