No se puede encontrar la ubicación de la OTAN UTM en Sentinel-2

10

Mirar coordenadas 31.96212, -103.004715

Los convertidores UTM dan sus coordenadas UTM 13/R/FR.

El convertidor de ejemplo está aquí: http://www.rcn.montana.edu/resources/converter.aspx

Pero hay muchos de ellos y dan respuestas similares para estas coordenadas.

Simultáneamente, en el conjunto de datos Sentinel-2 aquí http://sentinel-s2-l1c.s3-website.eu-central-1.amazonaws.com/#tiles/13/R/

No puedo encontrar el FRsubdirectorio.

En google esta ubicación está aquí:

ingrese la descripción de la imagen aquí

Y veo el mismo lugar en el navegador de imágenes Sentinel , ese mosaico es diferente

ingrese la descripción de la imagen aquí

que significa,13/S/FR es decir, la misma banda UTMcuadrada pero diferente.

¿Cómo es esto posible?

ACTUALIZAR

KML con mosaicos Sentinel-2 también informa Smosaico en una ubicación determinada

ingrese la descripción de la imagen aquí

ACTUALIZACIÓN 2

De acuerdo con esta imagen

ingrese la descripción de la imagen aquí

tomada desde aquí , la FRplaza se encuentra mitad en la Szona UTM y mitad en la Rzona. Obviamente, la mayoría de los convertidores automáticos asignan este cuadrado a la Rzona, mientras que Sentinel-2 lo representa por Szona.

¿Hay alguna verdad aquí?

ACTUALIZACIÓN 3

El código simple de Python, tomado de aquí https://gis.stackexchange.com/a/224994/32207

bandVals = "CDEFGHJKLMNPQRSTUVWXX"

lon = 31.96212
lat = -103.004715

zone = int(lat + 186.0) / 6

if (lon >= 84.0):
    band = 'Y' if (lat < 0.0) else 'Z'
elif (lon <= -80.0):
    band = 'A' if (lat < 0.0) else 'B'
else:
    band = bandVals[int(lon + 80.0) / 8]

print '{:02d}{:s}'.format(zone,band)

También vuelve 13R.

¿Es este error en los datos de Sentinel-2 o qué?

Dims
fuente
Es S/FR, mientras que los convertidores UTM dan R/FR. ¿Cómo calcular la ubicación si los convertidores UTM funcionan incorrectamente?
atenúa
El valor de latitud es un poco menos de 32 grados norte. Eso lo pone directamente en la "banda" de latitud R. Sentinel-2 puede haber sido en mosaico utilizando el punto central de la ficha que podría estar en la banda "S".
mkennedy
@mkennedy ¿cómo simular este algoritmo a partir de coordenadas?
atenúa
2
También puede considerar informar esto a [email protected], ya que de hecho parece un comportamiento inesperado.
Kersten

Respuestas:

1

En respuesta a su pregunta de comentario "cómo simular este algoritmo":

Esta es una solución bastante bruta, pero fácil de implementar y debería ofrecer un buen rendimiento:

  1. Utilice cualquiera de los convertidores UTM que funcionan "como se esperaba", colocando las coordenadas en 13R.
  2. Luego, verifique si la carpeta existe en la estructura de datos de Sentinel 2. Si es así, has terminado, ¡hurra!

  3. De lo contrario, verifique las cuadrículas UTM vecinas y vea si el mosaico / carpeta "FR" existe en ellas. Dado que hay superposiciones en todas partes, debería verificar todas las 8 cuadrículas circundantes.
    El orden más probable para verificar sería 13S, 13Q, 12R, 14R, 12S, 14S, 12Q, 14Q.
    Los últimos cuatro podrían ser relevantes si sus coordenadas se encuentran en las esquinas de una zona UTM, pero son altamente improbables.

Dada la forma en que Sentinel2 etiqueta los mosaicos, solo uno de los vecinos debería tener dicha carpeta, lo que garantiza que obtenga el archivo correcto.

Cualquier otra solución geográficamente más "correcta" implicaría una carga computacional mucho mayor de lo que creo se justifica aquí.

Y definitivamente, definitivamente informe esto al equipo de la ESA como lo sugiere Kersten en los comentarios. Realmente no entiendo por qué eligieron un sistema organizativo tan innecesariamente complicado.

Senshi
fuente
0

Publicación relacionada aquí

Lo que me ha funcionado es utilizar el KML S2 proporcionado por ESA para calcular todos los mosaicos que se cruzan con mi AOI, y luego buscar estos mosaicos en AWS.

Este KML parece funcionar como una definición de todos los identificadores de mosaico posibles generados por S2, eliminando muchas opciones superpuestas.

Al mirar el KML (solo inspección visual, no 100% seguro) me parece que en el peor de los casos tendrías que buscar 4 mosaicos.

Sería bueno tener el algoritmo que ESA usó para definir el KML para hacerlo más eficiente.

Fred Liporace
fuente