本文为原创文章，转载注明出处，欢迎关注网站https://hkvision.cn
缘起
一个项目需要用GDAL这个库，然后去稍微深入看了下这玩意，感觉还是有很多坑的。
简介
就不介绍了，需要用这玩意的肯定知道这是啥
读写
最重要的肯定是读写
读:
1
2
3
4
5
6
7
8
9
| from osgeo import gdal
ds = gdal.Open("{image_path}")
"""
有很多种Open的方法
OpenShared
OpenEx
还可以设置只读
flags=gdal.GA_ReadOnly
"""
|
写:
1
2
| drv = gdal.GetDriverByName("GTiff")
drv.Create("{filepath}", xsize, ysize, bds_num, eType. options=[])
|
看起来一切都很美好，但是这里面蕴含着一个很大的坑，读影像还好，但是写影像需要指定路径，并且，gdal没有额外的方式能创建一个空的DataSet，这意味着如果你想弄一些临时的DataSet（这非常普遍，稍微复杂一点的算法都会需要有这个），你必须要存在硬盘上一份文件，这简直太蠢了。
VRT
当然这也不是没法解决的，GDAL其实有解决方案，那就是VRT文件，这个是GDAL的一种虚拟文件，是对影像的一种抽象表述，这个设定是很有意思的，给大家看个VRT文件的例子
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
| <VRTDataset rasterXSize="4000" rasterYSize="4000" subClass="VRTPansharpenedDataset">
<Metadata domain="IMAGE_STRUCTURE">
<MDI key="INTERLEAVE">PIXEL</MDI>
</Metadata>
<VRTRasterBand dataType="Byte" band="1" subClass="VRTPansharpenedRasterBand">
<ColorInterp>Red</ColorInterp>
</VRTRasterBand>
<VRTRasterBand dataType="Byte" band="2" subClass="VRTPansharpenedRasterBand">
<ColorInterp>Green</ColorInterp>
</VRTRasterBand>
<VRTRasterBand dataType="Byte" band="3" subClass="VRTPansharpenedRasterBand">
<ColorInterp>Blue</ColorInterp>
</VRTRasterBand>
<BlockXSize>512</BlockXSize>
<BlockYSize>512</BlockYSize>
<PansharpeningOptions>
<Algorithm>WeightedBrovey</Algorithm>
<AlgorithmOptions>
<Weights>0.3333333333333333,0.3333333333333333,0.3333333333333333</Weights>
</AlgorithmOptions>
<Resampling>Cubic</Resampling>
<SpatialExtentAdjustment>Union</SpatialExtentAdjustment>
<PanchroBand>
<SourceFilename relativeToVRT="1">pan_mask_opencv.tif</SourceFilename>
<SourceBand>1</SourceBand>
</PanchroBand>
<SpectralBand dstBand="1">
<SourceFilename relativeToVRT="1">origin_mask_opencv.tif</SourceFilename>
<SourceBand>1</SourceBand>
</SpectralBand>
<SpectralBand dstBand="2">
<SourceFilename relativeToVRT="1">origin_mask_opencv.tif</SourceFilename>
<SourceBand>2</SourceBand>
</SpectralBand>
<SpectralBand dstBand="3">
<SourceFilename relativeToVRT="1">origin_mask_opencv.tif</SourceFilename>
<SourceBand>3</SourceBand>
</SpectralBand>
</PansharpeningOptions>
</VRTDataset>
|
这里我们定义了一个Pansharpened文件，也就是一个影像融合之后的影像，但是实际上在硬盘上的只是一个这样的描述文件，我如果想把他变成真正的影像文件，我可以这样做
1
2
3
4
5
| from osgeo import gdal
ds = gdal.Open("pansharpened.vrt")
drv = gdal.GetDriverByName("GTiff")
drv.CreateCopy("pansharpened.tif", ds)
|
这样我就将一个很小的vrt文件变成了tif文件。
聪明的你肯定看出来了，这个vrt文件在实际编程的时候和这个tif文件没啥区别，但是不一样的是，vrt文件只占用这么小的空间，这是很爽的事情。
那么我们可以更进一步，这个vrt文件我都不想看到，怎么办？也可以，你只需要这样做
1
2
3
4
5
| # doing something
# now we want to create a new DataSet
drv = gdal.GetDriverByName("VRT")
ds = drv.Create("", xoff, yoff, bds_num, eType. options=[])
|
是的，如果是vrt文件，大可以在filename这个参数上填空串，这样就不会出现这个文件了。至于这个文件是放在了内存中还是再怎么样，我也就不关心了（GDAL也有内存文件，可以参考这个）
一些命令
GDAL除了提供库以外，还提供了一些命令行程序，可以对影像进行一些简单的操作，具体的内容就不细说了，看这里
B站上还有个视频教程，不过是全英文无字幕的，看着操作还是能明白是什么意思，在这里
part1
part2
