Questions tagged «gdal-translate»

3
创建透明的山体阴影?
通过组合使用gdal,我得到了优雅的透明山体阴影convert。与通常使用的基于灰色的山体阴影相比,这种透明的山体阴影非常酷,因为它们可以放置在地图背景与其他上层(道路,建筑物)之间,以提供3D感觉,无论背景的类型和颜色如何。 怎么运行的 诀窍:从产生的灰度山丘阴影开始gdal hillshade,诀窍是获取灰度通道的值,将每个值求反,然后将结果流到新的不透明度通道中。黑色像素[0,0,0]变为[0,0,0,255](不透明度= 255),灰色像素[120,120,120]降低至[120,120,120,135](不透明度= 135,又称为255-120),白色像素[255,255,255]变为透明[255,255,255,0](不透明度= 0,又称为255-255),依此类推。丘陵的阴影是不透明的和黑色的,平原逐渐变得(白色)透明。从概念上讲,像素的波段方程是这样的: alpha_value([0]) = 255-0 = 255; // black pixels get full opacity alpha_value([Grey]) = 255-Grey; // grey pixels get 255-GreyValue opacity. alpha_value([255]) = 255-255 = 0 // white pixels get null opacity ! 有关此方法的悠闲视频(由Photoshop设计师解释),请参阅在Photoshop中添加阴影浮雕(16分钟)。 题 将ETOPO或SRTM派生的基于灰色的山体阴影(gdaldem hillshade crop_xl.tmp.tif shadedrelief.tmp.tif -s 111120 -z 5 -az …

1
PyQGIS将栅格另存为渲染图像,然后在其上使用GDAL工具
我尝试使用QGis的“将栅格保存为渲染图像”中的代码。i是栅格图层: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), pipelayer.height()) piperenderer = pipelayer.renderer() pipeprovider = pipelayer.dataProvider() crs = pipelayer.crs().toWkt() pipe = QgsRasterPipe() pipe.set(pipeprovider.clone()) pipe.set(piperenderer.clone()) pipedFile = os.path.join(tempfile.gettempdir(), safeLayerName + '_pipe.tif') print pipedFile file_writer = QgsRasterFileWriter(pipedFile) file_writer.writeRaster(pipe, pipewidth, pipeheight, pipeextent, pipelayer.crs()) in_raster = pipedFile prov_raster = os.path.join(tempfile.gettempdir(), 'json_' + …

1
用gdal拼接图像时能否保留褪色的alpha图层?
我正在用gdal镶嵌一些图像,并想通过朝每个图像的边缘使用淡入/渐变Alpha层来消除马赛克中间的尖锐边缘来改善最终结果。我遇到的问题是,每个带有渐变Alpha层的图像的一部分在最终的镶嵌图中都将其下方的图像掩盖了,而不是半透明的,如下所示: 理想情况下,我希望使用这种渐变透明度将一张图像逐渐淡入另一张图像。 我执行的生成镶嵌的步骤如下: 将gcps添加到原始图像以对其进行地理位置定位和正确定向(依次对每个图像进行处理): gdal_translate -of GTiff -a_srs EPSG:4326 -a_srs EPSG:4326 -gcp 1616 0 -88.2728612066 40.5175787437 -gcp <etc., etc.> <original_image_with_gradual_alpha>.tif <image_with_gradual_alpha_and_gcps>.tif 将图像扭曲到正确定向的新geotiff(依次对每个图像执行): gdalwarp -s_srs EPSG:4326 -t_srs EPSG:4326 -dstnodata 0 <image_with_gradual_alpha_and_gcps>.tif <warped_geotiff_with_alpha>.tif 将所有变形的图像合并到一个马赛克中: gdalbuildvrt -srcnodata 0 mosaic.vrt <warped_geotiff_with_alpha_root>*.tif gdal_translate mosaic.vrt mosaic.tif 我链接的图像是mosaic.tif。 gdalinfo提供样本输入文件: Driver: GTiff/GeoTIFF Files: dsc00562.tif Size is 1616, 1080 …

2
将巨大的XYZ CSV转换为GeoTIFF
我有大量CSV格式的数据,其中包含UTM坐标X和Y和高程值作为Z信息。我需要将这些数据转换为GeoTIFF的DEM以便进一步分析。在这种情况下,大量意味着16 m。线,具有在一个点X,Y和Z每行。这些点是均匀分布的,因此不需要插值;每个点仅需要转换为栅格像元。 原始数据不带分隔符,列宽固定。我已经弄清楚了如何使用流文本编辑器sed将文件语法转换为使用分隔符而不是固定宽度并消除所有空格字符。从这里开始,通常,我的工作流程是通过根据和数据创建要素类,将数据导入ArcGIS X,Y然后Z第二步,使用“ 指向栅格”工具将点shapefile转换为GeoTIFF 。但是,对于该过程,我当前拥有的文件太大。 除了寻找上述工作流程之外,我一直在寻找一种有效的替代方法,并发现了GDAL。但是,在gdal_translate,我可以在支持的文件类型列表中找到的最接近的受支持格式是ASCII网格,但没有逗号分隔的XYZ。另一个困难是,我有UTM坐标,而大多数示例似乎都使用十进制度坐标。但是,我需要停留在UTM系统内(或者至少,我的输出GeoTIFF必须位于UTM坐标系内)。 因此,我正在寻找一种使用GDAL将CSV XYZ转换为GeoTIFF的方法,但到目前为止,尚无法找到处理此确切问题的示例。对于某些提示甚至代码示例,我将感到非常高兴。

1
使用GDAL在范围内将16位缩放为8位
我正在寻找使用GDAL(OSGeo4W)将16位无符号tiff图像缩放为8位图像。但是,我想限制像素值,以使它们保持在8位图像的范围内。我的意思是,我想确保16位图像中的视觉比例在8位图像中表示(0-255)。使用gdal_translate将16位转换为8位会剪切像素值,并且不会生成RGB图像。我不确定从gdal_translate中调用什么函数和设置,或者是否有更好的选择。 GDALINFO C:\>gdalinfo C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Driver: GTiff/GeoTIFF Files: C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Size is 17310, 11310 Coordinate System is `' Metadata: TIFFTAG_DATETIME=2015:05:03 12:27:13 TIFFTAG_IMAGEDESCRIPTION=UltraCam-Lvl03 -------------- CAM_ID: UC-SXp-1-50215465 [3] IMG_NO: 4761 RECORD_GUID: 8fc3d7a9-d2e5-40e5-babf-939f1f803dcf IMG_GUID: C850D0A4-39CA-4BC8-ACD7-81A363D78C6E FILE_GUID: 86C4F8E4-A2FD-43FB-B3B1-4C4ED8DF0E42 LICENSE_ID: 820342819 SOFTWARE: UltraCam Aerial Radiometry Core 12.6.1408.2501 PIXEL_SIZE_WIDTH: 6 [micron] PIXEL_SIZE_HEIGHT: 6 [micron] APERTURE: F_8 EXPOSURE_TIME: 0.002000000000 …

1
如何使用GDAL将色表添加到一个波段的Tiff中?
我有1个波段的GeoTIFF文件。带的颜色解释为灰色。我想要一个颜色表,以将其转换Col.tiff为Col.jpeg或Col.png。 我尝试通过.vrt在此tiff中放置一个颜色表: gdal_translate -of VRT Col.tiff Col.vrt 然后添加.vrt: <ColorInterp>Palette</ColorInterp> <ColorTable> <Entry c1="0" c2="0" c3="255" c4="240"/> <Entry c1="0" c2="0" c3="254" c4="240"/> <Entry c1="0" c2="0" c3="253" c4="240"/> <Entry c1="0" c2="0" c3="252" c4="240"/> </ColorTable> 并再次翻译: gdal_translate Col.vrt Col.tiff 但它回答: Warning 1: Unable to export color table to GeoTIFF file. Color tables can only …

1
为什么多个栅格合并的结果这么大?[关闭]
已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗?更新问题,使其仅通过编辑此帖子来关注一个问题。 2年前关闭。 我尝试像这样合并14个geotiff: 每个geotiff约为50Mb。我在输出端需要一个Geotiff 我的工作流程: gdalbuildvrt -input_file_list list.txt test.vrt (我的列表中包含tif的名称) 然后 : gdal_translate -of Gtiff test.vrt test.tif Input file size is 79841, 59955 它可以工作,但是结果是13,3 Gb的geotiff!对于14个文件(每个50 Mb),我尝试了700 Mb(而不是13 Gb)的geotiff。 我知道gdal默认不会压缩,所以我尝试了以下命令: gdal_translate -of Gtiff -co COMPRESS=JPEG test.vrt test_compressed.tif 但是文件的“合并”对于JPEG压缩来说太大了: Input file size is 79841, 59955 0ERROR 1: JPEGPreEncode:Strip/tile too large for …
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.