回答此问题可获得 20 贡献值,回答如果被采纳可获得 50 分。
<p>我有一个GDAL数据集,它只存在于内存中(inmem\ds):</p>
<pre><code>In: gdal.Info(inmem_ds)
Out:
('Driver: NUMPY/Numeric Python Array\n'
'Files: none associated\n'
'Size is 500, 500\n'
'Coordinate System is:\n'
'PROJCS["WGS 84 / UTM zone 19N",\n'
' GEOGCS["WGS 84",\n'
' DATUM["WGS_1984",\n'
' SPHEROID["WGS 84",6378137,298.257223563,\n'
' AUTHORITY["EPSG","7030"]],\n'
' AUTHORITY["EPSG","6326"]],\n'
' PRIMEM["Greenwich",0,\n'
' AUTHORITY["EPSG","8901"]],\n'
' UNIT["degree",0.0174532925199433,\n'
' AUTHORITY["EPSG","9122"]],\n'
' AUTHORITY["EPSG","4326"]],\n'
' PROJECTION["Transverse_Mercator"],\n'
' PARAMETER["latitude_of_origin",0],\n'
' PARAMETER["central_meridian",-69],\n'
' PARAMETER["scale_factor",0.9996],\n'
' PARAMETER["false_easting",500000],\n'
' PARAMETER["false_northing",0],\n'
' UNIT["metre",1,\n'
' AUTHORITY["EPSG","9001"]],\n'
' AXIS["Easting",EAST],\n'
' AXIS["Northing",NORTH],\n'
' AUTHORITY["EPSG","32619"]]\n'
'Origin = (617085.000000000000000,-1512585.000000000000000)\n'
'Pixel Size = (30.000000000000000,-30.000000000000000)\n'
'Corner Coordinates:\n'
'Upper Left ( 617085.000,-1512585.000) ( 67d55\' 2.38"W, 13d40\'47.58"S)\n'
'Lower Left ( 617085.000,-1527585.000) ( 67d55\' 0.13"W, 13d48\'55.77"S)\n'
'Upper Right ( 632085.000,-1512585.000) ( 67d46\'43.12"W, 13d40\'45.26"S)\n'
'Lower Right ( 632085.000,-1527585.000) ( 67d46\'40.59"W, 13d48\'53.42"S)\n'
'Center ( 624585.000,-1520085.000) ( 67d50\'51.55"W, 13d44\'50.54"S)\n'
'Band 1 Block=500x1 Type=Int16, ColorInterp=Undefined\n'
' NoData Value=-9999\n'
'Band 2 Block=500x1 Type=Int16, ColorInterp=Undefined\n'
' NoData Value=-9999\n'
'Band 3 Block=500x1 Type=Int16, ColorInterp=Undefined\n'
' NoData Value=-9999\n'
'Band 4 Block=500x1 Type=Int16, ColorInterp=Undefined\n'
' NoData Value=-9999\n'
'Band 5 Block=500x1 Type=Int16, ColorInterp=Undefined\n'
' NoData Value=-9999\n'
'Band 6 Block=500x1 Type=Int16, ColorInterp=Undefined\n'
' NoData Value=-9999\n'
'Band 7 Block=500x1 Type=Int16, ColorInterp=Undefined\n'
' NoData Value=-9999\n'
'Band 8 Block=500x1 Type=Int16, ColorInterp=Undefined\n'
' NoData Value=-9999\n')
</code></pre>
<p>现在我想把这个数据集直接传递给一个使用gdal.打开()打开其输入数据集。所以我需要一个可以传递到的物理文件路径gdal.打开(). 但我不想为了节省IO时间而把内存中的GDAL数据集写到磁盘上。在</p>
<p><strong>是否有方法将内存中的数据集传递给gdal.打开()?</strong></p>
<p><strong>我尝试</strong>使用VRT数据集,但不起作用:</p>
^{pr2}$
<p>保存的VRT文件(mypath/测试.vrt)看起来像这样:</p>
<pre><code><VRT<a href="https://www.cnpython.com/pypi/dataset" class="inner-link">dataset</a> rasterXSize="500" rasterYSize="500">
<SRS>PROJCS["WGS 84 / UTM zone 19N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-69],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","32619"]]</SRS>
<GeoTransform> 6.1708500000000000e+05, 3.0000000000000000e+01, 0.0000000000000000e+00, -1.5125850000000000e+06, 0.0000000000000000e+00, -3.0000000000000000e+01</GeoTransform>
<VRTRasterBand dataType="Int16" band="1">
<NoDataValue>-9999</NoDataValue>
<SimpleSource>
<SourceFilename relativeToVRT="0"></SourceFilename>
<SourceBand>1</SourceBand>
<SourceProperties RasterXSize="500" RasterYSize="500" DataType="Int16" BlockXSize="500" BlockYSize="1" />
<SrcRect xOff="0" yOff="0" xSize="500" ySize="500" />
<DstRect xOff="0" yOff="0" xSize="500" ySize="500" />
</SimpleSource>
</VRTRasterBand>
<VRTRasterBand dataType="Int16" band="2">
<NoDataValue>-9999</NoDataValue>
<SimpleSource>
<SourceFilename relativeToVRT="0"></SourceFilename>
<SourceBand>2</SourceBand>
<SourceProperties RasterXSize="500" RasterYSize="500" DataType="Int16" BlockXSize="500" BlockYSize="1" />
<SrcRect xOff="0" yOff="0" xSize="500" ySize="500" />
<DstRect xOff="0" yOff="0" xSize="500" ySize="500" />
</SimpleSource>
</VRTRasterBand>
<VRTRasterBand dataType="Int16" band="3">
<NoDataValue>-9999</NoDataValue>
<SimpleSource>
<SourceFilename relativeToVRT="0"></SourceFilename>
<SourceBand>3</SourceBand>
<SourceProperties RasterXSize="500" RasterYSize="500" DataType="Int16" BlockXSize="500" BlockYSize="1" />
<SrcRect xOff="0" yOff="0" xSize="500" ySize="500" />
<DstRect xOff="0" yOff="0" xSize="500" ySize="500" />
</SimpleSource>
</VRTRasterBand>
<VRTRasterBand dataType="Int16" band="4">
<NoDataValue>-9999</NoDataValue>
<SimpleSource>
<SourceFilename relativeToVRT="0"></SourceFilename>
<SourceBand>4</SourceBand>
<SourceProperties RasterXSize="500" RasterYSize="500" DataType="Int16" BlockXSize="500" BlockYSize="1" />
<SrcRect xOff="0" yOff="0" xSize="500" ySize="500" />
<DstRect xOff="0" yOff="0" xSize="500" ySize="500" />
</SimpleSource>
</VRTRasterBand>
<VRTRasterBand dataType="Int16" band="5">
<NoDataValue>-9999</NoDataValue>
<SimpleSource>
<SourceFilename relativeToVRT="0"></SourceFilename>
<SourceBand>5</SourceBand>
<SourceProperties RasterXSize="500" RasterYSize="500" DataType="Int16" BlockXSize="500" BlockYSize="1" />
<SrcRect xOff="0" yOff="0" xSize="500" ySize="500" />
<DstRect xOff="0" yOff="0" xSize="500" ySize="500" />
</SimpleSource>
</VRTRasterBand>
<VRTRasterBand dataType="Int16" band="6">
<NoDataValue>-9999</NoDataValue>
<SimpleSource>
<SourceFilename relativeToVRT="0"></SourceFilename>
<SourceBand>6</SourceBand>
<SourceProperties RasterXSize="500" RasterYSize="500" DataType="Int16" BlockXSize="500" BlockYSize="1" />
<SrcRect xOff="0" yOff="0" xSize="500" ySize="500" />
<DstRect xOff="0" yOff="0" xSize="500" ySize="500" />
</SimpleSource>
</VRTRasterBand>
<VRTRasterBand dataType="Int16" band="7">
<NoDataValue>-9999</NoDataValue>
<SimpleSource>
<SourceFilename relativeToVRT="0"></SourceFilename>
<SourceBand>7</SourceBand>
<SourceProperties RasterXSize="500" RasterYSize="500" DataType="Int16" BlockXSize="500" BlockYSize="1" />
<SrcRect xOff="0" yOff="0" xSize="500" ySize="500" />
<DstRect xOff="0" yOff="0" xSize="500" ySize="500" />
</SimpleSource>
</VRTRasterBand>
<VRTRasterBand dataType="Int16" band="8">
<NoDataValue>-9999</NoDataValue>
<SimpleSource>
<SourceFilename relativeToVRT="0"></SourceFilename>
<SourceBand>8</SourceBand>
<SourceProperties RasterXSize="500" RasterYSize="500" DataType="Int16" BlockXSize="500" BlockYSize="1" />
<SrcRect xOff="0" yOff="0" xSize="500" ySize="500" />
<DstRect xOff="0" yOff="0" xSize="500" ySize="500" />
</SimpleSource>
</VRTRasterBand>
</VRTDataset>
</code></pre>