Skip to main content
Search type Search syntax
Tags [tag]
Exact "words here"
Author user:1234
user:me (yours)
Score score:3 (3+)
score:0 (none)
Answers answers:3 (3+)
answers:0 (none)
isaccepted:yes
hasaccepted:no
inquestion:1234
Views views:250
Code code:"if (foo != bar)"
Sections title:apples
body:"apples oranges"
URL url:"*.example.com"
Saves in:saves
Status closed:yes
duplicate:no
migrated:no
wiki:no
Types is:question
is:answer
Exclude -[tag]
-apples
For more details on advanced search visit our help page
Results tagged with
Search options answers only not deleted user 2856

projection description, using a well-known text representation of coordinate reference systems

3 votes
Accepted

Reading coordinates from shapefile using shapely

You can reproject your points from EPSG:4326 to EPSG:2953 (NAD_1983_CSRS_New_Brunswick_Stereographic) using pyproj like so: import shapefile import pyproj from shapely.geometry import (Point, shape) # …
user2856's user avatar
  • 69.6k
6 votes
Accepted

Why ogr doesn't read this kind of prj

The WKT string in that .prj is corrupt. I don't know what wrote it, but the PROJCS and UNITS elements are repeated. i.e. … from osgeo import osr srs = osr.SpatialReference() prj='PROJCS["WGS 84 / UTM zone 24N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG …
user2856's user avatar
  • 69.6k