3

I am trying to georeference a raster map from Nepal. What CRS should I use? PROJ.4 or WKT form. I will use the UTM grid to base my reference points on.

On the map there is reference to: Scale 1:50000 Everest 1830 Spheroid Modified UTM Projection Origin Long 87E Lat 0N False Coordinates at Origin Easting 500000 Northing 0 Scale Factor at Central Meridian 0.9999

Centre of map is about 460000E 3080000N

5
  • Are you using QGIS? Can you provie a link to the map? Commented Aug 22, 2017 at 9:10
  • Yes using QGIS. Maps are at link I am using 2786 03 Namche Bajar map
    – John Flint
    Commented Aug 22, 2017 at 9:56
  • Someone should double-check: +proj=tm +a=6377276.345 +b=6356075.41314 +x_0=500000.0 +y_0=0.0 +k_0=0.9999 +lon_0=87.0 +lat_0=0.0 +wgs84=296.207,731.545,273.001 +units=m
    – mkennedy
    Commented Aug 22, 2017 at 20:06
  • 1
    I ended up using +proj=tmerc +lat_0=0 +lon_0=87 +k=0.9999 +x_0=500000 +y_0=0 +a=6377276.345 +b=6356075.41314024 +towgs84=296.207,737.545,273.001,0,0,0,0 +units=m +no_defs Can't be 100% sure that is right for that map as it seems there are various Everest Spheroids and datums but this matches up nicely with other maps
    – John Flint
    Commented Aug 23, 2017 at 4:34
  • You can find a list of all modified Everest 1830 ellipsoids at georepository.com/search/by-name/?query=Everest+1830 . The 1937 or 1975 adjustment fit for India (where the map comes from).
    – AndreJ
    Commented Aug 27, 2017 at 6:04

2 Answers 2

4

After a bit of searching around and a tip from afalciano on needing the towgs84 transformation parameters I was able to use +proj=tmerc +lat_0=0 +lon_0=87 +k=0.9999 +x_0=500000 +y_0=0 +a=6377276.345 +b=6356075.41314024 +towgs84=296.207,737.545,273.001,0,0,0,0 +units=m +no_defs as the source coordinate system. This enabled me to use the UTM grid coordinates on the map to georeference.

There were a number of slightly different definitions that I found for the Everest 1830 datum as it has been revised over time. After some trial and error I ended up using the above proj4 definition from https://geonet.esri.com/thread/48571.

When the map was then transformed to be in the WGS84 coordinate system I found the roads and prominent features matched up very well with some GPS tracks I found elsewhere and OSM maps.

The method of using actual features on the map and matching up with the corresponding features on a known map such as OSM as detailed by Nitzan Matan would have also worked but I think at reduced accuracy

2

When you are Georeferencing in QGIS you need to adjust the raster to the CRS of your satellite map, the raster is without CRS it's just a photo at this point. Look at QGIS lesson for more explanation. I am using OSM maps, so when I am uploading the raster into the Georeferencer I set the CRS to WGS 84 / Pseudo Mercator projection. Now you need to mark location in the raster and then in the satellite map. Me advice is to find a location of a city like Chame in the raster and then use the GeoSearch plugin to find the same location in the satellite map.

From my experience in this kinds of map it is best to make at list ten points, to erase points with a long red line and scatter the point as much as passable. For example would erase all of this points . enter image description here

Now you have to set the transformation settings (the yellow wheel icon). At first I used this ones enter image description here

and the map position was not accurate. enter image description here Then I tried the Helmert transformation type and the map position looks much better.enter image description here

I think that your main concern are the lacks of details that make it hard to much the exact location in the raster and in the satellite maps. It is pretty passable that you cannot make a perfect georeferenc especially if the raster maps are not accurate. On the other hand, I used only for points and it the map was in the general area, I am pretty sure that more points will make it much better.

4
  • I was not trying to reference the overview map but one of the detail maps. 2786 03 Namche Bajar [link] (pahar.in/pahar/Maps--Primary/Nepal/Nepal Topo Maps/2786 03 Namche Bajar.jpg)
    – John Flint
    Commented Aug 22, 2017 at 14:01
  • The same georeference method works for all kinds of maps Commented Aug 22, 2017 at 14:04
  • Ok, well I thought to use the UTM grid as the reference points, I would think more accurate
    – John Flint
    Commented Aug 22, 2017 at 15:42
  • Theoretically you are right, but I am not sure it is possible. Commented Aug 22, 2017 at 15:49

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.