private IFunctionRasterDataset calCov()
{
double N = (SumRws * SumClms * 2);
IFunctionRasterDataset mvRst = getShiftedRaster(InRaster);
IFunctionRasterDataset pRs = rsUtil.calcArithmaticFunction(InRaster, mvRst, esriRasterArithmeticOperation.esriRasterPlus);
IFunctionRasterDataset pRs2T = rsUtil.calcArithmaticFunction(InRaster, mvRst, esriRasterArithmeticOperation.esriRasterMultiply);
IFunctionRasterDataset pRs2 = rsUtil.calcArithmaticFunction(pRs2T, 2, esriRasterArithmeticOperation.esriRasterMultiply);
IFunctionRasterDataset mRsT = rsUtil.calcFocalStatisticsFunction(pRs, SumClms, SumRws, rasterUtil.focalType.SUM);
IFunctionRasterDataset mRsT2 = rsUtil.calcArithmaticFunction(mRsT, 2, esriRasterArithmeticOperation.esriRasterPower);
IFunctionRasterDataset mRs = rsUtil.calcArithmaticFunction(mRsT2, N, esriRasterArithmeticOperation.esriRasterDivide);
IFunctionRasterDataset mRs2 = rsUtil.calcFocalStatisticsFunction(pRs2, SumClms, SumRws, rasterUtil.focalType.SUM);
IFunctionRasterDataset dif = rsUtil.calcArithmaticFunction(mRs2, mRs, esriRasterArithmeticOperation.esriRasterMinus);
return(rsUtil.calcArithmaticFunction(dif, N, esriRasterArithmeticOperation.esriRasterDivide));
}
private void btnExecute_Click(object sender, EventArgs e)
{
string inRst1Nm = cmbRaster.Text;
string inFocalStat = cmbFocalStat.Text;
string windowType = cmbWindowType.Text;
string outNmRst = txtOutName.Text;
int clms = System.Convert.ToInt32(nudColumns.Value);
int rws = System.Convert.ToInt32(nudRows.Value);
if (inRst1Nm == "" || inRst1Nm == null)
{
MessageBox.Show("You must specify an input raster for In Raster 1 or type in a number", "Error", MessageBoxButtons.OK, MessageBoxIcon.Error);
return;
}
if (inFocalStat == "" || inFocalStat == null)
{
MessageBox.Show("You must specify an Focal Statistic", "Error", MessageBoxButtons.OK, MessageBoxIcon.Error);
return;
}
if (windowType == "" || windowType == null)
{
MessageBox.Show("You must select a widow type from dropdown menu", "Error", MessageBoxButtons.OK, MessageBoxIcon.Error);
return;
}
if (outNmRst == "" || outNmRst == null)
{
MessageBox.Show("You must specify an output raster name", "Error", MessageBoxButtons.OK, MessageBoxIcon.Error);
return;
}
IRaster rs1 = rstDic[inRst1Nm];
IRaster outRs = null;
rasterUtil.windowType wdType = (rasterUtil.windowType)Enum.Parse(typeof(rasterUtil.windowType), windowType);
rasterUtil.focalType fcType = (rasterUtil.focalType)Enum.Parse(typeof(rasterUtil.focalType), inFocalStat);
esriUtil.Forms.RunningProcess.frmRunningProcessDialog rp = new RunningProcess.frmRunningProcessDialog(false);
DateTime dt = DateTime.Now;
rp.addMessage("Creating Raster. This may take a while...");
rp.stepPGBar(10);
rp.TopMost = true;
rp.Show();
try
{
if (wdType == rasterUtil.windowType.CIRCLE)
{
outRs = rsUtil.returnRaster(rsUtil.calcFocalStatisticsFunction(rs1, clms, fcType));
}
else
{
outRs = rsUtil.returnRaster(rsUtil.calcFocalStatisticsFunction(rs1, clms, rws, fcType));
}
outraster = outRs;
outrastername = outNmRst;
if (mp != null && addToMap)
{
rp.addMessage("Calculating Statistics...");
rp.Show();
rp.Refresh();
IRasterLayer rsLyr = new RasterLayerClass();
//rsUtil.calcStatsAndHist(((IRaster2)outraster).RasterDataset);
rsLyr.CreateFromRaster(outraster);
rsLyr.Name = OutRasterName;
rsLyr.Visible = false;
mp.AddLayer(rsLyr);
}
this.DialogResult = DialogResult.OK;
}
catch (Exception ex)
{
rp.addMessage(ex.ToString());
}
finally
{
DateTime dt2 = DateTime.Now;
TimeSpan ts = dt2.Subtract(dt);
string t = " in " + ts.Days.ToString() + " days " + ts.Hours.ToString() + " hours " + ts.Minutes.ToString() + " minutes and " + ts.Seconds.ToString() + " seconds .";
rp.stepPGBar(100);
rp.addMessage("Finished Focal Analysis" + t);
rp.enableClose();
this.Close();
}
}