Tag Info

Hot answers tagged

5

Assuming you have a reference to the feature class pFC: ERSI.ArcGIS.Geometry.ISpatialReference pSP = ((ESRI.ArcGIS.GeoDatabase.IGeoDataset)pFC).SpatialReference; coord=pSP.FactoryCode coord now contains the code for the coordinate system. Note that not all projected coordinate systems are defined in esriSRProjCSType, there are more constants ...


4

The best way I can think is to get two UTM points, convert them to Lat/Long, and compare their geodesic distances to their UTM pythagorean distance. E.g. Take a point from this example: The CN Tower is ... in UTM zone 17, and the grid position is 630084m east, 4833438m north. So if we take A (17n 630084 4833438) and move it 30 km east, we get B (17n ...


4

If you divide those "strange" latitudes and longitudes by 11930465, you get the North and East degree values you expect: But don't ask me why ;-)


4

Based on the conversation on the comments, ArcGIS online is making a guess at the coordinate system, while ArcMap is not. When it encounters an unknown coordinate system, ArcGIS Online appears to automatically assume WGS 1984. ArcGIS for Desktop treats undefined coordinate systems differently and does not make assumptions about the data's actual coordinate ...


4

The geographic coordinates systems (also called lat/long) are defined by 3 features : a datum, a prime meridian (most of the time, it is Greenwich) and a unit (most of the time, it is degree). The EPSG code help to unambiguously identify a geographic coordinate sytem (you can find the description on spatial reference.org) GEOGCS["WGS 84", ...


4

Trying to find a solution, I've changed OpenLayers.js file, provided by plugin authors (which had 'Release 2.11' version in his body and was placed at ~/.qgis2/python/plugins/openlayers_plugin/weblayers/html folder in my system) to last 2.13.1 version from official OpenLayers site (pay attention, 2.+ version!). Just download the archive, extract to ...


4

Georeferencing is a process where you define which coordinates your image covers. It'll write in the coordinate system of the data frame. Therefore, it does not specify which coordinate system you use, only which coordinates it covers. The tool Define projection (under Data Management --> Projections and transformations) can be used on both vector and ...


4

I'll use spherical coordinates as defined here on Wikipedia which uses phi and theta (which is probably your lambda). Phi is the angle from the north pole. Hence if the WGS84 point is 10.0.0N, phi will be 80 degrees. For a point in the southern hemisphere, say 12.30.00S, phi will be 90 + 12.5 = 102.5 degrees. Theta is just the longitude in degrees, if the ...


3

You can have your WFS server reproject the data by adding the srsName parameter to your URL (i.e. http://ec2-54-69-8-151.us-west-2.compute.amazonaws.com:8080/geoserver/WRIA9/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=WRIA9:2009BuildingsCOS&maxFeatures=50&outputFormat=application/json&srsName=epsg:4326). Then your ...


3

I'd suggest not writing this yourself, but instead using one of the existing coordinate transform libraries, such as proj4j or the CRS part of GeoTools. From GeoTools (with some JTS help): import com.vividsolutions.jts.geom.Geometry; import org.geotools.geometry.jts.JTS; import org.geotools.referencing.CRS; import org.geotools.geojson.geom.GeometryJSON; ...


3

For the first part of your question: GDAL can guess the format of the input file from the file extension. The output format is defined by the -f option. If it is missing, Geotiff is assumed, but you get that warning if the file extension is not .tif. For a .grd output, you can select between GS7BG (rw+v): Golden Software 7 Binary Grid (.grd) GSAG (rwv): ...


3

The spatial reference for a map is defined either by the extent passed to the map constructor or by the first layer added to the map. Checkout the any projection sample to see how to specify different spatial references for a map. If you're using the basemap or center and zoom options with the map constructor, the map's spatial reference will be web ...


3

You are probably looking for the Define Projection tool (under Data Management, Projections and Transformations). Unlike Project, which will transform from one coordinate system to another (e.g. -119.789° becomes some large number of meters), Define Projection changes the coordinate system for the data itself and is appropriate to use when ArcMap is ...


3

Your Gauss-Krueger projection uses +datum=potsdam. Up to 2012, this was hard coded in proj4 to a very unprecise value using a 3-parameter-transformation. You find more exact values for 7-parameter transformations in this topic: http://forum.openstreetmap.org/viewtopic.php?id=12723 There is an even better ntv2-grid transformation available here (take the ...


3

To determine whether an arcpy SpatialReference object is projected or geographic use the property type: geoSR = arcpy.SpatialReference(4326) print geoSR.type Geographic projSR = arcpy.SpatialReference(28356) print projSR.type Projected


3

I am unclear why you prefer the name to the factory code. You can use the srid and there is no fuss with spaces/underscores. For projected systems: http://resources.arcgis.com/en/help/main/10.2/018z/pdf/projected_coordinate_systems.pdf For geog. Systems: http://resources.arcgis.com/en/help/main/10.2/018z/pdf/geographic_coordinate_systems.pdf If you are ...


2

Your dataset should appear in QGIS with EPSG:5652. If not, use Set CRS for Layer to get it. The 6-digit CRS is EPSG:25832. Rightclick on the layer, Save As ... under a different name and that CRS, and add it to the canvas. It might be even enough if you just change the project CRS to EPSG:25832, leaving the dataset unchanged.


2

Taking the suggestion by mkennedy for EPSG:28992, the points are located this way: which does not fit well, unless the WGS84 coordinates are rough or wrong. There is also an Amersfoort RD Old projection, but that is far off. Looking closer, the RD New points perfectly match to adresses in Openstreetmap, while the WGS84 are just road junctions; the lower ...


2

Yes, you may use this site for batch processing: Batch Conversions of Latitude/Longitude to Address (Reverse Geocoding) 46.7270,2.5059 = Rue des Varennes, Saint-Amand-Montrond, Cher, Centre, Metropolitan France, 18200, France This site uses MapQuest reverse geocode service. You will have to apply some additional logic to parse out the city/town name ...


2

As commented by @jbchurchill and @user23715, I recommend checking your Output Coordinates setting under your Environments tab to see what it is set to. I/we do this because in the Point To Raster (Conversion) help it lists the Environments respected by that tool as being: Environments Compression, Current Workspace, Output Coordinate System, ...


2

Every set of coordinates stored inside SDE.ST_GEOMETRY is stored as a compressed array of 8-byte long integers. The IEEE floating-point to integer conversion is performed by subtracting the coordinate reference X_OFFSET (or Y_OFFSET, or Z_OFFSET, or M_OFFSET) and multiplying by the XYUNITS (or Z_SCALE or M_SCALE). The process is reversed when querying for ...


2

If you need high accuracy distances, or "ground" distances, you need to convert your UTM "grid" distances (which you do indeed calculate via pythagorous) using a combined scale factor. This removes the distortion introduced by the combination of (a) reducing the horizontal distance at its elevated (above the ellipsoid) position on the earth and (b) ...


2

The code looks like the USNG or MGRS coordinate system, see https://en.wikipedia.org/wiki/United_States_National_Grid and https://en.wikipedia.org/wiki/Military_grid_reference_system for further details. 12RVU indicates that you are in UTM zone 12, somehwere in Western Mexico. You might need to read the manual on how to change the setting (I don't have the ...


1

You have to build the proj.4 string for the CRS from the parameters given: +proj=lcc +lat_1=44.883333 +lat_2=45.133333 +lat_0=44.791111 +lon_0=-93.383333 +x_0=152400.000000 +y_0=30480.000000 +a=6378418.9409999996 +b=6357033.3098455509 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs Note that x_0 and y_0 have to be in meters, while false Easting and Northing ...


1

QGIS is not designed to handle J2000 data. Your data would have 4 coordinates (X,Y,Z and time), while QGIS only handles 2-dimensional coordinates (long and lat or X and Y). It is possible to reproject coordinates with cs2cs if you set the prime meridian according to the time. Every hour from 12 UTC is a shift of 15° to the greenwich meridian: put 7 51 in a ...


1

Not sure by your question if you want the map spatial reference or a layer. For a map you may use the IMap.SpatialReference Property.


1

You would do something like this: dataset = "c:/data/landbase.gdb/Wetlands" spatial_ref = arcpy.Describe(dataset).spatialReference as per the ESRI help page here


1

As per the response syntax and example in the documentation here, ArcGIS Server should definitely be returning the spatialReference to you explicitly. see this thread for more info about the confirmed bug.


1

After many attempts and even contacting the man, Denis Carrière , who created the geocoder module on Python, please find the solution below. I feel a bit drained after providing this first foray into developement , although my script was different initially, but I feel happy with the results and the helps. :) All done in Python, of course. import ...


1

What you need to do is a spatial join on the lat/long coordinates and the french postal zones Here is the tutorial http://www.qgistutorials.com/en/docs/performing_spatial_joins.html The GIS DATA... postal zones can be found on this website https://www.data.gouv.fr/en/datasets/correspondances-code-insee-code-postal-idf/ Hope this helps!!



Only top voted, non community-wiki answers of a minimum length are eligible