Dear Wangyal,
Well, with the Perl bindings you could do something like this (note: not
tested code!). I'm assuming many things: the rasters are similar size,
they are not too big so they fit into memory, etc.
use Modern::Perl;
use Geo::GDAL;
my ($admin,$lu) = @ARGV;
my $ds_admin =
I am looking for a script to combine two raster data that have unique values.
For example, one raster has administrative unit name and another has land cover
data. I want to combine these two raster data to find a relationship between
two data. I can do the same process in another GIS software