Au niveau du mémorandum, j'ai créé une fonction pour convertir le code de maillage du 3ème ordre en WKT avec python. C'est sobre, mais il est facile de perdre la trace du moment où vous en avez besoin ...
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