I've an ASCII file of three points. Each is specified by longitude,
latitude, and name. I cannot properly import these to a map so they can be
reprojected and used.
1) Tried creating a location, placesLL, using the wxPython UI but it would
not let me enter the principal latitude or central meridian. To overcome
this hurdle I copied a .prj file from another location (also in LL
coordinates) and specified that for the new location.
2) From the command line (within the placesLL location and rbs mapset) ran
"v.in.ascii input=/usr4/keypoints/sites.txt output=features format=point
columns='x double precision, y double precision, name varchar(20)' x=1 y=2"
but no primitives were seen or processed.
What might I have done incorrectly?
Rich
On Wed, 25 Nov 2009, Rich Shepard wrote:
columns='x double precision, y double precision, name varchar(20)' x=1 y=2"
Oops! I actually specified x=0 and y=1 since there is no cat column.
Rich
If it's just 3 points, why not reproject them first and them import
into Grass inside your final location?
Take a look at the proj comand and the cs2cs. It will project your lat
long points. The cs2cs will allow for datun shifts
link to proj man page:
http://trac.osgeo.org/proj/wiki/man_proj
Cheers
Daniel
On Wed, Nov 25, 2009 at 6:05 PM, Rich Shepard <rshepard@appl-ecosys.com> wrote:
On Wed, 25 Nov 2009, Rich Shepard wrote:
columns='x double precision, y double precision, name varchar(20)' x=1
y=2"
Oops! I actually specified x=0 and y=1 since there is no cat column.
Rich
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
On Thu, 26 Nov 2009, Daniel Victoria wrote:
If it's just 3 points, why not reproject them first and them import into
Grass inside your final location?
Daniel,
Good idea.
Take a look at the proj comand and the cs2cs. It will project your lat
long points. The cs2cs will allow for datun shifts
I'm not getting results because of a syntax error. Here's the script:
cs2cs +proj=latlong +datum=NAD83 +to +proj=lcc +datum=NAD83 +ellps=GRS80
+lat_1=43.0 +lat_2=45.5 +lat_0=41.75 +lon_0=-120.5 +x_0=1312336 +y_0=0
+nadgrids=WO 122d30'32.43" 45d19'19.49"
122d30'17.92" 45d18'52.45"
122d29'34.08" 45d18'47.16"
EOF
When I run the script nothing is returned; after a while I kill the
process with ^C. What have I missed here?
Thanks,
Rich
You forgot the <<EOF ??
Not a shell guru myself but I believe the <<EOF will tell the command
to recieve input from stdin until it reaches an EOF. Then it will
process everything. But you need to put the <<EOF after the command
and before you start typing the input.
I believe there is also a way for the cs2cs command to read input from a file...
Cheers
Daniel
On Fri, Nov 27, 2009 at 3:01 PM, Rich Shepard <rshepard@appl-ecosys.com> wrote:
On Thu, 26 Nov 2009, Daniel Victoria wrote:
If it's just 3 points, why not reproject them first and them import into
Grass inside your final location?
Daniel,
Good idea.
Take a look at the proj comand and the cs2cs. It will project your lat
long points. The cs2cs will allow for datun shifts
I'm not getting results because of a syntax error. Here's the script:
cs2cs +proj=latlong +datum=NAD83 +to +proj=lcc +datum=NAD83 +ellps=GRS80
+lat_1=43.0 +lat_2=45.5 +lat_0=41.75 +lon_0=-120.5 +x_0=1312336 +y_0=0
+nadgrids=WO 122d30'32.43" 45d19'19.49"
122d30'17.92" 45d18'52.45"
122d29'34.08" 45d18'47.16"
EOF
When I run the script nothing is returned; after a while I kill the
process with ^C. What have I missed here?
Thanks,
Rich
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
On Fri, 27 Nov 2009, Daniel Victoria wrote:
You forgot the <<EOF ??
Nope.
Rich
On Fri, 27 Nov 2009, Rich Shepard wrote:
Here's the script:
Lost something when cutting and pasting:
cs2cs +proj=latlong +datum=NAD83 +to +proj=lcc +datum=NAD83 +ellps=GRS80
+lat_1=43.0 +lat_2=45.5 +lat_0=41.75 +lon_0=-120.5 +x_0=1312336 +y_0=0
+nadgrids=WO <<EOF
122d30'32.43" 45d19'19.49"
122d30'17.92" 45d18'52.45"
122d29'34.08" 45d18'47.16"
EOF
Rich
On Fri, 27 Nov 2009, Rich Shepard wrote:
I'm not getting results because of a syntax error. Here's the script:
cs2cs +proj=latlong +datum=NAD83 +to +proj=lcc +datum=NAD83 +ellps=GRS80
+lat_1=43.0 +lat_2=45.5 +lat_0=41.75 +lon_0=-120.5 +x_0=1312336 +y_0=0
+nadgrids=WO 122d30'32.43" 45d19'19.49"
122d30'17.92" 45d18'52.45"
122d29'34.08" 45d18'47.16"
EOF
Got it. Have to enter it on the command line rather than running it as a
script.
Thanks, Daniel,
Rich