One of the problems viewing GIS data, e.g. Shape Files, is that it is hard to navigate around them as there are often no common landmarks. An obvious solution
to this is to add a transparent map as an overlay, referenced to the same coordinate scheme. The Open Street Map data is available, for free, as rendered tiles
in a number of resolutions, covering virtually the whole planet. MapWinGIS supports overlaying image files onto map as long as they are geo-referenced,
so it is relatively straightforward to add a dynamic map overlay which scales with the MapWinGIS window, to maintain a OpenStreetMap overlay at all resolutions:
Details on Open Street Map Tiles (Slippy Maps) can be found here: http://wiki.openstreetmap.org/wiki/Slippy_map_tilenames#Tile_servers.
Details on World Files, used to Geo-reference image files, can be found here: http://en.wikipedia.org/wiki/World_file.
Locating the Slippy Tiles
The two functions below, from the Slippy Tile Wiki, convert between Slippy Map Tile references and WGS-84 Latitude and Longitude. These are used to work out the tiles required
to cover the map area and then to Geo-reference each tile before adding it to the Map.
Details of the method and how Slippy Tiles are numbered can be found here: http://wiki.openstreetmap.org/wiki/Slippy_map_tilenames#Derivation_of_tile_names
Private Function CalcTileXY(ByVal LatDeg As Single, _
ByVal LonDeg As Single, _
ByVal ZoomLevel As Long) As PointF
CalcTileXY.X = CLng(Math.Floor((LonDeg + 180) / 360 * 2 ^ ZoomLevel))
CalcTileXY.Y = CLng(Math.Floor((1 - Math.Log(Math.Tan(LatDeg * Math.PI / 180)+ 1
/ Math.Cos(LatDeg * Math.PI / 180))
/ Math.PI) / 2 * 2 ^ ZoomLevel))
Private Function CalcTileTLCoord(ByVal TileCoord As PointF, _
ByVal ZoomLevel As Long) As PointF
Dim n As Double = Math.PI - ((2.0 * Math.PI * TileCoord.Y) / Math.Pow(2.0, ZoomLevel))
CalcTileTLCoord.X = (TileCoord.X / Math.Pow(2.0, ZoomLevel) * 360) - 180
CalcTileTLCoord.Y = 180 / Math.PI * Math.Atan(Math.Sinh(n))
Creating a Map Overlay
In order to create a map overlay, we need to work out the area to cover. This can be found by looking at the Extents of the AxMapWinGIS Map:
Dim MapExtents As MapWinGIS.Extents = MapMain.Extents
Dim LonMin As Double = MapExtents.xMin
Dim LonMax As Double = MapExtents.xMax
Dim LatMin As Double = MapExtents.yMin
Dim LatMax As Double = MapExtents.yMax
NB This method will only work if the AxMapWinGIS Map is using the same coordinate system as OSM Slippy Tiles i.e., WGS-84.
If not, you will have to convert WGS-84 into the Map's coordinate system and vica-versa.
Once we have the area define by Min and Max Lat and Lon, we can use one of the Slippy Map functions to convert from the coordinates to Tile references for the
Top Left corner and Bottom right corner of the Map:
TLTile = CalcTileXY(LatMax, LonMin, ZoomLevel) BRTile = CalcTileXY(LatMin, LonMax, ZoomLevel)
Then it is a simple process to loop through all the tiles in between the two corners, to create a mosaic of map tiles which cover the whole map area:
Dim H As Integer, V As Integer Dim TileRef As Point
For H = TLTile.X To BRTile.X For V = TLTile.Y To BRTile.Y TileRef.X = H
TileRef.Y = V
Call LoadOSMTile(MapMain, TileRef, ZoomLevel)
Downloading and Geo-referencing a Tile
Each tile is stored, on the Slippy Tile Server, in a directory structure:
Const ServerName As String = "http://a.tile.openstreetmap.org/"
Dim URL As String = ServerName & ZoomLevel & "/" & TileRef.X & "/" & TileRef.Y & ".png"
To download a tile, we use Net.HttpWebRequest:
Dim request As Net.HttpWebRequest = DirectCast(Net.HttpWebRequest.Create(URL), Net.HttpWebRequest)
Dim response As Net.HttpWebResponse = DirectCast(request.GetResponse, Net.HttpWebResponse)
Dim Image As System.Drawing.Image = Image.FromStream(response.GetResponseStream())
Once we have downloaded the tile, we save it in the local cache:
Image.Save(CacheTileDirectory & "\" & CacheFileName)
And Geo-reference it by creating a World File for the tile (using the .pgw extension for a World file associated with a .png file):
Dim TileTLCoord As PointF = CalcTileTLCoord(TileRef, ZoomLevel)
Dim Point2 As PointF
Point2.X = TileRef.X + 1
Point2.Y = TileRef.Y + 1
Dim TileWidthDeg As Double = CalcTileTLCoord(Point2, ZoomLevel).X - CalcTileTLCoord(TileRef, ZoomLevel).X
Dim TileHeightDeg As Double = CalcTileTLCoord(TileRef, ZoomLevel).Y - CalcTileTLCoord(Point2, ZoomLevel).Y
Debug.Print("Tile Width = " & TileWidthDeg & ", Height = " & TileHeightDeg)
Dim oWrite As System.IO.StreamWriter = IO.File.CreateText(CacheTileDirectory & _
"\" & Left(CacheFileName, CacheFileName.Length - 4) & ".pgw")
oWrite.WriteLine(FormatNumber(TileWidthDeg / 256, 10)) oWrite.WriteLine("0.0") oWrite.WriteLine("0.0") oWrite.WriteLine(FormatNumber(-1 * TileHeightDeg / 256, 10)) oWrite.WriteLine(FormatNumber(TileTLCoord.X, 10)) oWrite.WriteLine(FormatNumber(TileTLCoord.Y, 10)) oWrite.Close()
Note that the Height is negative, to indicate pixels going down vertically from the top left corner (the geo-reference point).
Once we have saved the World File to the same directory as the PNG image, we can add the PNG file to the Map. MapWinGIS will automatically
detect the world file, as long as it is in the same directory:
the World file and loads that at the smae time)
Dim ImageType As MapWinGIS.ImageType = MapWinGIS.ImageType.PNG_FILE
Dim MWImage As New MapWinGIS.Image
If Not MWImage.Open(CacheTileDirectory & "\" & CacheFileName, ImageType) Then
Debug.Print("Image load failed: " & MWImage.ErrorMsg(MWImage.LastErrorCode))
MWImage.TransparencyPercent = 0.5
Dim PNGLayer As Integer = MapMain.AddLayer(MWImage, True)
Using the code
The complete Visual Studio 2010 Project is included which comprises the code above and a Windows Form wrapper which lets the user navigate a Shape file,
zooming in and zooming out. The mouse position is displayed in a text window at the bottom of the screen, so you can check the overlay is correct by locating known features etc:
The only external reference required is a Shape File to load, and the code currently uses one of the MapWindow example project files, which installs with MapWindow.
"C:\Program Files\MapWindow\Sample Projects\World\Shapefiles\WORLD30.shp"
So, either make sure this is present, or choose another file to play with.
First version 19 Sep 2012.