Re: [GRASS-dev] [GRASS GIS] #957: v.voronoi has extra lines in output
#957: v.voronoi has extra lines in output ---+ Reporter: helena | Owner: grass-dev@… Type: defect | Status: new Priority: major | Milestone: 6.4.4 Component: Vector | Version: svn-releasebranch64 Keywords: v.voronoi |Platform: All Cpu: All| ---+ Comment(by mmetz): Replying to [comment:26 hamish]: Hi, are these empty functions still useful for debug or future work? https://trac.osgeo.org/grass/browser/grass/trunk/vector/v.voronoi/sw_output.c#L7 These functions can probably be removed. I guess they come from the original code of Steve J. Fortune (1987). -- Ticket URL: https://trac.osgeo.org/grass/ticket/957#comment:27 GRASS GIS http://grass.osgeo.org ___ grass-dev mailing list grass-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-dev
Re: [GRASS-dev] How to calculate mean coordinates from big point datasets?
Markus Metz wrote: Please try the new addon v.centerpoint [0]. It calculates various center points for point clouds, lines, and areas. Standard options are the geometric mean (center of gravity) That's the arithmetic mean. The geometric mean of a set of N values is the Nth root of the product of the values (i.e. the logarithm of the geometric mean is the arithmetic mean of the logarithms of the values). It wouldn't be meaningful to compute the geometric mean of coordinate data. -- Glynn Clements gl...@gclements.plus.com ___ grass-dev mailing list grass-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-dev
Re: [GRASS-dev] [GRASS GIS] #957: v.voronoi has extra lines in output
#957: v.voronoi has extra lines in output
---+
Reporter: helena | Owner: grass-dev@…
Type: defect | Status: new
Priority: major | Milestone: 6.4.4
Component: Vector | Version: svn-releasebranch64
Keywords: v.voronoi |Platform: All
Cpu: All|
---+
Comment(by mlennert):
In 64release, I'm still getting a subset of the holes reported by Helena
in a raster resulting from v.voronoi polygons. Both
{{{
v.voronoi input=elev_lid792_randpts@PERMANENT output=voro
v.to.rast input=voro@user1 layer=1 type=point,line,area output=voro
use=attr column=value value=1 rows=4096
}}}
and
{{{
v.voronoi input=elev_lid792_randpts@PERMANENT output=voro
v.clean in=voro out=voro_snap tool=snap thresh=0.001 type=boundary
v.to.rast input=voro_snap@user1 layer=1 type=point,line,area
output=voro_snap use=attr column=value value=1 rows=4096
}}}
result in the same holes, although the vector polygons all look ok.
Moritz
--
Ticket URL: https://trac.osgeo.org/grass/ticket/957#comment:28
GRASS GIS http://grass.osgeo.org
___
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev
Re: [GRASS-dev] [GRASS GIS] #957: v.voronoi has extra lines in output
#957: v.voronoi has extra lines in output
---+
Reporter: helena | Owner: grass-dev@…
Type: defect | Status: new
Priority: major | Milestone: 6.4.4
Component: Vector | Version: svn-releasebranch64
Keywords: v.voronoi |Platform: All
Cpu: All|
---+
Comment(by mlennert):
Replying to [comment:28 mlennert]:
In 64release, I'm still getting a subset of the holes reported by Helena
in a raster resulting from v.voronoi polygons.
Actually I also still see that issue in grass7. The region used in both
cases was:
{{{
g.region -p
projection: 99 (Lambert Conformal Conic)
zone: 0
datum: nad83
ellipsoid: a=6378137 es=0.006694380022900787
north: 220750
south: 22
west: 638300
east: 639000
nsres: 1
ewres: 1
rows: 750
cols: 700
cells: 525000
}}}
--
Ticket URL: https://trac.osgeo.org/grass/ticket/957#comment:29
GRASS GIS http://grass.osgeo.org
___
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev
Re: [GRASS-dev] [GRASS GIS] #2014: r.sun using EPSG:3031 projection gives strange results
#2014: r.sun using EPSG:3031 projection gives strange results
---+
Reporter: pierreroudier | Owner: grass-dev@…
Type: defect | Status: new
Priority: normal | Milestone: 7.0.0
Component: Raster | Version: svn-trunk
Keywords: r.sun |Platform: Linux
Cpu: x86-64 |
---+
Comment(by hamish):
Hi, fwiw after fixing a small bug in r.sunmask I was able to write a
script to show the sun's position through the day on a polar graph.
Plot is for Jan 1st 00:00-23:45 (gap in the plot is 23:45-24:00) around
Ross Island in Antarctica (McMurdo Scott Bases) in a polar stereographic
location (epsg:3286). Red + are the 15 minute marks. So the more sun in
the south effect seems to be a combination of midnight shade on the north
slope, sun peeking over the ridges onto the south slope at midday, and the
cumulative time the Sun spends in the south at this time of year- the
azimuth is not strictly going around at 15deg/hour.
d.histogram of a r.slope.aspect slope map shows most of Erebus between 4
and 17 deg, very little over 25 deg, and a max of 35deg. So midday sun at
35.5 deg above the horizon would generally get to the south slope, even if
the shallow angle meant that little of it would be hitting it with much
impact.
so it's all about the cumulative timing the tortoise beats the hare.
Hamish
-
The r.sunmask module can calculate the position of the Sun in the sky
throughout the day:
{{{
for HOUR in `seq -f'%02.0f' 0 23` ; do
for MINUTE in 00 15 30 45 ; do
r.sunmask -sg elev=ramp200dem_wgs_v2.ross_island out=dummy \
year=2013 month=1 day=1 hour=$HOUR minute=$MINUTE second=0 \
timezone=12 \
solar_pos.hour_$HOUR.$MINUTE.txt
done
done
grep sun solar_pos.hour_* | grep -v set | grep -v rise solar_pos.day
grep above solar_pos.day | cut -f2 -d= angleabove.txt
grep azim solar_pos.day | cut -f2 -d= azimuth.txt
% matlab or octave
load azimuth.txt
load angleabove.txt
polar(azimuth * pi/180, angleabove)
hold on
polar(azimuth * pi/180, angleabove, 'r+')
set(gca, 'View', [90 -90]) % rotate so 0 deg as north-up not +x axis
}}}
--
Ticket URL: https://trac.osgeo.org/grass/ticket/2014#comment:9
GRASS GIS http://grass.osgeo.org
___
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev
Re: [GRASS-dev] [GRASS GIS] #2014: r.sun using EPSG:3031 projection gives strange results
#2014: r.sun using EPSG:3031 projection gives strange results ---+ Reporter: pierreroudier | Owner: grass-dev@… Type: defect | Status: new Priority: normal | Milestone: 7.0.0 Component: Raster | Version: svn-trunk Keywords: r.sun |Platform: Linux Cpu: x86-64 | ---+ Comment(by hamish): (radius of the plot line is angle of the Sun above the horizon, theta is compass angle, north is 0 deg, east is 90, ..) -- Ticket URL: http://trac.osgeo.org/grass/ticket/2014#comment:10 GRASS GIS http://grass.osgeo.org ___ grass-dev mailing list grass-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-dev
