Dear Friends,
I am still using 4.1 and maybe I have missed it in g.manual. Can anyone tell
me an easy way to convert a site_lists file from x,y projection to utm
coordinates?
Thank you for your help,
Francisco
Dr. Francisco Estrada Belli
Research Fellow
Boston University
Department of Archaeology
675 Commonwealth Ave.
Boston, MA 02215
phone: (617) 353-3415
fax: (617) 353-6800
fax: 617-353-6800
The Santa Rosa Archaeological Project 1996 web page is at:
http://www.bu.edu/ARCHAEOLOGY/www/grads/franc/Sap.htm
Francisco,
If you have at least 4 points that you have know UTM coordinates for,
you can use v.transform, however this will only work on an ascii file.
If your site_lists file is not in vector format yet, do s.to.vect.
Convert the file to an ascii coverage using v.out.ascii. Make a file
called "points" or something in the format:
old-x1 old-y1 new-x1 new-y1
old-x2 old-y2 new-x2 new-y2
old-x3 old-y3 new-x3 new-y3
...and so on...you need at least 4 points. Run the v.transform script to
generate a new ascii coverage. (It's easy to run interactively.) Then,
finally, use v.in.ascii to convert back to a binary vector file. I just
did this two days ago to convert a geodimeter file to UTM and it is
quite simple! I only used 4 points and got an RMS of 2.9.
Good Luck!
--
Katie Alvin
FOCS/FGI/GIS Specialist
USDA NRCS, Bozeman, MT
(406)587-6929
kalvin@mt.nrcs.usda.gov