3
votes
0answers
165 views

Reprojecting WGS 1984 Web Mercator (EPSG:3857) in Python with GDAL

I am reprojecting rasters in python using GDAL. I need to project several tiffs from geographic WGS 84 coordinates to WGS 1984 Web Mercator (Auxiliary Sphere), in order to use them later in Openlayers ...
1
vote
2answers
65 views

Custom BatchProject Script

Here's my problem. My GPS collects data in WGS1984 and I have decided to start working in UTM instead of state plane. Needless to say, I would edit shapefiles, not paying attention to what ...
2
votes
3answers
156 views

How to find and apply reprojection parameters?

I have a data set in one projection, and a shapefile in another, and I am trying to combine the two. I am using Python. I apologize for posting a "help me transform coordinates" question on this site ...
0
votes
0answers
46 views

Combining two separate shapefiles with undefined projections out of Leica

So this is a mess. When the software outputs I have two shapefiles. The projections aren't defined, but I know one is WGS84 and the other is NAD 83. The NAD 83 includes northing easting and ...
1
vote
1answer
74 views

Different output from running algorithm in toolbox and from script

I am unable to get the same results when running the qgis:reproject tool from QGIS processing toolbox (QGIS version 2.0.1) as when running it from a Python script (console). In this example I am ...
2
votes
2answers
108 views

No prj file output from qgis:reproject

I'm trying to reproject a vector layer using the 'qgis:reproject' algorithm in the processing/sextante module of QGIS, from Python. This kind of works, but the .prj and .qpj files created after ...
0
votes
2answers
427 views

Tile script (Mapnik) with OpenLayers - Problem with projection (900913)

I use Mapnik python bindings in a geodjango application to render tiles in openlayers with datasource stored in postGIS in 900913 projection. The feature seems to be rendered correctly in the ...
0
votes
1answer
142 views

Orthographic reprojection with mapnik

I am trying to create a google earth like wep application. My problem is that I can't get to render a full view of the earth. The zoom_all() function generates the following error : RuntimeError: ...
3
votes
2answers
371 views

What's the best projection for rasterization of random latitude and longitude data in the northern Atlantic?

I have a dataset in the WGS84 geographic coordinate system that I would like to interpolate along a grid using Python. It appears gdal_grid is not available through the Python bindings, so I plan to ...
1
vote
1answer
264 views

How to re-project the EASE (Equal Area Scalable Earth) grid with a ~25 km cylindrical projection to WGS84 0.25 degree?

I have nc files ,here is one file https://echange-fichiers.inra.fr/get?k=LUfqiyO1Dl7az5bihnS. from the metadata ,the projection is cylindrical and the resolution is 25 km: ...
2
votes
1answer
1k views

Mapping x, y pixel coordinates to latitude and longitudes on the WGS84 datum

Apologies in advance for the long question, I am fairly confused at this point and I figure if I show my thought process any flaws early on in the piece can be more easily pointed out. I have a grid ...
6
votes
1answer
874 views

Python script to convert Lat Long to ITM (Irish Transverse Mercator)

I am currently working on a python script that converts Lat Long (WGS84) to ITM. The formulas used are based on a paper published by OSGB (Appendix B): ...
4
votes
1answer
714 views

Reproject script in Python with GDAL

I'm having a lot of trouble with GDAL. Besides just the at times lacking documentation, there seems to be little support in Python. Anyways this is just a batch reproject script that goes thusly: ...
2
votes
1answer
141 views

Any reason NOT to convert to a single projection when loading GIS data?

I have been given the task of: Loading a list of properties with lat/long coordinates into a SQL Server 2008 database. Loading points and areas (polygons) of interest from a variety of shapefiles ...
2
votes
1answer
162 views

Converting 2D shapefile to Equal Area Projection?

I have thousands of shapefiles with coordinates in lat/long (-180 to +180 x, -90 to +90 y), that contain polygons of which I wish to calculate the area. I'm using the Python GEOS/OGR/GDAL API. If I ...
0
votes
4answers
466 views

How to unproject a shapefile in Python?

I have one shapefile with XY coordinates in meters. Now, I want to develop one script in Python which converts meters to decimal degrees. Please help me to develop this. Thanks
5
votes
4answers
1k views

Changing image projection using Python

Image corner coordinates are known: Top left: (70.20315, 20.92749) Top right: (70.20315, 31.45240) Bottom left: (59.70571, 20.92749) Bottom right: (59.70571, 31.45240) (Basically map over Finland ...
5
votes
2answers
5k views

Converting X,Y coordinates to lat/long using pyproj and Proj.4 returns the wrong coordinates

I'm writing a python script that reads multiple XML files containing x and y coordinates and combines them all into a single csv file. Latitude and Longitude are required fields in the csv, but I am ...
9
votes
6answers
5k views

ArcGIS/ArcPy Python, how to get latitude and longitude of a projected point?

I have a point feature in a feature class that is being accessed by ArcPy. The point is projected but I need to find an efficient means to get the unprojected latitude and longitude for this point. ...
3
votes
1answer
443 views

With CTM how can I get the latitude and longitude bounding box of a geoPDF?

Embedded in a geoPDF is info on the CTM - Coordinate Transformation Matrix. I am unfamiliar with this CTM. I know how to use python-numpy. How (or is it even possible) for me to use the CTM ...