La première étape pour déterminer la projection correcte de n'importe quel calque consiste à trouver les informations de projection, le cas échéant, fournies avec votre calque. Dans le cas d'un Shapefile
, comme ce que vous avez téléchargé à partir de Census.gov , cette information est contenue dans un .prj
fichier, abréviation de Projection.
Voici le contenu du fichier de projection des données du recensement:
GEOGCS ["GCS_North_American_1983", DATUM ["D_North_American_1983", SPHEROID ["GRS_1980", 6378137,298.257222101]], PRIMEM ["Greenwich", 0], UNIT ["Degree", 0.01745329251994
Les éléments clés de ce sont les premiers mots GEOGCS moyens Geographic Coordinate System
, la chaîne DATUM : D_North_American_1983
, SPHERIOD : GRS_1980
, et la dernière partie UNITÉ : Degree
.
Ceux-ci correspondent facilement à la définition Proj4 que vous avez reçue:+proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs
Vous pouvez voir que le datum et l'ellipsoïde correspondent au datum et à la période dans le ".prj". Les longlat
moyens Longitude/Latitude
, qui est mesuré en unités de Degrees
, qui correspond au fichier « .prj ».
Il n'y a "aucune définition" car aucune information supplémentaire n'est nécessaire pour identifier le système de coordonnées géographiques.
Un site pour rechercher des références spatiales est: SpatialReference.org
Le système de coordonnées géographiques dans NAD 83 est difficile à trouver directement, car il fait partie de nombreuses autres projections.
Pour le rendre plus facile, celui dont vous avez besoin est EPSG 4269 .
Ce lien présente les informations de référence spatiale dans un certain nombre de formats différents. Celui spécifiquement pour PostGIS est:
INSÉRER dans les valeurs spatial_ref_sys (srid, auth_name, auth_srid, proj4text, srtext) (94269, 'epsg', 4269, '+ proj = longlat + ellps = GRS80 + datum = NAD83 + no_defs', '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É [" degré ", 0,01745329251994328, AUTORITÉ [" EPSG "," 9122 "]], AUTORITÉ [" EPSG "," 4269 "]] ');
Pour répondre à votre dernier point concernant le SRID = -1, je pense qu'il est important de tester vos données avant de supposer qu'elles sont inutiles. Essayez de le faire apparaître dans certains logiciels SIG, comme QGIS ou ArcGIS . Regardez les coordonnées qui s'affichent lorsque vous déplacez le curseur. Si vous voyez des plages de nombres familières, comme y = 0 - +90 et x = -180 - +180, vous devez les reconnaître comme Latitude / Longitude. Commencer à reconnaître les plages de coordonnées qui se produisent dans les projections et les systèmes de coordonnées communs dans votre région vous aidera grandement à mesure que vous avancerez dans votre carrière.
Si vous souhaitez en savoir plus sur les projections et les systèmes de coordonnées, voici une référence fantastique: un manuel de travail ( PDF ) - John P. Snyder