I must have a syntax error in either the ASCII standard input file or the
command line and I don't see my error.
The command and the results:
v.in.ascii --overwrite input=project_bounds.txt output=projbnd \
format=standard separator=space
WARNING: Vector map <projbnd> already exists and will be overwritten
Building topology for vector map <projbnd@topography>...
Registering primitives...
One primitive registered
4 vertices registered
Building areas...
100%
0 areas built
0 isles built
Attaching islands...
Attaching centroids...
100%
Number of nodes: 2
Number of primitives: 1
Number of points: 0
Number of lines: 0
Number of boundaries: 1
Number of centroids: 0
Number of areas: 0
Number of isles: 0
WARNING: Number of incorrect boundaries: 1
The input file:
ORGANIZATION: Applied Ecosystem Services, Inc.
DIGIT DATE: 2016-09-08
DIGIT NAME: -
MAP NAME: project_area
MAP DATE: 2016
MAP SCALE: -
OTHER INFO: -
ZONE: 0
VERTI:
B 4
-123.102155 45.215480
-123.098186 45.215480
-123.098186 45.211372
-123.102155 45.211372
The projection of the location is:
NAD 1983 HARN State Plane, Oregon North; FIPS 3601 Feet Intl.
Geographic Coordinate System: North American 1983 HARN
Datum; North American 1983 HARN
Spheroid: GRS 1980, 6378137.0,298.257222101
Prime Meridiam: Greenwich, 0.0, unit Degree, 0.0174532925199433
Projection: Lambert Conformal Conic,
Parameter: False Easting, 8202099.737532808
Parameter: False Northing, 0.0]
Parameter: Central Meridian, -120.5
Parameter: Standard Parallel 1, 44.33333333333334
Parameter: Standard Parallel 2, 46.0
Parameter: Latitude Of Origin 43.66666666666666
Unit: Foot, 0.3048m
What do I need to change?
TIA,
Rich
Hi Rich,
On Thu, Sep 8, 2016 at 6:37 PM, Rich Shepard <rshepard@appl-ecosys.com>
wrote:
VERTI:
B 4
-123.102155 45.215480
-123.098186 45.215480
-123.098186 45.211372
-123.102155 45.211372
What do I need to change?
I think you need to add centroid [1]:
B 4
5959010.9323622 3401338.36037757
5959096.7459483 3401370.54047235
5959091.38259917 3401450.99070932
5959010.9323622 3401338.36037757
C 1 1
5959063.08352122 3401386.98533277
Best,
Vaclav
[1]
https://grass.osgeo.org/grass70/manuals/v.in.ascii.html#example-1a)—standard-format-mode
On Thu, 8 Sep 2016, Vaclav Petras wrote:
I think you need to add centroid [1]:
B 4
5959010.9323622 3401338.36037757
5959096.7459483 3401370.54047235
5959091.38259917 3401450.99070932
5959010.9323622 3401338.36037757
C 1 1
5959063.08352122 3401386.98533277
https://grass.osgeo.org/grass70/manuals/v.in.ascii.html#example-1a)—standard-format-mode
Vaclav,
I looked at the example but did not realize the centroid was required
input.
Are the above points projected by v.in.ascii when reading the lon/lat data
file?
Thanks,
Rich
On Thu, 8 Sep 2016, Vaclav Petras wrote:
I think you need to add centroid [1]:
B 4
5959010.9323622 3401338.36037757
5959096.7459483 3401370.54047235
5959091.38259917 3401450.99070932
5959010.9323622 3401338.36037757
C 1 1
5959063.08352122 3401386.98533277
https://grass.osgeo.org/grass70/manuals/v.in.ascii.html#example-1a)—standard-format-mode
Vaclav,
Following the link on that page to
<https://grass.osgeo.org/grass70/manuals/vectorascii.html>
provides details on each portion of the input file. I modified my input file
by adding a centroid (whose coordinates are approximately in the center):
ORGANIZATION: Applied Ecosystem Services, Inc.
DIGIT DATE: 2016-09-08
DIGIT NAME: -
MAP NAME: projbnd
MAP DATE: 2016
MAP SCALE: -
OTHER INFO: -
ZONE: 0
VERTI:
B 4
-123.102155 45.215480
-123.098186 45.215480
-123.098186 45.211372
-123.102155 45.211372
C 1 1
-123.100170 45.213426
and grass is still not satisfied:
v.in.ascii --overwrite input=project_bounds.txt output=projbnd \
format=standard separator=space
WARNING: Vector map <projbnd> already exists and will be overwritten
Building topology for vector map <projbnd@topography>...
Registering primitives...
2 primitives registered
5 vertices registered
Building areas...
100%
0 areas built
0 isles built
Attaching islands...
Attaching centroids...
100%
Number of nodes: 2
Number of primitives: 2
Number of points: 0
Number of lines: 0
Number of boundaries: 1
Number of centroids: 1
Number of areas: 0
Number of isles: 0
WARNING: Number of centroids exceeds number of areas: 1 > 0
WARNING: Number of incorrect boundaries: 1
WARNING: Number of centroids outside area: 1
GRASS 7.3.svn
(Oregon-Dayton):~/projects/oregon/washington-eichler/data/topography >
Still no areas are built.
Rich
1 5
On Thu, 8 Sep 2016, Rich Shepard wrote:
ORGANIZATION: Applied Ecosystem Services, Inc.
DIGIT DATE: 2016-09-08
DIGIT NAME: -
MAP NAME: projbnd
MAP DATE: 2016
MAP SCALE: -
OTHER INFO: -
ZONE: 0
VERTI:
B 4
-123.102155 45.215480
-123.098186 45.215480
-123.098186 45.211372
-123.102155 45.211372
C 1 1
-123.100170 45.213426
1 5
Uh, oh! That last line got out of place.
Rich
On Thu, Sep 8, 2016 at 7:47 PM, Rich Shepard <rshepard@appl-ecosys.com> wrote:
VERTI:
B 4
-123.102155 45.215480
-123.098186 45.215480
-123.098186 45.211372
-123.102155 45.211372
C 1 1
-123.100170 45.213426
and grass is still not satisfied:
v.in.ascii --overwrite input=project_bounds.txt output=projbnd
format=standard separator=space
…
Number of isles: 0
WARNING: Number of centroids exceeds number of areas: 1 > 0
The example in the specification [1] shows that the last vertex must be repeated:
B 4
5959010.9323622 3401338.36037757
5959096.7459483 3401370.54047235
5959091.38259917 3401450.99070932
5959010.9323622 3401338.36037757
[1] https://grass.osgeo.org/grass70/manuals/vectorascii.html
On Thu, 8 Sep 2016, Vaclav Petras wrote:
The example in the specification [1] shows that the last vertex must be
repeated:
B 4
5959010.9323622 3401338.36037757
5959096.7459483 3401370.54047235
5959091.38259917 3401450.99070932
5959010.9323622 3401338.36037757
Thanks, Vaclav! I totally missed that and it makes sense that the
beginning and ending points must be both the same and provided. Otherwide
there's a three-sided box rather than a closed one.
It's been several years since I've used grass in a project and it's taking
a while to recall the details.
Much appreciated,
Rich