Форумы GIS-Lab.info

Код: Выделить всё

PROJCS["Local Conditional",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Local"],PARAMETER["False_Easting",27478.878],PARAMETER["False_Northing",19463.17],PARAMETER["Scale_Factor",1.0],PARAMETER["Azimuth",-1.272472],PARAMETER["Longitude_Of_Center",49.3933174],PARAMETER["Latitude_Of_Center",53.4992761],UNIT["Meter",1.0]] 

Код: Выделить всё

Imports System.Drawing

Public Class Projection

    Public Shared Function LatLonToTile(lat As Double, lon As Double, zoom As Integer) As Point
        Dim x As Single = ((lon + 180.0) / 360.0 * (1 << zoom))
        Dim y As Single = ((1.0 - Math.Log(Math.Tan(lat * Math.PI / 180.0) + 1.0 / Math.Cos(lat * Math.PI / 180.0)) / Math.PI) / 2.0 * (1 << zoom))
        Return New Point(x, y)
    End Function
    Public Shared Function LatLonToTile(lat As Double, lon As Double, zoom As Integer, fixed As Boolean) As Point
        Dim x As Single = ((lon + 180.0) / 360.0 * (1 << zoom))
        Dim y As Single = ((1.0 - Math.Log(Math.Tan(lat * Math.PI / 180.0) + 1.0 / Math.Cos(lat * Math.PI / 180.0)) / Math.PI) / 2.0 * (1 << zoom))
        Return New Point(Math.Floor(x), Math.Floor(y))
    End Function


    ''' <summary>
    ''' Расчет точки координат в пикселях
    ''' </summary>
    ''' <param name="lat"></param>
    ''' <param name="lon"></param>
    ''' <param name="zoom"></param>
    ''' <returns></returns>
    ''' <remarks></remarks>
    Shared Function LatLonToPixelF(lat As Double, lon As Double, zoom As Double) As PointF
        Dim p As New PointF
        Dim centrPoint As Double = Math.Pow(2, zoom + 7)
        Dim totalPixel As Double = 2 * centrPoint
        Dim pixelsPerLngDegree As Double = totalPixel / 360
        Dim pixelsPerLngRadian As Double = totalPixel / (2 * Math.PI)
        Dim siny As Double = Math.Min(Math.Max(Math.Sin(lat * (Math.PI / 180)), -0.9999), 0.9999)
        p.X = Math.Round(centrPoint + lon * pixelsPerLngDegree)
        p.Y = Math.Round(centrPoint - 0.5 * Math.Log((1 + siny) / (1 - siny)) * pixelsPerLngRadian)
        Return p
    End Function

    ''' <summary>
    ''' Пиксель карты в координаты Lat,Lon
    ''' </summary>
    ''' <param name="px"></param>
    ''' <param name="py"></param>
    ''' <param name="zoom"></param>
    ''' <returns>Lat,Lon</returns>
    ''' <remarks></remarks>
    Shared Function PixelsToLatLon(px As Double, py As Double, zoom As Integer) As PointF
        Dim r() As Double = PixelsToMeters(px, py, zoom)
        r = MetersToLatLon(r(0), r(1))
        Dim res As New PointF(r(0), r(1))
        Return res
    End Function

    Shared Function MetersToLatLon(x As Double, y As Double) As Double()
        Dim lon As Double = (x / originShift) * 180
        Dim lat As Double = (y / originShift) * 180
        lat = 180 / Math.PI * (2 * Math.Atan(Math.Exp(lat * Math.PI / 180)) - Math.PI / 2)
        Return {lat, lon}
    End Function


    Shared Function PixelsToMeters(px As Double, py As Double, zoom As Integer) As Double()
        Dim res = Resolution(zoom)
        Dim x As Double = px * res - originShift
        Dim y As Double = Math.Abs(py * res - originShift)
        Return {x, y}
    End Function

    Shared Function Resolution(zoom As Integer) As Double
        Return initialResolution / Math.Pow(2, zoom)
    End Function

    Const tileSize = 256.0
    Const initialResolution = 2 * Math.PI * 6378137 / tileSize
    Const originShift = 2 * Math.PI * 6378137 / 2

End Class

Код: Выделить всё

String.Format("https://{0}.tile.openstreetmap.org/{1}/{2}/{3}.png", serverType(severId), zoom, x, y)

Код: Выделить всё

Private Shared serverType As New List(Of Char) From {"a", "b", "c"}