def reverse_proj_from_dll(outdir, condor_path, landscape, nx, ny, flag_show):
navicon_dll = init_navicon_dll(condor_path, landscape)
max_x, max_y = navicon_dll.GetMaxX(), navicon_dll.GetMaxY()
P = {}
P[0] = (0, 0)
P[1] = (max_x, 0)
P[2] = (0, max_y)
P[3] = (max_x, max_y)
d_df = {}
s = "ref"
d_df[s] = pd.DataFrame(index=np.arange(4),
columns=["PosX", "PosY", "Lat", "Lon"])
for i in d_df[s].index:
x, y = P[i][0], P[i][1]
d_df[s]["PosX"][i] = x
d_df[s]["PosY"][i] = y
d_df[s]["Lat"][i] = navicon_dll.XYToLat(x, y)
d_df[s]["Lon"][i] = navicon_dll.XYToLon(x, y)
a_x = np.linspace(0, max_x, nx)
a_y = np.linspace(0, max_y, ny)
values = cartesian([a_x, a_y])
s = "measures"
d_df[s] = pd.DataFrame(values, columns=["PosX", "PosY"])
d_df[s]["Lat"] = np.nan
d_df[s]["Lon"] = np.nan
for i in d_df[s].index:
x = d_df[s]["PosX"][i]
y = d_df[s]["PosY"][i]
d_df[s]["Lat"][i] = navicon_dll.XYToLat(x, y)
d_df[s]["Lon"][i] = navicon_dll.XYToLon(x, y)
#print(d_df["ref"])
#print(d_df["measures"])
print("""ref:
%s
measures:
%s""" % (d_df["ref"], d_df["measures"]))
filename_out = os.path.join(outdir, "geodesic_%s.xlsx" % landscape)
print("Output '%s'" % filename_out)
with pd.ExcelWriter(filename_out) as writer:
key = "ref"
d_df[key].to_excel(writer, sheet_name=key)
key = "measures"
d_df[key].to_excel(writer, sheet_name=key)
plot_geodesic(outdir, landscape, d_df["measures"])
if flag_show:
plt.show()
def reverse_proj_from_dll(outdir, condor_path, landscape, nx, ny, flag_show):
navicon_dll = init_navicon_dll(condor_path, landscape)
max_x, max_y = navicon_dll.GetMaxX(), navicon_dll.GetMaxY()
P = {}
P[0] = (0, 0)
P[1] = (max_x, 0)
P[2] = (0, max_y)
P[3] = (max_x, max_y)
d_df = {}
s = "ref"
d_df[s] = pd.DataFrame(index=np.arange(4), columns=["PosX", "PosY", "Lat", "Lon"])
for i in d_df[s].index:
x,y = P[i][0], P[i][1]
d_df[s]["PosX"][i] = x
d_df[s]["PosY"][i] = y
d_df[s]["Lat"][i] = navicon_dll.XYToLat(x, y)
d_df[s]["Lon"][i] = navicon_dll.XYToLon(x, y)
a_x = np.linspace(0, max_x, nx)
a_y = np.linspace(0, max_y, ny)
values = cartesian([a_x, a_y])
s = "measures"
d_df[s] = pd.DataFrame(values, columns=["PosX", "PosY"])
d_df[s]["Lat"] = np.nan
d_df[s]["Lon"] = np.nan
for i in d_df[s].index:
x = d_df[s]["PosX"][i]
y = d_df[s]["PosY"][i]
d_df[s]["Lat"][i] = navicon_dll.XYToLat(x, y)
d_df[s]["Lon"][i] = navicon_dll.XYToLon(x, y)
#print(d_df["ref"])
#print(d_df["measures"])
print("""ref:
%s
measures:
%s""" % (d_df["ref"], d_df["measures"]))
filename_out = os.path.join(outdir, "geodesic_%s.xlsx" % landscape)
print("Output '%s'" % filename_out)
with pd.ExcelWriter(filename_out) as writer:
key = "ref"
d_df[key].to_excel(writer, sheet_name=key)
key = "measures"
d_df[key].to_excel(writer, sheet_name=key)
plot_geodesic(outdir, landscape, d_df["measures"])
if flag_show:
plt.show()
def main(xls_filename, outdir):
basepath = os.path.dirname(__file__)
#basepath = os.path.dirname(os.path.abspath(__file__))
if outdir=='':
outdir = os.path.join(basepath, 'out')
#xls_filename = os.path.join(outdir, "Provence-Oisans2.xlsx")
filename_base, filename_ext = os.path.splitext(os.path.basename(xls_filename))
d_df = {}
d_df = pd.read_excel(xls_filename, sheetname=None)
max_x, max_y = d_df['ref']['PosX'].max(), d_df['ref']['PosY'].max()
print "max_x=%s, max_y=%s" % (max_x, max_y)
p = pyproj.Proj(
proj='utm',
zone=32,
ellps='WGS84'
)
d_df['measures']['PosX2'] = np.nan
d_df['measures']['PosY2'] = np.nan
for i, point in d_df['measures'].iterrows():
xy2 = p(point['Lat'], point['Lon'])
d_df['measures']['PosX2'][i] = xy2[0]
d_df['measures']['PosY2'][i] = xy2[1]
#print(xy2)
d_df['measures']['Eps'] = np.sqrt(
(d_df['measures']['PosX2'] - d_df['measures']['PosX'])**2 + \
(d_df['measures']['PosY2'] - d_df['measures']['PosY'])**2
)
print(d_df)
print(d_df['measures']['Eps'].mean())
plot_geodesic(outdir, filename_base, d_df['measures'])
plt.show()
def main(xls_filename, outdir):
basepath = os.path.dirname(__file__)
#basepath = os.path.dirname(os.path.abspath(__file__))
if outdir == '':
outdir = os.path.join(basepath, 'out')
#xls_filename = os.path.join(outdir, "Provence-Oisans2.xlsx")
filename_base, filename_ext = os.path.splitext(
os.path.basename(xls_filename))
d_df = {}
d_df = pd.read_excel(xls_filename, sheetname=None)
max_x, max_y = d_df['ref']['PosX'].max(), d_df['ref']['PosY'].max()
print "max_x=%s, max_y=%s" % (max_x, max_y)
p = pyproj.Proj(proj='utm', zone=32, ellps='WGS84')
d_df['measures']['PosX2'] = np.nan
d_df['measures']['PosY2'] = np.nan
for i, point in d_df['measures'].iterrows():
xy2 = p(point['Lat'], point['Lon'])
d_df['measures']['PosX2'][i] = xy2[0]
d_df['measures']['PosY2'][i] = xy2[1]
#print(xy2)
d_df['measures']['Eps'] = np.sqrt(
(d_df['measures']['PosX2'] - d_df['measures']['PosX'])**2 + \
(d_df['measures']['PosY2'] - d_df['measures']['PosY'])**2
)
print(d_df)
print(d_df['measures']['Eps'].mean())
plot_geodesic(outdir, filename_base, d_df['measures'])
plt.show()
