the-ONE 模拟器的使用 osm转换wkt - ToB企服应用市场:ToB评测及商务社交产业平台

import pandas as pd
import os
import re
import shutil
def process_files(directory):
for filename in os.listdir(directory):
if filename.endswith(".csv"): # 只处理.csv文件
with open(os.path.join(directory, filename), 'r+') as file:
content = file.read()
content = content.replace('"', '') # 删除所有引号
file.seek(0) # 将文件指针移回文件开头
file.write(content) # 覆盖原文件
file.truncate() # 删除原文件中未覆盖的部分
def transori(directory):
for filename in os.listdir(directory):
if filename.endswith('.csv'):
file_path = os.path.join(directory, filename)
df = pd.read_csv(file_path, skiprows=[0])
df = df[df.columns[0]].to_frame()
df.to_csv(file_path, index=False, header=False)
def transwkt(directory):
for dirpath, dirnames, filenames in os.walk(directory):
for filename in filenames:
if filename.endswith('.csv'):
base = os.path.splitext(filename)[0]
shutil.move(os.path.join(dirpath, filename), os.path.join(dirpath, base + '.wkt'))
def trans100000(directory):
files = [f for f in os.listdir(directory) if f.endswith('.wkt')]
for file in files:
with open(os.path.join(directory, file), 'r') as f:
data = f.read()
coordinates = re.findall(r'(\d+\.\d+ \d+\.\d+)', data)
mapping = {}
for coordinate in coordinates:
x, y = map(float, coordinate.split())
mapping[coordinate] = f"{x * 100000} {y * 100000}"
for original, new in mapping.items():
data = data.replace(original, new)
with open(os.path.join(directory, file), 'w') as f:
f.write(data)
def out_put(directory):
output_list = []
for i, filename in enumerate(os.listdir(directory), start=1):
if filename.endswith('.wkt'):
filepath = os.path.join(directory, filename)
output_str = f'MapBasedMovement.mapFile{i} = {filepath}'
output_list.append(output_str)
for line in output_list:
print(line)
# 使用函数
#directory 设置为你存放数据文件的目录
directory = '../data/mygeodata'
transori(directory)
process_files(directory)
transwkt(directory)
trans100000(directory)
out_put(directory)

复制代码