[GRASS-user] Importing Converted Points File

[Markus Neteler]

On Sat, Nov 28, 2009 at 12:45 AM, Rich Shepard <rshepard at appl-ecosys.com> wrote:
>  I converted a points file from lon/lat geographic coordinates to lcc using
> cs2cs. The command line used was:
>
> 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

You may also want to try m.proj (GRASS wrapper for cs2cs):

http://grass.osgeo.org/grass64/manuals/html64_user/m.proj.html


> and the output is:
>
> -5056134.11|5898037.45
> -5057005.89|5898228.79
> -5057306.14|5899149.81
>
> after I stripped off the trailing 0.0 (which I presume is a z value).
>
>  Now, when I try to use these values in v.in.ascii I'm told that column
> numbers must not be negative:
>
> GRASS 6.4.0svn (sitesLCC):/usr4/grassbase/sitesLCC > v.in.ascii
> input=keypoints/sites_lcc.txt output=sites format=point x=0 y=1
> ERROR: Column numbers must not be negative
>
>  Do I just change the longitude values to positive?

No - it says *number*: You specified x as column 0 but it is the
first. So x=1 y=2 will do it.

Markus