QGIS "ERROR 6: no hay traducción para Lambert_Conformal_Conic al formato PROJ.4 es conocida."

Question

Estoy tratando de clip de un condado DEM archivo utilizando el PLSS forma de Sección archivos como capas de máscara... es decir, yo quiero que recorta las secciones, de modo que puedo cargar mi DEM en QGIS una sección a la vez, sin tener que cargar todo 4GB de trama para cualquier pequeño proyecto que quiero trabajar.

He intentado hacer esto de ir a Raster > Extracción > Clipper ... sin embargo, se devuelve un error:

ERROR 6: No hay traducción para Lambert_Conformal_Conic a PROJ.4 formato es conocido. ERROR 6: No hay traducción para Lambert_Conformal_Conic a PROJ.4 formato es conocido.

He encontrado una forma para trazar la forma de archivos y el uso de la calculadora ráster, sin embargo, esto es más de un dolor y quiero entender por qué esto no funciona para mí.

Supongo que tiene algo que ver con el hecho de que estoy usando un custom CRS para mi condado? CRS bases de datos no lo reconoce.

+proj=lcc +lat_1=43.57503293972223 +lat_0=43.57503293972223 +lon_0=-90.78333333333333 +k_0=1.0000408158 +x_0=222504.44500889 +y_0=47532.0603505207 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +unidades=us-ft +no_defs

Ejecución de QGIS 2.0.1 en el Macbook Air con Mavericks 10.9

Cualquier ayuda es apreciada. Gracias!

EDIT: me encontré gdalsrsinfo en el archivo raster y la línea de comandos devuelve el mismo error, además de la WKT:

ERROR 6: No translation for Lambert_Conformal_Conic to PROJ.4 format is known.
PROJ.4 : ''

OGC WKT :
PROJCS["NAD_1983_HARN_WISCRS_Vernon_County_Feet",
    GEOGCS["GCS_North_American_1983_HARN",
        DATUM["D_North_American_1983_HARN",
            SPHEROID["GRS_1980",6378137.0,298.257222101]],
        PRIMEM["Greenwich",0.0],
        UNIT["Degree",0.017453292519943295]],
    PROJECTION["Lambert_Conformal_Conic"],
    PARAMETER["False_Easting",730000.0],
    PARAMETER["False_Northing",155944.768],
    PARAMETER["Central_Meridian",-90.78333333333333],
    PARAMETER["Standard_Parallel_1",43.57503293972223],
    PARAMETER["Scale_Factor",1.0000408158],
    PARAMETER["Latitude_Of_Origin",43.57503293972223],
    UNIT["Foot_US",0.30480060960121924],
    VERTCS["NAVD_1988",
        VDATUM["North_American_Vertical_Datum_1988"],
        PARAMETER["Vertical_Shift",0.0],
        PARAMETER["Direction",1.0],
        UNIT["Foot_US",0.3048006096012192]]]

Answer 1

Esto parece ser un error conocido en GDAL de trabajo en el fondo de QGIS acerca de las diferentes formas de definición de LCC 1-SP y 2-SP proyecciones:

Ver

• https://github.com/mapbox/node-srs/issues/1
• http://lists.osgeo.org/pipermail/gdal-dev/2011-May/028728.html
• http://trac.osgeo.org/gdal/ticket/2072

Añadiendo ESRI:: en un gdalwarp de la línea de comandos como se describe debe resolver el problema. O cambiar los archivos de forma .prj a

PROJCS["NAD_1983_HARN_WISCRS_Vernon_County_Feet",
GEOGCS["GCS_North_American_1983_HARN",
    DATUM["D_North_American_1983_HARN",
            SPHEROID["GRS80",6378137,298.257222101],
            TOWGS84[0,0,0,0,0,0,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_1SP"],
    PARAMETER["latitude_of_origin",43.57503293972223],
    PARAMETER["central_meridian",-90.78333333333333],
    PARAMETER["scale_factor",1.0000408158],
    PARAMETER["false_easting",730000],
    PARAMETER["false_northing",155944.768],
    UNIT["Foot_US",0.3048006096012192]]

y ejecutar gdal_translate -a_srs new.prj src_dataset dst_dataset sobre el archivo raster.

Tenga en cuenta que los falsos Coordenada y coordenada y en los pies para el .prj archivo, y siempre en metros para el +proj cadena.

---

ACTUALIZACIÓN

He encontrado el DEM en http://relief.ersc.wisc.edu/wisconsinview/form.php y PLSS de datos en http://www.geocommunicator.gov/GeoComm/lsis_home/home/PLSS_download_WIcounty.htm (NAD 83)

El PLSS también está impreso en el USGS Topo mapas. Así que me tomó un municipio (T13N R02W), y se vuelve a proyectar el shapefile Vernon WCCS con los parámetros indicados anteriormente. El DEM viene en formato adf, por lo que he traducido que a Vernon WCCS así:

gdal_translate -a_srs "+proj=lcc +lat_1=43.57503293972223 +lat_0=43.57503293972223 +lon_0=-90.78333333333333 +k_0=1.0000408158 +x_0=222504.44500889 +y_0=47532.0603505207 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs" -of GTiff D:/Karten/shp/Laender/USA/vernon/w001000.adf D:/Karten/shp/Laender/USA/Vernon-DEM.tif

gdalwarp -q -cutline D:/Karten/shp/Laender/USA/T13NR02W.shp -crop_to_cutline -of GTiff D:/Karten/shp/Laender/USA/Vernon-DEM.tif D:/Karten/shp/Laender/USA/Vernon-DEM-clipped.tif

Estas son en realidad las líneas de comandos generada por QGIS. El resultado se ajusta para el municipio de medida:

De fondo es el de La Farge Topo de 1965.