输入资料(经纬度).doc

輸入資料(經緯度) 編號 經度 緯度 0 120.697632719 24.135416905 1 120.634248518 24.156401633 2 120.687778143 24.136219644 3 120.622795329 24.108526787 4 120.616951028 24.104954300 5 120.685055794 24.127922475 6 120.652497652 24.154471448 7 120.635832970 24.106824480 8 120.639950471 24.151686914 9 120.685258920 24.133322313 圖式 0~9 分別用 A~J 代替 經緯度轉XYZ n=a/Math.sqrt(1-e*Math.pow(Math.sin(lot[i]*3.14159/180),2)); x=(n+h)*Math.cos(lot[i]*3.14159/180)*Math.cos(lat[i]*3.14159/180); y=(n+h)*Math.cos(lot[i]*3.14159/180)*Math.sin(lat[i]*3.14159/180); z=(n*(1-e)+h)*Math.sin(lot[i]*3.14159/180); xx=(n+h)*Math.cos(dlot[i]*3.14159/180)*Math.cos(dlat[i]*3.14159/180); yy=(n+h)*Math.cos(dlot[i]*3.14159/180)*Math.sin(dlat[i]*3.14159/180); zz=(n*(1-e)+h)*Math.sin(dlot[i]*3.14159/180); 產生隨機中心點 st.k=3; for(int i=0;ist.k;i++){ st.round=(Math.random()*10); st.round2=(((int)st.round)%st.px.size())+1; for(int j=0;jst.round2;j++){ st.r=(Math.random()*10); st.r2=((int)st.r)%st.px.size(); st.center[i][0]=st.center[i][0]+st.px.get(st.r2); st.center[i][1]=st.center[i][1]+st.py.get(st.r2); st.center[i][2]=st.center[i][2]+st.pz.get(st.r2); //用亂數產生座標當作中心點的初始值 } st.center[i][0]=st.center[i][0]/st.round2; st.center[i][1]=st.center[i][1]/st.round2; st.center[i][2]=st.center[i][2]/st.round2; //計算出初始的中心點 st.newgroup[i]=1; //將新分群的值設為1，使第一次運算能正常運作 } 第0中心點的X座標為-2973800.6269722194第0中心點的Y座標為-1331694.443885626第0中心點的Z座標為5464656.430296873第1中心點的X座標為-2977887.5060579544第1中心點的Y座標為-1334157.328507759第1中心點的Z座標為5461848.634295342第2中心點的X座標為-2974032.0559696243第2中心點的Y座標為-1331679.8558799431第2中心點的Z座標為5464534.295048131 計算分群 while(Arrays.equals(st.group,st.newgroup)==false){//如果新舊分群解果不相等則進行下一次運算 for(int i=0;ist.px.size();i++){ st.group[i]=st.newgroup[i];//將舊分群解果由新分群解果替代