I might try this..is it the good way to do what I need? 1- Use gdal_polygonize.py to create a boundary (excluding pixels with no value) of each image and save the boundary as a .shp 2- Somehow find the intersected region of the 2 polygons (boundaries) created above and save the result to .shp 3- Use gdal_rasterize on one of the original image to burn a nodata value based on the polygon representing the intersection area
I didn't try yet because I don't have gdal_polygonize.py My installation is from MS4W and this script is not there. I just install GDAL with OSGEO4w and it seems that this script doesn't exist. If the method above is valid, - how can I create a boundary excluding pixels with nodata? - And, of course, where can I get gdal_polygonize.py Thanks! Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques biologiques, environnementaux et occupationnels Institut national de santé publique du Québec 945, avenue Wolfe Québec, Qc G1V 5B3 Tél.: (418) 650-5115 #5281 Fax.: (418) 654-3144 [email protected] http://www.inspq.qc.ca Steve Toutant/INSPQ/SSSS 17/07/2009 11:26 AM A [email protected] cc Objet clip intersected Hi, I have 2 Erdas (.img) images that overlap. See attached image. I need to clip one of the 2 images by removing the intersected (dark red) part. If it's possible, how should I accomplish this? thanks for your help, Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques biologiques, environnementaux et occupationnels Institut national de santé publique du Québec 945, avenue Wolfe Québec, Qc G1V 5B3 Tél.: (418) 650-5115 #5281 Fax.: (418) 654-3144 [email protected] http://www.inspq.qc.ca
<<image/gif>>
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
