/// <summary>
/// 创建图元
/// </summary>
/// <param name="kml">模型的kml</param>
/// <param name="layer">模型所在的图层</param>
/// <returns></returns>
public IMFElement CreateElement(Kml kml, ILayer layer)
{
KmlModel3d modelkml = kml.Placemark.Graph as KmlModel3d;
if (modelkml == null)
{
return(null);
}
if (modelkml.Position == null)
{
return(null);
}
if (!File.Exists(modelkml.ModelFilePath))
{
return(null);
}
int index = -1;
//图层
IGlobeGraphicsLayer graphicLayer = layer as IGlobeGraphicsLayer;
if (graphicLayer == null)
{
return(null);
}
//实例化图元
Model3d_ArcGlobe modelElement = null;
//添加图元
this.Dosomething((Action) delegate()
{
IImport3DFile import3Dfile = null;
if (!filePathDic.ContainsKey(modelkml.ModelFilePath))
{
import3Dfile = new Import3DFileClass();
import3Dfile.CreateFromFile(modelkml.ModelFilePath);
filePathDic.Add(modelkml.ModelFilePath, import3Dfile);
}
else//模型已创建
{
import3Dfile = filePathDic[modelkml.ModelFilePath];
}
modelElement = new Model3d_ArcGlobe(graphicLayer, modelkml, import3Dfile);
//设置属性
GlobeGraphicsElementPropertiesClass properties = new GlobeGraphicsElementPropertiesClass();
properties.Rasterize = modelkml.Rasterize;
graphicLayer.AddElement(modelElement, properties, out index);
modelElement.Index = index;
modelElement.ElementName = kml.Placemark.Name;
}, true);
return(modelElement);
}
public void LoadModel(string sFile)
{
IImport3DFile file = new Import3DFileClass();
file.CreateFromFile(sFile);
IGeometry geometry = file.Geometry;
this.m_pMarker3DSymbol.Shape = geometry;
if (System.IO.Path.GetExtension(sFile).ToLower() == ".wrl")
{
IVector3D pAxis = new Vector3DClass();
pAxis.SetComponents(1.0, 0.0, 0.0);
this.RotateGeometry(this.m_pMarker3DSymbol, pAxis, 90.0);
}
this.m_pMarker3DSymbol.UseMaterialDraping = true;
this.m_pMarker3DSymbol.RestrictAccessToShape();
}
public override void OnMouseDown(int Button, int Shift, int X, int Y)
{
// TODO: Add ToolAddCraterOnTIN.OnMouseDown implementation
if (pTinLayer != null)
{
ITin pTin = pTinLayer.Dataset;
ISurface pSurface = ((ITinAdvanced)pTin).Surface;
IMapControl2 pMapCtr = (((IToolbarControl)m_hookHelper.Hook).Buddy) as IMapControl2;
//在mapctr操作
if (pMapCtr != null)
{
IPoint mapPoint = pMapCtr.ToMapPoint(X, Y);
IZAware za = mapPoint as IZAware;
za.ZAware = true;
double zVal;
zVal = pSurface.GetElevation(mapPoint);
if (double.IsNaN(zVal))
{
MessageBox.Show("获取模型的高度失败！");
return;
}
if (m_listModels.Count == 0)
{
MessageBox.Show("请先生成模型！");
return;
}
mapPoint.Z = zVal;
try
{
//IMultiPatch pMP = new MultiPatchClass();
IFeatureClass pFC = pFeatureLayer.FeatureClass;
IFeature pF = pFC.CreateFeature();
IImport3DFile pI3D = new Import3DFileClass();
pI3D.CreateFromFile(m_listModels[m_listModels.Count - 1]);
IMultiPatch pMP = pI3D.Geometry as IMultiPatch;
ITransform3D pT3D = pMP as ITransform3D;
pT3D.Move3D(mapPoint.X, mapPoint.Y, mapPoint.Z);
pF.Shape = pMP as IGeometry;
pF.Store();
if (pMapCtr != null)
{
pMapCtr.ActiveView.PartialRefresh(esriViewDrawPhase.esriViewGeography, null, null);
}
}
catch (SystemException e)
{
MessageBox.Show(e.Message);
}
//IMultiPatch pMP = new MultiPatchClass();
//for (int i = 0; i < 5; i++)
//{
// IFeature pF = pFC.CreateFeature();
// IImport3DFile pI3D = new Import3DFileClass();
// pI3D.CreateFromFile(@"C:\Users\Administrator\Desktop\sampleAnalysis\1002.3ds");
// IMultiPatch pMP = pI3D.Geometry as IMultiPatch;
// // ITransform3D pT3D = pMP as ITransform3D;
// pF.Shape = pMP as IGeometry;
// ITransform3D pT3D = pF.Shape as ITransform3D;
// pT3D.Move3D(5293 + 1600, -20427 , 800 + i*400);
// // pF.Shape = pI3D.Geometry;
}
}
if (pRasterLayer != null)
{
//ITin pTin = pTinLayer.Dataset;
//ISurface pSurface = ((ITinAdvanced)pTin).Surface;
IMapControl2 pMapCtr = (((IToolbarControl)m_hookHelper.Hook).Buddy) as IMapControl2;
//在mapctr操作
if (pMapCtr != null)
{
IPoint mapPoint = pMapCtr.ToMapPoint(X, Y);
//添加手工编辑窗口
FrmManualSetModelPara pManualSetPara = null;
if (m_pModel == null)
{
Pt2d ptCurrent = new Pt2d();
ptCurrent.X = mapPoint.X;
ptCurrent.Y = mapPoint.Y;
pManualSetPara = new FrmManualSetModelPara(ptCurrent, "Crater");
}
else
{
m_pModel.x = mapPoint.X;
m_pModel.y = mapPoint.Y;
pManualSetPara = new FrmManualSetModelPara(m_pModel);
}
//FrmManualSetModelPara pManualSetPara = new FrmManualSetModelPara(ptCurrent, "Crater");
if (pManualSetPara.ShowDialog() == DialogResult.OK)
{
if (m_pModel == null)
{
m_pModel = new Model();
}
bool bFlag = m_pModel.compareModelPara(pManualSetPara.m_pModel);
if (!bFlag) //模型参数已经被修改，需要重新生成新的模型文件
{
m_pModel.copyFromModel(pManualSetPara.m_pModel);
//根据设置的参数生成相应的模型
TriType triType = TriType.TriForward;
LibModelGen.MappingType mappingType = LibModelGen.MappingType.Flat;
String szModelOutputFilename = System.IO.Path.GetTempFileName();
szModelOutputFilename = szModelOutputFilename.Substring(0, szModelOutputFilename.LastIndexOf('.')) + ".3ds";
ModelBase crater = new CraterGen(m_pModel.dbSize, m_pModel.dbDepth);
crater.OutputFilename = szModelOutputFilename;
crater.triype = triType;
crater.mappingType = mappingType;
//将撞击坑模型添加到相应图层
if (crater.generate())
{
m_listModels.Add(szModelOutputFilename);
}
}
//保存当前添加的模型参数
Model pTmpModel = new Model();
pTmpModel.copyFromModel(pManualSetPara.m_pModel);
m_manualAddModels.Add(pTmpModel);
}
else
{
return;
}
IZAware za = mapPoint as IZAware;
za.ZAware = true;
double zVal;
IRaster2 pRaster2 = pRasterLayer.Raster as IRaster2;
int col, row;
pRaster2.MapToPixel(mapPoint.X, mapPoint.Y, out col, out row);
zVal = Convert.ToDouble(pRaster2.GetPixelValue(0, col, row));
if (double.IsNaN(zVal))
{
MessageBox.Show("获取模型的高度失败！");
return;
}
if (m_listModels.Count == 0)
{
MessageBox.Show("请先生成模型！");
return;
}
mapPoint.Z = zVal;
IFeature pfeature = null;
try
{
//IMultiPatch pMP = new MultiPatchClass();
IFeatureClass pFC = pFeatureLayer.FeatureClass;
IFeature pF = pFC.CreateFeature();
IImport3DFile pI3D = new Import3DFileClass();
pI3D.CreateFromFile(m_listModels[m_listModels.Count - 1]);
IMultiPatch pMP = pI3D.Geometry as IMultiPatch;
ITransform3D pT3D = pMP as ITransform3D;
// pT3D.Move3D(mapPoint.X, mapPoint.Y, mapPoint.Z);
pT3D.Move3D(pManualSetPara.m_pModel.x, pManualSetPara.m_pModel.y, mapPoint.Z);
//IRasterProps pRasterProps = pRasterLayer.Raster as IRasterProps;
//double xmax = pRasterProps.Extent.XMax;
//double xmin = pRasterProps.Extent.XMin;
//double ymax = pRasterProps.Extent.YMax;
//double ymin = pRasterProps.Extent.YMin;
////生成地形的地理范围
//double dbGeoRangeX = xmax - xmin;
//double dbGeoRangeY = ymax - ymin;
//double dbModelRangeX = pMP.Envelope.Width;
//double dbModelRangeY = pMP.Envelope.Height;
//double dbModelRangeZ = pMP.Envelope.ZMax - pMP.Envelope.ZMin;
//根据地形大小改变撞击坑大小
//double dbRandomSize = Math.Min(dbGeoRangeX, dbGeoRangeY) * (0.15); //[0.03 0.06]
//double dbRandomSize = pRasterProps.MeanCellSize().X * 300;
//double dbScale = pManualSetPara.m_pModel.dbSize;// / Math.Max(dbModelRangeX, dbModelRangeY);
//pT3D.Scale3D(mapPoint, dbScale, dbScale, dbScale);
pF.Shape = pMP as IGeometry;
pfeature = pF as IFeature;
pF.Store();
if (pMapCtr != null)
{
pMapCtr.ActiveView.PartialRefresh(esriViewDrawPhase.esriViewGeography, null, null);
}
}
catch (SystemException e)
{
if (e.Message == "The spatial index grid size is invalid.")
{
IFeatureClassLoad pFCL = pFeatureLayer.FeatureClass as IFeatureClassLoad;
pFCL.LoadOnlyMode = true;
pfeature.Store();
pFCL.LoadOnlyMode = false;
}
else
{
MessageBox.Show(e.Message);
}
}
}
}
}
public override void OnMouseDown(int Button, int Shift, int X, int Y)
{
// TODO: Add ToolAddStone.OnMouseDown implementation
if (pTinLayer != null)
{
ITin pTin = pTinLayer.Dataset;
ISurface pSurface = ((ITinAdvanced)pTin).Surface;
IMapControl2 pMapCtr = (((IToolbarControl)m_hookHelper.Hook).Buddy) as IMapControl2;
//在mapctr操作
if (pMapCtr != null)
{
IPoint mapPoint = pMapCtr.ToMapPoint(X, Y);
IZAware za = mapPoint as IZAware;
za.ZAware = true;
double zVal;
zVal = pSurface.GetElevation(mapPoint);
if (double.IsNaN(zVal))
{
MessageBox.Show("获取模型的高度失败！");
return;
}
if (m_listModels.Count == 0)
{
MessageBox.Show("请先生成模型！");
return;
}
mapPoint.Z = zVal;
IFeature pfeature = null;
try
{
//IMultiPatch pMP = new MultiPatchClass();
IFeatureClass pFC = pFeatureLayer.FeatureClass;
IFeature pF = pFC.CreateFeature();
IImport3DFile pI3D = new Import3DFileClass();
pI3D.CreateFromFile(m_listModels[m_listModels.Count - 1]);
IMultiPatch pMP = pI3D.Geometry as IMultiPatch;
ITransform3D pT3D = pMP as ITransform3D;
pT3D.Move3D(mapPoint.X, mapPoint.Y, mapPoint.Z);
pF.Shape = pMP as IGeometry;
pfeature = pF as IFeature;
pF.Store();
if (pMapCtr != null)
{
pMapCtr.ActiveView.PartialRefresh(esriViewDrawPhase.esriViewGeography, null, null);
}
}
catch (SystemException e)
{
if (e.Message == "The spatial index grid size is invalid.")
{
IFeatureClassLoad pFCL = pFeatureLayer.FeatureClass as IFeatureClassLoad;
pFCL.LoadOnlyMode = true;
pfeature.Store();
pFCL.LoadOnlyMode = false;
}
else
{
MessageBox.Show(e.Message);
}
}
//IMultiPatch pMP = new MultiPatchClass();
//for (int i = 0; i < 5; i++)
//{
// IFeature pF = pFC.CreateFeature();
// IImport3DFile pI3D = new Import3DFileClass();
// pI3D.CreateFromFile(@"C:\Users\Administrator\Desktop\sampleAnalysis\1002.3ds");
// IMultiPatch pMP = pI3D.Geometry as IMultiPatch;
// // ITransform3D pT3D = pMP as ITransform3D;
// pF.Shape = pMP as IGeometry;
// ITransform3D pT3D = pF.Shape as ITransform3D;
// pT3D.Move3D(5293 + 1600, -20427 , 800 + i*400);
// // pF.Shape = pI3D.Geometry;
}
}
if (pRasterLayer != null)
{
//ITin pTin = pTinLayer.Dataset;
//ISurface pSurface = ((ITinAdvanced)pTin).Surface;
IMapControl2 pMapCtr = (((IToolbarControl)m_hookHelper.Hook).Buddy) as IMapControl2;
//在mapctr操作
if (pMapCtr != null)
{
IPoint mapPoint = pMapCtr.ToMapPoint(X, Y);
//添加手工编辑窗口
FrmManualSetModelPara pManualSetPara = null;
if (m_pModel == null)
{
Pt2d ptCurrent = new Pt2d();
ptCurrent.X = mapPoint.X;
ptCurrent.Y = mapPoint.Y;
pManualSetPara = new FrmManualSetModelPara(ptCurrent, "Rock");
}
else
{
m_pModel.x = mapPoint.X;
m_pModel.y = mapPoint.Y;
pManualSetPara = new FrmManualSetModelPara(m_pModel);
}
if (pManualSetPara.ShowDialog() == DialogResult.OK)
{
if (m_pModel == null)
{
m_pModel = new Model();
}
//bool bFlag = m_pModel.compareModelPara(pManualSetPara.m_pModel);
//if (!bFlag)//模型参数已经被修改，需要重新生成新的模型文件
//{
// m_pModel.copyFromModel(pManualSetPara.m_pModel);
// //获得当前文件路径下的石头模型文件,并随机获得石块模型
// String szAppPathName = GetParentPathofExe() + @"Resource\RockModel";
// String[] szFileList = Directory.GetFiles(szAppPathName);
// if (szFileList.Length > 0)
// {
// Random r = new Random(TerrainGen.Chaos_GetRandomSeed());
// int nRandomPos = r.Next(szFileList.Length);
// String szModelOutputFilename = szFileList[nRandomPos];
// m_listModels.Add(szModelOutputFilename);
// }
//}
m_pModel.copyFromModel(pManualSetPara.m_pModel);
//获得当前文件路径下的石头模型文件,并随机获得石块模型
String szAppPathName = ClsGlobal.GetParentPathofExe() + @"Resource\RockModel";
String[] szFileList = Directory.GetFiles(szAppPathName);
if (szFileList.Length > 0)
{
Random r = new Random(TerrainGen.Chaos_GetRandomSeed());
int nRandomPos = r.Next(szFileList.Length);
String szModelOutputFilename = szFileList[nRandomPos];
m_listModels.Add(szModelOutputFilename);
}
//保存当前添加的模型参数
Model pTmpModel = new Model();
pTmpModel.copyFromModel(pManualSetPara.m_pModel);
m_manualAddModels.Add(pTmpModel);
}
else
{
return;
}
IZAware za = mapPoint as IZAware;
za.ZAware = true;
double zVal;
IRaster2 pRaster2 = pRasterLayer.Raster as IRaster2;
int col, row;
pRaster2.MapToPixel(mapPoint.X, mapPoint.Y, out col, out row);
zVal = Convert.ToDouble(pRaster2.GetPixelValue(0, col, row));
if (double.IsNaN(zVal))
{
MessageBox.Show("获取模型的高度失败！");
return;
}
if (m_listModels.Count == 0)
{
MessageBox.Show("请先生成模型！");
return;
}
mapPoint.Z = zVal;
try
{
//IMultiPatch pMP = new MultiPatchClass();
IFeatureClass pFC = pFeatureLayer.FeatureClass;
IFeature pF = pFC.CreateFeature();
IImport3DFile pI3D = new Import3DFileClass();
pI3D.CreateFromFile(m_listModels[m_listModels.Count - 1]);
IMultiPatch pMP = pI3D.Geometry as IMultiPatch;
ITransform3D pT3D = pMP as ITransform3D;
pT3D.Move3D(mapPoint.X, mapPoint.Y, mapPoint.Z);
IRasterProps pRasterProps = pRasterLayer.Raster as IRasterProps;
double xmax = pRasterProps.Extent.XMax;
double xmin = pRasterProps.Extent.XMin;
double ymax = pRasterProps.Extent.YMax;
double ymin = pRasterProps.Extent.YMin;
//生成地形的地理范围
double dbGeoRangeX = xmax - xmin;
double dbGeoRangeY = ymax - ymin;
double dbModelRangeX = pMP.Envelope.Width;
double dbModelRangeY = pMP.Envelope.Height;
double dbModelRangeZ = pMP.Envelope.ZMax - pMP.Envelope.ZMin;
//根据地形大小改变石块大小
//double dbRandomSize = Math.Min(dbGeoRangeX, dbGeoRangeY) * (0.03 ); //[0.03 0.06]
//double dbScale = dbRandomSize / Math.Max(dbModelRangeX, dbModelRangeY);
double dbScale = m_pModel.dbSize / Math.Max(dbModelRangeX, Math.Max(dbModelRangeY, dbModelRangeZ));
pT3D.Scale3D(mapPoint, dbScale, dbScale, dbScale);
pF.Shape = pMP as IGeometry;
pF.Store();
if (pMapCtr != null)
{
pMapCtr.ActiveView.PartialRefresh(esriViewDrawPhase.esriViewGeography, null, null);
}
}
catch (SystemException e)
{
MessageBox.Show(e.Message);
}
}
}
}
