I just tried reinstalling with newest GRASS 6.4 binary snapshot (grass-6.4.svn-i686-pc-linux-gnu-11_10_2008.tar.gz), and trying r.proj with a different map. I am trying to reproject from a WGS84, UTM Z36N location (that I know is correct) to a wgs84 Longlat location I created "interactively" in the shell with the "projection values" option at GRASS start up. The first image linked below is the original map. You can see that it is of a single watershed with all values outside the watershed as NULL, and is the flowacc map produced by r.watershed. The second map shows the result of running r.proj to reproject the map from the UTM location into the latlong location. You can see that now not only is there a small one pixel border at what was the extents of the original map layer, bu there is a wierd "rainbow" of valued pixels forming around the edge of the watershed, and two "rainbow" bars shooting off towards the nw and se corners. FYI, I get the same effect whether I choose "nearest", "bilinear", or "bicubic" for the interpolation procedure. IMO, something is fundamentally broken in r.proj but I can't think of what could be causing such wierd errors other than an interpolation problem...<br>
<br>Original Map:<br><a href="http://www.public.asu.edu/~iullah/files/MEDLAND/original_map.png">http://www.public.asu.edu/~iullah/files/MEDLAND/original_map.png</a><br><br>Reprojected_map:<br><a href="http://www.public.asu.edu/~iullah/files/MEDLAND/more_wierd_rproj_effects.png">http://www.public.asu.edu/~iullah/files/MEDLAND/more_wierd_rproj_effects.png</a><br>
<br><br><br><br><div class="gmail_quote">On Tue, Nov 18, 2008 at 1:35 PM, Isaac Ullah <span dir="ltr"><<a href="mailto:isaac.ullah@asu.edu">isaac.ullah@asu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Sorry to have been unclear. It seems that the new "boundary" pixels are created at the what seems to be edge of the input map's regional extents. To clarify, there is no "boundary of pixels" (eg. "frame") in the original raster stream seg maps produced by r.watershed, but r.proj creates a one to two pixel wide "frame" that outlines the reprojected maps. to repeat, there is NO boundary in the origianl maps, but there IS a boundary after reprojection. Why is this happening?<div>
<div></div><div class="Wj3C7c"><br>
<br><div class="gmail_quote">On Tue, Nov 18, 2008 at 11:59 AM, Glynn Clements <span dir="ltr"><<a href="mailto:glynn@gclements.plus.com" target="_blank">glynn@gclements.plus.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div><br>
Isaac Ullah wrote:<br>
<br>
> Hi all, I'm attching to png's of a wierd phenomenon that just occured to me<br>
> when reprojecting a couple of stream segment maps (produced by r.watershed)<br>
> from a UTM location ot a longlat location (so I can export to kml for use in<br>
> GE). The maps reprojected okay, but there is a 1-2 pixel wide "border"<br>
> around the reprojected maps that is approximately the same size as the<br>
> original map's boundary.<br>
<br>
</div>What do you mean by "the original map's boundary"?<br>
<br>
If the source map has a border embedded in its data, it will get<br>
re-projected along with the rest of the map.<br>
<font color="#888888"><br>
--<br>
Glynn Clements <<a href="mailto:glynn@gclements.plus.com" target="_blank">glynn@gclements.plus.com</a>><br>
</font></blockquote></div><br><br clear="all"><br></div></div><div><div></div><div class="Wj3C7c">-- <br><br>Isaac I Ullah, M.A.<br><br>Archaeology PhD Student,<br>ASU School of Evolution and Social Change<br><br>Research Assistant,<br>
Mediterranean Landscape Dynamics Project<br>
***************************************************<br><a href="mailto:isaac.ullah@asu.edu" target="_blank">isaac.ullah@asu.edu</a><br><a href="mailto:ullah@archaeologist.com" target="_blank">ullah@archaeologist.com</a><br>
<br><a href="http://www.public.asu.edu/%7Eiullah" target="_blank">http://www.public.asu.edu/~iullah</a><br>
***************************************************<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><br>Isaac I Ullah, M.A.<br><br>Archaeology PhD Student,<br>ASU School of Evolution and Social Change<br><br>Research Assistant,<br>Mediterranean Landscape Dynamics Project<br>
***************************************************<br><a href="mailto:isaac.ullah@asu.edu">isaac.ullah@asu.edu</a><br><a href="mailto:ullah@archaeologist.com">ullah@archaeologist.com</a><br><br><a href="http://www.public.asu.edu/~iullah">http://www.public.asu.edu/~iullah</a><br>
***************************************************<br>