Create Bounding Box von der Mitte der Koordinate

Ich versuche, eine Bounding Box in einem bestimmten Maßstab von einer Mittelkoordinate zu erstellen. Ich versuche, es im Seitenverhältnis eines 8,5 x 11 Zoll Blatt Papier (612 x 792 Pixel @ 72 dpi) zu halten.Create Bounding Box von der Mitte der Koordinate

Der Code, den ich unten verwende, funktioniert meistens, aber die Höhe scheint für das Seitenverhältnis eines Buchstabens etwas zu groß zu sein. Betrachte ich nicht die Mercator-Projektion? Was fehlt mir hier?

def bounding_box_from_point(center:, size:, scale_denominator:) 
    dpi = 72 
    inches_per_unit = 4374754 
    resolution = 1/(scale_denominator * inches_per_unit * dpi) 
    half_width_deg = (size.width * resolution)/2 
    half_height_deg = (size.height * resolution)/2 

    BoundingBox.new(
    north: center.lat + half_height_deg, 
    south: center.lat - half_height_deg, 
    east: center.lon + half_width_deg, 
    west: center.lon - half_width_deg 
) 
end

Aufruf bounding_box_from_point(center: center, size: size, scale_denominator: scale) mit:

scale = 0.0008861342166177423 (i.e. 1/18055.955520) 
center = Geometry::Location.new(lat: 37.806336, lon: -122.270625) 
size.width = 612, 
size.height = 792

Es gibt:

west: -122.27172131608657, 
east: -122.26952868391342, 
south: 37.804917238005615 
north: 37.80775476199439

Wenn Sie http://www.openstreetmap.org/export gehen und geben Sie diese Zeichen-Box-Koordinaten, können Sie sehen, dass das Verhältnis nicht überein von einem 8,5x11in Blatt Papier ... es ist etwas zu groß. Was mache ich hier falsch oder nicht?

Quelle

2017-02-13 theartofbeing

Es ist nicht so sehr die Mercator-Projektion, wie es die Tatsache ist, dass Längengradlinien näher zusammen kommen, wenn Sie sich vom Äquator entfernen. Ihre Karte ist nicht zu groß, es ist zu dünn, –

Siehe http://StackOverflow.com/A/2913249/4504607 –

@RonJensen so, wie passen Sie für diese? – theartofbeing

Das war meine Implementierung, die es gelöst hat!

def self.bounding_box_from_location(location:, scale:) 
    scale_meters_per_pixel = 0.0003 

    # scale.zoom is an integer between 0-19 
    # See http://wiki.openstreetmap.org/wiki/Zoom_levels for a description on `zoom` 
    # location.lat_rad is the latitude in radians, same for lon_rad 

    meters_per_pixel = EARTH_CIR * Math.cos(location.lat_rad)/(2 ** (scale.zoom + 8)) 
    earth_meters_per_pictureMeters = meters_per_pixel/scale_meters_per_pixel 


    # height/width in meters is the height/width of the box you're applying this to 
    # In my case I'm using the heigh/width of a Letter of paper 
    # You can convert Pixels to meters (@ 72dpi) with 
    # pixels * 0.00035277777777777776 

    meters_north = earth_meters_per_pictureMeters * scale.height_in_meters/2 
    meters_east = earth_meters_per_pictureMeters * scale.width_in_meters/2 

    meters_per_degree_lat = 111111.0 
    meters_per_degree_long = 111111.0 * Math.cos(location.lat_rad) 

    degrees_north = meters_north/meters_per_degree_lat 
    degrees_east = meters_east/meters_per_degree_long 

    BoundingBox.new(
    north: (location.lat + degrees_north), 
    south: location.lat - degrees_north, 
    east: location.lon + degrees_east, 
    west: location.lon - degrees_east, 
    width: scale.width_in_pixels, 
    height: scale.height_in_pixels 
) 
end

Quelle

2017-03-02 00:39:22 theartofbeing

Create Bounding Box von der Mitte der Koordinate

Antwort

Verwandte Themen