Hello, I'm playing around with geodjango and this is a GIS question but I'm starting my quest for an answer by asking here.
So have a shp file I'm using from the city of chicago: http://egov.cityofchicago.org/webportal/COCWebPortal/COC_ATTACH/CitywideWardsMap2008.zip There is a wards.shp file bundled in there that you can successfully import into a django model. Using ogrinfo, for example, I can see the polygon information looks like: ogrinfo -ro -al wards.shp POLYGON ((1160473.005399998277426 1921203.864600002765656,1160474.219999998807907 ... And the layer description is: ogrinfo -ro -so apps/wards/data/wards.shp wards INFO: Open of `apps/wards/data/wards.shp' using driver `ESRI Shapefile' successful. Layer name: wards Geometry: Polygon Feature Count: 53 Extent: (1091130.772400, 1813891.890100) - (1205199.881775, 1951669.020100) Layer SRS WKT: PROJCS["NAD_1983_StatePlane_Illinois_East_FIPS_1201_Feet", GEOGCS["GCS_North_American_1983", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980",6378137.0,298.257222101]], PRIMEM["Greenwich",0.0], UNIT["Degree",0.0174532925199433]], PROJECTION["Transverse_Mercator"], PARAMETER["False_Easting",984250.0], PARAMETER["False_Northing",0.0], PARAMETER["Central_Meridian",-88.33333333333333], PARAMETER["Scale_Factor",0.999975], PARAMETER["Latitude_Of_Origin",36.66666666666666], UNIT["Foot_US",0.3048006096012192]] OBJECTID: Integer (10.0) DATA_ADMIN: Real (19.8) PERIMETER: Real (19.8) WARD: String (4.0) ALDERMAN: String (60.0) CLASS: String (2.0) WARD_PHONE: String (12.0) HALL_PHONE: String (12.0) HALL_OFFIC: String (45.0) ADDRESS: String (39.0) EDIT_DATE1: String (10.0) SHAPE_AREA: Real (19.11) SHAPE_LEN: Real (19.11) But when I use an example from the geodjango tutorial ( http://geodjango.org/docs/tutorial.html#id4 ) the numbers are in a different format and resemble latitude/longitude. ogrinfo -ro -so apps/world/data/TM_WORLD_BORDERS-0.3.shp TM_WORLD_BORDERS-0.3 INFO: Open of `apps/world/data/TM_WORLD_BORDERS-0.3.shp' using driver `ESRI Shapefile' successful. Layer name: TM_WORLD_BORDERS-0.3 Geometry: Polygon Feature Count: 246 Extent: (-180.000000, -90.000000) - (180.000000, 83.623596) Layer SRS WKT: GEOGCS["GCS_WGS_1984", DATUM["WGS_1984", SPHEROID["WGS_1984",6378137.0,298.257223563]], PRIMEM["Greenwich",0.0], UNIT["Degree",0.0174532925199433]] FIPS: String (2.0) ISO2: String (2.0) ISO3: String (3.0) UN: Integer (3.0) NAME: String (50.0) AREA: Integer (7.0) POP2005: Integer (10.0) REGION: Integer (3.0) SUBREGION: Integer (3.0) LON: Real (8.3) LAT: Real (7.3) So I think what I need to use is the Transverse_Mercator projection but I don't know how. Given a lat/long I would like to calculate what ward it is in. Any ideas on how to accomplish this calculate a lat/long against this specific projection or vise versa? Thanks, Milan --~--~---------~--~----~------------~-------~--~----~ You received this message because you are subscribed to the Google Groups "Django users" group. To post to this group, send email to [email protected] To unsubscribe from this group, send email to [email protected] For more options, visit this group at http://groups.google.com/group/django-users?hl=en -~----------~----~----~----~------~----~------~--~---
