private Feature[] ReadFeatures(string shpFileName)
{
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(shpFileName) as IVectorFeatureDataReader)
{
return(dr.FetchFeatures());
}
}
public Feature[] GetFeatures(string fname)
{
if (fname == null)
{
return(null);
}
IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(fname) as IVectorFeatureDataReader;
if (dr == null)
{
return(null);
}
try
{
Feature[] fets = dr.FetchFeatures();
if (fets == null || fets.Length == 0)
{
return(null);
}
List <Feature> features = new List <Feature>();
foreach (Feature fet in fets)
{
if (fet.Geometry is ShapePolygon)
{
features.Add(fet);
}
}
return(features != null?features.ToArray() : null);
}
finally
{
dr.Dispose();
}
}
public string[] GetDisplayFieldNames()
{
if (_displayFieldName == null || _shpFileName == null || !File.Exists(_shpFileName))
{
return(null);
}
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(_shpFileName) as IVectorFeatureDataReader)
{
Feature[] fets = dr.FetchFeatures();
if (fets == null || fets.Length == 0)
{
return(null);
}
if (fets[0].FieldNames == null || !fets[0].FieldNames.Contains(_displayFieldName))
{
return(null);
}
List <string> names = new List <string>();
foreach (Feature f in fets)
{
names.Add(f.GetFieldValue(_displayFieldName));
}
return(names.Count > 0 ? names.ToArray() : null);
}
}
public float[] GetEvelope(string filename, string area)
{
List <Feature> fets = new List <Feature>();
IVectorFeatureDataReader dr = null;
int index = 0;
if (area == "中国区")
{
url = Path.Combine(AppDomain.CurrentDomain.BaseDirectory, "数据引用\\基础矢量\\行政区划\\面\\中国边界.shp");
index = 0;
}
else
{
url = Path.Combine(AppDomain.CurrentDomain.BaseDirectory, "数据引用\\基础矢量\\行政区划\\面\\省级行政区域_面.shp");
index = GetFeatureIndx(area);
}
dr = VectorDataReaderFactory.GetUniversalDataReader(url) as IVectorFeatureDataReader;
Feature[] temp = null;
temp = dr.FetchFeatures();
fets.AddRange(temp);
fets.ToArray();
Feature[] features = null;
features = fets.ToArray();
_blockName = area;
_vectorFeature = features[index];//_selectFeature
//string fieldValue = _vectorFeature.GetFieldValue(fieldIndex);
Envelope evp1 = new Envelope();
evp1 = GetMaskEnvelopeHDF(_vectorFeature.Geometry as ShapePolygon);
float[] envelope = new float[] { Convert.ToSingle(evp1.MinX), Convert.ToSingle(evp1.MaxX), Convert.ToSingle(evp1.MinY), Convert.ToSingle(evp1.MaxY) };
return(envelope);
}
public Feature[] LoadFeatures()
{
if (_urls == null)
{
return(null);
}
List <Feature> fets = new List <Feature>();
IVectorFeatureDataReader dr = null;
Feature[] temp = null;
foreach (string url in _urls)
{
dr = VectorDataReaderFactory.GetUniversalDataReader(url) as IVectorFeatureDataReader;
if (dr == null)
{
continue;
}
try
{
temp = dr.FetchFeatures();
if (temp == null || temp.Length == 0)
{
continue;
}
fets.AddRange(temp);
}
finally
{
dr.Dispose();
}
}
return(fets.Count == 0 ? null : fets.ToArray());
}
/// <summary>
/// 读取海冰控制点，获取海冰控制单信息
/// 最多6列，每列最多显示如下行数22,22,14,9,5,5=77，总共最多显示77个点，这样组织显示
/// </summary>
/// <param name="shpFile"></param>
private void ReadIceControlPointInfo(string shpFile)
{
_iceCptInfos.Clear();
List <string> infos = new List <string>();
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(shpFile) as IVectorFeatureDataReader)
{
if (dr == null)
{
return;
}
if (!dr.Fields.Contains("NO") || !dr.Fields.Contains("TEMP"))
{
return;
}
Feature[] features = dr.FetchFeatures();
if (features == null || features.Length == 0)
{
return;
}
_iceCptInfos.AppendLine(string.Format("{0}\t{1}\t{2}\t{3}\t", "编号", "经度", "纬度", "温度"));
foreach (Feature fet in features)
{
string iceNo = fet.GetFieldValue("NO");
string iceTemp = fet.GetFieldValue("TEMP");
ShapePoint sp = fet.Geometry as ShapePoint;
infos.Add(string.Format("{0}\t{1}\t{2}\t{3}\t", iceNo, sp.X.ToString("f2"), sp.Y.ToString("f2"), iceTemp));
}
}
int maxCol = 22;
int[] eachColCount = new int[] { 22, 21, 13, 7, 4, 4 };
int colCount = eachColCount.Length;
int infoCount = infos.Count;
for (int i = 0; i < maxCol && i < infoCount; i++)
{
int colBeginIndex = i;
StringBuilder curLine = new StringBuilder(infos[i]);
for (int j = 1; j < colCount; j++)
{
colBeginIndex += eachColCount[j - 1];
if (colBeginIndex >= infoCount)
{
break;
}
if (i < eachColCount[j])
{
curLine.AppendFormat(" {0}", infos[colBeginIndex]);
}
else
{
curLine.AppendFormat(" {0}", " ".PadRight(infos[colBeginIndex].Length));
}
}
_iceCptInfos.AppendLine(curLine.ToString());
}
}
public IFileExtractResult ClipSNWResult(string filename, string area)
{
inputFileName = filename;
using (IRasterDataProvider raster = GeoDataDriver.Open(filename) as IRasterDataProvider)
{
_resolutionX = raster.ResolutionX;
_resolutionY = raster.ResolutionY;
_fileSize = new Size(raster.Width, raster.Height);
_fileEnvelope = raster.CoordEnvelope;
_activeSpatialRef = raster.SpatialRef;
_activeCoordType = raster.CoordType;
raster.Dispose();
}
_outsizeRegion = new CoordEnvelope(_fileEnvelope.MinX, _fileEnvelope.MaxX, _fileEnvelope.MinY, _fileEnvelope.MaxY);
OutDir = Path.GetDirectoryName(filename);
string hdrfile = Path.Combine(Path.GetDirectoryName(filename), Path.GetFileNameWithoutExtension(filename) + ".hdr");
List <Feature> fets = new List <Feature>();
IVectorFeatureDataReader dr = null;
if (area == "中国区")
{
url = Path.Combine(AppDomain.CurrentDomain.BaseDirectory, "数据引用\\基础矢量\\行政区划\\面\\中国边界.shp");
}
else
{
url = Path.Combine(AppDomain.CurrentDomain.BaseDirectory, "数据引用\\基础矢量\\行政区划\\面\\省级行政区域_面.shp");
}
dr = VectorDataReaderFactory.GetUniversalDataReader(url) as IVectorFeatureDataReader;
Feature[] temp = null;
temp = dr.FetchFeatures();
fets.AddRange(temp);
Feature[] features = null;
features = fets.ToArray();
_blockName = area;
int index = GetFeatureIndx(area);
_vectorFeature = features[index];//_selectFeature
//string fieldValue = _vectorFeature.GetFieldValue(fieldIndex);
evp = GetMaskEnvelope(_vectorFeature.Geometry as ShapePolygon);
BlockDefWithAOI outEnvelope;
Size size;
List <BlockDefWithAOI> blockList = new List <BlockDefWithAOI>();
GetArgs(out outEnvelope, out size, out OutDir, out filename);
envelopes = new BlockDefWithAOI[] { outEnvelope };
aoiIndex = GetFeatureAOIIndex();
blockList.AddRange(envelopes);
blockList[0].AOIIndexes = aoiIndex;
string[] put = RasterClipT(inputFileName, blockList.ToArray(), OutDir, "Cut");
//put[0] = OutDir + put[0];
string getpath = Path.GetDirectoryName(put[0]);
IFileExtractResult res = new FileExtractResult("MWS", put[0], true);
res.SetDispaly(false);
return(res);
}
public PixelFeatures Export(Envelope geoEnvelope, Size rasterSize, string shpFile,
string fieldName, Func <Feature, double, bool> filter)
{
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(shpFile) as IVectorFeatureDataReader)
{
Feature[] fets = dr.FetchFeatures();
if (fets == null || !(fets[0].Geometry is ShapePoint))
{
return(null);
}
ShapePoint pt;
int row = 0, col = 0;
double spanY = geoEnvelope.Height / rasterSize.Height;
double spanX = geoEnvelope.Width / rasterSize.Width;
Dictionary <int, double> pixelFets = new Dictionary <int, double>();
double fetValue;
foreach (Feature fet in fets)
{
pt = fet.Geometry as ShapePoint;
row = (int)((geoEnvelope.MaxY - pt.Y) / spanY);
col = (int)((pt.X - geoEnvelope.MinX) / spanX);
if (row < 0 || row >= rasterSize.Height)
{
continue;
}
if (col < 0 || col >= rasterSize.Width)
{
continue;
}
fetValue = ToDouble(fet.GetFieldValue(fieldName));
int key = row * rasterSize.Width + col;
if (filter != null)
{
if (filter(fet, fetValue))
{
if (!pixelFets.ContainsKey(key))
{
pixelFets.Add(key, fetValue);
}
}
}
else
{
pixelFets.Add(key, fetValue);
}
}
if (pixelFets.Count == 0)
{
return(null);
}
PixelFeatures ret = new PixelFeatures();
ret.FeatureValues = pixelFets.Values.ToArray();
ret.RasterIndexes = pixelFets.Keys.ToArray();
return(ret);
}
}
public FetClassReaderForShpFile(string shpfilename)
{
_vdr = VectorDataReaderFactory.GetUniversalDataReader(shpfilename) as IVectorFeatureDataReader;
string annf = shpfilename.ToUpper().Replace(".SHP", "_注记.DBF");
if (File.Exists(annf))
{
_annoTable = Path.GetFileNameWithoutExtension(shpfilename) + "_注记";
}
}
public ShapePolygon[] GetGeometry(List <string> retNames, string nameField)
{
string fname = FindVector();
if (fname == null)
{
return(null);
}
IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(fname) as IVectorFeatureDataReader;
if (dr == null)
{
return(null);
}
try
{
Feature[] fets;
if (_where != null)
{
fets = dr.FetchFeatures((f) => { return(_where(f)); });
}
else
{
Feature f = dr.FetchFirstFeature();
if (f != null)
{
fets = new Feature[] { f }
}
;
else
{
return(null);
}
}
if (fets == null || fets.Length == 0)
{
return(null);
}
List <ShapePolygon> geometrys = new List <ShapePolygon>();
foreach (Feature fet in fets)
{
geometrys.Add(fet.Geometry as ShapePolygon);
if (retNames != null)
{
retNames.Add(fet.GetFieldValue(nameField) ?? string.Empty);
}
}
return(geometrys.ToArray());
}
finally
{
dr.Dispose();
}
}
public Feature GetFeatureByName(string name)
{
if (string.IsNullOrEmpty(name))
{
return(null);
}
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(_shpFileName) as IVectorFeatureDataReader)
{
return(dr.FetchFeature((fet) => { return fet.GetFieldValue(_displayFieldName) == name; }));
}
}
/// <summary>
/// 读取shpfile，获取范围信息
/// </summary>
/// <param name="shpFile"></param>
private void ReadShpFileAttribute(string shpFile)
{
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(shpFile) as IVectorFeatureDataReader)
{
if (dr == null)
{
return;
}
_shpEnvelope = dr.Envelope;
}
}
private void button15_Click(object sender, EventArgs e)
{
using (OpenFileDialog dlg = new OpenFileDialog())
{
dlg.Filter = "ESRI Shape Files(*.shp)|*.shp";
if (dlg.ShowDialog() == System.Windows.Forms.DialogResult.OK)
{
IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(dlg.FileName) as IVectorFeatureDataReader;
Feature fet = dr.FetchFirstFeature();
fet = dr.FetchFeature((f) => { return(f.GetFieldValue("CNTRY_NAME") == "China"); });
Vector2Bitmap(fet.Geometry as ShapePolygon);
}
}
}
private void LoadMaskTemplateFromShapeFile(string filename)
{
IVectorFeatureDataReader vdr = VectorDataReaderFactory.GetUniversalDataReader(filename) as IVectorFeatureDataReader;
if (vdr == null)
{
MessageBox.Show("输入的矢量文件\"" + filename + "\"格式错误,无法作为模板装入!");
return;
}
if (vdr.ShapeType != enumShapeType.Polygon)
{
MessageBox.Show("不支持非多边形矢量文件作为模板装入!");
return;
}
try
{
int fetCount = vdr.FeatureCount;
if (fetCount == 1)
{
SetMaskTemplateMatrixFromFeature(filename, vdr.Features[0]);
}
else
{
using (frmFeatureSelector frm = new frmFeatureSelector(vdr.Features))
{
//frm.TopMost = true;
if (frm.ShowDialog() == DialogResult.OK)
{
txtFileName.Text = filename;
Feature[] fets = frm.GetSelectedVectorFeatures();
if (fets == null || fets.Length == 0)
{
return;
}
for (int i = 0; i < fets.Length; i++)
{
SetMaskTemplateMatrixFromFeature(filename, fets[i]);
}
}
}
}
}
finally
{
if (vdr != null)
{
vdr.Dispose();
}
}
}
private void button21_Click(object sender, EventArgs e)
{
using (OpenFileDialog dlg = new OpenFileDialog())
{
dlg.Filter = "ESRI Shape Files(*.shp)|*.shp";
if (dlg.ShowDialog() == System.Windows.Forms.DialogResult.OK)
{
IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(dlg.FileName) as IVectorFeatureDataReader;
Feature fet = dr.FetchFirstFeature();
//fet = dr.FetchFeature((f) => { return f.GetFieldValue("CNTRY_NAME") == "China"; });
_aoiLayer.AddAOI(fet);
_canvas.Refresh(enumRefreshType.All);
}
}
}
private void button1_Click(object sender, EventArgs e)
{
Color c = Color.Transparent;
using (OpenFileDialog dlg = new OpenFileDialog())
{
dlg.Filter = "ESRI Shape Files(*.shp)|*.shp";
if (dlg.ShowDialog() == System.Windows.Forms.DialogResult.OK)
{
IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(dlg.FileName) as IVectorFeatureDataReader;
//Feature fet = dr.FetchFirstFeature();
//fet = dr.FetchFeature((f) => { return f.GetFieldValue("CNTRY_NAME") == "China"; });
v = dr.FetchFeatures();
}
}
}
public override void EndRead()
{
try
{
if (_reader != null)
{
_reader.Dispose();
_reader = null;
}
}
finally
{
_readIsFinished = true;
_isReady = false;
}
}
private void button7_Click(object sender, EventArgs e)
{
using (OpenFileDialog dlg = new OpenFileDialog())
{
dlg.Filter = "全球火点数据(*.hdf)|*.hdf";
if (dlg.ShowDialog() == System.Windows.Forms.DialogResult.OK)
{
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(dlg.FileName) as IVectorFeatureDataReader)
{
//dr.FeatureCount;
Feature[] features = dr.Features;
Console.WriteLine(features.Length);
}
}
}
}
private void LoadMaskTemplateFromShapeFile(string filename)
{
IVectorFeatureDataReader vdr = VectorDataReaderFactory.GetUniversalDataReader(filename) as IVectorFeatureDataReader;
if (vdr == null)
{
//MsgBox.ShowInfo("输入的矢量文件\"" + filename + "\"格式错误,无法作为模板装入!");
return;
}
if (vdr.ShapeType != enumShapeType.Polygon)
{
//MsgBox.ShowInfo("不支持非多边形矢量文件作为模板装入!");
return;
}
try
{
_maskObj = vdr.Features[0];
//int fetCount = vdr.FeatureCount;
//if (fetCount == 1)
//{
// SetMaskTemplateMatrixFromFeature(filename, vdr.FetchFirstFeature());
//}
//else
//{
// using (frmFeatureSelector frm = new frmFeatureSelector(vdr.Features))
// {
// if (frm.ShowDialog() == DialogResult.OK)
// {
// Feature[] fets = frm.GetSelectedVectorFeatures();
// if (fets == null || fets.Length == 0)
// return;
// for (int i = 0; i < fets.Length; i++)
// {
// SetMaskTemplateMatrixFromFeature(filename, fets[i]);
// }
// }
// }
//}
}
finally
{
if (vdr != null)
{
vdr.Dispose();
}
}
}
private void ImportFeatures(OdbcConnection conn, OdbcCommand cmd, IVectorFeatureDataReader vdr)
{
Feature[] features = vdr.Features;
if (features == null)
{
return;
}
string sql = "insert into " + Path.GetFileNameWithoutExtension(_dbffile) + "( ";
foreach (FieldMap map in _fieldMaps)
{
sql += (map.DesField.FieldName.ToString() + ",");
}
sql = sql.Substring(0, sql.Length - 1) + ") values ";
string fetsql = null;
foreach (Feature fet in features)
{
string status = fet.GetFieldValue("status").ToUpper();
if (status != null && status == "1")
{
continue;
}
fetsql = sql;
string values = "(";
foreach (FieldMap map in _fieldMaps)
{
if (map.ResField.FieldName.ToUpper() == "SHAPE.X")
{
values += ("'" + fet.Geometry.Centroid.X.ToString() + "',");
}
else if (map.ResField.FieldName.ToUpper() == "SHAPE.Y")
{
values += ("'" + fet.Geometry.Centroid.Y.ToString() + "',");
}
else
{
values += ("'" + fet.GetFieldValue(map.ResField.FieldName) + "',");
}
}
values = values.Substring(0, values.Length - 1) + ")";
fetsql += values;
//
cmd.CommandText = fetsql;
int ret = cmd.ExecuteNonQuery();
}
}
ProductColorTable _colorTable; //ISOT
private void CreateLegendItems(string shpFile)
{
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(shpFile) as IVectorFeatureDataReader)
{
if (dr == null)
{
return;
}
Feature[] features = dr.FetchFeatures();
if (features == null || features.Length == 0)
{
return;
}
string field = dr.Fields[0];//Contour
GetFeatureValueType(features, field);
}
}
public void ProcessVectorToRaster(string shpFileName, enumDataType dataType, double resolution, string rasterFileName)
{
if (string.IsNullOrEmpty(rasterFileName))
{
return;
}
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(shpFileName)
as IVectorFeatureDataReader)
{
if (dr == null)
{
return;
}
CoordEnvelope envelope = new CoordEnvelope(-180, 180, -90, 90);
int height = (int)Math.Ceiling((envelope.MaxY - envelope.MinY) / resolution);
int width = (int)Math.Ceiling((envelope.MaxX - envelope.MinX) / resolution);
IRasterDataProvider dataPrd = null;
try
{
if (Path.GetExtension(rasterFileName).ToUpper() == ".DAT")
{
IRasterDataDriver driver = GeoDataDriver.GetDriverByName("MEM") as IRasterDataDriver;
string mapInfo = envelope.ToMapInfoString(new Size(width, height));
dataPrd = driver.Create(rasterFileName, width, height, 1, dataType, mapInfo);
}
else
{
return;
}
Feature[] features = dr.FetchFeatures();
if (features == null || features.Length < 1)
{
return;
}
ProcessVectorToRaster(features, dataPrd);
}
finally
{
if (dataPrd != null)
{
dataPrd.Dispose();
}
}
}
}
private bool CheckVectorFile(string fileName)
{
using (IVectorFeatureDataReader vdr = VectorDataReaderFactory.GetUniversalDataReader(fileName) as IVectorFeatureDataReader)
{
if (vdr == null)
{
return(false);
}
if (vdr.ShapeType != enumShapeType.Polygon)
{
return(false);
}
else
{
return(true);
}
}
}
public static CodeCell.AgileMap.Core.IFeatureLayer CreateAndLoadVectorLayerForMicaps(ISmartSession session, ICanvas canvas, string fname, string dataTypeId, params object[] args)
{
if (string.IsNullOrEmpty(fname) || !File.Exists(fname))
{
return(null);
}
using (IVectorFeatureDataReader reader = MicapsDataReaderFactory.GetVectorFeatureDataReader(fname, dataTypeId))
{
if (reader == null)
{
return(null);
}
if (reader.Features != null)
{
MemoryDataSource mds = new MemoryDataSource(fname, reader.ShapeType, enumCoordinateType.Geographic, reader.Envelope, reader.Fields);
IFeatureClass fetClass = new FeatureClass(mds);
mds.AddFeatures(reader.Features);
CodeCell.AgileMap.Core.IFeatureLayer fetLayer = new FeatureLayer(fname, fetClass);
TryApplyStyle(fetLayer, dataTypeId);
IVectorHostLayer host = canvas.LayerContainer.VectorHost as IVectorHostLayer;
if (host != null)
{
host.Set(canvas);
IMapRuntime mapRuntime = host.MapRuntime as IMapRuntime;
if (mapRuntime != null)
{
IMap map = mapRuntime.Map as IMap;
if (map != null)
{
map.LayerContainer.Append(fetLayer);
FeatureLayer fetL = map.LayerContainer.Layers[0] as FeatureLayer;
FeatureClass fetC = fetL.Class as FeatureClass;
Envelope evp = fetC.FullEnvelope.Clone() as Envelope;
CoordEnvelope cvEvp = new CoordEnvelope(evp.MinX, evp.MaxX, evp.MinY, evp.MaxY);
canvas.CurrentEnvelope = cvEvp;
canvas.Refresh(enumRefreshType.All);
}
}
return(fetLayer);
}
}
}
return(null);
}
public override void BeginRead()
{
try
{
IUniversalVectorDataReader udr = VectorDataReaderFactory.GetUniversalDataReader(_fileUrl);
_reader = udr as IVectorFeatureDataReader;
_reader.SetArgsOfLeveling(_argOfLeveling);
_readIsFinished = false;
}
finally
{
if (GlobalCacher.VectorDataGlobalCacher != null &&
GlobalCacher.VectorDataGlobalCacher.IsEnabled &&
GlobalCacher.VectorDataGlobalCacher.IsGllPrjOfActiveViewer)
{
GlobalCacher.VectorDataGlobalCacher.Request(_fileUrl);
}
}
}
private void LoadFeatures(string shpFileName)
{
if (!File.Exists(shpFileName))
{
throw new FileNotFoundException(shpFileName);
}
//
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(shpFileName) as IVectorFeatureDataReader)
{
Feature[] fets = dr.FetchFeatures();
if (fets != null || fets.Length > 0)
{
_features.AddRange(fets);
}
}
//
TryLoadInterestRegions();
_shpFileName = shpFileName;
}
public Feature[] GetFeatures()
{
string fname = FindVector();
if (fname == null)
{
return(null);
}
IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(fname) as IVectorFeatureDataReader;
if (dr == null)
{
return(null);
}
try
{
Feature[] fets;
if (_where != null)
{
fets = dr.FetchFeatures((f) => { return(_where(f)); });
}
else
{
Feature f = dr.FetchFirstFeature();
if (f != null)
{
fets = new Feature[] { f }
}
;
else
{
return(null);
}
}
return(fets.Length == 0 ? null : fets.ToArray());
}
finally
{
dr.Dispose();
}
}
private object ReadFeatures(string filename)
{
if (string.IsNullOrEmpty(filename))
{
return(null);
}
IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(filename) as IVectorFeatureDataReader;
if (dr == null)
{
return(null);
}
try
{
return(dr.Features);
}
finally
{
dr.Dispose();
}
}
private void TryLoadInterestRegions()
{
if (_interestRegions == null || _interestRegions.Length == 0)
{
return;
}
foreach (string region in _interestRegions)
{
string fname = GetShpFileName(region);
if (fname != null)
{
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(fname) as IVectorFeatureDataReader)
{
Feature[] fets = dr.FetchFeatures();
if (fets != null || fets.Length > 0)
{
_waterBodyFeatures.AddRange(fets);
}
}
}
}
}
private IExtractResult GenerateExtractResult(string gfrFile, IRasterDataProvider dataPrd)
{
IPixelIndexMapper result = PixelIndexMapperFactory.CreatePixelIndexMapper("FIR", dataPrd.Width, dataPrd.Height, dataPrd.CoordEnvelope, dataPrd.SpatialRef);
Dataset _dataset = Gdal.Open(gfrFile, Access.GA_ReadOnly);
if (_dataset.RasterCount == 0)
{
return(result);
}
else
{
CoordEnvelope envelope = dataPrd.CoordEnvelope.Clone();
double maxX = envelope.MaxX;
double minX = envelope.MinX;
double maxY = envelope.MaxY;
double minY = envelope.MinY;
using (IVectorFeatureDataReader dr = VectorDataReaderFactory.GetUniversalDataReader(gfrFile) as IVectorFeatureDataReader)
{
Feature[] features = dr.FetchFeatures();
for (int i = 0; i < features.Length; i++)
{
double x, y;
if (double.TryParse(features[i].FieldValues[4], out x) && double.TryParse(features[i].FieldValues[3], out y))
{
if (IsInRange(minX, maxX, x) && IsInRange(minY, maxY, y))
{
int index = GetIndex(x, y);
if (index >= result.Count)
{
break;
}
result.Put(index);
}
}
}
}
return(result);
}
}
