<div>All,<br></div><div><br></div><div>I am currently trying to calculate east (direction=0) and west (direction=180) horizons for a single location (lon,lat=-105.007,40.05).  I am using a 30m DEM (USGS NED) in geographic coordinates:</div>
<div><br></div><div>projection: 3 (Latitude-Longitude)<br>zone:       0<br>datum:      nad83<br>ellipsoid:  grs80<br>north:      42:00:06N<br>south:      38:59:54N<br>west:       107:00:06W<br>east:       103:59:54W<br>nsres:      0:00:01<br>
ewres:      0:00:01<br>rows:       10812<br>cols:       10812<br>cells:      116899344<br></div><div><br></div><div>When I run r.horizon in single point mode and degree output:</div><div>r.horizon -d elevin=elev horizonstep=180 coord=-105.007,40.05 dist=1.0 </div>
<div><br></div><div>I get the following result:</div>
<div>0.000000, 0.574081<br>180.000000, 0.645220</div><div><br></div><div>However, when I run r.horizon in raster map mode for the w horizon and extract the corresponding grid cell, I get a different result:</div><div><br>
</div><div>r.horizon -d elevin=elev direction=180 horizon=horizonW</div><div>r.what input=horizonW_0 east_north=-105.007,40.05</div><div><br></div><div>produces:</div><div>-105.007|40.05|2.247817</div><div><br></div><div>
Is there any reason why single point mode and raster map mode would produce different horizon results for the same point (0.645220 vs. 2.247817)?</div><div><br></div><div>Thanks,</div><div><br></div><div>Jared Oyler<br>Numerical Terradynamic Simulation Group (NTSG)<br>
College of Forestry &amp; Conservation<br>The University of Montana<br>Email: <a href="mailto:jared.oyler@ntsg.umt.edu">jared.oyler@ntsg.umt.edu</a><br></div><div><br></div><div><br></div><div><br></div><div><br><br></div>
<div><br></div>