6 votos

Volviendo a proyectar MODIS-Esfera de CUADRÍCULA de ESRI archivo WGS84 utilizando QGIS/GDAL

Tengo un archivo de CUADRÍCULA de ESRI proyectado en MODIS-Esfera (aquí el conjunto de datos http://urban.yale.edu/sites/default/files/data/setoguneralphutyra_uf.zip) que necesito para convertir a un WGS84 GeoTIFF.

Traté de Raster > Proyecciones > Deformar (Reproyectar) y me esta tontería: QGIS_Output

Pensé usar el MODIS de reproyección herramienta o HEG para convertir, pero sólo aceptan la OL de los archivos y utilización de QGIS' urdimbre o traducir las funciones para convertir este archivo a HDF throwns un error. Este chico estaba tratando de hacer lo mismo en el Arco, ni alegría: MODIS-Esfera a WGS84 sistema de coordenadas en ArcGIS

Miré en busca de pistas en los metadatos, y el archivo parece un poco tonto; gdalinfo rendimientos:

Driver: AIG/Arc/Info Binary Grid
Files: seto_uf
       seto_uf.aux
       seto_uf/.aux.xml
       seto_uf/.DS_Store
       seto_uf/dblbnd.adf
       seto_uf/hdr.adf
       seto_uf/log
       seto_uf/metadata.xml
       seto_uf/sta.adf
       seto_uf/vat.adf
       seto_uf/w001001.adf
       seto_uf/w001001.tif
       seto_uf/w001001.tif.aux.xml
       seto_uf/w001001x.adf
Size is 7468, 2911
Coordinate System is `'
Origin = (-17345341.010356999933720,8423784.764957003295422)
Pixel Size = (5000.000000000000000,-5000.000000000000000)
Corner Coordinates:
Upper Left  (-17345341.010, 8423784.765) 
Lower Left  (-17345341.010,-6131215.235) 
Upper Right (19994658.990, 8423784.765) 
Lower Right (19994658.990,-6131215.235) 
Center      ( 1324658.990, 1146284.765) 
Band 1 Block=256x4 Type=Byte, ColorInterp=Undefined
  Description = Seto_UF
  Min=0.000 Max=101.000 
  NoData Value=255

De modo que no puede determinar el sistema de coordenadas.

Hay una .tif incluido en el directorio, pero gdalinfo que muestra:

Driver: GTiff/GeoTIFF
Files: w001001.tif
       w001001.tif.aux.xml
Size is 7468, 2911
Coordinate System is:
GEOGCS["Unspecified datum based upon the GRS 1980 Authalic Sphere",
    DATUM["unknown",
        SPHEROID["unnamed",6371007,1]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4047"]]
Origin = (-17345341.010356999933720,8423784.764957003295422)
Pixel Size = (5000.000000000000000,-5000.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (-17345341.010, 8423784.765) 
Lower Left  (-17345341.010,-6131215.235) 
Upper Right (19994658.990, 8423784.765) 
Lower Right (19994658.990,-6131215.235) 
Center      ( 1324658.990, 1146284.765) 

Así se dice que las unidades son grados, pero los de la esquina enorme coordenadas y el origen de los valores no están de acuerdo.

Alguna idea?

Mi configuración:
OSX 10.9.4
QGIS 2.6.1 Brighton
GDAL 1.11.2 (al parecer hay nuevas versiones pero esta es la más reciente que usted puede conseguir en OSX sin tener que compilarlo usted mismo)

8voto

SpliFF Puntos 214

El Readme.docx dentro del archivo zip da algunas pistas sobre la proyección, pero tal vez engañosa.

La trama es en una interrumpido Goode Homolosine de proyección, basada en la MODIS de la esfera.

Así que he creado el siguiente personalizado CRS:

+proj=igh +R=6371007 +no_defs

se vuelve a proyectar una Tierra Natural shapefile en ella, y se ajusta a los datos ráster con el mismo CRS asignado (no se vuelve a proyectar!):

enter image description here

Así que una vez que asigne la costumbre de CRS a la trama, se puede establecer que el proyecto de CRS sobre la marcha de la Web Mercator o algo:

enter image description here

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X