Question

我正在用底图创建一个世界的SVG图像，在该底图上我将柏林市作为“控制点”绘制（因为我想手动在其中放置一个圆圈...所以我有个参考）。 / p>

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt

# Berlin & New York & Sydney
lats = [52.516667] #[52.516667, 40.730610]
lons = [13.388889] # [13.388889, -73.935242]

plt.figure(figsize=(15,15/2))
m = Basemap(projection='robin', lon_0=0, resolution='c')
m.drawcountries(color='#ffffff', linewidth=0.75)
m.fillcontinents(color='#c0c0c0', lake_color='#e6f5ff')                
m.drawmapboundary(fill_color='#e6f5ff', linewidth=1, color='#000000') # Ocean          

x,y = m(lons,lats)

plt.plot(x, y, 'bo', color='r', markersize=5)

plt.savefig("basemap.svg", figsize=(24,12))
plt.show()

在下一步中，我想通过编辑已创建的SVG文件的代码在SVG图像上手动放置一个圆圈。可以通过在SVG图片代码末尾</svg>之前引入以下代码来完成此操作。

<circle fill="blue" cx="250" cy="470" r="2"/>

如何用我的Python代码确定cx和cy的正确值，以将蓝点放置在柏林所在的位置？

我认为我有一个mapWidth = 1080和一个mapHeigth = 540以及xMax = 33973600和yMax = 17231000。

通过这种方式，我可以计算cx = mapWidth - x/xMax*mapWidth并类似地计算cy = mapHeigth - y/yMax*mapHeigth。但是，如果我分别考虑底部和左侧边距分别为72和152 pt，这也不会将蓝点放置在正确的位置。有什么想法吗？

Answer 1

似乎您确实需要解决方案，所以我将为您提供一种解决方法。

注意：我不建议您手动操作.svg文件。但是，如果别无选择，那就去做吧！

解决方案假定，非python进程不会更改表单中的文件。

要克服创建.svg文件（至少对我来说是黑盒）的麻烦，您可以创建另一个绘制期望坐标的图像，将图像另存为临时.svg文件，找到坐标点（在.svg文件中），最后将它们添加到您的初始地图文件中。

我定义了3种方法： -createMap绘制地图并将输出另存为.png文件 -get_svg_coordinates：创建一个临时地图（.svg文件），读取点坐标，删除临时文件，返回点坐标。 -add_circle：在现有.svg地图文件上画圆。

代码如下：（工作示例）

# Import modules
import os
import re
# os.environ['PROJ_LIB'] = r'C:\Users\...\Library\share'
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt


def createMap(lats, lons, color='r', figure_name="basemap.svg", show=False):
    """ Create a map from points coordinates + export to .csv file
    Arguments
        :param lats: latitudes
        :param lons: longitudes
        :param color='r': color points
        :param figure_name="basemap.svg": name output file
        :param show=False: show the map
    """
    # Same code as yours
    plt.figure(figsize=(15, 15/2))
    m = Basemap(projection='robin', lon_0=0, resolution='c')
    m.drawcountries(color='#ffffff', linewidth=0.75)
    m.fillcontinents(color='#c0c0c0', lake_color='#e6f5ff')
    m.drawmapboundary(fill_color='#e6f5ff', linewidth=1, color='#000000')  # Ocean

    x, y = m(lons, lats)

    plt.plot(x, y, 'bo', color=color, markersize=5)
    plt.savefig(figure_name, figsize=(24, 12))

    if show: plt.show()


def get_svg_coordinates(lat, lon, color='#ff0000', figure_temp_name="tmp_figure.svg"):
    """ Create a temporary file using the createMap function
        Find the point coordinates inside (using regex check)
        Remove temporary csv file
    Arguments
        :param lat: new point latitude
        :param lon: new point  longitude
        :param color='#ff0000': point color
        :param figure_temp_name="tmp_figure.svg": temp file name
    """
    createMap(lat, lon, color=color, figure_name=figure_temp_name)

    with open(figure_temp_name, "r") as f:
        # read file
        content = f.read()  
        # Find x - y values (pattern ' x=' is unique is you are using 1 point)
        x = re.findall(r'<use.*x=\"(\d*\.*\d*)\"', content)
        y = re.findall(r'<use.*y=\"(\d*\.*\d*)\"', content)

    # remove file
    os.remove(figure_temp_name)
    return x, y


def add_circle(map_file_name, x, y):
    """ Draw circle at the end of file
    Arguments:
        :param map_file_name: filename (adding circle)
        :param x: x coordinates (results of get_svg_coordinates method)
        :param y: y coordinates (results of get_svg_coordinates method)
    """
    with open(map_file_name, "r+") as f:
        content = f.readlines()
        # get number of lines in file
        for i, l in enumerate(content):
            pass
        # Add content
        content.insert(i, '<circle fill="blue" cx="{0}" cy="{1}" r="2"/>'.format(x[0], y[0]))
        f.seek(0)                   # file pointer locates at the beginning to write the whole file again
        f.writelines(content)       # rewrite file


# Berlin & New York & Sydney
lats = [52.516667]  # [52.516667, 40.730610]
lons = [13.388889]  # [13.388889, -73.935242]

# create your initial svg map
map_file_name = "basemap.svg"
createMap(lats, lons, figure_name=map_file_name)

# Find new position point on svg file
# Define coordinates points
NewYork_lat = 40.730610
NewYork_long = -73.935242
x, y = get_svg_coordinates(NewYork_lat, NewYork_long)
add_circle(map_file_name, x, y)

注意：

我对.svg文件不熟悉。为了回答这个问题，我按预期在文件末尾添加了<circle fill="blue" cx="???" cy="???" r="2"/>。但是，最好标识整个DOM <g id="line2d_1">并复制过去。
该代码段适用于一个图像，我让您概括出一组要点。

底图m（lats，lons）到svg cx和cy

1 个答案: