来自：http://www.pudn.com/downloads153/sourcecode/graph/detail672576.html
using System;
using System.Drawing;
using System.Collections;
using System.ComponentModel;
using System.Windows.Forms;
using System.Data;
using System.IO;
using System.Runtime.InteropServices;
using ESRI.ArcGIS.esriSystem;
using ESRI.ArcGIS.Carto;
using ESRI.ArcGIS.Controls;
using ESRI.ArcGIS.ADF;
using ESRI.ArcGIS.SystemUI;
using ESRI.ArcGIS.Geodatabase;
using ESRI.ArcGIS.Display;
using ESRI.ArcGIS.GeoAnalyst;
using ESRI.ArcGIS.DataSourcesRaster;
using ESRI.ArcGIS.Geometry;
using ESRI.ArcGIS.Output;//输出栅格
namespace RasterInterpolate
{
public sealed partial class MainForm : Form
{
#region class private members
private IMapControl3 m_mapControl = null;
private string m_mapDocumentName = string.Empty;
#endregion
#region class constructor
public MainForm()
{
InitializeComponent();
}
#endregion
string[] PtName;
//double[] Cx;
//double[] Cy;
double[] Ptx;
double[] Pty;
double[] Railfall;
double[] Ptz;
double[] lat;
double[] lon;
double[] di;//大圆距离
private void MainForm_Load(object sender, EventArgs e)
{
//get the MapControl
m_mapControl = (IMapControl3)axMapControl1.Object;
m_mapControl.AddShapeFile(Application.StartupPath + @"\shp", "sheet2.shp");
// m_mapControl.LoadMxFile(Application.StartupPath + @"\data\站点.mxd","","");
//disable the Save menu (since there is no document yet)
menuSaveDoc.Enabled = false;
dt=oledata.operate();
PtName = new string[dt.Rows.Count];
//Cx = new double[dt.Rows.Count];
//Cy = new double[dt.Rows.Count];
Ptx = new double[dt.Rows.Count];
Pty = new double[dt.Rows.Count];
Railfall = new double[dt.Rows.Count];
Ptz = new double[dt.Rows.Count];
lon = new double[dt.Rows.Count];
lat=new double[dt.Rows.Count];
di = new double[dt.Rows.Count];
Init();//数据赋值
}
//private int m_ColumnCount;
//private int m_RowCount;
//private int m_size;
#region Main Menu event handlers
private void menuNewDoc_Click(object sender, EventArgs e)
{
//execute New Document command
ICommand command = new CreateNewDocument();
command.OnCreate(m_mapControl.Object);
command.OnClick();
}
private void menuOpenDoc_Click(object sender, EventArgs e)
{
//execute Open Document command
ICommand command = new ControlsOpenDocCommandClass();
command.OnCreate(m_mapControl.Object);
command.OnClick();
}
private void menuSaveDoc_Click(object sender, EventArgs e)
{
//execute Save Document command
if (m_mapControl.CheckMxFile(m_mapDocumentName))
{
//create a new instance of a MapDocument
IMapDocument mapDoc = new MapDocumentClass();
mapDoc.Open(m_mapDocumentName, string.Empty);
//Make sure that the MapDocument is not readonly
if (mapDoc.get_IsReadOnly(m_mapDocumentName))
{
MessageBox.Show("Map document is read only!");
mapDoc.Close();
return;
}
//Replace its contents with the current map
mapDoc.ReplaceContents((IMxdContents)m_mapControl.Map);
//save the MapDocument in order to persist it
mapDoc.Save(mapDoc.UsesRelativePaths, false);
//close the MapDocument
mapDoc.Close();
}
}
private void menuSaveAs_Click(object sender, EventArgs e)
{
//execute SaveAs Document command
ICommand command = new ControlsSaveAsDocCommandClass();
command.OnCreate(m_mapControl.Object);
command.OnClick();
}
private void menuExitApp_Click(object sender, EventArgs e)
{
//exit the application
Application.Exit();
}
#endregion
//listen to MapReplaced evant in order to update the statusbar and the Save menu
private void axMapControl1_OnMapReplaced(object sender, IMapControlEvents2_OnMapReplacedEvent e)
{
//get the current document name from the MapControl
m_mapDocumentName = m_mapControl.DocumentFilename;
//if there is no MapDocument, diable the Save menu and clear the statusbar
if (m_mapDocumentName == string.Empty)
{
menuSaveDoc.Enabled = false;
statusBarXY.Text = string.Empty;
}
else
{
//enable the Save manu and write the doc name to the statusbar
menuSaveDoc.Enabled = true;
statusBarXY.Text = System.IO.Path.GetFileName(m_mapDocumentName);
}
}
private void axMapControl1_OnMouseMove(object sender, IMapControlEvents2_OnMouseMoveEvent e)
{
statusBarXY.Text = string.Format("{0}, {1} {2}", e.mapX.ToString("#######.##"), e.mapY.ToString("#######.##"), axMapControl1.MapUnits.ToString().Substring(4));
}
public void CreateContourListFromShapefile()//内插生成高程信息
{
//用反距离IDW插值生成的栅格图像
IInterpolationOp pInterpolationOp;
pInterpolationOp = new RasterInterpolationOpClass();
// Create the input point object
IFeatureClass pFeatureClass;
IFeatureLayer pFeaturelayer;
pFeaturelayer = this.m_mapControl.Map.get_Layer(0) as IFeatureLayer;//这里我的数据在第八层
pFeatureClass = pFeaturelayer.FeatureClass;
// Define the search radius
IRasterRadius pRadius;
pRadius = new RasterRadiusClass();
object maxDistance = Type.Missing;
pRadius.SetVariable(12, ref maxDistance);
//Create FeatureClassDescriptor using a value field
IFeatureClassDescriptor pFCDescriptor;
pFCDescriptor = new FeatureClassDescriptorClass();
pFCDescriptor.Create(pFeatureClass, null, "高程（米）");
//set { cellsize for output raster in the environment
object dCellSize = cellSize;
IRasterAnalysisEnvironment pEnv;
pEnv = new RasterAnalysis();
pEnv = pInterpolationOp as IRasterAnalysisEnvironment;
pEnv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref dCellSize);
object objectbarrier = Type.Missing;
//Perform the interpolation
IGeoDataset rasDataset;
rasDataset = pInterpolationOp.IDW((IGeoDataset)pFCDescriptor, 2, pRadius, ref objectbarrier);
//IRasterDataset rsDt;//这里就是调用Esri函数的地方 克里金 双线性 那些都可以的
//rsDt.CreateDefaultRaster();
IRaster pOutRaster;
pOutRaster = rasDataset as IRaster;
//Add output into ArcMap as a raster layer
IRasterLayer pOutRasLayer;
pOutRasLayer = new RasterLayerClass();
pOutRasLayer.CreateFromRaster(pOutRaster);
pOutRasLayer.Name = "Contour";
//m_RowCount = pOutRasLayer.RowCount;
//m_ColumnCount = pOutRasLayer.ColumnCount;
this.m_mapControl.AddLayer(pOutRasLayer, 0);//
UsingRasterClassifyColorRampRenderer();
IRawPixels rawPixels = null;
IRasterBandCollection rasterBandCollection;
IRasterProps rasterProps;
// QI for IRawPixels and IRasterProps
rasterBandCollection = (IRasterBandCollection)rasDataset;
rawPixels = (IRawPixels)rasterBandCollection.Item(0);
rasterProps = (IRasterProps)rawPixels;
// int dHeight = rasterProps.Height;//当前栅格数据集的行数
// int dWidth = rasterProps.Width; //当前栅格数据集的列数
// double dX = rasterProps.MeanCellSize().X; //栅格的宽度
// double dY = rasterProps.MeanCellSize().Y; //栅格的高度
IEnvelope extent = rasterProps.Extent; //当前栅格数据集的范围
//IPoint pt=new PointClass();
IPnt pixelBlockSize = null;
IPnt pixelBlockOrigin = null;
// Create pixelblock
pixelBlockOrigin = new DblPntClass();
pixelBlockOrigin.SetCoords(0, 0);
IPixelBlock3 pixelBlock3;
pixelBlockSize = new DblPntClass();
pixelBlockSize.SetCoords(rasterProps.Width, rasterProps.Height);
pixelBlock3 = rawPixels.CreatePixelBlock(pixelBlockSize) as IPixelBlock3;
// Read pixelblock
rawPixels.Read(pixelBlockOrigin, (IPixelBlock)pixelBlock3);
// Get pixeldata array
System.Array pixelData;
pixelData = (System.Array)pixelBlock3.get_PixelDataByRef(0);
// IRasterDataset RsDateset;
createFileRasterDataset("c:\\temp", RainfallValue, extent.LowerLeft, rasterProps.Width, rasterProps.Height, pixelData);
//IRasterLayer pRasterLayer = new RasterLayerClass();
//pRasterLayer.CreateFromDataset(RsDateset);
//ILayer pLayer;
//pLayer = pRasterLayer;
//pLayer.Name = "Result";
//axMapControl1.Map.AddLayer(pLayer);
//axMapControl1.ActiveView.Refresh();
UsingRasterClassifyColorRampRenderer();//这里是要创建刚另一个图层（结果） 在c：\temp
//IWorkspaceFactory pWSF;//保存栅格
//IWorkspace pWS;
//pWSF = new RasterWorkspaceFactory();
//pWS = pWSF.OpenFromFile("c:\\temp", 0);
//IRaster raster;
//raster = RsDateset.CreateDefaultRaster();
//ISaveAs isa;
//isa = raster as ISaveAs;
//isa.SaveAs("test2.tif",pWS, "TIFF");
}
public int cellSize;
public int classifyNum;
public string RainfallValue;
private void gIDS法ToolStripMenuItem_Click(object sender, EventArgs e)
{
FrmSet frmSet = new FrmSet();
if(DialogResult.OK==frmSet.ShowDialog())
{
classifyNum = frmSet.classifyNum;
cellSize = frmSet.cellSize;
RainfallValue=frmSet.RailfallValue;
if (RainfallValue=="降雨量50")
{
Cx = 0.0006; Cy = 0.0005; Cz = 0.4175;
}
else if (RainfallValue == "降雨量100")
{
Cx = 0.0007;
Cy = 0.0006;
Cz = 0.5199;
}
else if (RainfallValue == "降雨量500")
{
Cx = 0.001;
Cy = 0.0009;
Cz = 0.769;
}
else if (RainfallValue == "降雨量1000")
{
Cx = 0.0011;
Cy = 0.001;
Cz = 0.8789;
}
ResetRailfall();
CreateContourListFromShapefile();
}
}
public static IColor ConvertColorToIColor(Color color)
{
IColor pColor = new RgbColorClass();
pColor.RGB = (int)color.R | (int)(color.G) < 8 | (int)(color.B) < 16;
return pColor;
}
private void UsingRasterClassifyColorRampRenderer() //分级渲染
{
//get { Map
IMap pMap;
pMap = m_mapControl.Map;
//get { raster input from layer
IRasterLayer pRLayer;
pRLayer = (IRasterLayer)pMap.get_Layer(0);
IRaster pRaster;
pRaster = pRLayer.Raster;
//Create classfy renderer and QI RasterRenderer interface
IRasterClassifyColorRampRenderer pClassRen;
pClassRen = new RasterClassifyColorRampRenderer();
IRasterRenderer pRasRen;
pRasRen = pClassRen as IRasterRenderer;
//set { raster for the render and update
pRasRen.Raster = pRaster;
pClassRen.ClassCount = classifyNum;
pRasRen.Update();
//Create a color ramp to use
IAlgorithmicColorRamp pRamp;
pRamp = new AlgorithmicColorRamp();
Color clr = new Color();
clr = Color.FromArgb(255, 255, 128);
pRamp.FromColor = ConvertColorToIColor(clr);
clr = Color.FromArgb(135, 42, 8);
pRamp.ToColor = ConvertColorToIColor(clr);
pRamp.Size = 10;
bool flag = true;
pRamp.CreateRamp(out flag);
//Create symbol for the classes
IFillSymbol pFSymbol;
pFSymbol = new SimpleFillSymbol();
//loop through the classes and apply the color and label
int I;
for (I = 0; I = pClassRen.ClassCount - 1; I++)
{
pFSymbol.Color = pRamp.get_Color(I);// .Color(I);
pClassRen.set_Symbol(I, pFSymbol as ISymbol);
pClassRen.set_Label(I, "等级" + Convert.ToString(I + 1));
} // I
//Update the renderer and plug into layer
pRasRen.Update();
pRLayer.Renderer = (IRasterRenderer)pClassRen;
m_mapControl.ActiveView.Refresh();
// m_mapControl;
this.axTOCControl1.Update();
//Release memeory
pMap = null;
pRLayer = null;
pRaster = null;
pRasRen = null;
pClassRen = null;
pRamp = null;
pFSymbol = null;
}
private void rasterDataToolStripMenuItem_Click(object sender, EventArgs e)
{
ExportImage(this.m_mapControl.ActiveView);
}
private void ExportImage(IActiveView pActiveView)
{
SaveFileDialog pSaveFileDialog = new SaveFileDialog();
pSaveFileDialog.Filter = "JPEG(*.jpg)|*.jpg|AI(*.ai)|*.ai|BMP(*.BMP)|*.bmp|EMF(*.emf)|*.emf|GIF(*.gif)|*.gif|PDF(*.pdf)|*.pdf|PNG(*.png)|*.png|EPS(*.eps)|*.eps|SVG(*.svg)|*.svg|TIFF(*.tif)|*.tif";
pSaveFileDialog.Title = "输出地图";
pSaveFileDialog.RestoreDirectory = true;
pSaveFileDialog.FilterIndex = 1;
if (pSaveFileDialog.ShowDialog() != System.Windows.Forms.DialogResult.OK)
{
return;
}
string sFilePath = pSaveFileDialog.FileName;
string sPathName = System.IO.Path.GetDirectoryName(sFilePath);
string sFileName = System.IO.Path.GetFileNameWithoutExtension(sFilePath);
IExport pExporter = null;
switch (pSaveFileDialog.FilterIndex)
{
case 1:
pExporter = new ExportJPEGClass();
break;
case 2:
pExporter = new ExportBMPClass();
break;
case 3:
pExporter = new ExportEMFClass();
break;
case 4:
pExporter = new ExportGIFClass();
break;
case 5:
pExporter = new ExportAIClass();
break;
case 6:
pExporter = new ExportPDFClass();
break;
case 7:
pExporter = new ExportPNGClass();
break;
case 8:
pExporter = new ExportPSClass();
break;
case 9:
pExporter = new ExportSVGClass();
break;
case 10:
pExporter = new ExportTIFFClass();
break;
default:
MessageBox.Show("输出格式错误");
return;
}
IEnvelope pEnvelope = new EnvelopeClass();
ITrackCancel pTrackCancel = new CancelTrackerClass();
tagRECT ptagRECT;
ptagRECT.left = 0;
ptagRECT.top = 0;
ptagRECT.right = (int)pActiveView.Extent.Width;
ptagRECT.bottom = (int)pActiveView.Extent.Height;
int pResolution = (int)(pActiveView.ScreenDisplay.DisplayTransformation.Resolution);
pEnvelope.PutCoords(ptagRECT.left, ptagRECT.bottom, ptagRECT.right, ptagRECT.top);
pExporter.Resolution = pResolution;
pExporter.ExportFileName = sFilePath;
pExporter.PixelBounds = pEnvelope;
pActiveView.Output(pExporter.StartExporting(), pResolution, ref ptagRECT, pActiveView.Extent, pTrackCancel);
pExporter.FinishExporting();
}
// IRasterDataset
private void createFileRasterDataset(string directoryName, string fileName, IPoint originPoint, int rowCount, int columCount, System.Array CountourpixelVale)//创建输出结果图层
{
//The file must not be there. Otherwise, dataset can not be got.
if (File.Exists(directoryName + fileName))
{
MessageBox.Show("File Exist!");
IRasterLayer pRasterLayer = new RasterLayerClass();
pRasterLayer.CreateFromFilePath(directoryName + fileName);
ILayer pLayer;
pLayer = pRasterLayer;
pLayer.Name = "Result";
axMapControl1.Map.AddLayer(pLayer);
axMapControl1.ActiveView.Refresh();
return;// null;
}
System.IO.FileStream fs = File.Create("C:\\temp\\" + RainfallValue + "均方差.txt");
StreamWriter sw = new StreamWriter((System.IO.Stream)fs);
sw.WriteLine("站点名，x坐标，y坐标," + RainfallValue + "均方差");
// This function creates a new img file in the given workspace
// and then assigns pixel values
try
{
IRasterDataset rasterDataset = null;
// double cellSize = 500.0;//栅格大小xugy
// IPoint originPoint = new PointClass();
double minX = originPoint.X - cellSize / 2;
double minY = originPoint.Y - cellSize / 2;
//double minX=GetMinX();
//double minY=GetMinY();
//double maxX=GetMaxX();
//double maxY=GetMaxY();
double minLon = GetMinLon();
double minLat = GetMinLat();
double maxLon = GetMaxLon();
double maxLat = GetMaxLat();
// originPoint.PutCoords(minX - cellSize / 2, minY - cellSize / 2);
// Create the dataset
IRasterWorkspace2 rasterWorkspace2 = null;
// IGeographicCoordinateSystem m_GeographicCoordinateSystem;
ISpatialReference SpRef = new UnknownCoordinateSystem() as ISpatialReference;
// IProjectedCoordinateSystem m_ProjectedCoordinateSystem;
//ISpatialReferenceFactory2 spatRefFact = new SpatialReferenceEnvironmentClass();
//m_ProjectedCoordinateSystem = spatRefFact.CreateProjectedCoordinateSystem((int)esriSRProjCSType.esriSRProjCS_WGS1984UTM_49N);
rasterWorkspace2 = createRasterWorkspace(directoryName);
// int rowCount = (int)Math.Round((maxX - minX)/cellSize)+1;
// int columCount = (int)Math.Round((maxY - minY) / cellSize)+1;
//经纬度计算
double dLon = (maxLon-minLon)/rowCount;
double dlat=(maxLat-minLat)/columCount;
rasterDataset = rasterWorkspace2.CreateRasterDataset(fileName, "GRID", originPoint, rowCount, columCount, cellSize, cellSize, 1, rstPixelType.PT_FLOAT, SpRef, true);//"IMAGINE Image",rstPixelType.PT_UCHAR,new UnknownCoordinateSystemClass()
IRawPixels rawPixels = null;
IPixelBlock3 pixelBlock3 = null;
IPnt pixelBlockOrigin = null;
IPnt pixelBlockSize = null;
IRasterBandCollection rasterBandCollection;
IRasterProps rasterProps;
// QI for IRawPixels and IRasterProps
rasterBandCollection = (IRasterBandCollection)rasterDataset;
rawPixels = (IRawPixels)rasterBandCollection.Item(0);
rasterProps = (IRasterProps)rawPixels;
// Create pixelblock
pixelBlockOrigin = new DblPntClass();
pixelBlockOrigin.SetCoords(0, 0);
pixelBlockSize = new DblPntClass();
pixelBlockSize.SetCoords(rasterProps.Width, rasterProps.Height);
pixelBlock3 = rawPixels.CreatePixelBlock(pixelBlockSize) as IPixelBlock3;
// Read pixelblock
rawPixels.Read(pixelBlockOrigin, (IPixelBlock)pixelBlock3);
// Get pixeldata array
System.Array pixelData;
pixelData = (System.Array)pixelBlock3.get_PixelDataByRef(0);
// Loop through all the pixels and assign value
for (int i = 0; i rasterProps.Width; i++)//这个数组中就是输入的Value值
{
for (int j = 0; j rasterProps.Height; j++)
{
//tempPt.X = minX + i * cellSize;
// object oValue =CountourpixelVale.GetValue(0, 0);
// tempPt.Y = minY + j * cellSize;
Filldi(minLon + i * dLon, minLat + j * dlat);
float temp = (float)GIDS(minX + i * cellSize, minY + j * cellSize, Convert.ToDouble(CountourpixelVale.GetValue(i, j)));
////求均方差
//for (int k=0;k<dt.Rows.Count;k++)
//{
// MeanSquareDeviation += (temp - Railfall[k])*(temp - Railfall[k]);
//}
//MeanSquareDeviation = Math.Sqrt(MeanSquareDeviation / dt.Rows.Count);
////均方差写入文件
//sw.WriteLine(i + "," + j + "," + MeanSquareDeviation);
pixelData.SetValue(temp, i, j);//Convert.ToByte((i * j) % 255)
}
}
for (int k=0;k<dt.Rows.Count;k++)
{
int rr = Convert.ToInt32(Math.Ceiling((Ptx[k] - minX) / cellSize));
int cc=Convert.ToInt32(Math.Ceiling((Pty[k]-minY)/cellSize));
double pp;
if (rr rowCount && cc columCount)
{
pp = Convert.ToDouble(pixelData.GetValue(rr, cc));
sw.WriteLine(PtName[k]+","+Ptx[k] + "," + Pty[k] + "," + Math.Abs(pp - Railfall[k]));
}
}
sw.Close();
pixelBlock3.set_PixelData(0, (System.Object)pixelData);
// Write the pixeldata back
//System.Object cachePointer;
//cachePointer = rawPixels.AcquireCache();
rawPixels.Write(pixelBlockOrigin, (IPixelBlock)pixelBlock3);
//rawPixels.ReturnCache(cachePointer);
//IRasterLayer pRasterLayer = new RasterLayerClass();
//pRasterLayer.CreateFromDataset(rasterDataset);
//ILayer pLayer;
//pLayer = pRasterLayer;
//pLayer.Name = "Result";
//axMapControl1.Map.AddLayer(pLayer);
//axMapControl1.ActiveView.Refresh();
//解除引用
rawPixels = null;
rasterDataset= null;
pixelBlock3= null;
rasterBandCollection= null;
rasterProps = null;
GC.Collect();//强制回收
IRasterLayer m_raster = new RasterLayerClass();
m_raster.CreateFromFilePath("c:\\temp\\" + RainfallValue);
axMapControl1.AddLayer(m_raster);
//IRasterEdit rasterEdit; //保存
//rasterEdit = pRasterLayer.Raster as IRasterEdit;
//rasterEdit.Write(pixelBlockOrigin, pixelBlock3 as IPixelBlock);
//rasterEdit.Refresh();
//return rasterDataset;
}
catch (Exception ex)
{
sw.Close();
MessageBox.Show("My God, Error!" + ex.Message);
//System.Diagnostics.Debug.WriteLine(ex.Message);
return;// null;
}
}
double Cx = 0.0007;
double Cy = 0.0006;
double Cz = 0.5199;
private double GIDS(double x,double y,double z)//反梯度法求值
{
double diSum=0;
double tempSum=0;
for (int i = 0; i dt.Rows.Count; i++)
{
diSum += 1 / (di[i] * di[i]);
tempSum += (Railfall[i] + (x - Ptx[i]) * Cx + (y - Pty[i]) * Cy + Cz*(z - Ptz[i])) / (di[i] * di[i]);
}
tempSum = tempSum / diSum;
return tempSum;
}
private void Filldi(double Templon,double Templat)
{
for (int i = 0; i dt.Rows.Count; i++)
{
di[i] = GetDistance(Templon, Templat,lon[i],lat[i]);
}
}
public IRasterWorkspace2 createRasterWorkspace(string pathName)
{
// Create RasterWorkspace
IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactoryClass();
return workspaceFactory.OpenFromFile(pathName, 0) as IRasterWorkspace2;
}
OleData oledata = new OleData();
DataTable dt =new DataTable();
private void Init()//初始化复制，从数据库
{
for(int i = 0; i dt.Rows.Count; i++)
{
PtName[i]=dt.Rows[i][1].ToString();
//Cx[i] = Convert.ToDouble(dt.Rows[i][2]);//转换前X坐标
//Cy[i] = Convert.ToDouble(dt.Rows[i][3]);
Ptx[i] = Convert.ToDouble(dt.Rows[i][4]);
Pty[i] = Convert.ToDouble(dt.Rows[i][5]);
lon[i] = Convert.ToDouble(dt.Rows[i][6]);
lat[i] = Convert.ToDouble(dt.Rows[i][7]);
Ptz[i] = Convert.ToDouble(dt.Rows[i][8]);
Railfall[i] = Convert.ToDouble(dt.Rows[i][9]);//降雨量100
}
}
private void ResetRailfall()
{
for (int i = 0; i dt.Rows.Count; i++)
{
Railfall[i] = Convert.ToDouble(dt.Rows[i][RainfallValue]);//降雨量100
}
}
private double GetMinX()
{
double min = Ptx[0];
for (int i = 1; i dt.Rows.Count;i++ )
{
if (Ptx[i] min)
min = Ptx[i];
}
return min;
}
private double GetMinY()
{
double min = Pty[0];
for (int i = 1; i dt.Rows.Count; i++)
{
if (Pty[i] min)
min = Pty[i];
}
return min;
}
private double GetMaxX()
{
double max = Ptx[0];
for (int i = 1; i dt.Rows.Count; i++)
{
if (Ptx[i] > max)
max = Ptx[i];
}
return max;
}
private double GetMaxY()
{
double max = Pty[0];
for (int i = 1; i dt.Rows.Count; i++)
{
if (Pty[i] > max)
max = Pty[i];
}
return max;
}
private double GetMinLon()
{
double min = lon[0];
for (int i = 1; i dt.Rows.Count; i++)
{
if (lon[i] min)
min = lon[i];
}
return min;
}
private double GetMinLat()
{
double min = lat[0];
for (int i = 1; i dt.Rows.Count; i++)
{
if (lat[i] min)
min = lat[i];
}
return min;
}
private double GetMaxLon()
{
double max = lon[0];
for (int i = 1; i dt.Rows.Count; i++)
{
if (lon[i] > max)
max = lon[i];
}
return max;
}
private double GetMaxLat()
{
double max = lat[0];
for (int i = 1; i dt.Rows.Count; i++)
{
if (lat[i] > max)
max = lat[i];
}
return max;
}
const double PAI=3.1415926535;
const double EARTH_RADUS = 6371004;
double GetDistance(double a_dFirLon, double a_dFirLat, double a_dSecdLon, double a_dSecdLat)
{
double dRes = 0;
double dTemp = 0, dTemp2 = 0;
dTemp = Math.Cos(a_dSecdLat * PAI / 180) * Math.Cos(a_dFirLat * PAI / 180) * Math.Cos((a_dFirLon - a_dSecdLon) * PAI / 180) + Math.Sin(a_dSecdLat * PAI / 180) * Math.Sin(a_dFirLat * PAI / 180);
dTemp2 = EARTH_RADUS * Math.Acos(dTemp);
dRes = dTemp2;
return dRes;
}
}
}