How to format the labels of a UTM Grid in a map Layout/Print Composer


I am doing a map in for a small area (1.5km by 1.5km on UTM zone 14N). I want to know if it is possible (and how to do it) to format the labels of a UTM grid. This is exactly what I want to show

a) On the Southwest corner of the map I want to see the full value of the coordinates (north or east). I want that the east coordinate shows the first 3 digits in a bigger font than the las 3 digits. For north coordinates I want the labels to show the first four digits in a bigger font than the remaining three digits. I also want them to show the suffix E or N accordingly

b) For the remaining east and north coordinates I just want to show the last 3 digits of the coordinate

I can use either QGis or ArcMap.

I tried to use ArcMap’s grid styles, particularly the USN National Grid and I have no luck. I also tried to use corner style labels feature no luck also

Difference between UTM 10U and UTM 10N? [closed]


I there a difference between UTM 10U and UTM 10N?

ArcGIS has 10N, but my scan map has 10U?

When I attempt to georeference, my scan gets skewed

Can I use Redfearn's formula to convert between WGS84 and UTM?


I have GPS data at 1 point per second for several weeks in WGS84 latitude/longitude form. It was collected using a UBLOX GPS in 2015 in Australia (Zone 55s)

Would it be correct to use Redfearn’s formula to project WGS84 data to Northings/Eastings? Normally it seems that Redfearn’s formula is used with GDA94 rather than with WGS84:

I would like to get Northings/Eastings in meters. The data is only accurate to 2-5m at the moment, but I don’t want to make it worse so would need to use a conversion that is accurate to within a few cm.

I could just use the ArcMap (ArcPy) projection but since I have very large files (1 point per second) even ArcPy is very slow. Therefore I’d rather put the formula directly into Python. I’ll probably use ArcPy later to check some of the results.

The data is collected on animal collars (in Zone 55s). The collars could be used anywhere in Australia (or possibly overseas) so I’d like to use a formula that could be transferable to another part of Australia (although that could involve user input about approximate location).

I’ve also been investigating how to convert WGS84 (presumably = ITRF2008?) into GDA94 and then to MGA, but this seems very complicated and I think it may be simpler just to work in WGS84 UTM since I don’t really need it for an official map.

What is most accurate coordinate system for calculating areas of polygons?


I have multiple polygons, all located offshore around the UK. I am trying to calculate the areas of these polygons in square kilometres, without splitting them by UTM zone.

Is it possible to do this accurately? I have been using World Mercator, but there are big distortions at this scale. If the polygons were limited to the east coast, say, I would use UTM31N. However, as they are all over the UK waters, I am unsure what to use.

Is there a coordinate system that covers this area? British National Grid does not extend far enough out to sea to be of use.

Measuring distances when crossing UTM zones?


I like to calculate distance between two UTM coordinates which are in different zones. Could any one point me how I can do this.

I found that we can map these coordinates in to a new projection for solving the problem. Could any one let me know which is the best map projection where cartesian coordinate related calculations are applicable like UTM.

If some points me to extended UTM zones and how it helps, that will also be very helpful.

Which UTM zone to use for a large country?


I want to reproject OpenStreetMap roads data from the current projection (WGS 84) to UTM, since I read in another GIS.SE discussion (How to get $length in meters in QGIS?) that WGS 84 cannot be used to measure length in m or km. My end goal is to calculate the total roadway kilometres within each district, so that I can enter that variable into my regression model. The country in question is Indonesia, and as you can see from the map below, it covers UTM zones 46-54.

UTM zones map from Wikimedia

Image source: Wikimedia, http://upload.wikimedia.org/wikipedia/commons/e/ed/Utm-zones.jpg

Does it really matter which one I choose? Can I choose the one that’s roughly in the middle and call it good (that would be zone 50) or should I pick the zone that is more densely populated with roads (Jakarta/Java, zone 48)?

I am only barely functionally literate in GIS. The only formal training I had was a 3-day course back in 2008. Really, the only reason I need the roadway kms variable was because one of my thesis committee members wanted me to account for it in my model…

And to confirm, once I have reprojected to the UTM CRS, I can use my GIS software’s built-in length tool to directly calculate the length in metres (or divide by 1000 for km), right?

I am using QGIS.

Georeferencing raster imagery using UTM eastings and northings


enter image description hereI have a scanned cadastral map in raster form and UTM co-ordinates in northings and eastings for that area. Is there any way that I can directly georeference my map using E/N values without converting them? I tried doing it directly through the Georeferencing toolbar in ArcMap but RMS was too high.

How to change GK-coordinates to UTM with the field calculator in QGIS?


i have a qgis project in ETRS89 (UTM) with some vectorshapes. i already transformes them from GK into UTM. but now the attribute table still show me the GK-coordinates. HOW can i change it? with the field calculator?!

