I have a two shapefiles. One of neighborhoods and one of points. The point file was originally a spreadsheet of X,Y coordinates. but have been exported to shapefile. They display but not over the neighborhood shapefile like they should. I've tried reprojecting,  project tool and define projection. It still will not overlap. Everything is displaying in US (Feet). What is wrong? Using arcmap 10.4   Here are the projection data in the Layer source:  Point file:  Projected Coordinate System:   NAD_1983_StatePlane_Pennsylvania_South_FIPS_3702_Feet Projection: Lambert_Conformal_Conic False_Easting:  1968500.00000000 False_Northing: 0.00000000 Central_Meridian:   -77.75000000 Standard_Parallel_1:    39.93333333 Standard_Parallel_2:    40.96666667 Latitude_Of_Origin: 39.33333333 Linear Unit:    Foot_US  Geographic Coordinate System:  GCS_North_American_1983 Datum:  D_North_American_1983 Prime Meridian:     Greenwich Angular Unit:   Degree  Neighborhood file:  Projected Coordinate System:   NAD_1983_StatePlane_Pennsylvania_South_FIPS_3702_Feet Projection: Lambert_Conformal_Conic False_Easting:  1968500.00000000 False_Northing: 0.00000000 Central_Meridian:   -77.75000000 Standard_Parallel_1:    39.93333333 Standard_Parallel_2:    40.96666667 Latitude_Of_Origin: 39.33333333 Linear Unit:    Foot_US  Geographic Coordinate System:  GCS_North_American_1983 Datum:  D_North_American_1983 Prime Meridian:     Greenwich Angular Unit:   Degree  Here's a sample of the Lat Long Coordinates:  Lat1        Long1 39.93128729   -75.14691545 39.93130836   -75.14711798 39.93143481   -75.14845446 39.93144498   -75.14853968 39.93117173   -75.14714819 39.93099668   -75.14705857 39.9310071    -75.14715449 39.93101477   -75.14725051 39.93119391   -75.14848865
3562
I put 4 different layers into my ArcMap document, but only one of them will show up on the map! I should be seeing the other 3 layers as well. When I do a Zoom to layer, I can see one layer, but I cannot figure out how to make the other 3 show up at once on top of each other.  Before I opened ArcMap, I used ArcCatalog where I right-clicked 3/4 of the shapefiles and  clicked Properties &rarr; Coordinate system &rarr; Import and then clicked on my Trees layer, which has the coordinate system I want the other 3 to have. But when I opened ArcMap the 4 layers did not overlap, nor could I move them above or below each other.
79921
I have several rectangular polygons that overlap each other and I'd like to label them with the ID's (for atlas page numbers).  The image below shows my initial setup.  I've offset each label from the centroid using the data defined override which I can easily understand how this one works.  But I would like to add more offset to each individual label using a CASE WHEN expression, so here's the issue.  At the bottom of the image are two inputs for "Offset X,Y".  The data defined override for this option says:  Data defined override Active: yes (ctrl/right-click toggles) Expected input: double coord [X,Y] as &lt;= 0.0 =&gt; Valid input types: string, int, double Current definition (expression):   I've used the overide for other individual coordinate values (just X or just Y) and they work fine because they expect just a single number.  But what is the format for the expression to use BOTH X and Y values?  I've tried two comma sperated numbers in square and curly brackets, as well as just numbers and also without the comma separation but I keep getting Expression is invalid from the output preview.  Any answers that I come across only really refer to python which I can't use.
27865
I'm using QGIS 2.4.0 and I have a shapefile of points with elevation data and X and Y coordinates. I'm using windows 7 OS.  The idea is to create a contour line SHP using my point data.  My question is: How do I start working on that issue? Is there a function? Is there a step by step process?
115738
I am trying to import an Excel table with around 10,000 XY points into ArcMap 10.2. I am using State Plane NAD 1983 Alabama West (US Feet) as the Coordinate System in my map data frame. I then add the excel table to my Table of Contents in the map then select Display XY Data. Next, I choose the same coordinate system as my data frame then select OK. Once it creates the Events layer for my points they are not in the correct location.   Should I switch everything in my map to WGS84 or is there a solution for my points to be added onto my current map with the state plane coordinate system?
102903
I am new to geotools and i find bit difficult to plot a point using the latitude and longitude values. can any one help me out to solve my problem. I tried the following code but not helpful.By quickstart i can able to pass a shapefile and getting the map but unable to plot the lat long values on it.     SimpleFeatureTypeBuilder b = new SimpleFeatureTypeBuilder();           b.setName( "MyFeatureType" );           b.setCRS( DefaultGeographicCRS.WGS84 );           b.add( "location", Point.class );            //building  the type           final SimpleFeatureType TYPE = b.buildFeatureType();          SimpleFeatureBuilder featureBuilder = new SimpleFeatureBuilder(TYPE);         GeometryFactory geometryFactory = JTSFactoryFinder.getGeometryFactory();          Point point = geometryFactory.createPoint(new Coordinate(longitude, latitude));           featureBuilder.add(point);           SimpleFeature feature = featureBuilder.buildFeature( "fid.1" );           DefaultFeatureCollection featureCollection = new DefaultFeatureCollection("internal",TYPE);           featureCollection.add(feature);   Adding to map:    List a=new ArrayList(); a.add(-155.2792); a.add(67.3623);  double longitude=(Double) a.get(0); double latitude=(Double)a.get(0);  File file = JFileDataStoreChooser.showOpenFile("shp", null);  if (file == null) { return; }  FileDataStore store = FileDataStoreFinder.getDataStore(file);  SimpleFeatureSource featureSource = store.getFeatureSource();  MapContext map = new DefaultMapContext();  map.setTitle("Show Map");      Layer layer = new FeatureLayer(featureSource, null);         map.addLayer(layer);     map.addLayer(featureSource, null);      JMapFrame.showMap(map);
21090
I am running QGIS 1.7.0 on windows, and I donwloaded a shapefile of .  This file appears to be in the WGS 84 coordinate system, and when I open it in QGIS the unit of measure is degrees.  I tried re-projecting the file into the NAD 83 coordinate system, and saving it as a new shapefile, but when I open the new file it still displays a ruler in degrees.  If I try to make a buffer, it interprets buffer distance as degrees.  Assuming I'm starting with a shapefile from TIGER, what's the best way to convert my unit of measure from degrees to meters?    For posterity, it seems like the following workflow make sense when dealing with US census TIGER shapefiles:           Open it in QGIS     Click the Settings menu and select project properties     Select a projected coordinate system, like NAD83 / UTM zone 19N     Click the box Enable 'on the fly' CRS transformation     Click ok     Click the 'zoom full' button (one of the magnifying glasses)
55838
I am trying to perform a spatial join between point data and polygon data.  I have data that indicate the spatial coordinates of an event in my csv file A and have another file, shapefile B, that contains the boundaries of an area as polygons.  head(A)   month   longitude latitude lsoa_code                   crime_type 1 2014-09 -1.550626 53.59740 E01007359        Anti-social behaviour 2 2014-09 -1.550626 53.59740 E01007359                 Public order 3 2014-09 -1.865236 53.93678 E01010646        Anti-social behaviour  head(B@data)   code      name                                  altname 0 E05004934 Longfield, New Barn and Southfleet    &lt;NA&gt; 1 E05000448                   Lewisham Central    &lt;NA&gt; 2 E05003149                            Hawcoat    &lt;NA&gt;   I want to join the crime data A to my shapefile B to map the crime events that happen in my area A. Unfortunately I cannot perform an attribute join based code as the code in A refers to different units than the code in B.  I've read a number of tutorials and posts but could not find an answer. I tried:  joined = over(A, B)   and overlay, but did not accomplish what I wanted.  Is there a way to do this join directly or would an intermediate transformation from A to another format be needed?  Conceptually I want to select those points of A that fall into the code areas of B (similar to "join based on spatial location in ArcGIS").  Did someone have this issue and solved it?
73587
I am trying to translate coordinates. I have a set of coordinates for city A which I need to map to city B, given a point P1 in city A maps to P2 in city B. How do I do this? This problem would be similar to a translation transformation but in this case, it's on a sphere and not on a plane.
57233
I like to think I'm relatively well-versed in datums, projections, and coordinate reference systems.  OpenStreetMap data is stored in WGS84 Lat/Lon (EPSG:4326). This CRS is geographic and therefore intended to store locations on a 3D globe. When I see  what am I looking at?  Looking at Reykjavik it appears that the longitude lines are tending towards the North Pole at the Greenwhich Meridian (or maybe they planned the city in perspective view??) so I assume 0,90 is top-centre of the map.  Does this class as a projection, and if so which one?   Is there one projection that is commonly used for lat/lon data?  I have data stored in Oracle using EPSG:4326 and when rendered by GeoServer without a specified projection it also displays these characteristics.
62321
I've got about a dozen polygons in a feature class loaded in ArcMap 10, all in geographic WGS 1984.  How do I easily obtain the coordinates associated with each vertex of each polygon in that feature class?   Ideally I would like them to be nicely tabulated in spreadsheet format.
72125
Excuse me if this is a stupid question, but I've honestly tried to find the answer in a few threads here on SE, and I cannot seem to follow the suggested steps correctly.    Project is in  (WGS 84 / Pseudo-Mercator) BaseMap (OpenStreetMap) layer is set to EPSG:3857 Vector data has been imported from CSV, and has Lat/Long coordinates, set to EPSG:3857.    The layer appears right in the middle of the map, around 0,0 obviously.   Why is this, what am I doing wrong?
10699
Using QGIS v. 1.9 (1.7 does not seem to be available for Ubuntu Linux 12.04), do the following:   Create new project, specify coordinate system by clicking on button in bottom right hand corner, enable on the fly CRS transformation, set projection to WGS 84 / UTM Zone 18S. Pass mouse over project window and UTM coordinates appear in coordinates window in status bar at the bottom. Add delimited text layer in decimal degrees. When asked for projection, specify same projection as in (1). Points in the layer are plotted but, mousing over the map, the coordinates in the coordinates window are now in decimal degrees (even though the projection is set to UTM).   Is there something I am missing here? I have tried doing this in so many ways, so many times, but the result is always the same. I seem to be setting the projection correctly but only decimal degrees show in the coordinates window.  With thanks, Martin
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
62433
I have a point layer and a polygon layer that are all in UTM51N projection but do not align. Just want to know when should I use Assign, Reproject, or Define.  Experimented with a polygon in EPSG:4326 - WGS 84 - Geographic, which is the project CRS since this polygon is the first layer I loaded. Target projection is UTM51N.  Reproject - aligned  Assign - not aligned  Define - not aligned  Save layer as - aligned  Using 3.4 so there is no "on the fly projection" feature.
67497
I have a vector floorplan (DXF) that I've added to a blank project in  using the Dxf2Shp plugin.   I've been reading tutorials and StackExchange threads, but I can't seem to find out how to "georeference" my layer?  I'm assuming I can somehow select the corners of the building and set the lat/long.  My goal is to export a shapefile that I can use with  and .  I realized that the project coordinates were actually lat/long coordinates:    and I need to use the "Affine (Rotation, Translation, Scale)" plugin to "place" my CAD drawing.
79048
I am plotting witness tree data to calculate pre-settlement tree density. I have a shapefile with points along section lines (square miles) at each half mile. I am making another table that contains tree data - distance, bearing, and diameter - related by survey point to the existing shapefile attribute table. Is there a way to create points on the map for the trees based on the distance and bearing from each survey point?   I'm not necessarily asking for a detailed set of instructions (although if anyone wishes to, I'd be grateful) but rather just if it's possible and a couple key words to continue searching how to do it.
69450
Copied from    The problem  Suppose we have two shapefiles that should border seamlessly. Only, they don't. Is there a way to force them to stick to one another without gaps?      The specific case  I have two shapefiles: one for European regions -- REG, the other for the neighbouring countries -- NEI. Both shapefiles are taken from  and should fit together nicely; but there are small gaps. Also, I need to simplify the polygons, and then the gaps become really notable.       The best I can think of  I've tried several approaches but with no success. The only way to achieve the desired result that I see requires following steps:   create a line sf with just the border between my shapefiles; from this line create a buffer polygon just big enough to cover all gaps; join and dissolve this buffer to the shapefile for neighbours -- NEI; clip off the expanded NEI with the REG shapefile.   Obviously, this is a rather clumsy workaround.    Is there a better way to go?    Reproducible example in     A minimal example  # install dev version of ggplot2 devtools::dev_mode() devtools::install_github("tidyverse/ggplot2")  library(tidyverse) library(sf) library(rmapshaper)  library(ggthemes)   # load data source(file = url("https://gist.githubusercontent.com/ikashnitsky/4b92f6b9f4bcbd8b2190fb0796fd1ec0/raw/1e281b7bb8ec74c9c9989fe50a87b6021ddbad03/minimal-data.R"))  # test how good they fit together ggplot() +          geom_sf(data = REG, color = "black", size = .2, fill = NA) +         geom_sf(data = NEI, color = "red", size = .2, fill = NA)+         coord_sf(datum = NA)+         theme_map()  ggsave("test-1.pdf", width = 12, height = 10)  # simplify REGs &lt;- REG %&gt;% ms_simplify(keep = .5, keep_shapes = TRUE) NEIs &lt;- NEI %&gt;% ms_simplify(keep = .5, keep_shapes = TRUE)   ggplot() +          geom_sf(data = REGs, color = "black", size = .2, fill = NA) +         geom_sf(data = NEIs, color = "red", size = .2, fill = NA)+         coord_sf(datum = NA)+         theme_map()  ggsave("test-2.pdf", width = 12, height = 10)