#1864: r.horizon generates incorrect results in single point mode
------------------------+-------------------------
Reporter: rmyorston | Owner: grass-dev@…
Type: defect | Status: new
Priority: normal | Milestone: 6.4.3
Component: Raster | Version: 6.4.2
Resolution: | Keywords: r.horizon
CPU: x86-32 | Platform: Linux
------------------------+-------------------------
Comment (by annakrat):
A test case for UTM32 location, on artificial surface:
{{{
g.region n=5118356 s=5088351 e=666798 w=636793 res=10
r.mapcalc --overwrite expression="test_surface = 100*sin( 0.2 *row()) +
100*sin( 0.2 *col())"
r.horizon elevation=test_surface direction=0 step=1
coordinates=656096,5102603 -d fi=horizon.csv maxd=800 --o
}}}
{{{
import numpy as np
import matplotlib.pyplot as plt
my_data = np.genfromtxt('horizon.csv', delimiter=',')
my_data = my_data[1:, :]
ax = plt.subplot(111, polar=True)
bars = ax.plot(my_data[:, 0] / 180 * np.pi, ( my_data[:, 1]) / 180 *
np.pi)
plt.show()
}}}
--
Ticket URL: <https://trac.osgeo.org/grass/ticket/1864#comment:2>
GRASS GIS <http://grass.osgeo.org>
_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev
