Tengo unos 800 archivos geoTIFF que tienen una proyección no especificada, pero como conozco la proyección (WGS 84 UTM Zona 32N), he utilizado la función "definir proyección" para definir por lotes las proyecciones de estos archivos. Ahora necesito recortar todos estos archivos a una nueva extensión utilizando un shapefile con la misma proyección. Sin embargo, cuando he comprobado los archivos he visto que no se solapan. Mi shapefile aparece en la parte superior de la pantalla y los geoTIFFs son pequeños puntos en el centro de la pantalla. En realidad, son tan pequeños que no puedo verlos en la vista global, pero sé que están ahí. Los datos tienen casi la misma extensión, así que deberían tener el mismo tamaño.
Los datos de referencia espacial y los datos de extensión se muestran a continuación para cada tipo de archivo. La diferencia que veo está en las extensiones, pero no sé cómo solucionar este problema, o incluso si éste es el origen de mi problema. ¿Alguien sabe qué puede estar mal y cómo hacer que mis archivos se superpongan?
geoTIFF file
Spatial reference
XY Coordinate System
WGS_1984_UTM_Zone_32N
Linear Unit Meter (1.000000)
Angular Unit Degree (0.0174532925199433)
False_Easting 500000
False_Northing 0
Central Meridian 9
Scale_Factor 0.9996
Latitutde of Origin 0
Datum D_WGS_1984
Extent Top 6720.20271636
Left 370.036258324
Right 461.540519838
Bottom 6651.1691216
Shapefile
WGS_1984_UTM_Zone_32N
WKID: 32632
Authority: EPSG
Projection: Transverse_Mercator
False_Easting: 500000.0
False_Northing: 0.0
Central_Meridian: 9.0
Scale_Factor: 0.9996
Latitude_Of_Origin: 0.0
Linear Unit: Meter (1.0)
Geographic Coordinate System:
GCS_WGS_1984
Angular Unit: Degree (0.0174532925199433)
Prime Meridian: Greenwich (0.0)
Datum: D_WGS_1984
Spheroid: WGS_1984
Semimajor Axis: 6378137.0
Semiminor Axis: 6356752.314245179
Inverse Flattening: 298.257223563
Extent
Top 6733114.02831481
Left 364349.954284994
Right 501090.000365895
Bottom 6620131.02208354
3 votos
¿Estás seguro de que los Geotiff están en el sistema de referencia de coordenadas correcto? Si es así, tienen una extensión de ~69 m x 91 m... Parece que la proyección está mal definida
0 votos
Todo lo que he hecho a los tiffs originales es un lote "definir proyección" al sistema de coordenadas WGS_1984_UTM_Zone_32N. La información en el texto anterior es directamente de Arc Catalog, así que creo que es correcta?
2 votos
Sí, pero usted conoce los datos... ¿La extensión de sus rastreos es coherente con los datos? Parece que son demasiado pequeños. Para las ortofotos podría estar bien. Por otro lado, tu shapefile tiene una extensión de ~100.000 m x 136.000 m, completamente diferente de las extensiones de tus rasters. Por esta diferencia, creo que
WGS_1984_UTM_Zone_32Nno representa sus datos0 votos
Sé que los datos están en UTM 32N, así que ¿cómo puede estar mal definida la proyección si la defino como tal? No sé cómo decirle a ArcGIS cuál es la extensión - ¿no debería estar esta información incrustada en el archivo GeoTIFF?
0 votos
Sí, la extensión es demasiado pequeña, pero no sé por qué.
3 votos
Hay muchos mensajes de "las extensiones no coinciden" aquí, y casi todos ellos se refieren a una respuesta de "tu proyección no es lo que pensabas que era". No tenemos tus datos, así que no podemos hacer ninguna investigación, pero tienes que considerar la posibilidad de que al menos uno de tus conjuntos de datos esté en una referencia de coordenadas incorrecta. Tenga en cuenta que hay un número infinito de proyecciones UTM 32N debido a la variación del origen falso.
0 votos
Ya he comprobado dos y tres veces la proyección/coordenadas de referencia, pero lo haré de nuevo. A qué te refieres con que "hay un número infinito de proyecciones UTM 32N debido a la falsa variación del origen".
0 votos
No puedo entrar en la teoría de la proyección en quinientos caracteres. Basta con decir que 500000 no es el único posible False.Easting. Si tomas una capa geográfica conocida en WGS84, extraes un subconjunto de condados y lo reproyectas a lo que crees que deberían ser tus datos, entonces sabrás si son tus imágenes o tu shapefile los que están mal definidos (a menos que sean ambos)
2 votos
¡Ah! Creo que los geotiffs están en kilómetros. Vas a tener que hacer un UTM 32N personalizado usando kilómetros (asegúrate de que el valor de la falsa easting se convierte en 500,0) o [...]
0 votos
He conseguido solucionar el problema. Sólo tuve que multiplicar la extensión por 1000 y los tiffs se superpusieron perfectamente con el shapefile. Lo hice en R utilizando la función extent que forma parte del paquete raster. ¡Gracias!