¿Cuál es el SRID de census.gov shapfiles?

Question

He descargado el condado de archivos de:

http://www.census.gov/cgi-bin/geo/shapefiles2011/main

No puedo averiguar qué SRID a utilizar con ellos. Cuando estoy tratando de generar un archivo de clase de ./manage.py ogrinspect ..., me sale un SRID de -1 en el polígono de campo. No creo que esto es correcto.

GDAL me dio la Proj4 de +proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs.

Sé que un SRID de -1 significa desconocido. No se que hacer mi geometrías inútil? ¿Cómo puedo averiguar qué SRID a utilizar?

Soy realmente nuevo en el SIG. Yo ya olvida los días en que la Tierra era plana.

Esto también puede ser útil:

SELECT PostGIS_full_version();
                                         postgis_full_version                                          
-------------------------------------------------------------------------------------------------------
 POSTGIS="1.5.3" GEOS="3.2.2-CAPI-1.6.2" PROJ="Rel. 4.7.1, 23 September 2009" LIBXML="2.7.8" USE_STATS
(1 row)

---

>>> import django
>>> django.VERSION
(1, 4, 0, 'final', 0)

Answer 1

El primer paso para determinar la correcta proyección de cualquier capa, es encontrar la información de proyección, si los hubiere, que vino con su capa. En el caso de un Shapefile, como lo has descargado desde Census.govque la información está contenida en un .prj archivo, corto para la Proyección.

Aquí está el contenido de la proyección de archivo de los datos del censo:

GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]

Las partes clave de este, son la primera palabra GEOGCS significa Geographic Coordinate System, la cadena de DATO: D_North_American_1983, SPHERIOD: GRS_1980, y la última parte de la UNIDAD: Degree.

Estos coincidir fácilmente hasta el Proj4 definición que has recibido: +proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs

Se puede ver que el datum y elipsoide, partido hasta el punto de referencia y spheriod en el ".prj". El longlat medio Longitude/Latitude, la cual se mide en Unidades de Degrees, lo que coincide con el ".prj" archivo.

No hay "defs" porque no hay información adicional es necesaria para identificar el Sistema de Coordenadas Geográficas.

Un sitio para buscar referencias espaciales es: SpatialReference.org

El Sistema de Coordenadas Geográficas en NAD 83 es difícil llegar directamente, porque es una parte de muchas otras proyecciones.

Para hacerlo más fácil, lo que usted necesita es EPSG 4269.

Este enlace presenta la información de referencia espacial en un número de diferentes formatos. El uno específicamente para PostGIS es:

INSERT into spatial_ref_sys (srid, auth_name, auth_srid, proj4text, srtext) valores ( 94269, 'epsg', 4269, '+proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs ', 'GEOGCS["NAD83",DATO["North_American_Datum_1983",ESFEROIDE["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4269"]]');

Para abordar su último punto sobre el SRID = -1, creo que es importante comprobar los datos antes de que se supone que ellos son inútiles. Tratar de llevar hasta en algunos programas de SIG, como QGIS, o ArcGIS. Mira las coordenadas que se muestran como mover el cursor alrededor. Si ves familiar rangos de números, como y = 0 – +90 y x = -180 a +180, usted debe reconocer estos como Latitud/Longitud. Comenzando a reconocer las coordenadas de los rangos que se producen en el común de las proyecciones y sistemas de coordenadas en su área le ayudará enormemente a medida que avance en su carrera.

Si usted está interesado en aprender más acerca de las proyecciones y sistemas de coordenadas, esta es una fantástica referencia: Un Manual de Trabajo (PDF) - Juan P. Snyder

Answer 2

@Obtener Espacial que es correcto y que no podía haber sido mejor respuesta a esta pregunta. además de él, usted puede comprobar fuera de spatial_ref_sys antes de agregar cualquier referencia espacial.

SELECT * FROM
       spatial_ref_sys WHERE auth_srid = 4269;

probablemente usted va a tomar algo como esto:

SRID  | AUTH_NAME  | AUTH_SRID  | SRTEXT
------+------------+------------+-------
4269  |   EPSG     |   1624     | <big SRTEXT string>

SRTEXT = GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 
         1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],
         PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,
         AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4269"]]

PROJ4TEXT = +proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs

Spatial_ref_sys es simplemente una tabla de coordenadas conocido de los sistemas de llave por su SRID. usted cabn obtener más información aquí...

• srid: El número de SRID. Esta debe ser la clave principal de la tabla.
• auth_name: Un nombre de la autoridad como una cadena. Esta es establecer si este sistema de coordenadas especificado por una autoridad exterior, tales como la EPSG.
• auth_srid: El ID numérico del sistema de coordenadas en el por encima de la autoridad del catálogo.
• srtext: El Conocido Texto (WKT) la representación del sistema de coordenadas (como se describe en la sección 4).
• proj4text: El Proj4 representación del sistema de coordenadas.

y la última cosa es que para cargar a la base de datos, use este comando:

shp2pgsql -c -D -s 4269 -I tl_2011_02_anrc.shp myschema.mytable > mysql.sql
psql -d mydb -f mysql.sql

espero te sirva de ayuda....