recortar ráster a polígono con 3 mallas

Question

Tengo 70 imágenes landsat rasterizadas con srid: 32647,32648,32649 y quiero recortar raster a polígono (shapefile, srid: 4326) en postgis.

He importado todos los raster en 1 tabla con el parámetro -s : 32648, transforme el srid del shapefile a 32648 pero el resultado no es como esperaba ( solo se superpone el 40% del polígono)

¿Cómo puedo solucionar este problema?

Answer 1

Suponiendo que haya utilizado raster2pgsql para importar los datos, -s es NO la opción correcta para reproyectar los datos en un único SIR.

Ha asignado un nuevo srid, pero no ha recalculado las coordenadas: http://postgis.net/docs/manual-2.2/using_raster_dataman.html

En su lugar, utilice gdalwarp -s_srs old_srs -t_srs new_srs srcfile destfile para reproyectar dos de sus rasters al CRS del tercero. A continuación, puede importar todos ellos en una tabla.

Answer 2

Intente importar sus rásteres sin la opción -s. raster2pgsql debería asignarles automáticamente el SRID correcto. Puede comprobarlo después de la carga con ST_SRID(rast). Una vez importados con su propio SRID correcto, sólo tiene que ST_Transform(geom, ST_SRID(rast)) el polígono antes de recortarlo.

Si en algún caso el importador no asigna a cada trama su SRID adecuado, puede reasignarlos con ST_SetSRID().