using System;
using SuperMapLib;

namespace ProfileCtrl
{
	/// <summary>
	/// 地層線類
	/// </summary>
	public class StratumLine
	{
		public StratumLine(Bore[] bores)
		{

			//
			// TODO: 在此處添加構造函數邏輯
			//
			this.scaleX=PublicValue.scaleX;
			this.scaleY=PublicValue.scaleY;
			this.bores=bores;
			zkCount=bores.Length;
			this.space=PublicValue.space;
		}
		float scaleX, scaleY, scaleT=0.005f;//文本縮放比例因子
		private Bore[] bores;//鉆孔數據的集合
		private int zkCount;//鉆孔個數
		private int space;
		/// <summary>
		/// 文本縮放比例因子
		/// </summary>
		public float ScaleT
		{
			get{return scaleT;}
			set{scaleT=value;}
		}
		public void Draw(soRecordset objRs)
		{
			DrawTB(objRs);
			DrawSumLine(objRs);
		}

	
		
		#region "自動生成概化地質界線"
		/// <summary>
		/// 概化地質界線方法
		/// </summary>
		/// <param name="objRs">存放剖面的數據集</param>
		private void DrawSumLine(soRecordset objRs)
		{
			soPoint objPoint;//存放地質線的臨時點
			soPoints objPointsT,objPointsB;	//頂、底板點集			
			soGeoLine objGeoLine;//地層線

			bool[][] isPop = new bool[zkCount][];	//設一臨時變量，判斷此地層是否取出
			for(int i=0;i<zkCount;i++)
			{
				isPop[i] = new bool[bores[i].GetHSCDM().Length];
			}

			for(int i=0;i<zkCount;i++)		//循環所有鉆孔
			{	
				#region 
				for(int j=1;j<bores[i].GetHSCDM().Length;j++)			//循環所有含水層
				{							
					if(bores[i].GetHSCDM()[j] != "隔水層" && !isPop[i][j])	//畫現所有的含水層
					{											
						objPoint = new soPoint();		
						objPointsT = new soPoints();
						objPointsB = new soPoints();
						if(i==0)	//第一個鉆孔和高程標尺連線
						{
							objPoint.x = -(float)space/this.scaleX;
							objPoint.y = bores[i].GetCDBG()[j-1]/scaleY;
							objPointsT.Add(objPoint);
							objPoint.y = bores[i].GetCDBG()[j]/scaleY;
							objPointsB.Add(objPoint);					
						}
						//??????????
						if(i != 0)	//不是第一個鉆孔，前一鉆孔缺失此地層
						{
							float gene = 1f/3f;	//運用地質規則，算出尖滅位置
							float depth = bores[i].GetCDBG()[j-1] - bores[i].GetCDBG()[j];
							if(depth > 5)	gene = 0.5f;
							if(depth > 8)	gene = 2f/3f;	
							objPoint.x = (bores[i].X - (bores[i].X - bores[i-1].X) * gene)/scaleX;
							//							objPoint.x = (bores[i-1].X + bores[i].X) * 0.5/scaleX;	//暫時取中點
							objPoint.y = (bores[i].GetCDBG()[j-1] + bores[i].GetCDBG()[j]) * 0.5/scaleY;		//暫時取中點
							objPointsT.Add(objPoint);
							objPointsB.Add(objPoint);
						}					
						objPoint.x = bores[i].X/scaleX;
						objPoint.y = bores[i].GetCDBG()[j-1]/scaleY;
						objPointsT.Add(objPoint);
						objPoint.y = bores[i].GetCDBG()[j]/scaleY;
						objPointsB.Add(objPoint);

						

						#region
						int forNum = j;					
						for(int n = i+1;n<zkCount;n++)		//循環此鉆孔后所有鉆孔及地層
						{	
							#region
							bool isFind = false;						//是否找到地層							
							for(int k=1;k<bores[n].GetHSCDM().Length;k++)
							{	
								if(!isPop[n][k] && bores[i].GetHSCDM()[j] == bores[n].GetHSCDM()[k])
								{
									objPoint.x = bores[n].X/scaleX;
									objPoint.y = bores[n].GetCDBG()[k-1]/scaleY;
									objPointsT.Add(objPoint);
									objPoint.y = bores[n].GetCDBG()[k]/scaleY;
									objPointsB.Add(objPoint);
									isPop[n][k] = true;
									if( n == (zkCount-1))	//最后一個鉆孔和高程標尺連線
									{
										objPoint.x = (bores[zkCount-1].X + space)/this.scaleX;
										objPoint.y = bores[zkCount-1].GetCDBG()[k-1]/scaleY;
										objPointsT.Add(objPoint);
										objPoint.y = bores[zkCount-1].GetCDBG()[k]/scaleY;
										objPointsB.Add(objPoint);					
									}

									isFind = true;
									forNum = k;
									break;
								}								
							}
							#endregion
							if(!isFind)	//后一鉆孔缺失此地層
							{
								float gene = 1f/3f;	//運用地質規則，算出尖滅位置
								float depth = bores[i].GetCDBG()[j-1] - bores[i].GetCDBG()[j];
								if(depth > 5)	gene = 0.5f;
								if(depth > 8)	gene = 2f/3f;
								objPoint = new soPoint();
								objPoint.x = (bores[n-1].X + (bores[n].X - bores[n-1].X) * gene)/scaleX;
								//								objPoint.x = (bores[n-1].X + bores[n].X) * 0.5/scaleX;							//暫時取中點
								objPoint.y = (bores[n-1].GetCDBG()[forNum-1] + bores[n-1].GetCDBG()[forNum]) * 0.5/scaleY;	//暫時取中點
								objPointsT.Add(objPoint);
								objPointsB.Add(objPoint);
								break;
							}
						}
						#endregion
						objGeoLine = new soGeoLine();
						objGeoLine.AddPart(objPointsT);
						objRs.AddNew((soGeometry)objGeoLine);					
						objRs.Update();	
						objGeoLine = new soGeoLine();
						objGeoLine.AddPart(objPointsB);
						objRs.AddNew((soGeometry)objGeoLine);					
						objRs.Update();	
					}					
				}
				#endregion
			}
			objPoint = null;
			objPointsT = null;
			objPointsB = null;
			objGeoLine = null;
		}

