<DIV>hi all,</DIV>

<DIV>i am using grass54 compiled from source code in a Debian (Sarge) OS...</DIV>

<DIV>i was trying to use i.class to get signatures to a supervisioned classification...but, after drawing the needed regions&nbsp;i could not use the "analyse" tool or save that specific signature:</DIV>

<DIV>error messages:</DIV>

<DIV>"cannot analyse until region is completed (click mouse to continue)"</DIV>

<DIV>and</DIV>

<DIV>"warning:prepare_signature: data error</DIV>

<DIV>&nbsp;</DIV>

<DIV>all my steps:</DIV>

<DIV>&nbsp;</DIV>

<DIV>click zoom</DIV>

<DIV>click draw</DIV>

<DIV>click done</DIV>

<DIV>(after several samples)</DIV>

<DIV>click analyze -&gt; error message</DIV>

<DIV>&nbsp;</DIV>

<DIV>thanks in advance</DIV>

<DIV>samuel</DIV>

<DIV><BR><BR><B><I>Tom Russo &lt;russo@bogodyn.org&gt;</I></B> wrote:</DIV>

<BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid">I've recently been trying to georeference scanned images of old (i.e. 100 year old) US Geological Survey topo maps. I'm using GRASS 5.7. I've run into  some problems, naturally. The maps all have Lat/Lon grids, and are all in "polyconic projection" ---  with no specified central meridian, so it is pointless to try to georeference  them into a GRASS polyconic projection location. Trying to use a first order affine transformation in i.rectify to georeference the map into a UTM location  does not work, because we need more than a simple scaling, rotation, and  translation. The USGS itself doesn't even do that when it georeferences its DRGs: they create a set of control points based on the lat/lon grid intersections, and do a piecewise linear rubbersheeting into UTM.  Since GRASS doesn't have a piecewise linear rubbersheeting module!

 , I

 tried<BR>to use i.group,i.target,i.points, and i.rectify with a transformation of <BR>higher order than the simple linear affine transformation, in the hopes that<BR>it would be a reasonable approximation. I ran into some issues:<BR><BR>1) i.points has an "analyze" function that is supposed to tell the RMS error<BR>of the transformation based on the points you supply. However, it only<BR>uses the linear affine transformation, and cannot tell you the quality of<BR>any higher order transformation you might use in i.rectify<BR><BR>2) i.rectify's documentation implies that it reports RMS error, but in fact<BR>it does not --- in order to get any output at all that tells how well<BR>higher order transformations are doing, you must recompile i.rectify with<BR>a "#define BDEBUG" uncommented in imagery/i.rectify/crs.c to get any such<BR>output, and it reports error in the image coordinates, not in meters in<BR>projected coordinates.<BR><BR>3) there is a module called i.vpoints that

apparently does have code that would allow you to specify 1st, 2nd, or 3rd order transformations in the analys function, but when I attempt to use that menu it enters some sort of  infinte loop waiting for mouse input, and no matter what I do it doesn't respond to mouse clicks or keyboard input. It has to be killed from another window. By using the debugging output of i.rectify with the debug #define included I can have a tiny bit of confidence that choosing a 3rd order transformation using the 16 lat/lon grid intersections is somewhat reasonable, but it's not as good as having the simple "is the RMS error below the pixel size?" test that the first order affine transformation can have.  Does anybody use i.vpoints? Has anyone else seen this strange behavior when choosing the transformation order menu? Is anyone aware of free (or cheap) software that can perform general piecewise linear rubbersheeting of rasters to georefe!

 rence

 them into a <BR>projection that doesn't match the ones they were created in? USGS refers<BR>to some in-house software that is unavailable. There's apparently a<BR>set of Matlab image processing functions that can do this (the method <BR>creates a Delaunay triangulation of the net of control points, and does<BR>a local linear affine transformation of points within each triangle). But<BR>I can't find anything that could be used easily with GRASS. Does anyone have<BR>suggestions?<BR><BR>-- <BR>Tom Russo KM5VY SAR502 DM64ux http://www.swcp.com/~russo/<BR>Tijeras, NM QRPL#1592 K2#398 SOC#236 AHTB#1 http://www.qsl.net/~km5vy/<BR>"When life gives you lemons, find someone with a paper cut."<BR><BR></BLOCKQUOTE><p>

<hr size=1><a href="http://br.rd.yahoo.com/mail/taglines/*http://br.acesso.yahoo.com/">Yahoo! Acesso Grátis</a> - Internet rápida e grátis. Instale o discador do Yahoo! agora.