Thanks <span class="Apple-style-span" style="border-collapse: collapse; ">Micha </span>for his help!<div>Following his advise, I install the qgis,and the grass plugin and create new location(epsg4326),and reproj my vector map to this proj,then export it to kml successfully,however when add the kml to googlearth I found the layer doesnot match the tiff,the annex can give a clear explaination.</div>
<div>The reason I export the map to kml and add it to googlearth is to make a thecmatic map by their areas(I have five vector maps),also I found the d.vect.thecmatic is to make this kind of map ,however I did not sure the "column" parameter,I wonder this commond is useful for my goal?<br>
<div class="gmail_quote">2009/6/21 Micha Silver <span dir="ltr"><<a href="mailto:micha@arava.co.il">micha@arava.co.il</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
maven apache wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div></div><div class="h5">
I am using 6.3,my tiff info:<br>
---------------------------------------<br>
Coordinate System is:<br>
PROJCS["WGS 84 / UTM zone 20S",<br>
GEOGCS["WGS 84",<br>
DATUM["WGS_1984",<br>
SPHEROID["WGS 84",6378137,298.2572235630016,<br>
AUTHORITY["EPSG","7030"]],<br>
AUTHORITY["EPSG","6326"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4326"]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",0],<br>
PARAMETER["central_meridian",-63],<br>
PARAMETER["scale_factor",0.9996],<br>
PARAMETER["false_easting",500000],<br>
PARAMETER["false_northing",10000000],<br>
UNIT["metre",1,<br>
AUTHORITY["EPSG","9001"]],<br>
AUTHORITY["EPSG","32720"]]<br>
--------------------------------------------------------------<br>
First I create location :g.proj -c georef=d:/swilAlphaTIFF.tif datumtrans=0 location=newlocation<br>
then import it to the new location: get four raster maps(it has three band),I composite them and then transfer the composted raster to vector,then extract each category to a vector,one of them is named bule.<br>
Then I start to covert the projection.<br>
I creation new epsg4326 location: then "v.proj input=bule location=newlocation mapset=PERMANENT dbase=D:\GrassDataBase output=newblue<br>
But when I display the newblue it is blank,just like the projection is not matched......<br>
I ...... have no idea.<br>
<br>
=============================<br>
Actually I have exported the vector data to shapefile ,then use arcgis to transfer the shapefile to kml,but I think the accurcy is low...<br>
<br></div></div>
2009/6/19 Hamish <<a href="mailto:hamish_b@yahoo.com" target="_blank">hamish_b@yahoo.com</a> <mailto:<a href="mailto:hamish_b@yahoo.com" target="_blank">hamish_b@yahoo.com</a>>><div class="im"><br>
<br>
<br>
Hamish:<br>
> create a new lat/lon WGS84 location with EPSG code #<br>
> 4326, then pull in maps from source mapset using r.proj, v.proj.<br>
<br>
maven:<br>
> I have lose my heart so.... So many<br>
> people have gave me help,.......I can not do it.<br>
<br>
<br>
ok probably that I broke g.proj in 6.4.0rc5 makes it my fault that you<br>
can't create any new locations :) sorry.<br>
<br>
<br>
<br>
try this:<br>
<br>
start grass from the command line with<br>
<br>
grass64 -text<br>
<br>
then type a new location name, esc-enter, and create new Lat/Lon<br>
location.<br>
follow prompts, when it asks for a datum type "wgs84".<br>
<br>
<br>
<br>
or try to install another version besides 6.4.0rc5<br>
(not 6.5svn or 7svn either).<br>
<br>
</div></blockquote>
Maven:<br>
Don't lose heart! It seems you've already made some good progress.<br>
Based on Hamish's suggestion (older version of GRASS) maybe you should try QGIS with it's GRASS plugin.<br>
On windows it's fairly easy to install the whole package with the OSGeo4W installer. Follow carefully the instructions here:<br>
<a href="http://sites.google.com/site/eospansite/qsig-for-windows" target="_blank">http://sites.google.com/site/eospansite/qsig-for-windows</a><br>
THen start QGIS, and activate the GRASS plugin, and there's a button on the GRASS toolbar to create a new location. You select the GISDBASE (D:\GrassDataBase) then use the EPSG ID (4326) to define projection (Be sure to click the "Find" button) , and you should be up and running in the new Lat/Lon location.<br>
Now use v.proj (from the QGIS GRASS toolbox) to project your vector from UTM to the new LatLon location, and then find v.out.ogr (again in the QGIS GRASS toolbox ) to make the kml file.<br>
BTW, If you prefer, you can close QGIS right after creating the location, and then open GRASS in this new Lat/Lon location and do these last two steps from within GRASS.<br>
-- <br>
Micha<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
<br>
good luck,<br>
Hamish<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
This mail was received via Mail-SeCure System.<br>
------------------------------------------------------------------------<div class="im"><br>
<br>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
<br></div>
This mail was received via Mail-SeCure System.<br>
<br>
<br>
<br>
</blockquote>
<br>
</blockquote></div><br></div>