¿Existe un estándar para las coordenadas de píxeles en rásteres georreferenciados?

Estamos trabajando con diferentes datasets ráster, y las coordenadas de los píxeles se almacenan en el centro del píxel o en una esquina. La información generalmente se proporciona en los metadatos, pero a veces tenemos que adivinar y esta publicación ilustra que no siempre es fácil de encontrar.

Como también producimos algunas capas, me preguntaba si había una regla estándar y / o mejores prácticas para asignar las coordenadas de un píxel, idealmente con una referencia o con pro / contras.

Esto es un poco desordenado. Debe leer al menos este documento GeoTIFF

http://www.remotesensing.org/geotiff/spec/geotiff2.5.html#2.5.2.2 y algunas consideraciones de GDAL

http://trac.osgeo.org/gdal/wiki/rfc33_gtiff_pixelispoint .

Como regla general, todos los rásteres (aéreos, imágenes satelitales) usan píxeles por área y los datos de medición, como los DEM, usan píxeles por punto. Pero, lamentablemente, cómo se refleja en la georreferenciación y las extensiones de imagen no es tan simple. Para GDAL, el punto de anclaje del píxel siempre está en la esquina superior izquierda de un píxel; si almacena un nuevo Punto de control del suelo (GPC), las coordenadas dadas siempre están vinculadas a la esquina superior izquierda del píxel correspondiente y es importante si los píxeles se definen como áreas o puntos. GDAL también informa los mismos alcances para las imágenes pixel-is-area y pixel-is-point. Depende del usuario o software que usa GDAL tratarlos de manera diferente si es necesario.

La regla actual de OGC para la georreferenciación de coberturas es utilizar una cuadrícula que se define en la especificación GML v.3.2.1 de 437 páginas. Eso hace una diferencia de medio píxel en la ubicación de origen para píxel es punto y píxel es área porque en el caso de píxeles de área, el origen se establecerá en el medio de un píxel. La siguiente cita es de GML en el estándar JPEG2000:

Nota: GMLJP2 sigue la definición de cuadrículas en GML 3.2.1 [OGC 07-036] cláusula 19.2.2: “Cuando un punto de cuadrícula se utiliza para representar un espacio de muestra (por ejemplo, píxel de imagen), el punto de cuadrícula representa el centro de la cuadrícula espacio muestral (ver ISO 19123: 2005, 8.2.2) ". Esto corresponde con el valor pixelInCell de ImageCRS establecido en CellCenter como se especifica en ISO 19111. Esto se puede interpretar como el origen de RectifiedGrid es el punto central del píxel de la esquina.

Para hacer que las cosas no parezcan demasiado simples, algunos programas como ERDAS imagine (AFAIK) siguen internamente el esquema de píxel es el punto, así como MapServer hasta la próxima versión 7.0, con el intento de cambiar sobre la marcha en El esquema de píxel es área. Leyendo sobre esto: http://mapserver.org/fr/development/rfc/ms-rfc-107.html .

Para una mejor interoperabilidad, recomendaría utilizar píxeles por área y GeoTIFF para todos los rásteres naturales (imágenes aéreas y satelitales, mapas topográficos). Esta es una combinación bien conocida y por lo tanto confiable. Si juega con DEM y otros datos de medición / sensor y el desplazamiento de medio píxel puede ser significativo, tenga en cuenta que pueden etiquetarse como píxel por punto y que diferentes programas pueden introducir o no los desplazamientos de medio píxel.

Ambos son comunes, y ninguno puede considerarse completamente estándar.

Para GeoTIFF, ambos son posibles; consulte la Sección 2.5.2.2 para conocer el GTRasterTypeGeoKeymétodo de interpretación. Las preguntas frecuentes de GeoTIFF sugieren usar el valor predeterminado ( PixelIsArea) de esa etiqueta para la compatibilidad con versiones anteriores de GDAL.

El formato de archivo mundial utiliza el centro del píxel.

El recorte de imagen de NITF (usando la extensión ICHIPB) codifica explícitamente la posición del primer píxel en la imagen fragmentada (subconjunto) como coordenada (0.5. 0.5), y la asigna de nuevo a una coordenada específica en la imagen principal (superconjunto), por ejemplo ( 2.5, 1.5).

No sé si existe una convención, pero las coordenadas en las imágenes normalmente se refieren a centros de píxeles. Esto puede diferir en el caso de rásteres producidos a partir de, por ejemplo, modelos de computadora.

Si no está seguro, puede usar GMT , las herramientas de mapeo genérico, para probar sus rásteres y convertirlos a algún formato apropiado, ya que ofrece explícitamente la opción de especificar el registro de cuadrícula o nodo. Consulte esta sección del tutorial y también el manual de la función grdedit .