type file;
type imagefile;
type landuse;

app (landuse output) getLandUse (imagefile input, int sortfield)
{
  getlanduse @input sortfield stdout=@output ;
}

app (file output, file tilelist) analyzeLandUse (landuse input[], int usetype, int maxnum)
{
  analyzelanduse @output @tilelist usetype maxnum @filenames(input);
}

app (imagefile output) colormodis (imagefile input)
{
  colormodis @input @output;
}

imagefile geos[]<filesys_mapper; suffix=".tif">;
landuse   land[]<structured_regexp_mapper; source=geos,match="(h..v..)", transform="\\1.landuse.byfreq">;

# Find the land use of each modis tile

foreach g,i in geos {
  land[i] = getLandUse(g,1);
}

# Find the top 10 most urban tiles (by area)

int UsageTypeURBAN=13;
file bigurban<"topurban.txt">;
file urbantiles<"urbantiles.txt">;
(bigurban, urbantiles) = analyzeLandUse(land, UsageTypeURBAN, 10);

# Map the files to an array

string urbanfilenames[] = readData(urbantiles);
imagefile urbanfiles[] <array_mapper;files=urbanfilenames>;

# Create a set of recolored images for just the urban tiles

foreach uf, i in urbanfiles {
  imagefile recoloredImage <single_file_mapper;
                       file=@strcat(@strcut(urbanfilenames[i],"(h..v..)"),".recolored.tif")>;
  recoloredImage = colormodis(uf);
}