[GRASS-dev] [GRASS GIS] #1864: r.horizon generates incorrect results in single point mode

[GRASS GIS]

#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>