<HTML>
<pre wrap=""><BR>
</pre><div class="moz-text-plain" wrap="true" graphical-quote="true" style="font-family: -moz-fixed; font-size: 12px;" lang="x-western"><blockquote type="cite"><pre wrap=""><span class="moz-txt-citetags">></span><BR>
<span class="moz-txt-citetags">> </span>Message: 4<BR>
<span class="moz-txt-citetags">> </span>Date: Sat, 3 Jan 2009 13:25:31 -0800 (PST)<BR>
<span class="moz-txt-citetags">> </span>From: Hamish <a class="moz-txt-link-rfc2396E" href="mailto:hamish_b@yahoo.com"><hamish_b@yahoo.com></a><BR>
<span class="moz-txt-citetags">> </span>Subject: Re: [GRASS-user] Re: GRASS-user] Help with reprojection<BR>
<span class="moz-txt-citetags">> </span>To: <a class="moz-txt-link-abbreviated" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><BR>
<span class="moz-txt-citetags">> </span>Message-ID: <a class="moz-txt-link-rfc2396E" href="mailto:62626.56685.qm@web110014.mail.gq1.yahoo.com"><62626.56685.qm@web110014.mail.gq1.yahoo.com></a><BR>
<span class="moz-txt-citetags">> </span>Content-Type: text/plain; charset=us-ascii<BR>
<span class="moz-txt-citetags">></span><BR>
<span class="moz-txt-citetags">> </span>Hamish wrote:<BR>
<span class="moz-txt-citetags">> </span> <BR>
</pre><blockquote type="cite"><blockquote type="cite"><pre wrap=""><span class="moz-txt-citetags">>> > </span>- create a simple xy location<BR>
<span class="moz-txt-citetags">>> > </span> (lat/lon location will not allow north > 90, and your image while<BR>
<span class="moz-txt-citetags">>> > </span> still not geo-referenced will go to 8100)<BR>
<span class="moz-txt-citetags">>> > </span>- run r.region to set n,s,e,w bounds to 90,-90,180,-180<BR>
<span class="moz-txt-citetags">>> > </span>- run g.setproj to rejig the location into a lat/lon one.<BR>
</pre></blockquote><pre wrap=""><span class="moz-txt-citetags">>> </span> <BR>
</pre></blockquote><pre wrap=""><span class="moz-txt-citetags">></span><BR>
<span class="moz-txt-citetags">> </span>umm, that might not work -- r.info will still know the map is XY even<BR>
<span class="moz-txt-citetags">> </span>if the location is changed to lat/lon.<BR>
<span class="moz-txt-citetags">></span><BR>
<span class="moz-txt-citetags">> </span> <BR>
</pre><blockquote type="cite"><blockquote type="cite"><pre wrap=""><span class="moz-txt-citetags">>> > </span>- check resolution is correct (nicely "0:01:20") with r.info.<BR>
<span class="moz-txt-citetags">>> > </span>- zoom to area of interest. you probably do not want to reproject<BR>
<span class="moz-txt-citetags">>> > </span> entire planet. after zooming check resolution is preserved,<BR>
<span class="moz-txt-citetags">>> > </span> (g.region -p, maybe with "g.region res=0:01:20 -a" after zoom to fix)<BR>
<span class="moz-txt-citetags">>> > </span> and run "r.mapcalc cropmap=fullmap" to perform the crop.<BR>
<span class="moz-txt-citetags">>> > </span><BR>
<span class="moz-txt-citetags">>> > </span>then from the lambert location run r.proj to pull the<BR>
<span class="moz-txt-citetags">>> > </span>cropped image across.<BR>
<span class="moz-txt-citetags">>> > </span><BR>
<span class="moz-txt-citetags">>> > </span>there are some examples of this process in the GRASS wiki,<BR>
<span class="moz-txt-citetags">>> > </span>look at the "Global datasets" page.<BR>
</pre></blockquote><pre wrap=""><span class="moz-txt-citetags">>> </span> <BR>
</pre></blockquote><pre wrap=""><span class="moz-txt-citetags">></span><BR>
<span class="moz-txt-citetags">> </span>the above issue should be covered there; also check the mailing list<BR>
<span class="moz-txt-citetags">> </span>archives.<BR>
<span class="moz-txt-citetags">></span><BR>
<span class="moz-txt-citetags">> </span>so your "easiest" solution is to create a "world file". see the GDAL<BR>
<span class="moz-txt-citetags">> </span>JPG or GeoTiff format import page, or do a web search for instructions.<BR>
<span class="moz-txt-citetags">></span><BR>
<span class="moz-txt-citetags">></span><BR>
<span class="moz-txt-citetags">> </span>Hamish<BR>
<span class="moz-txt-citetags">> </span> <BR>
</pre></blockquote><pre wrap=""><!---->Hamish,<BR>
<BR>
Grass-6.3.0 under Mac doesn't seem to be checking for the world file.<BR>
<BR>
My test:<BR>
<BR>
1) Export a raster using r.out.tiff with the "create world file" box checked<BR>
2) Check the output folder:<BR>
test.tiff<BR>
test.tfw<BR>
<BR>
3) Import using r.in.gdal:<BR>
r.in.gdal input=test.tiff output=test_2<BR>
<BR>
In mine, this produces a projection mismatch error - which seems to tell me that Grass isn't noticing the tfw file.<BR>
<BR>
Richard<BR>
<BR>
</pre></div> </HTML>
<BR>