Questions tagged [proj]

A library for performing conversions between cartographic projections.

Filter by
Sorted by
Tagged with
36
votes
5answers
10k views

Reprojecting raster from 0 360 to -180 180 with cutting 180 meridian using gdalwarp

I have a geotiff raster image that has a coordinate system with longitudes from 0 to 360. The horizontal center of the image is 180 longitude. See image below: I want to transform it to EPSG:4326 SRS ...
6
votes
1answer
15k views

Assigning CRS to shapefile when it doesn't have one, in R?

I think that by using readOGR to import a shapefile renders worthless to put proj4string to the script since if any coords are available in the shapefile, it will be transported through the readOGR. ...
10
votes
3answers
21k 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 ...
15
votes
1answer
3k views

Defining Winkel Tripel in proj 4.8.0?

I have QGIS with proj 4.8.0 installed. I'd like to add Winkel Tripel projection which seems to be defined as +proj=wintri But that definition string does not work. Is this a bug in my proj version, ...
6
votes
1answer
5k views

R-4.0.1 Not sure I understand this message: Warning message: In proj4string(x) : CRS object has comment, which is lost in output [closed]

Context: I want to store the TIFF data to use it for running a linear regression model later on. Ultimate goal: Create a carbon map using field data as a basis. Libraries: library(sp) library(raster) ...
4
votes
2answers
9k views

How to get map units to find current scale in OpenLayers?

I need my OpenLayers application to be able to determine the current scale. This is not for display (there is already the scale line for that). This is because I have a scale-dependent style ...
4
votes
3answers
4k views

Leaflet proj4leaflet with multiple projections

I am looking at having both EPSG4326 tiles and EPSG3857(cloudmade) tiles on a single map (but separate layers). I have been looking at the proj4leaflet plugin, and it "appears" to me that this might ...
15
votes
2answers
33k views

Choosing the correct value for proj4string for shapefile reading in R?

I am having a shapefile of polygons and another CSV file which contains a list of points as (Lat, Lng) pairs.. I want to check for each (lat, lng) pair from the CSV file which polygon does it fall ...
14
votes
4answers
8k views

What is the best way to programmatically convert between WKT and Proj4 string?

Some shapefiles have a .prj file associated with it, and the .prj file contains the projection info of the shapefile in the format of WKT. Sometimes I need to convert WKT to proj4 string, and ...
7
votes
2answers
4k views

Transforming geostationary satellite image to lon/lat

