Queeno11 / Satellite-Imagery-Income-Prediction

Scripts for: Abbate (2024). Mapeando el ingreso de AMBA en alta resolución
MIT License
2 stars 0 forks source link

Remover secciones negras de la grilla #48

Closed Queeno11 closed 4 months ago

Queeno11 commented 6 months ago

Tengo que sacar las lineas negras de acá y resolver el problema con el bache rojo... Image

Queeno11 commented 6 months ago

La imagen que hay que reprocesar es la PNEO4_202211031357171_PAN_ORT R5C2, intentar conseguir los parametros con los que hice el processing asi no hay que correr todo de nuevo

Queeno11 commented 6 months ago

Revisar por que hay bordes negros en las imagenes de 2022. Son bordes negros o hay otro error?

Idea Para remover los bordes negros tengo que agregar una función dentro de utils.get_dataset_extent(ds) que corte las áreas negras del dataset.

def crop_black_borders_from_dataset(ds):
  # Genero binarias para cada celda identificando si tienen data
  iszero = (ds.band_data == 0)
  # Calculo el agregado para filas y columnas
  valid_y_values = iszero.all(dim="x").all(dim="band")
  valid_x_values = iszero.all(dim="y").all(dim="band")
  # Identifico el índice a partir de cual todos son True
  ## Falta get_min_valid_index. Capaz calcular la media movil de 5 y cuando sea 1 sacar todo ese rango desde t-2 hasta t+2?
  valid_idx_x_min, valid_idx_x_max= get_min_valid_index(ds, valid_x_values) 
  valid_idx_y_min, valid_idx_y_max= get_min_valid_index(ds, valid_y_values) 

  ds = ds.isel(x=slice(valid_idx_x_min, valid_idx_x_max), y=slice(valid_idx_y_min, valid_idx_y_max))
  return ds
Queeno11 commented 4 months ago

Solved here! #49