I downloaded a mile marker geojson from the Idaho DOT and the CRS is listed as 102605, which doesn't seem to exist outside of the DOT office.   This is the extent of the layer metadata:  XMin: 2254295.3529000003 YMin: 1199867.0130000003 XMax: 2743212.4619000005 YMax: 1980420.2782000005 Spatial Reference: 102605  (102605)   link to mapserver info:
123095
Using Geoserver, I'm trying to get the correct projection for a map layer that is provided by the State of Idaho. I've tried looking up the EPSG code from . However, the search results throws an error.  Here is the location of the data -   Here are the results of the .prj file  PROJCS["NAD_1983_Transverse_Mercator",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",2500000.0],PARAMETER["False_Northing",1200000.0],PARAMETER["Central_Meridian",-114.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",42.0],UNIT["Meter",1.0]]   I found this page that should help, but I don't know what it all means.
60528
I have made a raster elevation map of Greenland, and have information of the region, cities and more as vectors. The information is not in the same place as my elevation map.  I have selected EPSG:5938 on the bottom right in QGIS and on my informations and elevation map I have set the CRS to the same like the picture    How do I get it all on the same position?
9430
I am having an issue with buffer. My data is in Geographic coordinate system. I am trying to do a buffer of 600m and i am finding different results for each case, have tried thse so far:    Map in web maercator map in UTM zone 12 map in geogaphic coordinate   I am missing some data points due to this and wanted to check how I can resolve this issue. I know that UTM is the most accurate projection for distance but I wanted to know if by using WebMercator map if there is a way to get correct buffer results , any ideas? Thanks jay
113120
I have a road network and want to test connectivity between the features.  I am having trouble with two specific features, the rest of the dataset is functioning just fine.  When I get those 2 links' start &amp; end locations using ST_StartPoint and ST_EndPoint, the result is that they are slightly different.  However, when I convert the result to text by wrapping it like so:  ST_AsText(ST_StartPoint(geom)) and ST_AsText(ST_EndPoint(geom))  then the results are exactly the same.  Is it possible they are actually off by a tiny amount like .000000001 and the ST_AsText loses a decimal point?
82784
I'm looking for a web service providing a full address based on coordinates (and reverse).  I'm ready to pay for that service.  I don't want to use Google's API.
33644
[EDIT: In my original post, I had received an error that was based on a typo, so I scrapped that part of the question, but left the questions that still apply.]  In the  it says:  -nlt type: Define the geometry type for the created layer. One of NONE, GEOMETRY, POINT, LINESTRING, POLYGON, GEOMETRYCOLLECTION, MULTIPOINT, MULTIPOLYGON or MULTILINESTRING. Add "25D" to the name to get 2.5D versions.   Does this part: Add "25D" to the name to get 2.5D versions apply to both loading things into postgres as well as exporting data from postgres to .shp files?  Additionally, I am assuming 25D to meant that one can have a z value that corresponds to each pair of xy coordinates (as is the case with PolygonZ shapefile types), but that these coordinates may not overlap. Is this correct? What is the intended distinction in this case between 2.5D and 3D?  Thanks
32326
How can one change a coordinate system for a .lasd (Lidar) dataset?  The file is currently in UTM_Zone_15N and I am trying to convert it to IL_East.
88671
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 proj4string&lt;-(*tmp*, value = ) :     Geographical CRS given to non-conformant data: 153.53332931   Here are relevant lines of the code,  # Loading suburb list  qs &lt;- read.csv(file = '../data/Australian_Post_Codes_Lat_Lon.csv',                 header = TRUE,                stringsAsFactors = FALSE)  suburb &lt;- qs %&gt;%    filter(state == 'QLD') %&gt;%    select(postcode,          suburb,          lat,          lon)   coordinates(suburb) &lt;- c('lat', 'lon') # Assigning proper Coordinate Reference System (CRS) to Lat/Long. proj4string(suburb) &lt;- CRS("+init=epsg:4326")   I am using the csv file provided
37637
I see that the Point type is made to store geometric coordinates, however I don't see what problem it is trying to solve, or what advantages it offers over simply using two numerical columns. It's 16 bytes, so is it just an abstraction for two 8 byte floats? Or is it a string containing some kind of geohash that optimizes the performance of certain queries?  Let's say I want to get all the points for a specific viewport on a Leaflet map... with numerical columns representing latitude and longitude, I'd just filter for all the rows that have a latitude number between the N/S bounds and a longitude number between the E/W bounds, easy. I looked into using Point here, but the syntax is cumbersome and it seems to me like unnecessary overhead for development and for database performance... is it?
60006
I have looked at solutions to what seem to be somewhat similar problems but they don't seem to solve my problem.  A georeferenced raster layer loads into the correct location in relation to other layers, so CRS is obviously OK, but has this black background which I don't know how to remove.
37219
I have an ASCII file with latitude, longitude, and data_val in the following format.  35-13.643782N, 080-57.190157W, 118.6 ...   I have a GeoTiff image file, and I can easily view it.  I want to place a "pin" (can be a dot/flag/star or whatever is easiest) on the image at the specific latitude/longitude position found in the ASCII file.  Here is what I've managed to do so far:  My source image looks like this:  Driver: GTiff/GeoTIFF Files: /tmp/Charlotte SEC 100.tif Size is 16867, 12358 Coordinate System is: PROJCS["Lambert Conformal Conic",     GEOGCS["NAD83",         DATUM["North_American_Datum_1983",             SPHEROID["GRS 1980",6378137,298.2572221010042,                 AUTHORITY["EPSG","7019"]],             AUTHORITY["EPSG","6269"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4269"]],     PROJECTION["Lambert_Conformal_Conic_2SP"],     PARAMETER["standard_parallel_1",38.66666666666666],     PARAMETER["standard_parallel_2",33.33333333333334],     PARAMETER["latitude_of_origin",34.11666666666667],     PARAMETER["central_meridian",-78.75],     PARAMETER["false_easting",0],     PARAMETER["false_northing",0],     UNIT["metre",1,         AUTHORITY["EPSG","9001"]]] Origin = (-365041.822331817995291,240536.419747152860509) Pixel Size = (42.334586069440391,-42.334898968590878) Metadata:   AREA_OR_POINT=Area   TIFFTAG_DATETIME=2016:06:24 12:46:45   TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)   TIFFTAG_SOFTWARE=Adobe Photoshop CS5 Windows   TIFFTAG_XRESOLUTION=300   TIFFTAG_YRESOLUTION=300 Image Structure Metadata:   COMPRESSION=LZW   INTERLEAVE=BAND Corner Coordinates: Upper Left  ( -365041.822,  240536.420) ( 82d48'55.43"W, 36d13' 4.92"N) Lower Left  ( -365041.822, -282638.262) ( 82d35'10.11"W, 31d30'17.00"N) Upper Right (  349015.641,  240536.420) ( 74d51'46.40"W, 36d13'26.16"N) Lower Right (  349015.641, -282638.262) ( 75d 4'55.60"W, 31d30'36.99"N) Center      (   -8013.091,  -21050.921) ( 78d50'12.11"W, 33d55'36.35"N) Band 1 Block=16867x1 Type=Byte, ColorInterp=Palette   Color Table (RGB with 256 entries)     0: 255,255,255,255 ...   Here is what I've managed to cobble together in Python:  from osgeo import gdal, osr  src_filename = '/tmp/Charlotte SEC 100.tif' dst_filename = '/tmp/foo.tiff'  # Opens source dataset src_ds = gdal.Open(src_filename) format = "GTiff" driver = gdal.GetDriverByName(format)  # Open destination dataset dst_ds = driver.CreateCopy(dst_filename, src_ds, 0)  # Specify raster location through geotransform array # (upperleftx, scalex, skewx, upperlefty, skewy, scaley) # Scale = size of one pixel in units of raster projection # this example below assumes 100x100 gt = [-365041.822, 100, 0, 240536.420, 0, -100]  # Set location dst_ds.SetGeoTransform(gt)  # Get raster projection epsg = 4269            # http://spatialreference.org/ref/sr-org/lambert_conformal_conic_2sp/ srs = osr.SpatialReference() srs.ImportFromEPSG(epsg) dest_wkt = srs.ExportToWkt()  # Set projection dst_ds.SetProjection(dest_wkt)  # Close files dst_ds = None src_ds = None   But, I can't quite figure out how to place a "red dot" at 35-13.643782N, 080-57.190157W  I'm having to learn some new details here (nomenclature about GIS).
27274
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 with -180 180 longtitude range.  And I want the center of the image to be at Greenwich meridian (0). I guess this srs is very widely used. I expect the result look like this:    So I use a gdalwarp command to reproject:  gdalwarp -s_srs '+proj=latlong +datum=WGS84 +pm=180dW' -t_srs EPSG:4326 test_col.tif test_4326.tif   But I only get a tiff with bigger dimensions (more pixels) and EPSG:4326 metadata. The image itself looks the same, as the initial one. But I expect it to swap the hemispheres.  How do I gdalwarp an image to be strictly -180 180 EPSG:4326 with the center in 0 longitude?  This is gdalinfo of my initial file:  Origin = (-0.102272598067084,89.946211604095552) Pixel Size = (0.204545196134167,-0.204423208191126) Metadata:   AREA_OR_POINT=Area Image Structure Metadata:   INTERLEAVE=BAND Corner Coordinates: Upper Left  (  -0.1022726,  89.9462116) (  0d 6' 8.18"W, 89d56'46.36"N) Lower Left  (  -0.1022726, -89.9462116) (  0d 6' 8.18"W, 89d56'46.36"S) Upper Right (     359.897,      89.946) (359d53'50.18"E, 89d56'46.36"N) Lower Right (     359.897,     -89.946) (359d53'50.18"E, 89d56'46.36"S) Center      ( 179.8975000,  -0.0000000) (179d53'51.00"E,  0d 0' 0.00"S)   This is gdalinfo after gdalwarp     Origin = (-180.102727401932952,89.946211604095552) Pixel Size = (0.091397622896436,-0.091420837939082) Metadata:   AREA_OR_POINT=Area Image Structure Metadata:   INTERLEAVE=PIXEL Corner Coordinates: Upper Left  (-180.1027274,  89.9462116) (180d 6' 9.82"W, 89d56'46.36"N) Lower Left  (-180.1027274, -89.9699975) (180d 6' 9.82"W, 89d58'11.99"S) Upper Right ( 179.8211116,  89.9462116) (179d49'16.00"E, 89d56'46.36"N) Lower Right ( 179.8211116, -89.9699975) (179d49'16.00"E, 89d58'11.99"S) Center      (  -0.1408079,  -0.0118929) (  0d 8'26.91"W,  0d 0'42.81"S)
61738
I just started to play with leflet/geojson a little. But my coordinates are not rendered properly and I have no clue what is going on.  My coordinates are: 52.23943, 4.97599. They work correct with the setView function.  var map = L.map('leaflet_map').setView([52.23943, 4.97599], 15);   But using a geojasonFeature they are, hmmm, 'projected', somewhere east of Somalia.  var geojsonFeature = {             "type": "Feature",             "properties": {             "name": "Coors Field",             "amenity": "Baseball Stadium",             "popupContent": "This is where the Rockies play!"             },             "geometry": {                 "type": "Point",                 "coordinates": [52.23943, 4.97599]             }         };         var myLayer = L.geoJson().addTo(map);         myLayer.addData(geojsonFeature).bindPopup("I am a gjson point.");   Anyone who knows what is happening here?  EDIT  Out of pure curiosity I changed the coordinates around [4.976143930893815,52.23925499011473] and the point appears at its correct location. A known bug!?
118235
I copied the example from geotools, which shows how to transform a shapefile into a different projection. () But when transforming from Gauss Krueger 4 into Pseudo Mercator (3857) the shapes are not at the correct position, they are displaced several meters.  What can I do to fix this?  The sample code contains :   CoordinateReferenceSystem dataCRS = schema.getCoordinateReferenceSystem();             CoordinateReferenceSystem worldCRS = map.getCoordinateReferenceSystem();             boolean lenient = true; // allow for some error due to different datums             MathTransform transform = CRS.findMathTransform(dataCRS, worldCRS, lenient);   ... and then ...  Geometry geometry2 = JTS.transform(geometry, transform);   additional information:  The sample data of the source as the following CRS  PROJCS["PD_83_3_degree_Gauss_Kruger_zone_4",    GEOGCS["GCS_PD/83",      DATUM["D_Potsdam_Datum_83",        SPHEROID["Bessel_1841", 6377397.155, 299.1528128]],      PRIMEM["Greenwich", 0.0],      UNIT["degree", 0.017453292519943295],      AXIS["Longitude", EAST],      AXIS["Latitude", NORTH]],    PROJECTION["Transverse_Mercator"],    PARAMETER["central_meridian", 12.0],    PARAMETER["latitude_of_origin", 0.0],    PARAMETER["scale_factor", 1.0],    PARAMETER["false_easting", 4500000.0],    PARAMETER["false_northing", 0.0],    UNIT["m", 1.0],    AXIS["x", EAST],    AXIS["y", NORTH]]   and the target has the following CRS:  PROJCS["WGS 84 / Pseudo-Mercator",    GEOGCS["WGS 84",      DATUM["World Geodetic System 1984",        SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],        AUTHORITY["EPSG","6326"]],      PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],      UNIT["degree", 0.017453292519943295],      AXIS["Geodetic latitude", NORTH],      AXIS["Geodetic longitude", EAST],      AUTHORITY["EPSG","4326"]],    PROJECTION["Popular Visualisation Pseudo Mercator"],    PARAMETER["semi_minor", 6378137.0],    PARAMETER["latitude_of_origin", 0.0],    PARAMETER["central_meridian", 0.0],    PARAMETER["scale_factor", 1.0],    PARAMETER["false_easting", 0.0],    PARAMETER["false_northing", 0.0],    UNIT["m", 1.0],    AXIS["Easting", EAST],    AXIS["Northing", NORTH],    AUTHORITY["EPSG","3857"]]   Any idea? D3
110668
First, I'll summarize what I understand so far so please confirm or correct me.   When On-the-Fly is selected a layer is reprojected from that of the layer's CRS to that of the project CRS, if the project CRS has been defined    and the layer CRS has been recognized   (which apparently is not always the case); The project CRS can be defined three ways: choose one from a menu, set it to a particular layer's CRS, let it default to a past CRS or a layer CRS. The layer CRS may be recognized on import (or not).  If it's not (or is) the way to configure the the CRS is to save the layer as a shape file and select a CRS when prompted.  However, it sounds like if a .qrj file has already been generated this will fail somehow, so we must assure ourselves a .qrj file is not present.   How am I doing so far?  Second:  I'm using QGIS 1.8 on a Windows 7 OS.  Let's say two different layers have been saved each with CRS designated as EPSG:2272 (NAD83 / Pennsylvania South (ftUS)) and we start a new project.  We select On-the-Fly reprojection and set the project CRS to EPSG:2272.  We open both .shp files as layers.  They are both files of Philadelphia, (PhilaCopIncsJune302012, PhilaCensusBlocks2012) and the the .prj files respectively are:     PROJCS["NAD83 / Pennsylvania South  (ftUS)",  GEOGCS["NAD83",DATUM["North_American_Datum_1983",  SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],  TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6269"]],  PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],  UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],  AUTHORITY["EPSG","4269"]],  PROJECTION["Lambert_Conformal_Conic_2SP"],  PARAMETER["standard_parallel_1",40.96666666666667],  PARAMETER["standard_parallel_2",39.93333333333333],  PARAMETER["latitude_of_origin",39.33333333333334],  PARAMETER["central_meridian",-77.75],  PARAMETER["false_easting",1968500],  PARAMETER["false_northing",0],UNIT["US survey  foot",0.3048006096012192,AUTHORITY["EPSG","9003"]],  AXIS["X",EAST],AXIS["Y",NORTH],AUTHORITY["EPSG","2272"]]   and    PROJCS["NAD83 / Pennsylvania South  (ftUS)",  GEOGCS["NAD83",  DATUM["North_American_Datum_1983",  SPHEROID["GRS  1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],  TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6269"]],  PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],  UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],  AUTHORITY["EPSG","4269"]],  PROJECTION["Lambert_Conformal_Conic_2SP"],  PARAMETER["standard_parallel_1",40.96666666666667],  PARAMETER["standard_parallel_2",39.93333333333333],  PARAMETER["latitude_of_origin",39.33333333333334],  PARAMETER["central_meridian",-77.75],  PARAMETER["false_easting",1968500],  PARAMETER["false_northing",0],  UNIT["US survey foot",0.3048006096012192,AUTHORITY["EPSG","9003"]],  AXIS["X",EAST],AXIS["Y",NORTH],AUTHORITY["EPSG","2272"]]   The CRS found in the PhilaCopIncsJune302012 info document is WGS84 decimal degrees.  The CRS found in the PhilaCensusBlocks2012 XML document is      NAD_1983_StatePlane_Pennsylvania_South_FIPS_3702_Feet   The two layers are far far away from each other.  Any ideas why?
120503
To start off with, I'm very new to QGIS. I'm attempting to place a shapefile with locations for 3 substations on a map of Peniche, Portugal. The locations for the substations are given in latitude and longitude coordinates.   When I start off by importing the substation shapefile, I make sure that both the layer CRS and project CRS are set to EPSG:32629, which gives each substation point the correct coordinates, however the scale is at 1:1, which is not correct. As soon as I then add a raster map of the region or a shapefile of Portugal using the same CRS the position of the substations also becomes incorrect.   I'm sure there's a simple solution I'm not getting, but I've been trying for a while and have had no success. Here's a dropbox link to the files
111538
I have slight experience with ArcGIS, but am new to QGIS.  I'm trying to set the projection of my map, and I understand that this is done by altering the CRS (coordinate reference system) of each layer and/or the entire project.  When I right click a layer and select "Change Layer CRS", however, the visual representation of the map (the distortion) does not reflect the selected projection. My world map looks like it has a Robinson (or similar) projection, even if I select a Mercator CRS. The EPSG code in the lower right corner changes, as do the coordinates.    Setting the project CRS to the same projection does not help, nor does enabling on-the-fly projection.  My goal is to export the map to SVG, so the visual distortions are important.  Am I doing something wrong?