Please help me :)
Lena

Erroneous location when converting from a GCS_WGS_1984 to UTM


I have a shapefile map of Belgium divided by cantons that I want to re-project in order to have units in meters and to be able to calculate distances. I know how to do it with ArcGIS Projection tool, but the result always gets wrong map location.

I will detail my steps. First of all, my original map had unknown projection and coordinate system. However, it is pretty clear that it is in GCS_WGS_1984, with decimal degrees, since the Left/Top/Right/Bottom limits are precisely equal to the real ones of Belgium, when a map of this country is projected that way. So far, so good. I defined the GCS_WGS_1984 projection with DD units in ArcGIS, using the Define Projection tool.

The problem is in the next step. I then want to convert from that to another projection. One from which I can have units in meters. I tried a bunch of UTM options, using the Project tool after finding out that Belgium is around UTM 31N. However, no matter which UTM projection I choose, the location of the country after re-projected becomes wrong.

For instance, re-projecting from thw GCS_WGS_1984 to WGS_1984_UTM_Zone_31N gives me the following extent limits: Left:468861,233828 m; Top:5707355,258593 m; Right:742547,667605 m; Bottom:5485665,513781 m. That’s very off the target.

Could anyone point me what am I doing wrong or what is missing?

PS: Another detail is the following. I tried opening the original map among other country maps with the GCS_WGS_1984 projection in decimal degrees. All countries are positioned correctly. Then, I do the same operation described above to convert some countries besides Belgium to check what happens. All of them get located wrong – what means, my problem doest not seem to be in my Belgium files, but rather in me doing something wrong.

Convert shapefile to LonLat (R)


Just a small disclaimer: I don’t have any background in cartography, hence I might have missed the answer to my question in search just because I didn’t even know the keywords to look for.

So, I am using R to plot statistical data on maps. It was all fine and understandable when all of the materials came in same units of measurments (lon/lat in form of decimal degrees or degrees/minutes/seconds which have obvious and easy conversion).

So what I have encountered – doing routine fortify(shapeMap) to plot data with ggplot2 yields dataframe with values:

    long    lat order   hole    piece   group   id
1   563656.6    373826.0    1   FALSE   1   0.1 0
2   563864.6    373664.3    2   FALSE   1   0.1 0
3   564272.9    373895.4    3   FALSE   1   0.1 0
4   564491.6    373805.5    4   FALSE   1   0.1 0
5   565151.1    372709.0    5   FALSE   1   0.1 0
6   565328.3    372347.0    6   FALSE   1   0.1 0

...

As far as I understand that this is UTM coordinate system. And hence the question – is there an easy way to convert these coordinates to decimal degrees? I have dug through quite a few materials and online convertors but I can’t really seem to get the expected result.

Just for reference – shapefile coordinates above should be somewhere in Latvia. And Latvia is split among 3 UTM grid cells: 34V; 35V and 35U (if that is relevant at all).

Update: I have found out that coordinates are in LKS92 system. Hence now there is a clear objective. Convert LKS92 to WGS84. Guess now it is up to figuring out spTransform() function from rgdal package. Unfortunately it is not an obvious function with “from” and “to” arguments. Or at least at first glance.

Update2: Reference for those that will encounter the same problem:

  1. You have to be sure what coordinate system you want to conver from/into. In my case I found out it was LKS92.

  2. Find specific parameters for the coordinate system on spatialreference.org. In my case parameters look like “+proj=tmerc +lat_0=0 +lon_0=24 +k=0.9996 +x_0=500000 +y_0=-6000000 +ellps=GRS80 +units=m +no_defs”. Not that these are parameters for coordinate system FROM which I want to convert. If this is not WGS84 then you will want another parameter string for target coordinate system.

  3. R package proj4 has function project(). An example of LKS92 -> WGS84 conversion:

xy <- cbind(c(508148.442), c(319423.294))
project(xy, "+proj=tmerc +lat_0=0 +lon_0=24 +k=0.9996 +x_0=500000 +y_0=-6000000 +ellps=GRS80 +units=m +no_defs", inverse = T)


[,1] [,2]
[1,] 24.1342 57.01823

Note the inverse=T which means convert to WGS84. Making inverse = F suggests that input coordinates are in WGS84 format.

Another small update:

To convert whole shapefile you would still need rgdal package:

proj4string(shapeMap) <- CRS("+proj=tmerc +lat_0=0 +lon_0=24 +k=0.9996 +x_0=500000 +y_0=-6000000 +ellps=GRS80 +units=m +no_defs")
shapeMap <- spTransform(shapeMap, CRS("+proj=longlat"))

First line defines the CRS of shapefile, second line tells what CRS to transform into.

Question and Answer is proudly powered by WordPress.
Theme "The Fundamentals of Graphic Design" by Arjuna
Icons by FamFamFam