I have a shapefile with following projection information:  PROJCS["kl2",GEOGCS["kl3",DATUM["D_Everest_1830",SPHEROID["Everest_1830",6377299.36,300.8017]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Polyconic"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",76.5],PARAMETER["Latitude_Of_Origin",8.5],UNIT["Meter",1.0]]  How can I convert this datum to WGS 1984?
122735
I got sample shapefiles with the following *prj file content:     PROJCS["RGF93_Lambert_93",GEOGCS["GCS_RGF_1993",DATUM["D_RGF_1993",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",700000.0],PARAMETER["False_Northing",6600000.0],PARAMETER["Central_Meridian",3.0],PARAMETER["Standard_Parallel_1",44.0],PARAMETER["Standard_Parallel_2",49.0],PARAMETER["Latitude_Of_Origin",46.5],UNIT["Meter",1.0]]   Trying to convert the projection using the shp2pgsql command shp2pgsql -s 4326 kept the original projection and didn't convert it as requested.  I also check spatial_ref_sys table srtext field for similar known projection unsuccessfully.  How can I convert it to WGS 84 projection?
46179
I have a world map in correct Longitudes from -180° to 180°. Additionally, I have a shapefile of data from South America, which shows at 290° instead of -70°. How can I reproject it using QGIS?
56584
I've downloaded a map from  (limited to San Francisco census tracts), and I'd like to convert it to KML format. However, after unzipping the shapefile .zip file and running  $ ogr2ogr -f KML pop_by_census_tract.kml 140_00.shp   I get the following error message:  ERROR 6: No translation for Mercator_Auxiliary_Sphere to PROJ.4 format is known. Warning 1: Failed to create coordinate transformation between the input coordinate system and WGS84.  This may be because they are not transformable, or because projection services (PROJ.4 DLL/.so) could not be loaded. KML geometries may not render correctly. This message will not be issued any more.   Source: PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",     GEOGCS["GCS_WGS_1984",         DATUM["WGS_1984",             SPHEROID["WGS_1984",6378137.0,298.257223563]],         PRIMEM["Greenwich",0.0],         UNIT["Degree",0.0174532925199433]],     PROJECTION["Mercator_Auxiliary_Sphere"],     PARAMETER["False_Easting",0.0],     PARAMETER["False_Northing",0.0],     PARAMETER["Central_Meridian",0.0],     PARAMETER["Standard_Parallel_1",0.0],     PARAMETER["Auxiliary_Sphere_Type",0.0],     UNIT["Meter",1.0]]  ERROR 1: Latitude 4553026.920500 is invalid. Valid range is [-90,90]. This warning will not be issued any more Warning 1: Longitude -13627851.842700 has been modified to fit into range [-180,180]. This warning will not be issued any more   When I open the resulting KML file in Google Earth, I get randomly scattered lines across the entire planet. The example at code.google.com/apis/kml/articles/vector.html works fine, so the installation of GDAL (1.8.0 from macports) isn't totally broken.   seems related, but I'm not sure if it's saying that Mercator_Auxiliary_Sphere is just not supported by GDAL, or if it's implying I should pass another command line option to make it work.    I wound up accomplishing my goal by downloading , which contains shapefiles that ogr2ogr and shp2pgsql can understand, and joining it with the pure table data from factfinder2. I'm still curious about the format factfinder2 is using, but it's no longer blocking progress for me.
33403
I want to do something relatively simply but it seems to be very difficult to do, though I am very new to using R for spatial data so apologies if I've missed something obvious.  I have a list containing 161 shapefiles of bird distributions in the world, they are each in the format of a SpatialPolygonsDataFrame. I also have world climate data in a raster stack format for 8 variables. Both sets of files are in the same projection (lat and long WGS84). I want to create new columns for the climate data with 1's where the species polygon overlaps with the climate data and 0 where it doesn't. I tried various functions to do this, quite a few initial options didn't work, the best option appeared to be:  pnt.in.poly(XY coordinates of climate data points, XY coordinates of species polygon)   I changed my raster of climate data into points to get the XY coordinates. I unified the polygons for each species shapefile and then used:  species@polygons[[1]]@Polygons[[1]]@cords   to get the coordinates for my species polygons as a matrix. However, I discovered that:  species@polygons[[1]]@Polygons[[1]]@coords    wasn't actually getting all the coordinates of the polygons only some of them.  I tried using fortify from the ggplot package, and it gives me a dataframe with coordinates but some of them are for holes in the polygons. I ran that with the pnt.in.poly function excluding the holes but it took an hour to run and when I looked at the output maps after it matched for some of the range but not all of it- there were holes in the distribution where there shouldn't be holes.   I also tried this method that I found somewhere else:  lonlat &lt;- lapply(slot(species, "polygons"), function(x) lapply(slot(x,"Polygons"), function(y) slot(y, "coords")))    which gives me a list with 436 items- but I'm still not sure how to get the coordinates into a matrix form that will work in pnt.in.poly.  I'm sure this should be easier to do! Does anyone know a better way?
115605
I am attempting to create a 24x36 map of New York City, and overlay it with a shapefile from City Planning. I'm using the Openlayers plugin, and I've tried all of the different types of base street map layers. I have been able to successfully convert the shapefile from a NAD 83 projection (EPSG 2263) to a WGS 84 projection (EPSG 3857), but I'm still suffering from some type of shift whenever I draw the map out on the print composer on a large scale piece of paper. The larger I try to print, the more dramatic the shift the misalignment.
4203
While importing my shapefile data to PostGIS, I did not select the proper Projection.  How do I now change the SRID of the data, without transforming the Coordinates?
12074
I have a shapefile of size 224 MB. I am using ArcGIS to convert it to a KML file. The process involves converting the shapefile to a feature layer and then converting the layer to KML.   An alternative method that I have tried is using ogr2ogr.  In both the cases, the KML file produced is around 500 MBs. I need it to be less than 250 MBs, so I can import it as a Google Fusion table.  I was informed that one of the best ways to cut size is to reduce precision of the coordinates to 5 decimal places (from its current 10 or so).  How do I go about changing the decimal places using either ogr2ogr or ArcGIS? Is there some other software that may help in this and also, are there any other strategies I may use to reduce the size of the KML file?  I tried using Notepad++ to open the KML file and then use regular expressions to change coordinate precision but for some reason this method is not working. If I make the change and try and save the file, all the data gets deleted and the KML file size changes to 0 bytes.  An alternative strategy that I was thinking of was to split my shapefiles into two or three files and convert them to KML separately, upload them to Google fusion tables and then append those. However, I am not sure of two things.   How to split shapefiles into two or three smaller shapefiles? Whether I can easily append fusion tables together to create a larger fusion table?
83757
I have a feature class with x/y values populated on the native projection of the feature class (state plane).  In the same feature class I also have Lat/Long fields that I need to populate (WGS 84 in decimal degree).  The problem I have is the calculate geometry function doesn't allow me to choose a different projection other than State Plane when I try to auto populate the Lat/Long values.  Does anyone have any python code or ideas that might help with this request?
83537
I need to print a map with coordinates in the format of Degree, Minute and second. CRS is set to WGS 84 (EPSG 3857) for both layers and on the project properties it is set to enable transformation on the fly. But then, when in map composer I get the following coordinates, which make no sense to me.
58124
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 having difficulty converting the x,y coordinates in Ohio North State Plane usFt to WGS84.  &gt;&gt;&gt; p = Proj(r'+proj=lcc +lat_1=41.7 +lat_2=40.43333333333333 +lat_0=39.66666666666666 +lon_0=-82.5 +x_0=600000 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=us-ft +no_defs') #Nad83 State Plane Ohio North US Feet Proj object using parameters &gt;&gt;&gt; p(739400.91,2339327.3,inverse=True) (-80.138057868777224, 60.278230707978487) &gt;&gt;&gt; p1 = Proj(init="epsg:3734") #Nad83 State Plane Ohio North US Feet Proj object using EPSG code &gt;&gt;&gt; p1(739400.91,2339327.3,inverse=True) (-80.138057868777224, 60.278230707978487)   Both methods above return the same result, however this lat long is somewhere in Hudson Bay.  When I plot the coordinates in ArcMap, the correct lat long is: -81.142311,41.688205.  *Notice all lat longs are provided long,lat as this is the order Proj uses  Does anyone know why I would be getting the wrong coordinates from Proj.4 and pyproj?
39912
I am trying to figure out how I can transform my unit of length from decimal degree to km when I import the following simple ArcGIS polygon shape.  Ultimately, I would like to have this information available when I assign this spatial polygon to my window in spatstat [i.e., Window area = xx square unit &amp; Unit of length: xx km].   I did came across this post:   This is what I've done so far, but I don't understand why/if I need to do so. Doesn't my div0A SpatialPolygonsDataFrame already has my CRS argument defined? Why do I need to spTransform then?  library(spatstat) library(maptools) library(lattice) library(sp) library(foreign) require("rgdal") div0A &lt;- readOGR(dsn="Archive", layer="Projections") &gt; summary(div0A) Object of class SpatialPolygonsDataFrame Coordinates:      min       max x -80.00 -57.63430 y  66.25  78.16666 Is projected: FALSE proj4string : [+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0] Data attributes:        Id      NPAzimutha         UTM20            UTM19           AlberEA        Min.   :0   Min.   :348877   Min.   :349232   Min.   :349162   Min.   :348656    1st Qu.:0   1st Qu.:348877   1st Qu.:349232   1st Qu.:349162   1st Qu.:348656    Median :0   Median :348877   Median :349232   Median :349162   Median :348656    Mean   :0   Mean   :348877   Mean   :349232   Mean   :349162   Mean   :348656    3rd Qu.:0   3rd Qu.:348877   3rd Qu.:349232   3rd Qu.:349162   3rd Qu.:348656    Max.   :0   Max.   :348877   Max.   :349232   Max.   :349162   Max.   :348656   proj4string(div0A) &lt;- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs") xxlcc &lt;- spTransform(div0A, CRS("+init=epsg:4326")) res1 &lt;- lapply(slot(xxlcc, "polygons"), function(x)   sapply(slot(x, "Polygons"), slot, "area"))  &gt; summary(unlist(res1))    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.   89.72   89.72   89.72   89.72   89.72   89.72  newW&lt;-as(xxlcc, "owin") &gt; summary(newW) Window: polygonal boundary single connected closed polygon with 29 vertices enclosing rectangle: [-80, -57.6343]x[66.25, 78.16666]units Window area =  89.7204 square units   Does this look correct? How can I specify/know my unit of length in km? Moreover, this window is ~Baffin Bay, west of Greenland and I am wondering if I should be using a different datum? E.g. EPSG Projection 4747 - GR96? Although, I will have to create ppp object, for which the coordinates have been taken with GPS satellite navigation system.  I would appreciate if someone can walk me through this or suggest an online example or further readings. I read chapter 4: Spatial Data Import and Export in Applied Spatial Data Analysis with R (2008), but I'm obviously very new to all of this.
30642
How can I convert an excel file with x, y columns to a point shapefile?  There are some somewhat optional requirements in addition to the correct creation of a shapefile:   Column types (as per Excel's format specifiers) should be retained (especially date types) Column names should be taken from the header I would like to do this from the command-line If I can include heterogeneous spatial references for the points in a third column I would be very happy :)
71486
I'm not sure how to remove the line that represents the 180th meridian (international dateline) on my polygon. I have the "Russia Country Boundary 2017" feature layer added from the Living Atlas. My Projected Coordinate System is set to "Asia North Albers Equal Area Conic". I've tried to google this question, but I haven't come across an answer I can fully understand. I'm a student and very new to all this!  The line is actually a gap between the polygons.   Using ArcGIS Pro 2.3.1  ANSWER:  It turns out I was not projecting correctly, which was creating issues down the line when editing polygons. Here's what I did...  1. Start a blank project in ArcGIS. Insert new map, go into its     properties, change the Coordinate System to "Asia North Albers Equal     Area Conic".  2. Add "World Countries (Generalized)" from the Living     Atlas. Select Russia and Export Features to new feature class layer.  3. Search and select the Project tool from Analysis>Geoprocessing>Tools.   4. Select the Russia layer for input. Select "Current Map" for Output Coordinate System. Run.  5.Now edit the polygons to close the gap and merge to remove any lines. All done!  The only thing that bothers me with this fix is the addition of shape area (484,914.49674 m²) and shape length (4.59 m) to the original Russia polygon. These numbers could be different depending on how you close the gap, and possibly affect data projections depending on your use of this map?
96543
I have a file with these information:    lat_min -60.0 lat_max 80.0 lon_min -180.0 lon_max 179.9 lon_sampling 0.1 degree lat_sampling 0.1 degree samples = 3600 lines   = 1401   I want to re-project the file to   I tried    gdalwarp -of "ENVI" -t_srs `EPSG:3410`  -ts 1383 586 but give error  gdalwarp -of "ENVI" -t_srs "+proj=cea +lon_0=0 +lat_ts=30 +x_0=0 +y_0=0 +a=6371228 +b=6371228 +units=m +no_defs" -ts 1383 586 -ot Float32  file.envi file_1383586.envi  ERROR 1: No PROJ.4 translation for source SRS, coordinate transformation initialization has failed
17985
I'm trying to use those parameters and I receive a difference of ~0.5m    ALON and BSHM APN Points:  coords = [   (dms2rad([31,42,28.5694]),  dms2rad([34,36,23.92196]),  55.671),   (dms2rad([32,46,44.34472]), dms2rad([35,1,22.74061]), 225.046), ]   Transform:  def lla_to_ig05(lon, lat, alt):     import pyproj     lla = pyproj.Proj(proj='latlong', ellps='WGS84', datum='WGS84')     ig05 = pyproj.Proj("+proj=tmerc +lat_0=31.7343936111111 +lon_0=35.2045169444445 +k=1.0000067 +x_0=219529.584 +y_0=626907.39 +ellps=GRS80 +towgs84=-24.002400,-17.103200,-17.844400,-0.33077,-1.852690,1.669690,5.424800 +units=m +no_defs")     x, y, h = pyproj.transform(lla, ig05, lon, lat, alt, radians=True)     return x, y, h for lat, lon, alt in coords:   print(pd.DataFrame(lla_to_ig05(lon, lat, alt))) # pd for pandas   print("")   Output:                 0 0  162790.443392 1  624087.953905 2      55.671000                 0 0  202456.903152 1  742717.095893 2     225.046000   Real numbers:                 0 0  162789.564 1  624088.355 2      55.671000                 0 0  202456.453 1  742718.779 2     225.046000   Any hint will be very helpful!
297
How do I change the projection of my shapefile using QGIS? The default projection is set to EPSG:3003, I want to change it to WGS84 EPSG:32632. When I do a reprojection or I change projection in the properties mask I don't have the correct result. On the left, I have the plan coordinates but they are not correct for the EPSG shown on the right:
31731
I have a project and a topomap in WGS84 (EPSG4326). On the fly CRS transformation is activated. I open Google satellite layer from the OpenLayers plugin. The Google layer is WGS84 / Pseudo Mercator (EPSG3857). But then the whole project has also become WGS84 / Pseudo Mercator (EPSG3857).  I thought on the fly transformation would reproject to EPSG4326. Why is that? Is it better to work directly in WGS84 / Pseudo Mercator (EPSG3857) systematically when we use Google Earth?