Michael, >>> On 7/6/2011 at 5:35 PM, in message <4e14ff42.50...@cironline.org>, Michael Corey <mco...@cironline.org> wrote: > Sure, I've uploaded samples here. > > http://www.mikejcorey.com/spatial/diablo-box-sample.tif > http://www.mikejcorey.com/spatial/diablo-cutout-sample.tif
I don't notice the semi-transparency in these scaled down images. Perhaps it is the way your viewer reads the mask? > > These are the same as the images created by the process I described (but > scaled down). > > To your point about specifying size in the first step -- will that make > the process run faster, or does it do the scaling down after it builds > the full-resolution image? > > Also, I notice that my filesize always gets significantly bigger when I > do the cutout step, which seems counter-intuitive to me since in theory > shouldn't there be less information present once the cutout is done? -cutline does not 'discard' any data. The extent of the data remains the same unless you reset those extents. You can do that with -crop_to_cutline. Here are some details from the gdalwarp page, http://gdal.org/gdalwarp.html : -crop_to_cutline: (GDAL >= 1.8.0) Crop the extent of the target dataset to the extent of the cutline. Polygon cutlines may be used as a mask to restrict the area of the destination file that may be updated, including blending. If the OGR layer containing the cutline features has no explicit SRS, the cutline features must be in the georeferenced units of the destination file. When outputing to a not yet existing target dataset, its extent will be the one of the original raster unless -te or -crop_to_cutline are specified. Best Regards, Eli > > Thanks for your help! > > Michael Corey > > > On 7/6/11 5:01 PM, Chaitanya kumar CH wrote: >> Michael, >> >> Can you provide screenshots of >> diablo-combined-center-utm10-70pct-box.tif and >> diablo-combined-center-utm10-70pct-cutout.tif for comparison? >> >> By the way, you can perform the actions of the two gdal_translate >> commands in the first step with the gdal_merge.py script itself unless >> you want to use a specific resampling algorithm. >> >> On Thu, Jul 7, 2011 at 4:28 AM, Michael Corey <mco...@cironline.org >> <mailto:mco...@cironline.org>> wrote: >> >> Hi all: >> >> I'm using a shapefile as a clipping mask to cut out the shoreline >> from some DOQ files that I have merged together. But when I do the >> clipping step, I end up with unwanted semitransparency in the >> non-clipped areas. >> >> I'm pretty sure the problem is only with my gdalwarp step at the end. >> >> Here's my process: >> >> gdal_merge.py -init "255" -o diablo-combined-center-utm10.tif file >> file file file >> >> gdal_translate -outsize 70% 70% diablo-combined-center-utm10.tif >> diablo-combined-center-utm10-70pct.tif >> >> ogrinfo -al ./diablo_canyon_detail_clipper.shp >> //Extent: (XXXX, YYYY) - (XXXX, YYYY) >> >> gdal_translate -projwin XXXX YYYY XXXX YYYY >> diablo-combined-center-utm10-70pct.tif >> diablo-combined-center-utm10-70pct-box.tif >> >> gdalwarp -co COMPRESS=DEFLATE -cutline >> ./diablo_canyon_detail_clipper.shp >> diablo-combined-center-utm10-70pct-box.tif >> diablo-combined-center-utm10-70pct-cutout.tif >> >> Can anyone help? >> >> Thanks! >> >> -- >> Michael Corey >> >> _______________________________________________ >> gdal-dev mailing list >> gdal-dev@lists.osgeo.org <mailto:gdal-dev@lists.osgeo.org> >> http://lists.osgeo.org/mailman/listinfo/gdal-dev >> >> >> >> >> -- >> Best regards, >> Chaitanya kumar CH. >> /t?a???nj?/ /k?m?r/ >> +91-9494447584 >> 17.2416N 80.1426E _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