The ICARE Data and Services Center distributes the SEV_AERUS-AEROSOL-D3 product in HDF5 format (http://www.icare.univ-lille1.fr/archive/?dir=GEO/MSG+0000/SEV_AERUS-AEROSOL-D3.v1.03/). The main ...
6
votes
2answers
7k views

How to perform transformations using ProjNet?

While writing some unit tests I thought it would be convenient/possible to use ProjNet to convert from EPSG:4326 into several other coordinate systems (3875,3395). I figured it would be as easy as ...
15
votes
4answers
20k views

Repairing broken GDAL and PROJ.4 on Ubuntu?

I solved the library problem, but others remain, only vaguely related to that of PROJ.4 library not loaded. I asked a follow-up question for these. My GDAL package and PROJ.4 library are damaged, ...
10
votes
4answers
4k views

Using PROJ.4 library to transform from local coordinate system coordinates to global coordinate system using ground control points?

I have a point cloud whose coordinates are with respect to a local coordinate system. I also have ground control points with GPS values. Can I convert these local coordinates to a global coordinate ...
7
votes
5answers
3k views

WKT for EPSG:4326 with LON 0 to +360 instead of -180 to +180

I have GRIB2 files that seem to be in a projection similar to EPSG: 4326, except that instead of longitude between -180 and +180 it is between 0 and 360. My intention is to get these files into ...
6
votes
1answer
2k views

Georeferencing Himawari-8 in GDAL (or other)

Does anyone know the appropriate project information to georeference (and hence reproject) Himawari-8 AHI data in the IR? The images are 5500x5500 pixels. Currently I try these: ulx=-5570248.832537 ...
3
votes
1answer
3k views

Creating custom projections with Proj4?

I'm reading into how to customise projections with Proj4 - specifically to achieve tilted projections such as the example below (from vis4.net) The Proj4 literature is pretty heavy going and having ...
6
votes
1answer
466 views

Explaining pyproj to_epsg min_confidence parameter?

Following on from this question How can I get epsg code in pyproj with version==2.1.3, I discovered that sometimes when intialising a pyproj Proj instance with a string ESPG code eventually results in ...
3
votes
3answers
1k views

Proj4 String for NAD83(2011) / Louisiana South (ftUS)

I am trying to convert NAD83(2011) / Louisiana South (ftUS) coordinates on survey maps to lat/lon using Proj4 in Python. Could someone point me in the right direction? I have verified the lat/...
1
vote
2answers
104 views

“non finite transformation detected” using sp R

This is a part of a broader question already posted (Non finite transformation detected - lidR spTransform) but I am creating a new question for the sp related part of that. I'm trying to reproject ...
14
votes
2answers
7k views

use proj4 to specify Robinson projection with R ggmap and ggplot2 packages?

I want to project this map in robinson projection: library(ggmap) world <- map_data("world") ggplot() + geom_path(data = world, aes(long, lat, group = group)) And I ...
22
votes
3answers
35k views

How can I get the proj4 string or EPSG code from a shapefile .prj file? [duplicate]

I have to work with a shapefile defined in an unusual projection, which my GIS software does not know about. How can I get the proj4 definition or the EPSG code for the projection?
2
votes
1answer
1k views

Choosing CRS / PROJ4 string for EURO-Cordex rotated pole projection?

I have to analyze some CORDEX models which have a rotated pole projection. The projection is described in this document, see page 10, Table 1, row "Europe high res". I have tried to describe the ...
8
votes
3answers
1k views

Adding Equal Earth Projection in QGIS?

There is a new accurate projection for world maps ► Equal Earth. I know its proj4 but I can't add it in QGIS (see below). Does it need some kind of QGIS update?
4
votes
0answers
2k views

How to set GDAL/cs2cs/proj to use NTv2 grid shift file when using EPSG codes?

I am using the GDAL API for Python for a script that takes only EPSG codes for inputs to specify source and target coordinate systems. While the cs2cs command seems to handle quite well ...
4
votes
1answer
320 views

proj.4 definition for (obsolete) UK War Office “Italy South” grid

I'm helping a friend convert his late father's war diary coordinates to WGS84. I have a proj.4 definition that appears to work in QGIS, but if I transform a scanned map with gdalwarp, it's not even ...
6
votes
1answer
778 views

Custom transverse Mercator projection with fiona and shapely

AS a build off from this question I posted over a year ago Buffer miles in decimal degrees with fiona and shapely I am attempting to create 50-mile buffers around each earthquake. I have taken the ...
5
votes
1answer
4k views

Assigning CRS to Lat/Lon in R giving error: “Geographical CRS given to non-conformant data”

I have a list of suburbs in Queensland, Australia along with their lat/longs. However, when I try to read in the file in R and assign a CRS, R is giving the following error message: Error in ...
3
votes
1answer
1k views

Rotating ~90° using Two-Point Equidistant projection with Proj4

I'm trying to rotate a tif raster in R using Proj4's Two-Point Equidistant projection - I think the only way to rotate spatial geometries to custom angles. But I'm finding the operation gets very ...
3
votes
1answer
5k views

Openlayers3 using custom EPSG:25832

I want to use EPSG:25832 as the default projection for all my layers. The projection is not working correctly because the Position (center point) should be in Germany but I end up in France. When I ...
2
votes
1answer
112 views

Reproject using WKT for just the horizontal CRS

I'm trying to reproject LiDAR point clouds from California zone 2 (ftUS) using spTransform from sp R package to a metric CRS like EPSG 3857 (initial post here) and got non finite transformation ...
2
votes
2answers
188 views

Multiply or Divide certain numbers on a projection file(?) of Japan

the title seems weird though, i am dealing with Japanese shapefile and it has its own "bizarre" projection. I have two hints to solve this issue. it's based on espg 4301 (http://pacificprojections....
1
vote
1answer
269 views

How to update CRS r object from PROJ4 to PROJ6

I have a script which contains several r objects which are CRSs that I use for transforming spatial objects. After updating R and many packages, there is a warning which wasn't there before. LCC <- ...
1
vote
3answers
286 views

Can someone recognize the CRS?

I'm using QGIS 1.8.0 and have some shapefiles with a "User Defined CRS" which I don't recognize. Is anyone able to recognize it? Here is an image of the settings I see on QGIS.
6
votes
1answer
7k views

Converting Lambert Conformal Conic to WGS84 using Proj.4

I am attempting to transform data from Lambert Conformal Conic lat/long and I am having some issues with the results. The shapefile I am attemtping to transform can be found here: http://tfr.faa.gov/...
5
votes
2answers
5k views

Stereographic projection of WGS84 ellipsoid on a plane[python]

I would like to ask you for help in following task. What I need to do is to derive local x,y coordinates from latitude, longitude, applying stereographic projection of WGS84 ellipsoid on a plane, when ...
9
votes
2answers
200 views

Setting proj parameters of old map

I would like to georeference these fantastic old public domains maps https://legacy.lib.utexas.edu/maps/ams/italy_50k/ I started from this one, and it has these projection parameters. Then I have ...
6
votes
3answers
3k views

Does anyone know the proj4 Coordinte Reference System command line for Rectangular Polyconic Gores?

I am a GIS beginner using QGIS 1.8.0 I am trying to transform a world map from the standard WGS84 Coordinate Projection System to a set of rectangular polyconic gores, with a gore count of 12. ...
5
votes
2answers
2k views

LSA-SAF data reprojection with GDAL and R

The Land Surface Analysis Satellite Applications Facility (LSA SAF) distributes several products derived from the EUMETSAT satellites data. I am working with R to read and process the data (the file ...
5
votes
1answer
2k views

Slight translation with coordinate-transformation from Gauss-Kruger to WGS 84. Possible Reasons?

I have a number of points in Gauss-Kruger, Zone 3 that I want to convert to display them within Google-Earth. Assuming WGS84 for Google-Earth the transformed coordinated differ from the original ones ...
3
votes
1answer
1k views

Get the geotransform in Mercator given a pixel size away from the origin

I have a radar image with the following information: The projection is Mercator (Mercator_1SP in GDAL) The projection origin is 0º lat 0º lon At a given geographic coordinates (e.g. lon=10º lat=40º), ...
10
votes
2answers
3k views

PROJ.4 - General Parameters: +no_defs

Within a QGIS-course I would like to tell a bit about the proj.4 parameters. However, I'm not sure about the "+no_defs" parameter. The proj.4 website (http://proj4.org/parameters.html#parameter-list) ...
8
votes
4answers
2k views

Proj.4/GDAL/QGIS Transformation between CRSs that are defined the same

I'm helping to ensure that open source software can handle Australia's new datum appropriately, see ICSM's website for details on the GDA2020 project. Now, QGIS has already got the definitions of ...
7
votes
1answer
2k views

projection/pyproj puzzle, and understanding SRS format

I have a shapefile that has a .prj file, with the following wkt contents: PROJCS["NAD_1927_StatePlane_New_York_West_FIPS_3103",GEOGCS["GCS_North_American_1927",DATUM["D_North_American_1927",SPHEROID["...
6
votes
1answer
1k views

How should I handle CRS properly after the major change in PROJ library?

As explained here, the PROJ library has changed a lot and one should stop using proj.4 strings. To me it is not clear, what I should do instead. This answer explains, how I could deal with the problem ...
6
votes
2answers
13k views

Basic proj4 conversion UTM grid -> longlat and back

I'm a GIS novice, and am struggling to understand how to drive proj4 to accomplish what seems like a simple thing. I've googled for a couple of days, but everything I find seems to presuppose quite ...
4
votes
1answer
990 views

PROJ.4 custom projection that is Transverse Mercator with Affine post-process

I am trying to define some of my local TMAF grids in PROJ.4 so I can reproject points between the local grid and MGA95 and/or WGS84. For these grids I have the MapInfo CoordSys clause which defines a ...
4
votes
1answer
2k views

Can we use proj to offset coordinates?

Can someone please help me construct a proj string to offset coordinates? The detail we have - one coordinate point in a local construction grid, and the corresponding WGS1984 LatLong which we ...
3
votes
1answer
860 views

Creating a local projection system in QGIS using Proj.4

I'm trying to recreate a local mine grid for the QGIS users within our group. This grid is a modification of the Northing (-8000000). I imported the original UTM zone (35S). What is the function in ...
3
votes
2answers
744 views

How can I convert points in EPSG:2956 to EPSG:4326 with pyproj?

I have some points in EPSG:4326 -110.28879750 64.46310187 405.49938863 -110.28879359 64.46310086 405.45967110 -110.28880175 64.46310446 405.93823604 -110.28880245 64.46310258 405.48060850 -110....
3
votes
2answers
2k views

What is the SRS code for Lambert Conformal Conic projection and GCS Everest Definition 1962?

I have a shapefile with Geographic Coordinate System is "GCS_Everest Definition 1962" and the Projected Coordinate System is "Lambert_Conformal_Conic". I want to add this shapefile in Geoserver. But I ...