User Tools

Site Tools


georef

This is an old revision of the document!


Georeferencing Historic Maps

Intro

This project aims to georeference historic maps which are of interest to bushwalkers in NSW, and make them available in a format which allows for easy viewing.

General parts involved are:

  • Obtain digitised map through either
    • Obtaining physical map and digitising using scanner or camera
    • Downloading / acquiring a map that has been digitised by someone else (eg from the NLA via Trove)
  • Georeference map
    • By features
      • Manually, by picking common points in both old map and one in the desired coordinate system
      • Automated, through some algorithm that can identify common points and perform the process
    • By coordinates
      • Making use of the coordinates listed on the map (if available) and using a known coordinate transformation
  • Display online

Examples

Georeferencing by Coordinates

Here, twelve 1:63,360 scale military topographic maps were georeferenced and overlayed on the NSW LPI basemap. These digitised maps were sourced as low-quality JPGs (large TIFs also available) from the NLA (National Library of Australia) via Trove.

The maps were then georeferenced using the lat,long coordinates of corners of each map. These maps were generally published to cover 0.5 deg of long by 0.25 deg of lat, resulting in a regular grid of maps. Finally, they were clipped to cover only their mapped extent, removing the boards and information contained within.

No attempt was made to correct for the difference between the ANG and current GDA datums, so although the maps generally align with each other well, there can be discrepancies of 100s metres between the georeferenced maps and the LPI basemap.

Further work should be done in exploring variations of this process, particularly if a datum transformation can be applied.

Resources

Coordinate Transformation

When using old maps which have a CRS (coordinate reference system) significantly different from that which we use today, it can be useful to be able to transform coordinates from one CRS to another.

A brief overview of 'grids & datums' history for Australia can be found in this 2003 publication

Transformations listed before are quite well known, and many GIS packages can be used:

  • AGD66 ↔ AGD84
  • AGD66 ↔ GDA94
  • GDA94 ↔ GDA2020
  • GDA94 ↔ WGS84

Note that these refer to geographic datums, rather than map grids of a projection

Geographic datum (short) Map projection

Some details for pre AGD66 transformations can be found in the AGD technical manual (pp 52-54) which have been included below for reference.

 AGD Technical Manual - p52

 AGD Technical Manual - p53

 AGD Technical Manual - p54

Coordinate Conversion

PROJ is a useful program here.

'Clarke Coordinate' have the following parameters

Parameter Value Units
ellipsoid Clarke 1858
a 20,926,348 british feet
a 6,378,339.78 metres
b 20,855,233 british feet
b 6,356,663.92 metres
1/f 294.26
projection transverse mercator
false easting 400,000 british feet
false easting 365,759.36 metres
false northing 800,000 british feet
false northing 731,518.73 metres
zones 1 to 8
zone width 5 degrees
zone 1 central meridian 116E degrees
zone latitude origin 34S degrees
scale factor 1

Using these details, the below PROJ code will convert lat,long to E,N for a given zone.

+proj=tmerc +south +lat_0=-34 +lon_0=151 +k_0=1 +a=6378339.78 +b=6356663.92 +units=yd +x_0=365760 +y_0=731520

georef.1539104689.txt.gz · Last modified: 2018/10/10 04:04 by allchin09

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki