Then, you can run a separate command to process the non-alpha bands.
You need two processes because there may not be one-to-one mapping of the
non-alpha and alpha band polygons.
Finally, you need to make intersections of the polygons in the two layers
and store the values of the alpha band polygon
John,
gdal_polygonize script stores the common pixel value in a field for every
polygon. By default it is named 'DN', which can be changed by specifying the
fieldname option.
You must use the alpha channel to polygonize. So you have to use the -b
option to specify that and may have to use the -no
Thanks Jay,
Ideally I would like the data to be in a column of the shapefile output so I
can load it in the database and use it to query. With gdal_polygonize.py
this doesn't seem possible so I was wondering what other options are out
there.
--
View this message in context:
http://osgeo-org.180
