public void TestSetBitmap()
{
using (GribFile file = new GribFile(Settings.REG_LATLON_GRB1))
{
var msg = file.First();
int numVals = 10;
double[] vals = new double[numVals];
double val = 42;
for(int i =0; i < numVals; i++)
{
vals[i] = val;
}
msg.HasBitmap = true;
msg.SetValues(vals);
double[] vals2;
msg.Values(out vals2);
for (int i = 0; i < numVals; i++)
{
Assert.AreEqual(vals[i], vals2[i]);
Assert.AreEqual(val, vals2[i]);
}
}
}
public void TestIterateKeyValuePairs()
{
using (var file = new GribFile(Settings.SPHERICAL_PRESS_LVL)) {
Assert.IsTrue(file.MessageCount > 0);
Assert.IsTrue(file.First().Any());
}
}
public void TestOpenPng()
{
using (GribFile file = new GribFile(Settings.PNG_COMPRESSION))
{
Assert.IsTrue(file.MessageCount > 0);
var msg = file.First();
try
{
Assert.IsTrue(msg["packingType"].AsString().ToLower().EndsWith("_png"));
Assert.IsTrue(msg.ValuesCount > 0);
Assert.IsTrue(msg.GeoSpatialValues.Any());
int i = 0;
foreach (var v in msg.GeoSpatialValues)
{
Assert.AreNotEqual(Double.NaN, v.Value);
if (i++ > 1000)
{
break;
}
}
} catch (GribApiException e)
{
Console.WriteLine(e.Message);
Console.WriteLine(msg.ShortName);
Console.WriteLine(msg.ToString());
Assert.IsTrue(false);
}
}
}
public void TestSetBitmap()
{
using (GribFile file = new GribFile(Settings.REG_LATLON_GRB1))
{
var msg = file.First();
int numVals = 10;
double[] vals = new double[numVals];
double val = 42;
for (int i = 0; i < numVals; i++)
{
vals[i] = val;
}
msg.HasBitmap = true;
msg.SetValues(vals);
double[] vals2;
msg.Values(out vals2);
for (int i = 0; i < numVals; i++)
{
Assert.AreEqual(vals[i], vals2[i]);
Assert.AreEqual(val, vals2[i]);
}
}
}
public void TestPrecision ()
{
var file = new GribFile(Settings.GAUSS);
var msg = file.First();
Assert.AreEqual(msg["bitsPerValue"].AsInt(), 13);
msg.DecimalPrecision = 4;
Assert.AreEqual(msg["bitsPerValue"].AsInt(), 20);
}
public void TestIterateKeyValuePairs ()
{
using(var file = new GribFile(Settings.BIN))
{
Assert.IsTrue(file.MessageCount > 0);
Assert.IsTrue(file.First().Any());
}
}
public void TestDumpNoError()
{
using (GribFile file = new GribFile(Settings.REG_LATLON_GRB1))
{
var msg = file.First();
msg.Dump();
}
}
public void TestPrecision()
{
var file = new GribFile(Settings.GAUSS);
var msg = file.First();
Assert.AreEqual(msg["bitsPerValue"].AsInt(), 13);
msg.DecimalPrecision = 4;
Assert.AreEqual(msg["bitsPerValue"].AsInt(), 20);
}
public void TestPrecision()
{
using (var file = new GribFile(Settings.REG_GAUSSIAN_SURFACE_GRB2)) {
var msg = file.First();
Assert.AreEqual(msg["bitsPerValue"].AsInt(), 14);
msg.DecimalPrecision = 4;
Assert.AreEqual(msg["bitsPerValue"].AsInt(), 20);
}
}
public void TestSteographic()
{
using (GribFile file = new GribFile(Settings.STEREO)) {
GribMessage msg = file.First();
Assert.Greater(msg.ValuesCount, 1);
foreach (GeoSpatialValue gs in msg.GeoSpatialValues)
{
Assert.Greater(gs.Latitude, 15);
}
}
}
public void TestGetNativeType()
{
using (GribFile file = new GribFile(Settings.REG_LATLON_GRB1)) {
var msg = file.First();
Assert.AreEqual(msg["packingType"].NativeType, GribValueType.String);
Assert.AreEqual(msg["longitudeOfFirstGridPointInDegrees"].NativeType, GribValueType.Double);
Assert.AreEqual(msg["numberOfPointsAlongAParallel"].NativeType, GribValueType.Int);
Assert.AreEqual(msg["values"].NativeType, GribValueType.DoubleArray);
// TODO: test other types
}
}
public void TestGetGrib2()
{
double delta = 0.1d;
using (GribFile file = new GribFile(Settings.REDUCED_LATLON_GRB2))
{
var msg = file.First();
// "InDegrees" is a magic token that converts coordinate double values to degrees
// explicit degree conversion via key name
double latitudeOfFirstGridPointInDegrees = msg["latitudeOfFirstGridPoint"].AsDouble();
Assert.AreEqual(latitudeOfFirstGridPointInDegrees, 90, delta);
// degree conversion via accessor
double longitudeOfFirstGridPointInDegrees = msg["longitudeOfFirstGridPoint"].AsDouble();
Assert.AreEqual(longitudeOfFirstGridPointInDegrees, 0, delta);
// degree conversion via accessor
double latitudeOfLastGridPointInDegrees = msg["latitudeOfLastGridPoint"].AsDouble();
Assert.AreEqual(latitudeOfLastGridPointInDegrees, -90, delta);
// degree conversion via accessor
double longitudeOfLastGridPointInDegrees = msg["longitudeOfLastGridPoint"].AsDouble();
Assert.AreEqual(longitudeOfLastGridPointInDegrees, 360, .5);
// degree conversion via accessor
double jDirectionIncrementInDegrees = msg["jDirectionIncrement"].AsDouble();
Assert.AreEqual(jDirectionIncrementInDegrees, 0.36, delta);
// degree conversion via accessor
double iDirectionIncrementInDegrees = msg["iDirectionIncrement"].AsDouble();
Assert.AreEqual(iDirectionIncrementInDegrees, -1.0E+100, delta);
int numberOfPointsAlongAParallel = msg["numberOfPointsAlongAParallel"].AsInt();
Assert.AreEqual(numberOfPointsAlongAParallel, 0x7FFFFFFF);
int numberOfPointsAlongAMeridian = msg["numberOfPointsAlongAMeridian"].AsInt();
Assert.AreEqual(numberOfPointsAlongAMeridian, 501);
string packingType = msg["packingType"].AsString();
Assert.AreEqual("grid_simple", packingType);
}
}
public void TestEnumDisposal()
{
using (var file1 = new GribFile(Settings.COMPLEX_GRID))
using (var file2 = new GribFile(Settings.TIME)) {
var msg1 = file1.First();
var msg2 = file2.First();
var msgEnumerator1 = msg1.GetEnumerator();
var msgEnumerator2 = msg2.GetEnumerator();
int i = 0;
while (msgEnumerator1.MoveNext() && msgEnumerator2.MoveNext())
{
i++;
Assert.IsNotEmpty(msgEnumerator1.Current.Key);
Assert.IsNotEmpty(msgEnumerator2.Current.Key);
}
Assert.IsTrue(i > 3);
}
}
public void TestWrite()
{
if (File.Exists(Settings.OUT_GRIB))
{
File.Delete(Settings.OUT_GRIB);
}
int count = 0;
int valCount = 0;
using (var readFile = new GribFile(Settings.PACIFIC_WIND))
{
var msg = readFile.First();
Assert.AreNotEqual(33, msg["latitudeOfFirstGridPoint"].AsDouble());
msg["latitudeOfFirstGridPoint"].AsDouble(33);
valCount = msg.ValuesCount;
GribFile.Write(Settings.OUT_GRIB, msg);
}
using (var readFile = new GribFile(Settings.OUT_GRIB))
{
var msg = readFile.First();
count = readFile.MessageCount;
Assert.AreEqual(valCount, msg.ValuesCount);
Assert.AreEqual(count, readFile.MessageCount);
Assert.AreEqual(33, msg["latitudeOfFirstGridPoint"].AsDouble());
}
using (var readFile = new GribFile(Settings.PACIFIC_WIND))
{
GribFile.Write(Settings.OUT_GRIB, readFile as IEnumerable<GribMessage>, FileMode.Append);
count += readFile.MessageCount;
}
using (var readFile = new GribFile(Settings.OUT_GRIB))
{
Assert.AreEqual(count, readFile.MessageCount);
Assert.AreEqual(33, readFile.First()["latitudeOfFirstGridPoint"].AsDouble());
}
}
public void TestSetBitmapMissing ()
{
using (GribFile file = new GribFile(Settings.REG_LATLON_GRB1))
{
var msg = file.First();
int missing = 3333;
// set the value used to represent missing data
msg.MissingValue = missing;
Assert.AreEqual(msg.MissingValue, missing);
Assert.IsFalse(msg.HasBitmap);
int numVals = 10;
double[] vals = new double[numVals];
for (int i = 0; i < numVals; i++)
{
vals[i] = missing;
}
msg.HasBitmap = true;
msg.SetValues(vals);
double[] vals2;
msg.Values(out vals2);
for (int i = 0; i < numVals; i++)
{
Assert.AreEqual(vals[i], vals2[i]);
Assert.AreEqual(missing, vals2[i]);
}
missing = 9898;
msg.MissingValue = missing;
msg.Values(out vals);
for (int i = 0; i < numVals; i++)
{
Assert.AreEqual(missing, vals[i]);
}
}
}
public void TestWrite()
{
if (File.Exists(Settings.OUT_GRIB))
{
File.Delete(Settings.OUT_GRIB);
}
int count = 0;
int valCount = 0;
using (var readFile = new GribFile(Settings.PACIFIC_WIND))
{
var msg = readFile.First();
Assert.AreNotEqual(33, msg["latitudeOfFirstGridPoint"].AsDouble());
msg["latitudeOfFirstGridPoint"].AsDouble(33);
valCount = msg.ValuesCount;
GribFile.Write(Settings.OUT_GRIB, msg);
}
using (var readFile = new GribFile(Settings.OUT_GRIB))
{
var msg = readFile.First();
count = readFile.MessageCount;
Assert.AreEqual(valCount, msg.ValuesCount);
Assert.AreEqual(count, readFile.MessageCount);
Assert.AreEqual(33, msg["latitudeOfFirstGridPoint"].AsDouble());
}
using (var readFile = new GribFile(Settings.PACIFIC_WIND))
{
GribFile.Write(Settings.OUT_GRIB, readFile as IEnumerable <GribMessage>, FileMode.Append);
count += readFile.MessageCount;
}
using (var readFile = new GribFile(Settings.OUT_GRIB))
{
Assert.AreEqual(count, readFile.MessageCount);
Assert.AreEqual(33, readFile.First()["latitudeOfFirstGridPoint"].AsDouble());
}
}
public void TestSetBitmapMissing()
{
using (GribFile file = new GribFile(Settings.REG_LATLON_GRB1))
{
var msg = file.First();
int missing = 3333;
// set the value used to represent missing data
msg.MissingValue = missing;
Assert.AreEqual(msg.MissingValue, missing);
Assert.IsFalse(msg.HasBitmap);
int numVals = 10;
double[] vals = new double[numVals];
for (int i = 0; i < numVals; i++)
{
vals[i] = missing;
}
msg.HasBitmap = true;
msg.SetValues(vals);
double[] vals2;
msg.Values(out vals2);
for (int i = 0; i < numVals; i++)
{
Assert.AreEqual(vals[i], vals2[i]);
Assert.AreEqual(missing, vals2[i]);
}
missing = 9898;
msg.MissingValue = missing;
msg.Values(out vals);
for (int i = 0; i < numVals; i++)
{
Assert.AreEqual(missing, vals[i]);
}
}
}
public void TestCanConvertToDegress()
{
using (GribFile file = new GribFile(Settings.REDUCED_LATLON_GRB2)) {
var msg = file.First();
// true
// Assert.IsTrue(msg["latitudeOfFirstGridPointInDegrees"].CanConvertToDegrees);
Assert.IsTrue(msg["latitudeOfFirstGridPoint"].CanConvertToDegrees);
Assert.IsTrue(msg["longitudeOfFirstGridPointInDegrees"].CanConvertToDegrees);
Assert.IsTrue(msg["longitudeOfFirstGridPoint"].CanConvertToDegrees);
Assert.IsTrue(msg["latitudeOfLastGridPointInDegrees"].CanConvertToDegrees);
Assert.IsTrue(msg["latitudeOfLastGridPoint"].CanConvertToDegrees);
Assert.IsTrue(msg["jDirectionIncrement"].CanConvertToDegrees);
Assert.IsTrue(msg["iDirectionIncrement"].CanConvertToDegrees);
// false
Assert.IsFalse(msg["numberOfPointsAlongAParallel"].CanConvertToDegrees);
Assert.IsFalse(msg["numberOfPointsAlongAParallelInDegrees"].CanConvertToDegrees);
Assert.IsFalse(msg["numberOfPointsAlongAMeridian"].CanConvertToDegrees);
Assert.IsFalse(msg["numberOfPointsAlongAMeridianInDegrees"].CanConvertToDegrees);
Assert.IsFalse(msg["packingType"].CanConvertToDegrees);
}
}
public void TestOpenPng()
{
using (GribFile file = new GribFile(Settings.PNG_COMPRESSION)) {
Assert.IsTrue(file.MessageCount > 0);
var msg = file.First();
try {
Assert.IsTrue(msg["packingType"].AsString().ToLower().EndsWith("_png"));
Assert.IsTrue(msg.ValuesCount > 0);
Assert.IsTrue(msg.GeoSpatialValues.Any());
int i = 0;
foreach (var v in msg.GeoSpatialValues) {
Assert.AreNotEqual(Double.NaN, v.Value);
if (i++ > 1000) break;
}
} catch (GribApiException e) {
Console.WriteLine(e.Message);
Console.WriteLine(msg.ShortName);
Console.WriteLine(msg.ToString());
Assert.IsTrue(false);
}
}
}
public void TestCanConvertToDegress ()
{
using (GribFile file = new GribFile(Settings.REDUCED_LATLON_GRB2))
{
var msg = file.First();
// true
Assert.IsTrue(msg["latitudeOfFirstGridPointInDegrees"].CanConvertToDegrees);
Assert.IsTrue(msg["latitudeOfFirstGridPoint"].CanConvertToDegrees);
Assert.IsTrue(msg["longitudeOfFirstGridPointInDegrees"].CanConvertToDegrees);
Assert.IsTrue(msg["longitudeOfFirstGridPoint"].CanConvertToDegrees);
Assert.IsTrue(msg["latitudeOfLastGridPointInDegrees"].CanConvertToDegrees);
Assert.IsTrue(msg["latitudeOfLastGridPoint"].CanConvertToDegrees);
Assert.IsTrue(msg["jDirectionIncrement"].CanConvertToDegrees);
Assert.IsTrue(msg["iDirectionIncrement"].CanConvertToDegrees);
// false
Assert.IsFalse(msg["numberOfPointsAlongAParallel"].CanConvertToDegrees);
Assert.IsFalse(msg["numberOfPointsAlongAParallelInDegrees"].CanConvertToDegrees);
Assert.IsFalse(msg["numberOfPointsAlongAMeridian"].CanConvertToDegrees);
Assert.IsFalse(msg["numberOfPointsAlongAMeridianInDegrees"].CanConvertToDegrees);
Assert.IsFalse(msg["packingType"].CanConvertToDegrees);
}
}
public void TestCompression(string grid)
{
if (File.Exists(Settings.OUT_GRIB))
{
File.Delete(Settings.OUT_GRIB);
}
int count = 0;
double val = Double.NaN;
Dictionary <int, GeoSpatialValue> orig = new Dictionary <int, GeoSpatialValue>();
using (var readFile = new GribFile(Settings.REDUCED_LATLON_GRB2))
{
var msg = readFile.First();
count = msg.ValuesCount;
val = msg["latitudeOfFirstGridPoint"].AsDouble();
Assert.AreNotEqual(val, Double.NaN);
Assert.AreNotEqual(msg["packingType"].AsString(), grid);
{
for (int j = 0; j < 50; j++)
{
int k = 0;
do
{
// over 200k values, so only pick from a subset
k = GetRandomNumber(0, (count) / 20);
}while (orig.ContainsKey(k));
orig.Add(k, new GeoSpatialValue());
}
int x = 0;
int y = 0;
foreach (var gsv in msg.GeoSpatialValues)
{
if (orig.ContainsKey(x))
{
orig[x] = gsv;
y++;
if (y == orig.Count)
{
break;
}
}
x++;
}
Assert.IsTrue(y == orig.Count);
Assert.IsTrue(x > 0);
}
msg["packingType"].AsString(grid);
Assert.AreEqual(msg["packingType"].AsString(), grid);
Assert.AreEqual(count, msg.ValuesCount);
Assert.AreEqual(val, msg["latitudeOfFirstGridPoint"].AsDouble());
GribFile.Write(Settings.OUT_GRIB, msg);
var newGsv = msg.GeoSpatialValues.ToArray();
Assert.IsTrue(newGsv.Any());
foreach (var kvp in orig)
{
var nv = newGsv[kvp.Key];
Assert.AreEqual(kvp.Value.Value, nv.Value, 0.01);
}
}
using (var readFile = new GribFile(Settings.OUT_GRIB))
{
var msg = readFile.First();
Assert.AreEqual(count, msg.ValuesCount);
Assert.AreEqual(val, msg["latitudeOfFirstGridPoint"].AsDouble());
Assert.AreEqual(msg["packingType"].AsString(), grid);
var newGsv = msg.GeoSpatialValues.ToArray();
Assert.IsTrue(newGsv.Any());
foreach (var kvp in orig)
{
var nv = newGsv[kvp.Key];
Assert.AreEqual(kvp.Value.Value, nv.Value, 0.01);
}
}
}
public void TestPrecision()
{
var file = new GribFile(Settings.REG_GAUSSIAN_SURFACE_GRB2);
var msg = file.First();
Assert.AreEqual(msg["bitsPerValue"].AsInt(), 14);
msg.DecimalPrecision = 4;
Assert.AreEqual(msg["bitsPerValue"].AsInt(), 20);
}
public void TestEnumDisposal ()
{
using (var file1 = new GribFile(Settings.COMPLEX_GRID))
using (var file2 = new GribFile(Settings.TIME)) {
var msg1 = file1.First();
var msg2 = file2.First();
var msgEnumerator1 = msg1.GetEnumerator();
var msgEnumerator2 = msg2.GetEnumerator();
int i = 0;
while (msgEnumerator1.MoveNext() && msgEnumerator2.MoveNext()) {
i++;
Assert.IsNotEmpty(msgEnumerator1.Current.Key);
Assert.IsNotEmpty(msgEnumerator2.Current.Key);
}
Assert.IsTrue(i > 3);
}
}
public void TestGetNativeType()
{
using (GribFile file = new GribFile(Settings.REG_LATLON_GRB1))
{
var msg = file.First();
Assert.AreEqual(msg["packingType"].NativeType, GribValueType.String);
Assert.AreEqual(msg["longitudeOfFirstGridPointInDegrees"].NativeType, GribValueType.Double);
Assert.AreEqual(msg["numberOfPointsAlongAParallel"].NativeType, GribValueType.Int);
Assert.AreEqual(msg["values"].NativeType, GribValueType.DoubleArray);
// TODO: test other types
}
}
public void TestGetGrib2 ()
{
double delta = 0.1d;
using (GribFile file = new GribFile(Settings.REDUCED_LATLON_GRB2))
{
var msg = file.First();
// "InDegrees" is a magic token that converts coordinate double values to degrees
// explicit degree conversion via key name
double latitudeOfFirstGridPointInDegrees = msg["latitudeOfFirstGridPoint"].AsDouble();
Assert.AreEqual(latitudeOfFirstGridPointInDegrees, 90, delta);
// degree conversion via accessor
double longitudeOfFirstGridPointInDegrees = msg["longitudeOfFirstGridPoint"].AsDouble(/* inDegrees == true */);
Assert.AreEqual(longitudeOfFirstGridPointInDegrees, 0, delta);
// degree conversion via accessor
double latitudeOfLastGridPointInDegrees = msg["latitudeOfLastGridPoint"].AsDouble(/* inDegrees == true */);
Assert.AreEqual(latitudeOfLastGridPointInDegrees, -90, delta);
// degree conversion via accessor
double longitudeOfLastGridPointInDegrees = msg["longitudeOfLastGridPoint"].AsDouble(/* inDegrees == true */);
Assert.AreEqual(longitudeOfLastGridPointInDegrees, 360, .5);
// degree conversion via accessor
double jDirectionIncrementInDegrees = msg["jDirectionIncrement"].AsDouble(/* inDegrees == true */);
Assert.AreEqual(jDirectionIncrementInDegrees, 0.36, delta);
// degree conversion via accessor
double iDirectionIncrementInDegrees = msg["iDirectionIncrement"].AsDouble(/* inDegrees == true */);
Assert.AreEqual(iDirectionIncrementInDegrees, -1.0E+100, delta);
int numberOfPointsAlongAParallel = msg["numberOfPointsAlongAParallel"].AsInt();
Assert.AreEqual(numberOfPointsAlongAParallel, -1);
int numberOfPointsAlongAMeridian = msg["numberOfPointsAlongAMeridian"].AsInt();
Assert.AreEqual(numberOfPointsAlongAMeridian, 501);
string packingType = msg["packingType"].AsString();
Assert.AreEqual("grid_simple", packingType);
}
}
public void TestCompression (string grid)
{
if (File.Exists(Settings.OUT_GRIB))
{
File.Delete(Settings.OUT_GRIB);
}
int count = 0;
double val = Double.NaN;
Dictionary<int, GeoSpatialValue> orig = new Dictionary<int, GeoSpatialValue>();
using (var readFile = new GribFile(Settings.REDUCED_LATLON_GRB2))
{
var msg = readFile.First();
count = msg.ValuesCount;
val = msg["latitudeOfFirstGridPoint"].AsDouble();
Assert.AreNotEqual(val, Double.NaN);
Assert.AreNotEqual(msg["packingType"].AsString(), grid);
{
for (int j = 0; j < 50; j++)
{
int k = 0;
do
{
// over 200k values, so only pick from a subset
k = GetRandomNumber(0, (count) / 20);
}
while (orig.ContainsKey(k));
orig.Add(k, new GeoSpatialValue());
}
int x = 0;
int y = 0;
foreach(var gsv in msg.GeoSpatialValues)
{
if (orig.ContainsKey(x))
{
orig[x] = gsv;
y++;
if (y == orig.Count) break;
}
x++;
}
Assert.IsTrue(y == orig.Count);
Assert.IsTrue(x > 0);
}
msg["packingType"].AsString(grid);
Assert.AreEqual(msg["packingType"].AsString(), grid);
Assert.AreEqual(count, msg.ValuesCount);
Assert.AreEqual(val, msg["latitudeOfFirstGridPoint"].AsDouble());
GribFile.Write(Settings.OUT_GRIB, msg);
var newGsv = msg.GeoSpatialValues.ToArray();
Assert.IsTrue(newGsv.Any());
foreach (var kvp in orig)
{
var nv = newGsv[kvp.Key];
Assert.AreEqual(kvp.Value.Value, nv.Value, 0.01);
}
}
using (var readFile = new GribFile(Settings.OUT_GRIB))
{
var msg = readFile.First();
Assert.AreEqual(count, msg.ValuesCount);
Assert.AreEqual(val, msg["latitudeOfFirstGridPoint"].AsDouble());
Assert.AreEqual(msg["packingType"].AsString(), grid);
var newGsv = msg.GeoSpatialValues.ToArray();
Assert.IsTrue(newGsv.Any());
foreach (var kvp in orig)
{
var nv = newGsv[kvp.Key];
Assert.AreEqual(kvp.Value.Value, nv.Value, 0.01);
}
}
}