我正在使用OGR/GDAL的C#管理API来写入GeoPackage数据库,但当我尝试将数据库同步到磁盘时,我总是遇到“数据库已锁定”的错误。经过多次试验,这似乎只有在我首先打开GeoPackage进行读取、关闭数据库(处置DataSource)、重新打开数据库并尝试写入时才会发生。Web搜索表明,如果在尝试写入之前未处理读取的GDAL和SQLite对象,则会发生此错误。但是,我认为我正在正确地处理所有托管对象,所以我想知道是否有我错过的调用。任何帮助都将不胜感激。我已经包含了一个示例控制台应用程序,它显示了我收到的错误。它创建了一个新的层,但如果我尝试更新一个现有的层,我会得到相同的错误。我在Windows上使用的是GDALV1.2、GDALV2.3和SQLite v3.21.0。
示例控制台应用程序:
using System;
using System.Collections.Generic;
using OSGeo.OGR;
namespace GeoPackageWriteTest
{
class GeoPackageTest
{
public static void usage()
{
Console.WriteLine( "usage (to create a new layer):" );
Console.WriteLine( "gpkgWriteTest {gpkg filename} {layername}" );
System.Environment.Exit( 0 );
}
static void Main( string[] args )
{
if ( args.Length < 2 )
usage();
Ogr.RegisterAll();
// simulate an object opening a geopackage db to do a simple read
Console.WriteLine( "Opening database readonly: " + args[0] );
DataSource dataSource1 = Ogr.OpenShared( args[0], 0 ); // open readonly
if ( dataSource1 == null )
throw new ArgumentException( string.Format( "Error opening GeoPackage database: {0}", args[0] ) );
using ( dataSource1 )
{
Driver driver = dataSource1.GetDriver();
if ( driver == null )
throw new ArgumentException( "Error getting GDAL driver!" );
Console.WriteLine( "Using driver: " + driver.name );
string pfile = GetParameterTableName( dataSource1 );
Console.WriteLine( "Parameter table: " + pfile );
} // DataSource gets disposed and geopackage is closed
dataSource1 = null;
// make sure all Disposed objects are garbage collected
// (just trying to see if it makes a difference, it doesn't)
GC.Collect( GC.MaxGeneration, GCCollectionMode.Forced, true, true );
GC.WaitForPendingFinalizers();
// simulate another object reopening a geopackage db to do a simple write
Console.WriteLine( "Opening database readwrite: " + args[0] );
DataSource dataSource2 = Ogr.OpenShared( args[0], 1 ); // reopen readwrite
if ( dataSource2 == null )
throw new ArgumentException( string.Format( "Error opening GeoPackage database: {0}", args[0] ) );
using ( dataSource2 )
{
List<string> options = new List<string>();
options.Add( "SPATIAL_INDEX=YES" );
Console.WriteLine( "Creating layer: " + args[1] );
Layer tableLayer = dataSource2.CreateLayer( args[1], null, wkbGeometryType.wkbNone, options.ToArray() );
if ( tableLayer == null )
throw new ArgumentException( string.Format( "Error creating layer: {0}", args[2] ) );
int rc = Ogr.OGRERR_NONE;
int approx_ok = 1;
int fieldCount = 3;
for ( int i = 0; i < fieldCount; i++ )
{
//FieldDefn fieldDefn = GetFieldDefn( args[n], args[n + 1] );
string fname = string.Format( "TestField{0}", i );
FieldDefn fieldDefn = GetFieldDefn( fname, "int32" );
rc = tableLayer.CreateField( fieldDefn, approx_ok );
if ( rc != Ogr.OGRERR_NONE )
throw new ArgumentException( string.Format( "Error creating field: {0}", fname ) );
}
// write changes to disk
Console.WriteLine( "Writing changes to disk..." );
dataSource2.FlushCache(); // this throws database locked error
}
}
static string GetParameterTableName( DataSource ogrDataSource )
{
try
{
string s = string.Empty;
string sql = "select * from sqlite_master where type='table' and name like '%Parameter%'";
OSGeo.OGR.Layer layer = ogrDataSource.ExecuteSQL( sql, null, "SQLITE" );
using ( layer )
{
layer.ResetReading();
OSGeo.OGR.Feature fea = layer.GetNextFeature();
if ( fea != null )
{
s = fea.GetFieldAsString( "name" );
}
}
return s;
}
catch ( Exception )
{
return string.Empty;
}
}
private static FieldDefn GetFieldDefn( string fldName, string fldType )
{
//OSGeo.OGR.FieldType fieldType = GetOGRFieldType( fldType );
OSGeo.OGR.FieldType fieldType = FieldType.OFTInteger;
OSGeo.OGR.FieldSubType subType = FieldSubType.OFSTNone;
FieldDefn fieldDefn = new FieldDefn( fldName, fieldType );
fieldDefn.SetSubType( subType );
return fieldDefn;
}
}
}发布于 2019-06-11 01:12:27
以防有人好奇。这里的问题是我第一次读取结果集层的方式。我只读了一个特写(因为那是我唯一感兴趣的)。显然,这会使结果集处于繁忙状态。我需要读取整个层,即
OSGeo.OGR.Feature fea = layer.GetNextFeature();
while ( fea != null )
{
s = fea.GetFieldAsString( "name" );
layer.GetNextFeature();
}或者我需要在读完后重置图层,即
OSGeo.OGR.Feature fea = layer.GetNextFeature();
if ( fea != null )
{
s = fea.GetFieldAsString( "name" );
}
layer.ResetReading();https://stackoverflow.com/questions/56452006
复制相似问题