public MrSidFileClass(IRasterDataset dataset, string filename, IPolygon polygon)
{
_dataset = dataset;
_filename = filename;
FileInfo fi = new FileInfo(filename);
switch (fi.Extension.ToLower())
{
case ".sid":
_type = RasterType.sid;
//_reader = MrSidWrapper.LoadMrSIDReader(filename, ref _geoCoord);
break;
case ".jp2":
_type = RasterType.jp2;
//_reader = MrSidWrapper.LoadJP2Reader(filename, ref _geoCoord);
break;
default:
_type = RasterType.unknown;
break;
}
FileInfo fiPrj = new FileInfo(fi.FullName.Substring(0, fi.FullName.Length - fi.Extension.Length) + ".prj");
if (fiPrj.Exists)
{
StreamReader sr = new StreamReader(fiPrj.FullName);
string wkt = sr.ReadToEnd();
sr.Close();
_sRef = gView.Framework.Geometry.SpatialReference.FromWKT(wkt);
}
if (polygon == null)
{
if (!calcPolygon())
{
return;
}
}
else
{
_polygon = polygon;
}
_isValid = true;
}
public IRasterSet Create(string fileName, int xsize, int ysize, int bands, RasterType eType)
{
IRasterSet rasterSet = null;
var driverCount = OSGeo.GDAL.Gdal.GetDriverCount();
for (int i = 0; i < driverCount; i++)
{
using (var driver = OSGeo.GDAL.Gdal.GetDriver(i))
{
var dataset = driver.Create(fileName, xsize, ysize, bands, eType.ToRasterType(), null);
if (dataset != null)
{
rasterSet = GetRasterSet(fileName, dataset);
}
}
}
return(rasterSet);
}
public RasterDataDetails(RasterType type)
{
Type = type;
}
public RasterClass(IRasterDataset dataset, string filename, IPolygon polygon)
{
try
{
FileInfo fi = new FileInfo(filename);
_title = fi.Name;
_filename = filename;
_dataset = dataset;
OSGeo.GDAL.Gdal.AllRegister();
_gDS = OSGeo.GDAL.Gdal.Open(fi.FullName, 0);
if (_gDS == null && _gDS.RasterCount == 0)
{
_valid = false;
return;
}
_iWidth = _gDS.RasterXSize;
_iHeight = _gDS.RasterYSize;
switch (fi.Extension.ToLower())
{
case ".adf":
case ".gsd":
_type = RasterType.grid;
break;
//case ".jp2":
// _type = RasterType.wavelet;
// break;
}
using (OSGeo.GDAL.Band band = _gDS.GetRasterBand(1))
{
if (_gDS.RasterCount == 1)
{
if (band.DataType != OSGeo.GDAL.DataType.GDT_Byte)
{
_type = RasterType.grid;
}
}
band.GetMinimum(out _min, out _hasNoDataVal);
band.GetMaximum(out _max, out _hasNoDataVal);
band.GetNoDataValue(out _nodata, out _hasNoDataVal);
}
OSGeo.GDAL.Driver driver = _gDS.GetDriver();
double[] tfw = new double[6];
_gDS.GetGeoTransform(tfw);
string tfwFilename = fi.FullName.Substring(0, fi.FullName.Length - fi.Extension.Length);
switch (fi.Extension.ToLower())
{
case ".jpg":
case ".jpeg":
tfwFilename += ".jgw";
break;
case ".jp2":
tfwFilename += ".j2w";
break;
case ".tif":
case ".tiff":
tfwFilename += ".tfw";
break;
case ".ecw":
tfwFilename += ".eww";
break;
default:
break;
}
FileInfo tfwInfo = new FileInfo(tfwFilename);
_tfw = new TFWFile(tfw[0], tfw[3], tfw[1], tfw[2], tfw[4], tfw[5]);
if (tfwInfo.Exists)
{
_tfw.Filename = tfwFilename;
}
if (_tfw.X == 0.0 && _tfw.Y == 0.0 &&
Math.Abs(_tfw.dx_X) == 1.0 && _tfw.dx_Y == 0.0 &&
Math.Abs(_tfw.dy_Y) == 1.0 && _tfw.dy_X == 0.0 && driver != null)
{
if (tfwInfo.Exists)
{
_tfw = new TFWFile(tfwFilename);
}
else
{
_tfw.isValid = false;
}
}
else
{
// Bei dem Driver schein es nicht Pixelmitte sein, oder ist das bei GDAL generell
//if (driver.ShortName.ToLower() == "jp2openjpeg")
{
_tfw.X += (_tfw.dx_X / 2.0D + _tfw.dx_Y / 2.0D);
_tfw.Y += (_tfw.dy_X / 2.0D + _tfw.dy_Y / 2.0D);
}
}
FileInfo fiPrj = new FileInfo(fi.FullName.Substring(0, fi.FullName.Length - fi.Extension.Length) + ".prj");
if (fiPrj.Exists)
{
StreamReader sr = new StreamReader(fiPrj.FullName);
string wkt = sr.ReadToEnd();
sr.Close();
_sRef = gView.Framework.Geometry.SpatialReference.FromWKT(wkt);
}
else
{
fiPrj = new FileInfo(fi.FullName.Substring(0, fi.FullName.Length - fi.Extension.Length) + ".wkt");
if (fiPrj.Exists)
{
StreamReader sr = new StreamReader(fiPrj.FullName);
string wkt = sr.ReadToEnd();
sr.Close();
_sRef = gView.Framework.Geometry.SpatialReference.FromWKT(wkt);
}
}
if (polygon != null)
{
_polygon = polygon;
}
else
{
calcPolygon();
}
}
catch (Exception ex)
{
string errMsg = ex.Message;
_valid = false;
}
}
public RasterDataDetails(RasterType type)
{
Type = type;
}
