public int[] GetBitmapIndexes()
{
int[] retAOI = null;
foreach (object obj in _aoiContainer.AOIs)
{
int[] aoi = null;
if (obj is GeometryOfDrawed)
{
aoi = GetBitmapIndexes(obj as GeometryOfDrawed);
}
else if (obj is Feature)
{
aoi = GetBitmapIndexes(obj as Feature);
}
//
if (aoi == null)
{
continue;
}
if (retAOI == null)
{
retAOI = aoi;
}
else
{
retAOI = AOIHelper.Merge(new int[][] { retAOI, aoi });
}
}
return(retAOI);
}
void canvasHost1_Load(object sender, EventArgs e)
{
canvasHost1.Canvas.Container.MouseMove += new MouseEventHandler(Container_MouseMove);
canvasHost1.Canvas.Container.MouseDown += new MouseEventHandler(Container_MouseDown);
canvasHost1.Canvas.OnEnvelopeChanged += new EventHandler((senser, arg) =>
{
if (_onCoordEnvelopeChanged != null)
{
_onCoordEnvelopeChanged(this, null);
}
});
canvasHost1.Canvas.Container.Paint += new PaintEventHandler(Container_Paint);
//
_initializer = new CanvasViewerInitializer(canvasHost1.Canvas, _session);
_aoiProvider = new AOIProvider(canvasHost1.Canvas);
_rasterInteractiver = new CurrentRasterInteractiver(this);
//
canvasHost1.Canvas.AOIGetter = () =>
{
return(_aoiProvider.GetBitmapIndexes());
};
canvasHost1.Canvas.MaskGetter = () =>
{
return(AOIHelper.Reverse(_aoiProvider.GetBitmapIndexes(), _aoiProvider.BitmapSize));
};
//
CreateVectorHost();
}
/// <summary>
/// 计算总覆盖度
/// 总覆盖度=(有蓝藻的像元面积)/(指定区域面积,例如太湖);
/// </summary>
/// <returns>覆盖度(0-1)</returns>
public float StatTotalConvertDegree(IRasterDataProvider dataProvider, int[] aoi)
{
int convertedCount = 0;
bool[] convertedPixels = new bool[aoi.Length];
int i = 0;
Size size = new Size(dataProvider.Width, dataProvider.Height);
Rectangle rect = AOIHelper.ComputeAOIRect(aoi, size);
ArgumentProvider ap = new ArgumentProvider(dataProvider, null);
RasterPixelsVisitor <float> visitor = new RasterPixelsVisitor <float>(ap);
visitor.VisitPixel(rect, aoi, new int[] { 1 }, (index, values) =>
{
if (values[0] >= 0f && values[0] <= 1f)
{
convertedPixels[i] = true;
i++;
}
});
foreach (bool isConverted in convertedPixels)
{
if (isConverted)
{
convertedCount++;
}
}
return((float)convertedCount / (float)convertedPixels.Length);
}
public void btnGetInfos_Click(object sender, EventArgs e)
{
if (_arp != null)
{
double lstBandRoom = (double)_arp.GetArg("lstband_Zoom");
IRasterDataProvider prd = _arp.DataProvider;
if (prd == null)
{
return;
}
IBandNameRaster bandNameRaster = prd as IBandNameRaster;
int lstbandNo = TryGetBandNo(bandNameRaster, "lstband");
if (lstbandNo == -1 || lstBandRoom == -1)
{
return;
}
ArgumentProvider ap = new ArgumentProvider(prd, null);
Size size = new Size(prd.Width, prd.Height);
if (_arp.AOIs != null && _arp.AOIs.Length != 0)
{
txtIdenfiy.Text = GetInfoExt(_arp.AOIs);
}
Rectangle rect = AOIHelper.ComputeAOIRect(_arp.AOI, size);
Dictionary <string, Dictionary <string, float> > anlysis = new Dictionary <string, Dictionary <string, float> >();
InitAnlysisInfo(anlysis);
int count = 0;
Int16[] cloudValues = GetNanValues("CloudyValue");
Int16[] waterValues = GetNanValues("WaterValue");
Int16[] invaildValues = GetNanValues("InvailValue");
using (RasterPixelsVisitor <Int16> visitor = new RasterPixelsVisitor <Int16>(ap))
{
visitor.VisitPixel(rect, _arp.AOI, new int[] { lstbandNo },
(index, values) =>
{
if (IsNanValue(values[0], cloudValues) || IsNanValue(values[0], waterValues) || IsNanValue(values[0], invaildValues))
{
return;
}
AnlysisValues(anlysis, "地表温度", values[0]);
count++;
});
}
anlysis["地表温度"]["avg"] = (float)Math.Round(anlysis["地表温度"]["avg"] / count, 4);
using (RasterPixelsVisitor <Int16> visitor = new RasterPixelsVisitor <Int16>(ap))
{
visitor.VisitPixel(rect, _arp.AOI, new int[] { lstbandNo },
(index, values) =>
{
anlysis["地表温度"]["pc"] += (float)Math.Pow(values[0] - anlysis["地表温度"]["avg"], 2);
});
}
anlysis["地表温度"]["pc"] = (float)Math.Round(Math.Sqrt(anlysis["地表温度"]["pc"]) / count, 4);
StringBuilder sb = new StringBuilder();
sb.Append(string.Format("地表温度:\n 最小值：{0}\n 最大值：{1}\n 平均值：{2}\n 偏差值：{3}\n", anlysis["地表温度"]["min"], anlysis["地表温度"]["max"], anlysis["地表温度"]["avg"], anlysis["地表温度"]["pc"]));
txtInfos.Text = sb.ToString();
}
}
private int[] GetAOI(IRasterDataProvider prd, bool isLand)
{
using (VectorAOITemplate v = VectorAOITemplateFactory.GetAOITemplate("海陆模版"))
{
Size size;
Envelope evp = GetEnvelope(out size, prd);
int[] landAoi = v.GetAOI(evp, size); //陆地
int[] reverseAOI = AOIHelper.Reverse(landAoi, size); //海洋
return(isLand ? landAoi : reverseAOI);
}
}
private int[] GetAOI()
{
VectorAOITemplate template = VectorAOITemplateFactory.GetAOITemplate("洞庭湖");
Size size;
IRasterDataProvider prd;
Envelope evp = GetEnvelope(out size, out prd);
aoi = template.GetAOI(evp, size);
//
int[] reverseAOI = AOIHelper.Reverse(aoi, size);
return(reverseAOI);
}
private IExtractResult FRILAlgorithm()
{
int VisibleCH = Obj2Int(_argumentProvider.GetArg("Visible"));
int NearInfraredCH = Obj2Int(_argumentProvider.GetArg("NearInfrared"));
int MiddleInfraredCH = Obj2Int(_argumentProvider.GetArg("MiddleInfrared"));
int FarInfraredCH = Obj2Int(_argumentProvider.GetArg("FarInfrared"));
string currRaster = Obj2String(_argumentProvider.GetArg("CurrentRasterFile"));
string currDBLV = Obj2String(_argumentProvider.GetArg("DBLVFile"));
if (VisibleCH == -1 || FarInfraredCH == -1 || NearInfraredCH == -1 || MiddleInfraredCH == -1)
{
PrintInfo("获取波段序号失败,可能是波段映射表配置错误。");
return(null);
}
if (string.IsNullOrEmpty(currRaster) || string.IsNullOrEmpty(currDBLV))
{
PrintInfo("获取算法所用文件失败。");
return(null);
}
List <int> vertifyIndexiex = null;
IPixelIndexMapper result = CreataPixelMapper(currDBLV, out vertifyIndexiex);
if (result == null || result.Indexes.ToArray().Length == 0)
{
PrintInfo("当前判识结果中无火点信息。");
return(null);
}
int[] filteredAOI = result.Indexes.ToArray();
Size size = new Size(_argumentProvider.DataProvider.Width, _argumentProvider.DataProvider.Height);
DoubtFirPixelFilter doubtFilter = CreateArgument.CreateDoubtFilter(_argumentProvider, _contextMessage);
if (doubtFilter == null)
{
return(null);
}
Rectangle aoiRect = AOIHelper.ComputeAOIRect(filteredAOI, size);
//背景温度计算
PrintInfo("[开始]背景亮温计算...");
BackTmpComputer backTmpComputer = CreateArgument.CreateBackTmpComputer(_argumentProvider, doubtFilter as IBackTmpComputerHelper, _contextMessage);
if (backTmpComputer == null)
{
return(null);
}
aoiRect = AOIHelper.ComputeAOIRect(filteredAOI, size);
Dictionary <int, PixelFeature> curFeatures = backTmpComputer.Compute(_argumentProvider, aoiRect, filteredAOI);
VertifyFirPixel(ref curFeatures, vertifyIndexiex, ref result);
return(GetOtherExtractResult.GetExtractResult(_argumentProvider, curFeatures, result, _contextMessage, _progressTracker));
}
private void button7_Click(object sender, EventArgs e)
{
PercentPixelFilterUInt16 u = new PercentPixelFilterUInt16(1, 0.04f, false);
string fname = @"f:\\FY3A_Mersi_2010_06_24_10_00_1000M_L1B_PRJ_Whole.LDF";
//fname = @"f:\\FY3A_MERSI_2010_06_24_10_00_1000M_L1B - 副本.HDF";
IRasterDataProvider prd = GeoDataDriver.Open(fname) as IRasterDataProvider;
int[] aoi = null;
Stopwatch sw = new Stopwatch();
sw.Start();
aoi = u.Filter(prd, AOIHelper.ComputeAOIRect(aoi, new Size(prd.Width, prd.Height)), aoi);
sw.Stop();
Text = sw.ElapsedMilliseconds.ToString();
}
private void button8_Click(object sender, EventArgs e)
{
string fname = "f:\\4_大昭寺_IMG_GE.tif";
//fname = "f:\\4600_3400.jpg";
//fname = "f:\\Penguins.jpg";
IRasterDataProvider prd = GeoDataDriver.Open(fname) as IRasterDataProvider;
int[] aoi = null;
PercentPixelFilterByte u = new PercentPixelFilterByte(1, 0.04f, false);
Stopwatch sw = new Stopwatch();
sw.Start();
aoi = u.Filter(prd, AOIHelper.ComputeAOIRect(aoi, new Size(prd.Width, prd.Height)), aoi);
sw.Stop();
Text = sw.ElapsedMilliseconds.ToString();
}
public void btnGetAOIIndex(object sender, EventArgs e)
{
this.AOIIndexs = _arp.AOI;
//根据选定的AOI区域进行自动计算水体NDVI最小值
IRasterDataProvider prd = _arp.DataProvider;
int bandNI = (int)_arp.GetArg("NearInfrared");
int bandVI = (int)_arp.GetArg("Visible");
ArgumentProvider ap = new ArgumentProvider(prd, null);
Size size = new Size(prd.Width, prd.Height);
Rectangle rect = AOIHelper.ComputeAOIRect(this.AOIIndexs, size);
using (RasterPixelsVisitor <UInt16> visitor = new RasterPixelsVisitor <UInt16>(ap))
{
List <double> ndvis = new List <double>();
visitor.VisitPixel(rect, this.AOIIndexs, new int[] { bandVI, bandNI },
(index, values) =>
{
if (values[1] + values[0] != 0) //是判断两个值都不等0的意思？
{
ndvis.Add(Math.Round((float)(values[1] - values[0]) / (values[1] + values[0]), 4));
}
});
if (ndvis.Count >= 0)
{
//设置界面值
float a = float.Parse(txta.Text);
float b = float.Parse(txtb.Text);
ndvis.Sort(); //从小到大排列
ndvis.RemoveRange(0, (int)(0.01 * ndvis.Count)); //去除百分之一最小值
float avgndvi = (float)ndvis.Min();
float MinNDVI = a * avgndvi + b;
this.txtNDVI.Text = avgndvi.ToString();
this.txtminndvi.Text = MinNDVI.ToString();
this.txtndvimin.Text = MinNDVI.ToString();
this.NDVIMultiBar.SetValues(new double[] { this.txtndvimin.Value, this.txtndvimax.Value });
}
}
this.ckbaoi.Checked = false;//重置状态
//设置完成之后自动生成
if (_handler != null)
{
_handler(GetArgumentValue());
return;
}
}
private int[] GetAOITemplates(IMonitoringSubProduct subProduct)
{
string[] aoiTemplates = _extractPanel.AOIs;
if (aoiTemplates == null || aoiTemplates.Length == 0)
{
return(null);
}
IRasterDataProvider dataprovider = GetCurrentRasterDataProvider();
if (dataprovider == null)
{
return(null);
}
CoordEnvelope evp = dataprovider.CoordEnvelope;
int[] aoi;
int[] retAOI = null;
foreach (string aoiname in aoiTemplates)
{
if (string.IsNullOrEmpty(aoiname))
{
continue;
}
aoi = AOITemplateFactory.MakeAOI(aoiname, evp.MinX, evp.MaxX, evp.MinY, evp.MaxY, new Size(dataprovider.Width, dataprovider.Height));
if (aoi == null)
{
continue;
}
if (subProduct != null)
{
if (aoiname.Contains("海陆模版") && subProduct.Identify == "ICE")
{
aoi = AOIHelper.Reverse(aoi, new Size(dataprovider.Width, dataprovider.Height));
}
}
if (retAOI == null)
{
retAOI = aoi;
}
else
{
retAOI = AOIHelper.Merge(new int[][] { aoi, retAOI });
}
}
return(retAOI);
}
public void btnGetAOIIndex(object sender, EventArgs e)
{
this.AOIIndexs = _arp.AOI;
//根据选定的AOI区域进行自动计算水体近红外反射率最小值
IRasterDataProvider prd = _arp.DataProvider;
int bandNI = (int)_arp.GetArg("NearInfrared");
ArgumentProvider ap = new ArgumentProvider(prd, null);
Size size = new Size(prd.Width, prd.Height);
Rectangle rect = AOIHelper.ComputeAOIRect(this.AOIIndexs, size);
using (RasterPixelsVisitor <UInt16> visitor = new RasterPixelsVisitor <UInt16>(ap))
{
List <double> nears = new List <double>();
visitor.VisitPixel(rect, this.AOIIndexs, new int[] { bandNI },
(index, values) =>
{
if (values[0] > 0) //近红外反射率>0
{
nears.Add(Math.Round(values[0] / 1000f, 4));
}
});
if (nears.Count >= 0)
{
//设置界面值
float pa = float.Parse(txtpa.Text);
float pb = float.Parse(txtpb.Text);
nears.Sort(); //从小到大排列
nears.RemoveRange(0, (int)(0.05 * nears.Count)); //去除百分之一最小值
float minnear = (float)nears.Min();
float MinN = minnear + (pa + pb) / 10;
this.txtNear.Text = minnear.ToString();
this.txtnearmin.Text = MinN.ToString();
this.NDVIMultiBar.SetValues(new double[] { MinN });
}
}
this.ckbone.Checked = false;//重置状态
//设置完成之后自动生成
if (_handler != null)
{
_handler(GetArgumentValue());
return;
}
}
internal Dictionary <int, FireAreaFeature> GetFireArea(IArgumentProvider argProvider, IPixelIndexMapper pixelMapper, Dictionary <int, PixelFeature> features)
{
RasterIdentify rid = new RasterIdentify();
Size size = new Size(argProvider.DataProvider.Width, argProvider.DataProvider.Height);
string iir = string.Empty;
using (InterestedRaster <int> result = new InterestedRaster <int>(rid, size, argProvider.DataProvider.CoordEnvelope))
{
foreach (int index in pixelMapper.Indexes)
{
result.Put(index, index);
}
iir = result.FileName;
}
List <int> aoiIndexs = new List <int>();
foreach (int keys in features.Keys)
{
if (features[keys].IsVertified)
{
aoiIndexs.Add(features[keys].PixelIndex);
}
}
int[] aoi = aoiIndexs.ToArray();
Rectangle aoiRect = AOIHelper.ComputeAOIRect(aoi, size);
IRasterDataProvider dataProvider = GeoDataDriver.Open(iir) as IRasterDataProvider;
PixelWndAccessor <int> wndVisitor = new PixelWndAccessor <int>(dataProvider);
wndVisitor.Reset(new int[] { 1 }, aoiRect);
int fAreaNum = -1;
Dictionary <int, FireAreaFeature> fireAreaFeatures = new Dictionary <int, FireAreaFeature>();
foreach (int index in aoi)
{
if (features[index].FireAreaNum == -1)
{
fAreaNum++;
features[index].FireAreaNum = fAreaNum;
fireAreaFeatures.Add(fAreaNum, GetFireAreaInfo(features[index], fAreaNum, argProvider.DataProvider));
GetFireAreaNum(wndVisitor, index, features, fAreaNum, fireAreaFeatures, argProvider.DataProvider);
}
}
return(fireAreaFeatures);
}
/// <summary>
/// 计算实际覆盖度
/// 实际覆盖度=(求和（每像元覆盖度*像元面积）)/(指定区域面积,例如太湖)
/// </summary>
/// <param name="dataProvider"></param>
/// <param name="aoi"></param>
/// <returns></returns>
public float CalcActualConvertDegree(IRasterDataProvider dataProvider, int[] aoi)
{
double actualArea = 0;
double totalArea = aoi.Count() * dataProvider.ResolutionX;
Size size = new Size(dataProvider.Width, dataProvider.Height);
Rectangle rect = AOIHelper.ComputeAOIRect(aoi, size);
ArgumentProvider ap = new ArgumentProvider(dataProvider, null);
RasterPixelsVisitor <float> visitor = new RasterPixelsVisitor <float>(ap);
//solution 使用横向？
visitor.VisitPixel(rect, aoi, new int[] { 1 }, (index, values) =>
{
if (values[0] <= 1f && values[0] >= 0f)
{
actualArea += values[0] * dataProvider.ResolutionX;
}
});
return((float)(actualArea / totalArea));
}
private void button16_Click(object sender, EventArgs e)
{
using (VectorAOITemplate v = VectorAOITemplateFactory.GetAOITemplate("贝尔湖")) //贝尔湖
{
Size size;
IRasterDataProvider prd;
Envelope evp = GetEnvelope(out size, out prd);
Stopwatch sw = new Stopwatch();
sw.Start();
aoi = v.GetAOI(evp, size);
//
int[] reverseAOI = AOIHelper.Reverse(aoi, size);
IBinaryBitmapBuilder b = new BinaryBitmapBuilder();
Bitmap bm = b.CreateBinaryBitmap(size, Color.Red, Color.Black);
b.Fill(reverseAOI, size, ref bm);
//
string express = "NDVI(band1,band2)";
int[] bandNos = new int[] { 1, 2 };
//构造栅格计算判识器
ArgumentProvider argprd = new ArgumentProvider(prd, null);
argprd.AOI = aoi;
IRasterExtracter <UInt16, float> extracter = new SimpleRasterExtracter <UInt16, float>();
extracter.Reset(argprd, bandNos, express);
//判识
resultNDVI = new MemPixelFeatureMapper <float>("NDVI", 1000, new Size(prd.Width, prd.Height), prd.CoordEnvelope, prd.SpatialRef);
extracter.Extract(resultNDVI);
//
RasterIdentify id = new RasterIdentify();
id.ThemeIdentify = "CMA";
id.ProductIdentify = "FIR";
id.SubProductIdentify = "NDVI";
id.Satellite = "FY3A";
id.Sensor = "MERSI";
id.Resolution = "250M";
id.OrbitDateTime = DateTime.Now.Subtract(new TimeSpan(1, 0, 0, 0, 0));
id.GenerateDateTime = DateTime.Now;
IInterestedRaster <float> iir = new InterestedRaster <float>(id, new Size(prd.Width, prd.Height), prd.CoordEnvelope.Clone());
iir.Put(resultNDVI);
iir.Dispose();
}
}
private int[] GetIndexes()
{
int[] retAOI = null;
foreach (object obj in _aoiContainer.AOIs)
{
int[] aoi = null;
aoi = GetIndexes(obj as Feature);
if (aoi == null)
{
continue;
}
if (retAOI == null)
{
retAOI = aoi;
}
else
{
retAOI = AOIHelper.Merge(new int[][] { retAOI, aoi });
}
}
return(retAOI);
}
public Feature Export(Envelope geoEnvelope, Size size, int[] aoi)
{
Dictionary <int, AOIHelper.AOIRowRange> rowRanges = AOIHelper.ComputeRowIndexRange(aoi, size);
if (rowRanges == null || rowRanges.Count == 0)
{
return(null);
}
int bCol = 0, eCol = 0;
int width = size.Width;
double top = geoEnvelope.MaxY;
double ySpan = geoEnvelope.Height / size.Height;
double left = geoEnvelope.MinX;
double xSpan = geoEnvelope.Width / size.Width;
List <ShapePoint> leftPoints = new List <ShapePoint>();
List <ShapePoint> retPoints = new List <ShapePoint>();
Rectangle aoiRect = AOIHelper.ComputeAOIRect(aoi, size);
double x = 0, y = 0;
int iRow = 0;
foreach (int row in rowRanges.Keys)
{
bCol = aoi[rowRanges[row].BeginIndex] % width;
eCol = aoi[rowRanges[row].EndIndex - 1] % width;
x = left + bCol * xSpan;
y = top - row * ySpan;
leftPoints.Add(new ShapePoint(x, y));
x = left + eCol * xSpan;
retPoints.Add(new ShapePoint(x, y)); //right point
iRow++;
}
leftPoints.Reverse();
retPoints.AddRange(leftPoints);
//
ShapePolygon ply = new ShapePolygon(new ShapeRing[] { new ShapeRing(retPoints.ToArray()) });
Feature fet = new Feature(0, ply, new string[] { "NAME" }, new string[] { "AOI" }, null);
return(fet);
}
/// <summary>
/// 计算指定等级的总实际覆盖面积
/// 比较条件为：minConvertDegree< x && x<= maxConvertDegree
/// </summary>
/// <param name="dataProvider">像元覆盖度文件DataProvider</param>
/// <param name="minConvertDegree"></param>
/// <param name="maxConvertDegree"></param>
/// <returns></returns>
public double CalcActualConvertArea(IRasterDataProvider dataProvider, int[] aoi, float minConvertDegree, float maxConvertDegree)
{
if (dataProvider == null || minConvertDegree > maxConvertDegree)
{
return(0);
}
double convertArea = 0;
double pixelArea = BAGStatisticHelper.CalPixelArea(dataProvider.ResolutionX);
ArgumentProvider ap = new ArgumentProvider(dataProvider, null);
Size size = new Size(dataProvider.Width, dataProvider.Height);
Rectangle rect = AOIHelper.ComputeAOIRect(aoi, size);
RasterPixelsVisitor <float> visitor = new RasterPixelsVisitor <float>(ap);
visitor.VisitPixel(rect, aoi, new int[] { 1 },
(index, values) =>
{
if (values[0] <= maxConvertDegree && values[0] > minConvertDegree)
{
convertArea += (pixelArea * values[0]);
}
});
return(convertArea);
}
private void SetAOIArgument(IMonitoringSubProduct subProduct)
{
int[] drawedAOI = null;
if (_session.SmartWindowManager.ActiveCanvasViewer != null && subProduct.Definition.IsNeedCurrentRaster)
{
drawedAOI = _session.SmartWindowManager.ActiveCanvasViewer.AOIProvider.GetIndexes();
}
int[] aoiFromTemplates = GetAOITemplates(subProduct);
//
ICanvasViewer activeCanvas = _session.SmartWindowManager.ActiveCanvasViewer;
if (activeCanvas != null)
{
IRasterDrawing drawing = _session.SmartWindowManager.ActiveCanvasViewer.ActiveObject as IRasterDrawing;
if (drawing != null)
{
IRasterDataProvider prd = drawing.DataProvider;
Size size = new System.Drawing.Size(prd.Width, prd.Height);
//
subProduct.ArgumentProvider.AOI = AOIHelper.Intersect(drawedAOI, aoiFromTemplates, size);
if (drawedAOI != null && drawedAOI.Length != 0)
{
List <Feature> fets = new List <Feature>();
AOI2ShapeFile aoi2ShpFile = new AOI2ShapeFile();
Feature fet = aoi2ShpFile.Export(GetEnvelope(drawing.DataProviderCopy),
new System.Drawing.Size(drawing.DataProviderCopy.Width, drawing.DataProviderCopy.Height),
drawedAOI);
fets.Add(fet);
subProduct.ArgumentProvider.AOIs = fets.ToArray();
}
else
{
subProduct.ArgumentProvider.AOIs = null;
}
}
}
}
private IPixelIndexMapper CalcBackTmp(IPixelIndexMapper result)
{
int[] filteredAOI = result.Indexes.ToArray();
if (filteredAOI == null || filteredAOI.Length == 0)
{
return(result);
}
int[] firePoints = GetPreExtractIndex();
Size size = new Size(_argumentProvider.DataProvider.Width, _argumentProvider.DataProvider.Height);
DoubtFirPixelFilter doubtFilter = CreateArgument.CreateDoubtFilter(_argumentProvider, _contextMessage);
if (doubtFilter == null)
{
return(result);
}
filteredAOI = MergeFirePoints(result, firePoints, filteredAOI);
Rectangle aoiRect = AOIHelper.ComputeAOIRect(filteredAOI, size);
//背景温度计算
PrintInfo("[开始]背景亮温计算...");
BackTmpComputer backTmpComputer = CreateArgument.CreateBackTmpComputer(_argumentProvider, doubtFilter as IBackTmpComputerHelper, _contextMessage);
if (backTmpComputer == null)
{
return(null);
}
aoiRect = AOIHelper.ComputeAOIRect(filteredAOI, size);
_curFeatures = backTmpComputer.ComputeNight(_argumentProvider, aoiRect, filteredAOI);
VertifyFirPixel(ref _curFeatures);
if (_curFeatures.Count != 0)
{
GetFirePointIndeiex.WriteFireIndexiexFilename(_argumentProvider, _curFeatures.Keys.ToArray());
}
return(result);
}
private void 海陆模板_Click(object sender, EventArgs e)
{
using (VectorAOITemplate v = VectorAOITemplateFactory.GetAOITemplate("海陆模版")) //贝尔湖
{
Size size;
IRasterDataProvider prd;
Envelope evp = GetEnvelope(out size, out prd);
Stopwatch sw = new Stopwatch();
sw.Start();
landAoi = v.GetAOI(evp, size);
//
int[] reverseAOI = AOIHelper.Reverse(landAoi, size);
//IBinaryBitmapBuilder b = new BinaryBitmapBuilder();
//Bitmap bm = b.CreateBinaryBitmap(size, Color.Red, Color.Black);
//b.Fill(reverseAOI, size, ref bm);
//对陆地区域使用陆地判识算法
Dictionary <string, object> args = new Dictionary <string, object>();
//args.Add("a", 28);
//args.Add("b", 78);
//args.Add("c", 245);
//args.Add("d", 293);
//args.Add("f", 28);
//args.Add("g", 0);
//args.Add("h", 15);
//args.Add("i", 250);
//构造判识表达式
//string express = "((band2/10f) > var_a) && (band2/10f < var_b) && (band5/10f > var_c) && (band5/10f < var_d) && "
// + "(band6/10f > var_f) && ((band6 - band2)>var_g) && ((band6/10f - band5/10f + var_i)>var_h)";
////band2:可见光，0.525~0.575（波长范围）
////band5:远红外，10.3~11.55
////band6:短波红外，1.60~1.69
//int[] bandNos = new int[] { 2, 5, 6 };
////构造栅格计算判识器
//ArgumentsProvider argprd = new ArgumentsProvider(prd, args);
//argprd.AOI = reverseAOI;
//IThresholdExtracter<UInt16> extracter = new SimpleThresholdExtracter<UInt16>();
//extracter.Reset(argprd, bandNos, express);
////判识
//landResult = new MemPixelIndexMapper("SAND", 1000);
//extracter.Extract(landResult);
//landInd = landResult.Indexes.ToArray();
//对海洋区域使用海洋判识算法
args.Clear();
args.Add("a", 8);
args.Add("b", 26);
args.Add("c", 265);
args.Add("d", 283);
args.Add("e", -10);
args.Add("f", 0);
string express = "(band2/10f > var_a) && (band2/10f < var_b) && (band5/10f > var_c) && (band5/10f < var_d) && "
+ "(band6/10f > var_a)&& ((band6/10f - band4/10f)>var_e) &&((band2/10f - band4/10f )>var_f)";
int[] bandNos = new int[] { 2, 5, 6, 4 };
ArgumentProvider argprd = new ArgumentProvider(prd, args);
argprd.AOI = reverseAOI;
IThresholdExtracter <UInt16> extracter = new SimpleThresholdExtracter <UInt16>();
extracter.Reset(argprd, bandNos, express);
//判识
seaResult = PixelIndexMapperFactory.CreatePixelIndexMapper("DST", prd.Width, prd.Height, prd.CoordEnvelope, prd.SpatialRef);
extracter.Extract(seaResult);
seaIdx = seaResult.Indexes.ToArray();
//
IBinaryBitmapBuilder builder = new BinaryBitmapBuilder();
Size bmSize = new Size(prd.Width / 2, prd.Height / 2);
Bitmap bitmap = builder.CreateBinaryBitmap(bmSize, Color.Red, Color.White);
builder.Fill(landInd, new Size(prd.Width, prd.Height), ref bitmap);
bitmap.Save("h:\\陆地沙尘.png", ImageFormat.Png);
builder.Fill(seaIdx, new Size(prd.Width, prd.Height), ref bitmap);
bitmap.Save("h:\\海洋沙尘.png", ImageFormat.Png);
//
RasterIdentify id = new RasterIdentify();
id.ThemeIdentify = "CMA";
id.ProductIdentify = "SAND";
id.SubProductIdentify = "2VAL";
id.Satellite = "FY3A";
id.Sensor = "MERSI";
id.Resolution = "1000M";
id.OrbitDateTime = DateTime.Now.Subtract(new TimeSpan(1, 0, 0, 0, 0));
id.GenerateDateTime = DateTime.Now;
IInterestedRaster <UInt16> iir = new InterestedRaster <UInt16>(id, new Size(prd.Width, prd.Height), prd.CoordEnvelope.Clone());
iir.Put(landInd, 1);
iir.Put(seaIdx, 1);
sw.Stop();
Text = sw.ElapsedMilliseconds.ToString();
iir.Dispose();
}
////取海洋矢量作为AOI，使用海洋的判识算法
//Dictionary<string, object> seaargs = new Dictionary<string, object>();
//seaargs.Add("sa", 8);
//seaargs.Add("sb", 26);
//seaargs.Add("sc", 265);
//seaargs.Add("sd", 283);
//seaargs.Add("se", -10);
//seaargs.Add("sf", 0);
//string seaExpress = "(band2/10f > var_sa) && (band2/10f < var_sb) && (band5/10f > var_sc) && (band5/10f < var_sd) && "
// + "(band6/10f > var_sa)&& ((band6/10f - band4/10f)>var_se) &&((band2/10f - band4/10f )>var_sf)";
//int[] seabandNos = new int[] { 2, 5, 6, 4 };
//ArgumentsProvider seaargprd = new ArgumentsProvider(prd, seaargs);
////构造判识表达式
////构造基于阈值的判识器
//IThresholdExtracter<UInt16> seaextracter = new SimpleThresholdExtracter<UInt16>();
//seaextracter.Reset(seaargprd, seabandNos, seaExpress);
////判识
//IPixelIndexMapper searesult = new MemPixelIndexMapper("SAND", 1000);
//seaextracter.Extract(searesult);
//int[] seaidxs = searesult.Indexes.ToArray();
//将两者判识的idex合并
}
public int[] Filter(IArgumentProvider argProvider, Dictionary <int, PixelFeature> features)
{
using (IRasterPixelsVisitor <UInt16> vistor = new RasterPixelsVisitor <UInt16>(argProvider))
{
int[] bandNos = new int[] { _visibleBandNo, _midInfraredBandNo, _farInfraredBandNo };
_visibleBandIndex = 0;
_midInfraredBandIndex = 1;
_farInfraredBandIndex = 2;
//
int[] aoi = GetAOI(features, (fet) => { return(fet.IsDoubtFirPixel); });
if (aoi == null)
{
return(null);
}
Rectangle aoiRect = AOIHelper.ComputeAOIRect(aoi, new System.Drawing.Size(argProvider.DataProvider.Width, argProvider.DataProvider.Height));
vistor.VisitPixel(aoiRect, aoi, bandNos, (idx, values) =>
{
float pvFarIfr = values[_farInfraredBandIndex] / _bandZoom;
float pvMidIfr = values[_midInfraredBandIndex] / _bandZoom;
float pvVis = values[_visibleBandIndex] / _bandZoom;
float pvMidIfr_FarIfr = values[_midInfraredBandIndex] - values[_farInfraredBandIndex];
//by chennan 20120821 一期火点确认错误代码
if (pvFarIfr < _maxFarIfrValue_fir || pvVis > _minVisValue_fir)
{
return;
}
PixelFeature fet = features[idx];
if (pvMidIfr > _secondMidIfrValue_fir)
{
//if (pvMidIfr > fet.MidIfrAvgValue + fet.MidIfr_StdDev)
if (pvMidIfr > fet.MidIfrAvgValue + _backTmpFactor_fir * fet.MidIfr_StdDev)
{
fet.IsVertified = true;
}
}
else
{
if (
(pvMidIfr > fet.MidIfrAvgValue + _backTmpFactor_fir * fet.MidIfr_StdDev)
&&
(pvMidIfr_FarIfr > fet.MidIfr_FarIfr_Diff_AvgValue + _backTmpFactor_fir * fet.MidIfr_FarIfr_Diff_StdDev)
)
{
fet.IsVertified = true;
}
}
//by chennan 20120821 增加中红外偏移量与远红外偏移量比值火点确认条件
if (fet.IsVertified == true && _isUsedMidFar)
{
float bzCH4 = (pvFarIfr - fet.FarIfrAvgValue);
if (bzCH4 <= 0)
{
bzCH4 = _defaultPercent;
}
//bzCH4 = 0.2f;
if ((pvMidIfr - fet.MidIfrAvgValue) / bzCH4 < _midFarRatio)
{
fet.IsVertified = false;
return;
}
}
});
}
return(GetAOI(features, (fet) => { return fet.IsVertified; }));
}
private IExtractResult CompareSHPFile(string backWaterPath, string binWater)
{
string shpPrimaryField = null;
Feature[] features = GetFeatures(backWaterPath, out shpPrimaryField);
if (features == null || features.Length == 0)
{
return(null);
}
//矢量栅格化
Dictionary <string, Color> nameColors = new Dictionary <string, Color>();
IRasterDataProvider dataPrd = GeoDataDriver.Open(binWater) as IRasterDataProvider;
IPixelFeatureMapper <Int16> result = new MemPixelFeatureMapper <Int16>("FLOD", 1000, new Size(dataPrd.Width, dataPrd.Height), dataPrd.CoordEnvelope, dataPrd.SpatialRef);
using (Bitmap bitmap = VectorsToBitmap(dataPrd, features, shpPrimaryField, out nameColors))
{
int[] aoi;
Color color;
string name;
List <int> items = new List <int>();
foreach (Feature fea in features)
{
name = fea.GetFieldValue(shpPrimaryField);
if (String.IsNullOrEmpty(name))
{
continue;
}
color = nameColors[name];
aoi = GetAOIByFeature(bitmap, color);
if (aoi != null)
{
items.AddRange(aoi);
}
}
//分情况处理
//无aoi区域
if (items.Count == 0)
{
ArgumentProvider ap = new ArgumentProvider(dataPrd, null);
RasterPixelsVisitor <Int16> rpVisitor = new RasterPixelsVisitor <Int16>(ap);
rpVisitor.VisitPixel(new int[] { 1 },
(idx, values) =>
{
if (values[0] == 1)
{
result.Put(idx, 1);
}
});
}
//有aoi区域
else
{
ArgumentProvider ap = new ArgumentProvider(dataPrd, null);
RasterPixelsVisitor <Int16> rpVisitor = new RasterPixelsVisitor <Int16>(ap);
rpVisitor.VisitPixel(new int[] { 1 },
(idx, values) =>
{
if (values[0] == 1)
{
result.Put(idx, 4);
}
});
Size size = new Size(dataPrd.Width, dataPrd.Height);
Rectangle aoiRect = AOIHelper.ComputeAOIRect(items.ToArray(), size);
rpVisitor.VisitPixel(aoiRect, items.ToArray(), new int[] { 1 },
(idx, values) =>
{
if (values[0] == 1)
{
result.Put(idx, 1);
}
else if (values[0] == 0)
{
result.Put(idx, 5);
}
});
}
}
return(result);
}
public void btnGetInfos_Click(object sender, EventArgs e)
{
if (_arp != null)
{
double visiBandRoom = (double)_arp.GetArg("Visible_Zoom");
double niBandRoom = (double)_arp.GetArg("NearInfrared_Zoom");
double shortInfraredRoom = (double)_arp.GetArg("ShortInfrared_Zoom");
IRasterDataProvider prd = _arp.DataProvider;
if (prd == null)
{
return;
}
IBandNameRaster bandNameRaster = prd as IBandNameRaster;
int visiBandNo = TryGetBandNo(bandNameRaster, "Visible");
int niBandNo = TryGetBandNo(bandNameRaster, "NearInfrared");
int shortInfraredBandNo = TryGetBandNo(bandNameRaster, "ShortInfrared");
if (visiBandNo == -1 || niBandNo == -1 || shortInfraredBandNo == -1 || shortInfraredRoom == -1 || visiBandRoom == -1 || niBandRoom == -1)
{
return;
}
ArgumentProvider ap = new ArgumentProvider(prd, null);
Size size = new Size(prd.Width, prd.Height);
if (_arp.AOIs != null && _arp.AOIs.Length != 0)
{
cmbType.Text = GetInfoExt(_arp.AOIs);
}
Rectangle rect = AOIHelper.ComputeAOIRect(_arp.AOI, size);
Dictionary <string, Dictionary <string, float> > anlysis = new Dictionary <string, Dictionary <string, float> >();
InitAnlysisInfo(anlysis);
int count = 0;
int ndviCount = 0;
using (RasterPixelsVisitor <UInt16> visitor = new RasterPixelsVisitor <UInt16>(ap))
{
visitor.VisitPixel(rect, _arp.AOI, new int[] { visiBandNo, niBandNo, shortInfraredBandNo },
(index, values) =>
{
AnlysisValues(anlysis, "可见光", values[0]);
AnlysisValues(anlysis, "近红外", values[1]);
AnlysisValues(anlysis, "短波红外", values[2]);
if (values[1] + values[0] != 0)
{
AnlysisValues(anlysis, "NDVI", (float)Math.Round((float)(values[1] - values[0]) / (values[1] + values[0]), 4));
ndviCount++;
}
count++;
});
}
anlysis["可见光"]["avg"] = (float)Math.Round(anlysis["可见光"]["avg"] / count, 4);
anlysis["近红外"]["avg"] = (float)Math.Round(anlysis["近红外"]["avg"] / count, 4);
anlysis["短波红外"]["avg"] = (float)Math.Round(anlysis["短波红外"]["avg"] / ndviCount, 4);
anlysis["NDVI"]["avg"] = (float)Math.Round(anlysis["NDVI"]["avg"] / ndviCount, 4);
using (RasterPixelsVisitor <UInt16> visitor = new RasterPixelsVisitor <UInt16>(ap))
{
visitor.VisitPixel(rect, _arp.AOI, new int[] { visiBandNo, niBandNo, shortInfraredBandNo },
(index, values) =>
{
anlysis["可见光"]["pc"] += (float)Math.Pow(values[0] - anlysis["可见光"]["avg"], 2);
anlysis["近红外"]["pc"] += (float)Math.Pow(values[1] - anlysis["近红外"]["avg"], 2);
anlysis["短波红外"]["pc"] += (float)Math.Pow(values[2] - anlysis["短波红外"]["avg"], 2);
if (values[1] + values[0] != 0)
{
anlysis["NDVI"]["pc"] += (float)Math.Pow((float)(values[1] - values[0]) / (values[1] + values[0]) - anlysis["NDVI"]["avg"], 2);
}
});
}
anlysis["可见光"]["pc"] = (float)Math.Round(Math.Sqrt(anlysis["可见光"]["pc"]) / count, 4);
anlysis["近红外"]["pc"] = (float)Math.Round(Math.Sqrt(anlysis["近红外"]["pc"]) / count, 4);
anlysis["短波红外"]["pc"] = (float)Math.Round(Math.Sqrt(anlysis["短波红外"]["pc"]) / count, 4);
anlysis["NDVI"]["pc"] = (float)Math.Round(Math.Sqrt(anlysis["NDVI"]["pc"]) / ndviCount, 4);
StringBuilder sb = new StringBuilder();
sb.Append(string.Format("可见光:\n 最小值：{0}\n 最大值：{1}\n 平均值：{2}\n 偏差值：{3}\n", anlysis["可见光"]["min"], anlysis["可见光"]["max"], anlysis["可见光"]["avg"], anlysis["可见光"]["pc"]));
sb.Append(string.Format("近红外:\n 最小值：{0}\n 最大值：{1}\n 平均值：{2}\n 偏差值：{3}\n", anlysis["近红外"]["min"], anlysis["近红外"]["max"], anlysis["近红外"]["avg"], anlysis["近红外"]["pc"]));
sb.Append(string.Format("NDVI:\n 最小值：{0}\n 最大值：{1}\n 平均值：{2}\n 偏差值：{3}\n", anlysis["NDVI"]["min"], anlysis["NDVI"]["max"], anlysis["NDVI"]["avg"], anlysis["NDVI"]["pc"]));
sb.Append(string.Format("短波红外:\n 最小值：{0}\n 最大值：{1}\n 平均值：{2}\n 偏差值：{3}\n", anlysis["短波红外"]["min"], anlysis["短波红外"]["max"], anlysis["短波红外"]["avg"], anlysis["短波红外"]["pc"]));
txtInfos.Text = sb.ToString();
InitWriteTxt(anlysis);
}
}
private IExtractResult ExtractLST()
{
float curNDVI = 0f;
IBandNameRaster bandNameRaster = _argumentProvider.DataProvider as IBandNameRaster; //
int VisibleCH = TryGetBandNo(bandNameRaster, "Visible"); //(int)_argumentProvider.GetArg("Visible");
int NearInfraredCH = TryGetBandNo(bandNameRaster, "NearInfrared"); //(int)_argumentProvider.GetArg("NearInfrared");
int MiddInfraredCH = TryGetBandNo(bandNameRaster, "MiddInfrared"); //(int)_argumentProvider.GetArg("MiddInfrared");
int FarInfrared11CH = TryGetBandNo(bandNameRaster, "FarInfrared11"); //(int)_argumentProvider.GetArg("FarInfrared11");
int FarInfrared12CH = TryGetBandNo(bandNameRaster, "FarInfrared12"); //(int)_argumentProvider.GetArg("FarInfrared12");
double VisibleZoom = (double)_argumentProvider.GetArg("Visible_Zoom");
double NearInfraredZoom = (double)_argumentProvider.GetArg("NearInfrared_Zoom");
double MiddInfraredZoom = (double)_argumentProvider.GetArg("MiddInfrared_Zoom");
double FarInfrared11Zoom = (double)_argumentProvider.GetArg("FarInfrared11_Zoom");
double FarInfrared12Zoom = (double)_argumentProvider.GetArg("FarInfrared12_Zoom");
if (VisibleCH == -1 || NearInfraredCH == -1 || MiddInfraredCH == -1 || FarInfrared11CH == -1 || FarInfrared12CH == -1)
{
PrintInfo("获取波段序号失败,可能是波段映射表配置错误或判识算法波段参数配置错误。");
return(null);
}
IRasterDataProvider currPrd = _argumentProvider.DataProvider;
if (currPrd == null)
{
PrintInfo("计算城市热岛的主文件未提供,请检查!");
return(null);
}
bool useCLMFile = true;
string clmRasterFile = _argumentProvider.GetArg("clmFile").ToString();
if (string.IsNullOrEmpty(clmRasterFile) || !File.Exists(clmRasterFile))
{
useCLMFile = false;
}
double lstZoom = (double)_argumentProvider.GetArg("LSTZoom");
UInt16 defCloudy = (UInt16)_argumentProvider.GetArg("defCloudy");
UInt16 defWater = (UInt16)_argumentProvider.GetArg("defWater");
float a0 = (float)_argumentProvider.GetArg("A0");
float alfa = (float)_argumentProvider.GetArg("alfa");
float beta = (float)_argumentProvider.GetArg("beta");
float gama = (float)_argumentProvider.GetArg("gama");
float alfa2 = (float)_argumentProvider.GetArg("alfa2");
float beta2 = (float)_argumentProvider.GetArg("beta2");
float b4m = (float)_argumentProvider.GetArg("b4m");
float b4n = (float)_argumentProvider.GetArg("b4n");
float b5m = (float)_argumentProvider.GetArg("b5m");
float b5n = (float)_argumentProvider.GetArg("b5n");
float b4Water = (float)_argumentProvider.GetArg("b4Water");
float b5Water = (float)_argumentProvider.GetArg("b5Water");
float b4Soil = (float)_argumentProvider.GetArg("b4Soil");
float b5Soil = (float)_argumentProvider.GetArg("b5Soil");
float b4VGT = (float)_argumentProvider.GetArg("b4VGT");
float b5VGT = (float)_argumentProvider.GetArg("b5VGT");
float NearInfraredCLMMin = (float)_argumentProvider.GetArg("NearInfraredCLMMin");
float FarInfrared11CLMMax = (float)_argumentProvider.GetArg("FarInfrared11CLMMax");
float FarInfrared1112CLMMin = (float)_argumentProvider.GetArg("FarInfrared1112CLMMin");
float FarInfrared11WaterMin = (float)_argumentProvider.GetArg("FarInfrared11WaterMin");
float NDVIWaterMax = (float)_argumentProvider.GetArg("NDVIWaterMax");
ArgumentProvider ap = new ArgumentProvider(currPrd, null);
RasterPixelsVisitor <UInt16> rpVisitor = new RasterPixelsVisitor <UInt16>(ap);
float b4Emiss = 0;
float b5Emiss = 0;
double pv = 0;
//4 5 通道比辐射率均值 差值
float EmissMean = 0;
float DeltaEmiss = 0;
//LST计算系数 P M
double P = 0;
double M = 0;
int[] bands = new int[] { VisibleCH, NearInfraredCH, MiddInfraredCH, FarInfrared11CH, FarInfrared12CH };
for (int i = 0; i < bands.Length; i++)
{
if (bands[i] > currPrd.BandCount)
{
PrintInfo("波段序号超过数据的波段范围，可能是影像数据选择不正确。");
return(null);
}
}
int[] aoi = _argumentProvider.AOI;
Size size = new Size(currPrd.Width, currPrd.Height);
Rectangle aoiRect = AOIHelper.ComputeAOIRect(aoi, size);
MemPixelFeatureMapper <UInt16> result = new MemPixelFeatureMapper <UInt16>("0LST", 1000, size, currPrd.CoordEnvelope, currPrd.SpatialRef);
rpVisitor.VisitPixel(aoiRect, aoi, bands,
(index, values) =>
{
if (values[0] == 0)
{
result.Put(index, 0);
}
else if (values[1] / NearInfraredZoom > NearInfraredCLMMin && values[3] / FarInfrared11Zoom <FarInfrared11CLMMax && Math.Abs(values[2] / MiddInfraredZoom - values[3] / FarInfrared11Zoom)> FarInfrared1112CLMMin)
{
result.Put(index, defCloudy);
}
else
{
curNDVI = GetNDVI(values, 1, 0);
if (values[3] / FarInfrared11Zoom > FarInfrared11WaterMin && curNDVI < NDVIWaterMax)
{
result.Put(index, defWater);
b4Emiss = b4Water;
b5Emiss = b5Water;
}
else
{
if (curNDVI < 0.2)
{
b4Emiss = b4Soil;
b5Emiss = b5Soil;
}
else if (curNDVI > 0.5)
{
b4Emiss = b4VGT;
b5Emiss = b5VGT;
}
else if (curNDVI >= 0.2 && curNDVI <= 0.5)
{
pv = Math.Pow((curNDVI - 0.2) / (0.5 - 0.2), 2);
b4Emiss = (float)(b4m * pv + b4n);
b5Emiss = (float)(b5m * pv + b5n);
}
EmissMean = (b4Emiss + b5Emiss) / 2;
DeltaEmiss = (b4Emiss - b5Emiss);
P = 1 + alfa * (1 - EmissMean) / EmissMean + beta * DeltaEmiss / Math.Pow(EmissMean, 2);
M = gama + alfa2 * (1 - EmissMean) / EmissMean + beta2 * DeltaEmiss / Math.Pow(EmissMean, 2);
result.Put(index, (UInt16)((a0 + P * (values[3] / FarInfrared11Zoom + values[4] / FarInfrared12Zoom) / 2
+ M * (values[3] / FarInfrared11Zoom - values[4] / FarInfrared12Zoom) / 2) * lstZoom));
}
}
});
return(GenrateIInterested(result, currPrd));
}
private void GetDBLVAOIIndex(string dblvFile, string cloudFile, List <int> dblvAOI)
{
Dictionary <string, FilePrdMap> filePrdMap = new Dictionary <string, FilePrdMap>();
filePrdMap.Add("dblvFile", new FilePrdMap(dblvFile, 1, new VaildPra(float.MinValue, float.MaxValue), new int[] { 1 }));
filePrdMap.Add("cloudFile", new FilePrdMap(cloudFile, 1, new VaildPra(float.MinValue, float.MaxValue), new int[] { 1 }));
ITryCreateVirtualPrd tryVPrd = new TryCreateVirtualPrdByMultiFile();
IVirtualRasterDataProvider vrd = tryVPrd.CreateVirtualRasterPRD(ref filePrdMap);
int minWnd = 3, maxWnd = 3;
if (vrd == null)
{
if (filePrdMap != null && filePrdMap.Count > 0)
{
foreach (FilePrdMap value in filePrdMap.Values)
{
if (value.Prd != null)
{
value.Prd.Dispose();
}
}
}
return;
}
try
{
List <int> temp = new List <int>();
temp.AddRange(dblvAOI);
dblvAOI.Clear();
ArgumentProvider ap = new ArgumentProvider(vrd, null);
Size size = new Size(vrd.Width, vrd.Height);
int dblvBand = filePrdMap["dblvFile"].StartBand;
int cloudBand = filePrdMap["cloudFile"].StartBand;
int[] tempAOI = temp.ToArray();
Rectangle aoiRect = AOIHelper.ComputeAOIRect(tempAOI, size);
int[] wndBandNos = new int[] { dblvBand, cloudBand };
bool isCloudEgde = false;
using (IRasterPixelsVisitor <float> visitor = new RasterPixelsVisitor <float>(ap))
{
visitor.VisitPixelWnd(aoiRect, tempAOI, wndBandNos, wndBandNos, minWnd, maxWnd,
(pixelIdx, crtWndSize, bandValues, wndValues) =>
{
return(false);
},
(pixelIdx, crtWndSize, bandValues, wndValues) =>
{
isCloudEgde = IsCloudEdgeMethod(bandValues, wndValues, 1);
if (isCloudEgde)
{
temp.Remove(pixelIdx);
}
});
}
dblvAOI.AddRange(temp);
temp = null;
}
finally
{
if (vrd != null)
{
vrd.Dispose();
}
}
}
private IExtractResult CalcMix(string dblvFile, string ldfFile, double bandZoom, int bandNo, double resultZoom, int[] dblvAOI, float waterNir, float LandNir, double zoom)
{
Dictionary <string, FilePrdMap> filePrdMap = new Dictionary <string, FilePrdMap>();
filePrdMap.Add("dblvFile", new FilePrdMap(dblvFile, 1, new VaildPra(float.MinValue, float.MaxValue), new int[] { 1 }));
filePrdMap.Add("ldfFile", new FilePrdMap(ldfFile, bandZoom, new VaildPra(float.MinValue, float.MaxValue), new int[] { bandNo }));
ITryCreateVirtualPrd tryVPrd = new TryCreateVirtualPrdByMultiFile();
IVirtualRasterDataProvider vrd = tryVPrd.CreateVirtualRasterPRD(ref filePrdMap);
int minWnd = 3, maxWnd = 3;
if (vrd == null)
{
if (filePrdMap != null && filePrdMap.Count > 0)
{
foreach (FilePrdMap value in filePrdMap.Values)
{
if (value.Prd != null)
{
value.Prd.Dispose();
}
}
}
return(null);
}
try
{
ArgumentProvider ap = new ArgumentProvider(vrd, null);
RasterPixelsVisitor <float> rpVisitor = new RasterPixelsVisitor <float>(ap);
Size size = new Size(vrd.Width, vrd.Height);
MemPixelFeatureMapper <Int16> result = new MemPixelFeatureMapper <Int16>("0MIX", 1000, size, vrd.CoordEnvelope, vrd.SpatialRef);
int dblvBand = filePrdMap["dblvFile"].StartBand;
int ldfBand = filePrdMap["ldfFile"].StartBand;
int[] wndBandNos = new int[] { dblvBand, ldfBand };
Rectangle aoiRect = AOIHelper.ComputeAOIRect(dblvAOI, size);
bool isEdge = false;
using (IRasterPixelsVisitor <float> visitor = new RasterPixelsVisitor <float>(ap))
{
visitor.VisitPixelWnd(aoiRect, dblvAOI, wndBandNos, wndBandNos, minWnd, maxWnd,
(pixelIdx, crtWndSize, bandValues, wndValues) =>
{
return(false);
},
(pixelIdx, crtWndSize, bandValues, wndValues) =>
{
isEdge = IsEdgeMethod(bandValues, wndValues, 0);
if (isEdge)
{
result.Put(pixelIdx, GetResult(bandValues, wndValues, waterNir, LandNir, zoom));
}
});
}
return(result);
}
finally
{
if (vrd != null)
{
vrd.Dispose();
}
if (filePrdMap != null && filePrdMap.Count > 0)
{
foreach (FilePrdMap value in filePrdMap.Values)
{
if (value.Prd != null)
{
value.Prd.Dispose();
}
}
}
}
}
private void DoFilters(IArgumentProvider argProvider, out Dictionary <int, PixelFeature> features, Action <int, string> progressTracker)
{
features = null;
ICandidatePixelFilter filter = null;
int[] aoi = argProvider.AOI;
Size size = new Size(argProvider.DataProvider.Width, argProvider.DataProvider.Height);
Rectangle aoiRect = AOIHelper.ComputeAOIRect(aoi, size);
int[] filteredAOI = null;
if (progressTracker != null)
{
progressTracker.Invoke(5, "正在过滤高温象元,请稍候...");
}
//4% 高温像元(中红外通道)
CreateArgument.PrintInfo(_contextMessage, "[开始]x%高温像元过滤...");
filter = CreateArgument.CreateHighTmpXPercentFilter(argProvider, _contextMessage);
if (filter == null)
{
return;
}
filteredAOI = filter.Filter(argProvider.DataProvider, aoiRect, aoi);
//by chennan 20130402 感兴趣内计算
UpdateFilteredAOIByAOI(ref filteredAOI, argProvider.AOI);
if (filteredAOI != null)
{
_waitingFirPixels.Put(filteredAOI);
CreateArgument.PrintInfo(_contextMessage, " x%高温点个数:" + filteredAOI.Length.ToString());
}
CreateArgument.PrintInfo(_contextMessage, "[完成]x%高温像元过滤。");
if (progressTracker != null)
{
progressTracker.Invoke(10, "正在提取陆地高温点,请稍候...");
}
//陆地高温点判识
CreateArgument.PrintInfo(_contextMessage, "[开始]陆地高温点过滤...");
filter = CreateArgument.CreateLandHighTmpFilter(argProvider, _contextMessage);
//by chennan 20130402 感兴趣内计算
UpdateFilteredAOIByAOI(ref filteredAOI, argProvider.AOI);
if (filter != null)
{
aoiRect = AOIHelper.ComputeAOIRect(filteredAOI, size);
filteredAOI = filter.Filter(argProvider.DataProvider, aoiRect, filteredAOI);
if (filteredAOI != null)
{
_candidateFirPixels.Put(filteredAOI);
_waitingFirPixels.Remove(filteredAOI);
CreateArgument.PrintInfo(_contextMessage, " 陆地高温点个数:" + filteredAOI.Length.ToString() + ",候选像元个数:" + _candidateFirPixels.Indexes.Count().ToString());
}
}
CreateArgument.PrintInfo(_contextMessage, "[完成]陆地高温点过滤。");
if (progressTracker != null)
{
progressTracker.Invoke(20, "正在提取异常高温点,请稍候...");
}
//异常高温点判识
CreateArgument.PrintInfo(_contextMessage, "[开始]邻域内异常高温点过滤...");
filter = CreateArgument.CreateAbnormalHighTmpFile(argProvider, _contextMessage);
//by chennan 20130402 感兴趣内计算
UpdateFilteredAOIByAOI(ref filteredAOI, argProvider.AOI);
if (filter != null)
{
filteredAOI = _waitingFirPixels.Indexes.ToArray();
aoiRect = AOIHelper.ComputeAOIRect(filteredAOI, size);
filteredAOI = filter.Filter(argProvider.DataProvider, aoiRect, filteredAOI);
//by chennan 20130402 感兴趣内计算
UpdateFilteredAOIByAOI(ref filteredAOI, argProvider.AOI);
if (filteredAOI != null)
{
_candidateFirPixels.Put(filteredAOI);
_waitingFirPixels.Remove(filteredAOI);
CreateArgument.PrintInfo(_contextMessage, " 异常高温点个数:" + filteredAOI.Length.ToString() + ",候选像元个数:" + _candidateFirPixels.Indexes.Count().ToString());
}
}
CreateArgument.PrintInfo(_contextMessage, "[完成]邻域内异常高温点过滤。");
if (progressTracker != null)
{
progressTracker.Invoke(30, "正在提取疑似高温点,请稍候...");
}
//疑似火点判识
CreateArgument.PrintInfo(_contextMessage, "[开始]疑似火点判识...");
DoubtFirPixelFilter doubtFilter = CreateArgument.CreateDoubtFilter(argProvider, _contextMessage);
if (doubtFilter == null)
{
return;
}
filteredAOI = _candidateFirPixels.Indexes.ToArray();
aoiRect = AOIHelper.ComputeAOIRect(filteredAOI, size);
filteredAOI = doubtFilter.Filter(argProvider, aoiRect, filteredAOI);
//by chennan 20130402 感兴趣内计算
UpdateFilteredAOIByAOI(ref filteredAOI, argProvider.AOI);
if (filteredAOI == null)
{
CreateArgument.PrintInfo(_contextMessage, " 没有监测出疑似火点");
_candidateFirPixels.Reset();
return;
}
_candidateFirPixels.Reset();
_candidateFirPixels.Put(filteredAOI);
CreateArgument.PrintInfo(_contextMessage, "[完成]疑似火点判识。");
CreateArgument.PrintInfo(_contextMessage, " 监测出疑似火点:" + filteredAOI.Length.ToString());
if (progressTracker != null)
{
progressTracker.Invoke(50, "正在计算背景温度,请稍候...");
}
//背景温度计算
CreateArgument.PrintInfo(_contextMessage, "[开始]背景亮温计算...");
BackTmpComputer backTmpComputer = CreateArgument.CreateBackTmpComputer(argProvider, doubtFilter as IBackTmpComputerHelper, _contextMessage);
if (backTmpComputer == null)
{
return;
}
aoiRect = AOIHelper.ComputeAOIRect(filteredAOI, size);
features = backTmpComputer.Compute(argProvider, aoiRect, filteredAOI);
//记录火点特征集合
_candidateFirPixels.Tag = new FirFeatureCollection("火点像元特征", features);
CreateArgument.PrintInfo(_contextMessage, " 为" + features.Count.ToString() + "个疑似火点计算了背景亮温");
CreateArgument.PrintInfo(_contextMessage, "[完成]背景亮温计算。");
if (progressTracker != null)
{
progressTracker.Invoke(70, "正在进行火点确认,请稍候...");
}
//火点确认
CreateArgument.PrintInfo(_contextMessage, "[开始]火点确认...");
VertifyFirPixelFiter vertifyFilter = CreateArgument.CreateVertifyFilter(argProvider, _contextMessage);
if (vertifyFilter == null)
{
return;
}
filteredAOI = vertifyFilter.Filter(argProvider, features);
_candidateFirPixels.Reset();
if (filteredAOI != null && filteredAOI.Length > 0)
{
_candidateFirPixels.Put(filteredAOI);
CreateArgument.PrintInfo(_contextMessage, " 确认火点数:" + filteredAOI.Length.ToString());
GetFirePointIndeiex.WriteFireIndexiexFilename(argProvider, features.Keys.ToArray());
}
CreateArgument.PrintInfo(_contextMessage, "[完成]火点确认。");
}
