[UW-GIS-L] (no subject)

[Joowon Park]

I hope somebody can help me.

My goal is making a AML code for i) selecting cells from a raster map, ii) getting the (x,y) coordinates of the selected cells, and iii) putting those coordinates into selectcircle funtion which can clip only the circle parts of the original raster map.

How can I do it? 
So far, I'm trying to use $$xmap, $$ymap. However, I don't know how to put each center location into a SELECTCIRCLE(INGRID, X, Y, R), especially iteratively..

FYI, I add the code for the only one cell case. It does not work, though.

x_ctr = scalar(0)
y_ctr = scalar(0)

DOCELL
if (RasterCell = WhatIwant)
  begin
   x_ctr = $$xmap
   y_ctr = $$ymap
  end
END

&sv x_input = [show scalar x_ctr]
&sv y_input = [show scalar y_ctr]

OUTGRID = SELECTCIRCLE(RasterCell, x_input, y_input, %radius%)

I'll deeply appreciate your any help. 
-Joowon