/// <summary>
/// Get a feature of the <see cref="VectorTileLayer"/>
/// </summary>
/// <param name="layer"><see cref="VectorTileLayer"/> containing the feature</param>
/// <param name="data">Raw byte data of the feature</param>
/// <param name="validate">If true, run checks if the tile contains valid data. Decreases decoding speed.</param>
/// <param name="clippBuffer">
/// <para>'null': returns the geometries unaltered as they are in the vector tile. </para>
/// <para>Any value >=0 clips a border with the size around the tile. </para>
/// <para>These are not pixels but the same units as the 'extent' of the layer. </para>
/// </param>
/// <returns></returns>
public static VectorTileFeature GetFeature(
VectorTileLayer layer
, byte[] data
, bool validate = true
, uint?clippBuffer = null
)
{
PbfReader featureReader = new PbfReader(data);
VectorTileFeature feat = new VectorTileFeature(layer);
bool geomTypeSet = false;
while (featureReader.NextByte())
{
int featureType = featureReader.Tag;
if (validate)
{
if (!ConstantsAsDictionary.FeatureType.ContainsKey(featureType))
{
throw new System.Exception(string.Format("Layer [{0}] has unknown feature type: {1}", layer.Name, featureType));
}
}
switch ((FeatureType)featureType)
{
case FeatureType.Id:
feat.Id = (ulong)featureReader.Varint();
break;
case FeatureType.Tags:
#if NET20
List <int> tags = featureReader.GetPackedUnit32().ConvertAll <int>(ui => (int)ui);
#else
List <int> tags = featureReader.GetPackedUnit32().Select(t => (int)t).ToList();
#endif
feat.Tags = tags;
break;
case FeatureType.Type:
int geomType = (int)featureReader.Varint();
if (validate)
{
if (!ConstantsAsDictionary.GeomType.ContainsKey(geomType))
{
throw new System.Exception(string.Format("Layer [{0}] has unknown geometry type tag: {1}", layer.Name, geomType));
}
}
feat.GeometryType = (GeomType)geomType;
geomTypeSet = true;
break;
case FeatureType.Geometry:
if (null != feat.Geometry)
{
throw new System.Exception(string.Format("Layer [{0}], feature already has a geometry", layer.Name));
}
//get raw array of commands and coordinates
List <uint> geometryCommands = featureReader.GetPackedUnit32();
//decode commands and coordinates
List <List <Point2d> > geom = DecodeGeometry.GetGeometry(
layer.Extent
, feat.GeometryType
, geometryCommands
);
if (clippBuffer.HasValue)
{
geom = clipGeometries(geom, feat.GeometryType, (long)layer.Extent, clippBuffer.Value);
}
feat.Geometry = geom;
break;
default:
featureReader.Skip();
break;
}
}
if (validate)
{
if (!geomTypeSet)
{
throw new System.Exception(string.Format("Layer [{0}]: feature missing geometry type", layer.Name));
}
if (null == feat.Geometry)
{
throw new System.Exception(string.Format("Layer [{0}]: feature has no geometry", layer.Name));
}
if (0 != feat.Tags.Count % 2)
{
throw new System.Exception(string.Format("Layer [{0}]: uneven number of feature tag ids", layer.Name));
}
if (feat.Tags.Count > 0)
{
#if NET20
int maxKeyIndex = -9999;
for (int i = 0; i < feat.Tags.Count; i += 2)
{
if (feat.Tags[i] > maxKeyIndex)
{
maxKeyIndex = feat.Tags[i];
}
}
int maxValueIndex = -9999;
for (int i = 1; i < feat.Tags.Count; i += 2)
{
if (feat.Tags[i] > maxValueIndex)
{
maxValueIndex = feat.Tags[i];
}
}
#else
int maxKeyIndex = feat.Tags.Where((key, idx) => idx % 2 == 0).Max();
int maxValueIndex = feat.Tags.Where((key, idx) => (idx + 1) % 2 == 0).Max();
#endif
if (maxKeyIndex >= layer.Keys.Count)
{
throw new System.Exception(string.Format("Layer [{0}]: maximum key index equal or greater number of key elements", layer.Name));
}
if (maxValueIndex >= layer.Values.Count)
{
throw new System.Exception(string.Format("Layer [{0}]: maximum value index equal or greater number of value elements", layer.Name));
}
}
}
return(feat);
}
public static int Main(string[] args)
{
string vtIn = string.Empty;
uint? clipBuffer = null;
bool outGeoJson = false;
ulong? zoom = null;
ulong? tileCol = null;
ulong? tileRow = null;
for (int i = 0; i < args.Length; i++)
{
string argLow = args[i].ToLower();
if (argLow.Contains("vt:"))
{
vtIn = argLow.Replace("vt:", "");
}
else if (argLow.Contains("clip:"))
{
clipBuffer = Convert.ToUInt32(argLow.Replace("clip:", ""));
}
else if (argLow.Contains("out:"))
{
outGeoJson = argLow.Replace("out:", "").Equals("geojson");
}
else if (argLow.Contains("tileid:"))
{
parseArg(argLow.Replace("tileid:", ""), out zoom, out tileCol, out tileRow);
}
}
if (!File.Exists(vtIn))
{
Console.WriteLine($"file [{vtIn}] not found");
usage();
return(1);
}
// z-x-y weren't passed via parameters, try to get them from file name
if (!zoom.HasValue || !tileCol.HasValue || !tileRow.HasValue)
{
if (!parseArg(Path.GetFileName(vtIn), out zoom, out tileCol, out tileRow))
{
usage();
return(1);
}
}
var bufferedData = File.ReadAllBytes(vtIn);
VectorTile tile = new VectorTile(bufferedData);
if (outGeoJson)
{
Console.WriteLine(tile.ToGeoJson(zoom.Value, tileCol.Value, tileRow.Value, clipBuffer));
}
else
{
foreach (string lyrName in tile.LayerNames())
{
VectorTileLayer lyr = tile.GetLayer(lyrName);
Console.WriteLine(string.Format("------------ {0} ---------", lyrName));
//if (lyrName != "building") { continue; }
int featCnt = lyr.FeatureCount();
for (int i = 0; i < featCnt; i++)
{
VectorTileFeature feat = lyr.GetFeature(i, clipBuffer);
Console.WriteLine(string.Format("feature {0}: {1}", i, feat.GeometryType));
Dictionary <string, object> props = feat.GetProperties();
foreach (var prop in props)
{
Console.WriteLine(string.Format(" {0}\t : ({1}) {2}", prop.Key, prop.Value.GetType(), prop.Value));
}
}
}
}
return(0);
}
