Auf Memorandum-Ebene habe ich eine Funktion zum Konvertieren von Mesh-Code 3. Ordnung in WKT mit Python erstellt. Es ist nüchtern, aber es ist leicht, den Überblick zu verlieren, wenn Sie es brauchen ...
def toWKTfromMeshCode3(mcode3: 'mesh code(str)'):-> 'wkt(str)'
south = (int(mcode3[0:2]) * 80 + int(mcode3[4:5]) * 10 + (int(mcode3[6:7]))) * 30 / 3600
west = (int(mcode3[2:4]) * 80 + int(mcode3[5:6]) * 10 + (int(mcode3[7:8]))) * 45 / 3600 + 100
north = south + 30 / 3600
east = west + 45 / 3600
wkt = 'POLYGON(({e} {n},{w} {n},{w} {s},{e} {s},{e} {n}))'.format(s=south, w=west, n=north, e=east)
return wkt
Recommended Posts