		#endregion

		#region 畫地面和層底連線
		/// <summary>
		/// 連接每個探孔的地表和底部，生產地表起伏線和層底連線
		/// </summary>
		/// <param name="objRs">存放剖面線的數據集</param>
		private void DrawTB(soRecordset objRs)
		{
			//地面線、孔底線
			soPoints objPointsT = new soPoints();
			soPoints objPointsB = new soPoints();
			soPoint objPoint = new soPoint();
			objPoint.x = -(float)space/scaleX;			
			objPoint.y = bores[0].GetCDBG()[0]/scaleY;
			objPointsT.Add(objPoint);
			objPoint.y = bores[0].GetCDBG()[bores[0].GetCDBG().Length-1]/scaleY;
			objPointsB.Add(objPoint);
			for(int i=0;i<bores.Length;i++)
			{				
				objPoint.x = bores[i].X/scaleX;
				objPoint.y = bores[i].GetCDBG()[0]/scaleY;
				objPointsT.Add(objPoint);
				objPoint.y = bores[i].GetCDBG()[bores[i].GetCDBG().Length-1]/scaleY;
				objPointsB.Add(objPoint);
			}
			objPoint.x = (bores[zkCount-1].X + space)/this.scaleX;		
			objPoint.y = bores[zkCount-1].GetCDBG()[0]/scaleY;
			objPointsT.Add(objPoint);
			objPoint.y = bores[zkCount-1].GetCDBG()[bores[zkCount-1].GetCDBG().Length-1]/scaleY;
			objPointsB.Add(objPoint);
			soGeoLine objGeoLine = new soGeoLine();
			objGeoLine.AddPart(objPointsT);
			objRs.AddNew((soGeometry)objGeoLine);
			objRs.Update();
			objGeoLine = new soGeoLine();
			objGeoLine.AddPart(objPointsB);
			objRs.AddNew((soGeometry)objGeoLine);
			objRs.Update();

			objPoint = null;
			objPointsT = null;
			objPointsB = null;
			objGeoLine = null;
		}
		#endregion

	}
}