四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法
【專利摘要】本發(fā)明公開了一種基于四參數(shù)和坐標(biāo)文件正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法。包括以下步驟：A、讀取舊坐標(biāo)文件的坐標(biāo)參數(shù)信息，包括北方向坐標(biāo)、東方向坐標(biāo)、北方向像元分辨率、東方向像元分辨率；B、計(jì)算不同坐標(biāo)系下正射影像的平面坐標(biāo)轉(zhuǎn)換四參數(shù)即豎軸北方向X坐標(biāo)的偏移量ΔX、橫軸東方向Y坐標(biāo)的偏移量ΔY、縮放參數(shù)K、豎軸北方向旋轉(zhuǎn)角α；C、通過以上參數(shù)將舊的坐標(biāo)文件的坐標(biāo)參數(shù)進(jìn)行計(jì)算，得到新坐標(biāo)系下的對應(yīng)參數(shù)，將新信息寫入坐標(biāo)文件；D、讀取新坐標(biāo)系下正射影像的坐標(biāo)數(shù)信息，糾正新坐標(biāo)系下正射影像數(shù)據(jù)。本發(fā)明方法可以節(jié)約大量數(shù)據(jù)處理時(shí)間和硬件資源，非常適用于一般作業(yè)部門或小型單位的正射影像坐標(biāo)系轉(zhuǎn)換。
【專利說明】
四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法
技術(shù)領(lǐng)域
[0001]本發(fā)明涉及正射影像處理技術(shù)領(lǐng)域，具體涉及一種基于四參數(shù)和坐標(biāo)文件正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法。
【背景技術(shù)】
[0002]正射影像的坐標(biāo)在不同坐標(biāo)系下的相互轉(zhuǎn)換處理，常常需要借助專業(yè)軟件來完成，形成多套不同坐標(biāo)系下的正射影像成果;這種模式將耗費(fèi)大量時(shí)間和存儲空間，造成資源的巨大浪費(fèi)。
[0003]現(xiàn)有的遙感或GIS平臺軟件，如ArcGIS，Erdas等已經(jīng)具備比較完備的地理信息處理功能，包括不同格式影像圖的配準(zhǔn)以及融合處理，對配準(zhǔn)影像的裁剪和正確定位顯示。但工程應(yīng)用中，這些軟件的使用難度大，購置成本高，不適宜小型生產(chǎn)單位或一線作業(yè)部門的實(shí)際生產(chǎn)條件。
[0004]坐標(biāo)轉(zhuǎn)換雖然可以通過專業(yè)軟件實(shí)現(xiàn)，但存在如下缺點(diǎn)和不足:
[0005]1、正射影像坐標(biāo)轉(zhuǎn)換必須采用專業(yè)平臺軟件，軟件成本高；
[0006]2、處理過程復(fù)雜，必須采用高性能計(jì)算機(jī)，硬件成本高，一般人員難于掌握；
[0007]3、通常需要進(jìn)行不同坐標(biāo)系下數(shù)據(jù)的重采樣，浪費(fèi)大量計(jì)算機(jī)存儲空間；
[0008]4、通用性和易用性較差，一般作業(yè)部門或小型單位難以實(shí)現(xiàn)。

【發(fā)明內(nèi)容】

[0009]本發(fā)明解決的技術(shù)問題在于提供一種基于四參數(shù)和坐標(biāo)文件正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，通過坐標(biāo)轉(zhuǎn)換四參數(shù)和坐標(biāo)文件，結(jié)合常見的開發(fā)平臺，實(shí)現(xiàn)了高效的正射影像坐標(biāo)轉(zhuǎn)換。
[0010]本發(fā)明解決上述技術(shù)問題的技術(shù)方案是:
[0011]所述的方法包括以下步驟:
[0012]A、讀取舊坐標(biāo)文件的坐標(biāo)參數(shù)信息，包括北方向坐標(biāo)、東方向坐標(biāo)、北方向像元分辨率、東方向像元分辨率；
[0013]B、計(jì)算不同坐標(biāo)系下正射影像的平面坐標(biāo)轉(zhuǎn)換四參數(shù)即豎軸北方向X坐標(biāo)的偏移量Δ X、橫軸東方向Y坐標(biāo)的偏移量△ Y、縮放參數(shù)K、豎軸北方向旋轉(zhuǎn)角α ；
[0014]C、通過以上參數(shù)將I日的坐標(biāo)文件的坐標(biāo)參數(shù)進(jìn)行計(jì)算，得到新坐標(biāo)系下的對應(yīng)參數(shù)，將新信息寫入坐標(biāo)文件；
[0015]D、讀取新坐標(biāo)系下正射影像的坐標(biāo)數(shù)信息，糾正新坐標(biāo)系下正射影像數(shù)據(jù)。
[0016]更為具體而言，
[0017]所述步驟A包括以下步驟:
[0018]Al、檢查正射影像坐標(biāo)文件完整性，確保正射影像文件名稱與坐標(biāo)文件名稱一致；
[0019]Α2、讀取坐標(biāo)文件信息，與文件名稱做一致性檢查，確保坐標(biāo)文件的北方向坐標(biāo)、東方向坐標(biāo)、北方向像元分辨率、東方向像元分辨率、坐標(biāo)走向等信息內(nèi)容準(zhǔn)確；
[0020]A3、判定坐標(biāo)信息無誤，讀取坐標(biāo)文件信息。
[0021]所述的步驟Al中，坐標(biāo)信息錯(cuò)誤包括坐標(biāo)文件信息與文件名稱隱含的坐標(biāo)信息不一致的情形。
[0022]所述步驟B包括以下步驟:
[0023 ] B1、通過坐標(biāo)轉(zhuǎn)換軟件，讀取不同坐標(biāo)系下同名點(diǎn)的兩套坐標(biāo)；
[0024]B2、計(jì)算不同坐標(biāo)系下正射影像的平面坐標(biāo)轉(zhuǎn)換四參數(shù)即豎軸北方向X坐標(biāo)的偏移量△ X、橫軸東方向Y坐標(biāo)的偏移量△ Y、縮放參數(shù)K、豎軸北方向旋轉(zhuǎn)角α，具體公式如下:
[0025]設(shè)(乂0、￥0)，(乂1、￥1)為舊坐標(biāo)系下兩點(diǎn)坐標(biāo)，
[0026](XO'、Y0')，(XI'、Y1')為新坐標(biāo)系下的坐標(biāo)；
[0027]DX'=Xl'-XO'，
[0028]DY'=Yl'-YO'，
[0029]DX=Xl-XO,
[0030]DY = Yl-YO,
[0031]Kl=K*cosa，
[0032]K2=K*sina，
[0033]貝lj，
[0034]a = arcsin((DY'*DX_DX'*DY)/(DX'2+DY'2))，
[0035]K = s qrt ((DY' 2+DX'2) / (DY2+DX2))，
[0036]ΔΧ=Χ0'+Χ0*Κ1_Υ0*Κ2，
[0037]ΔΥ = Υ0' +Y0*K2+X0*KI ；
[0038]其中，cos是余弦函數(shù)，sin是正弦函數(shù)，arcsin是反正弦函數(shù)，sqrt平方根計(jì)算函數(shù)。
[0039]所述步驟C包括以下步驟:
[0040]Cl、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)北方向偏移量Δ X和東方向偏移量Δ Y，將坐標(biāo)文件的起始點(diǎn)平移到新的坐標(biāo)系下；
[0041 ] C2、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)中的縮放參數(shù)K，將坐標(biāo)文件的像元大小按照縮放參數(shù)進(jìn)行縮放；
[0042]C3、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)豎軸北方向旋轉(zhuǎn)角α，將坐標(biāo)文件以起始點(diǎn)為圓心進(jìn)行旋轉(zhuǎn)，記錄轉(zhuǎn)換角度；
[0043]C4、將新的坐標(biāo)文件的起始點(diǎn)、北方向像元大小、東方向像元大小、旋轉(zhuǎn)角α等信息寫入新的坐標(biāo)文件，具體計(jì)算公式如下:
[0044]設(shè)(Χ0、Υ0)為就坐標(biāo)，(XO'、Υ0')為新坐標(biāo)，則
[0045]XO'= AX+X0*K*cosa-Y0*K*sina，
[0046]YO'= AY+Y0*K*sina+X0*K*cosa，
[0047]其中，cos是余弦函數(shù)，sin是正弦函數(shù)。
[0048]所述步驟D包括以下步驟:
[0049]Dl、讀取新坐標(biāo)系下正射影像的坐標(biāo)數(shù)信息；
[0050]D2、將正射影像數(shù)據(jù)按照新的坐標(biāo)信息縮放、旋轉(zhuǎn)、移動(dòng)到新的坐標(biāo)文件指定的位置及范圍。
[0051]本發(fā)明通過將坐標(biāo)文件名稱轉(zhuǎn)換四參數(shù)與坐標(biāo)文件結(jié)合，與常見的影像顯示平臺及開發(fā)平臺進(jìn)行結(jié)合，如基于PC的Windows 7操作系統(tǒng)上的常用影像顯示軟件及開發(fā)平臺為例，Microsta1n V8、Auto CAD完成影像數(shù)據(jù)顯不,Microsoft VB、AutoCAD 2005平臺自帶的VBA、Microsta1n V8平臺自帶的VBA完成新的坐標(biāo)轉(zhuǎn)換中參數(shù)求解的開發(fā)，S卩可實(shí)現(xiàn)正射影像的坐標(biāo)轉(zhuǎn)換問題。完成了以往必須借助專業(yè)化遙感軟硬件才能完成的任務(wù)，省去了專業(yè)化遙感軟硬件的使用費(fèi)用。該方法具有原理簡單，易于實(shí)現(xiàn)，高效穩(wěn)定的特點(diǎn)，可以大大節(jié)約數(shù)據(jù)處理的難度、成本和時(shí)間。
【附圖說明】
[0052]下面結(jié)合附圖對本發(fā)明進(jìn)一步說明:
[0053]圖1是本發(fā)明方法的流程示意圖；
[0054]圖2是本發(fā)明不同坐標(biāo)系下正射影像的平面坐標(biāo)轉(zhuǎn)換四參數(shù)示意圖。
【具體實(shí)施方式】
[0055]如圖1所示，是本發(fā)明基于四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，本實(shí)施例坐標(biāo)轉(zhuǎn)換1:2000正射影像在北京54坐標(biāo)系到西安80坐標(biāo)系轉(zhuǎn)換過程，運(yùn)行平臺為PC上的Windows 7操作系統(tǒng)。具體包括以下步驟:
[0056]A、讀取北京54坐標(biāo)文件的坐標(biāo)參數(shù)信息，包括北方向坐標(biāo)X、東方向坐標(biāo)Y、北方向像元分辨率Ce IIX、東方向像元分辨率Ce 11Y;
[0057]B、計(jì)算不同坐標(biāo)系下正射影像的平面坐標(biāo)轉(zhuǎn)換四參數(shù)，即豎軸北方向X坐標(biāo)的偏移量△ X、橫軸東方向Y坐標(biāo)的偏移量△ Y、縮放參數(shù)K及在兩個(gè)方向的份量Kl及K2、豎軸北方向旋轉(zhuǎn)角a; (X'，Y')為新坐標(biāo)系，(Χ，Υ)為舊坐標(biāo);則P點(diǎn)新坐標(biāo)為X' = A X+X*K*cosa-Y*K*
sina，Υ'= ΔY+Y*K*sina+X*K*cosa；
[0058]C、通過以上參數(shù)將I日的坐標(biāo)文件的坐標(biāo)參數(shù)進(jìn)行計(jì)算，得到新坐標(biāo)系下的對應(yīng)參數(shù)，將新信息寫入坐標(biāo)文件，將舊的坐標(biāo)文件名稱文件備份；
[0059]D、讀取新坐標(biāo)系下正射影像的坐標(biāo)數(shù)信息，糾正新坐標(biāo)系下正射影像數(shù)據(jù)。
[0060]更為具體而言，
[0061 ] 所述步驟A包括以下步驟:
[0062]Al、檢查正射影像坐標(biāo)文件完整性、文件格式的正確性，確保正射影像文件名稱與坐標(biāo)文件名稱一致，僅文件格式不同；
[0063]A2、讀取坐標(biāo)文件信息，與文件名稱做一致性檢查，確保坐標(biāo)文件內(nèi)容準(zhǔn)確，如北方向坐標(biāo)、東方向坐標(biāo)、北方向像元分辨率、東方向像元分辨率，坐標(biāo)走向信息等；
[0064]A3、判定坐標(biāo)信息無誤，讀取坐標(biāo)文件信息。
[0065]步驟Al中，坐標(biāo)信息錯(cuò)誤包括坐標(biāo)文件信息與文件名稱隱含的坐標(biāo)信息不一致的情況。
[0066]所述步驟B包括以下步驟:
[0067]B1、通過坐標(biāo)轉(zhuǎn)換軟件，讀取不同坐標(biāo)系下同名點(diǎn)的兩套坐標(biāo)；
[0068]B2、計(jì)算西安80坐標(biāo)系下正射影像的平面坐標(biāo)轉(zhuǎn)換四參數(shù)即豎軸北方向X坐標(biāo)的偏移量△ X、橫軸東方向Y坐標(biāo)的偏移量△ Y、縮放參數(shù)K、豎軸北方向旋轉(zhuǎn)角a，具體公式如下:
[0069]設(shè)(乂0、￥0)，(乂1、￥1)為舊坐標(biāo)系下兩點(diǎn)坐標(biāo)，
[0070](XO'、Y0')，(XI'、Y1')為新坐標(biāo)系下的坐標(biāo)，
[0071]DX'=Xl'-XO'，
[0072]DY'=Yl'-YO'，
[0073]DX=Xl-XO,
[0074]DY = Yl-YO,
[0075]α為旋轉(zhuǎn)角，
[0076]K為縮放因子，
[0077]Kl=K*cosa，
[0078]K2=K*sina，
[0079]Δ X為X方向偏移量，
[0080]Δ Y為Y方向偏移量，
[0081 ]則由步驟C4四參數(shù)公式推理可知
[0082]a = arcsin((DY'*DX_DX'*DY)/(DX'2+DY'2))，
[0083]K = sqrt ((DY' 2+DX'2)/ (DY2+DX2))，
[0084]ΔΧ=Χ0'+Χ0*Κ1_Υ0*Κ2，
[0085]Δ Y = YO' +Υ0*Κ2+Χ0*ΚI ；
[0086]其中，cos是余弦函數(shù)，sin是正弦函數(shù)，arcsin是反正弦函數(shù)，sqrt平方根計(jì)算函數(shù)；
[0087]注:此處X，Y坐標(biāo)是指北方向和東方向即上方向和右方向，與計(jì)算機(jī)坐標(biāo)的X，Y坐標(biāo)方向互換。
[0088]所述步驟C包括以下步驟:
[0089]Cl、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)北方向偏移量Δ X和東方向偏移量Δ Y，將坐標(biāo)文件的起始點(diǎn)平移到新的坐標(biāo)系下；
[°09°]C2、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)中的縮放參數(shù)K，將坐標(biāo)文件的像元大小按照縮放參數(shù)進(jìn)行縮放；
[0091]C3、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)豎軸北方向旋轉(zhuǎn)角α，將坐標(biāo)文件以起始點(diǎn)為圓心進(jìn)行旋轉(zhuǎn)，記錄轉(zhuǎn)換角度；
[0092]C4、將新的坐標(biāo)文件的起始點(diǎn)、北方向像元大小、東方向像元大小、旋轉(zhuǎn)角α等信息寫入新的坐標(biāo)文件，具體計(jì)算公式如下:
[0093]設(shè)(Χ0、Υ0)為就坐標(biāo)，(XO'、Υ0')為新坐標(biāo)，則
[0094]XO'= Δ X+X0*K*cosa-Y0*K*sina，
[0095]YO'= AY+Y0*K*sina+X0*K*cosa，
[0096]其中K為尺度因子即縮放參數(shù)、a為旋轉(zhuǎn)角，ΔΧ，△Y為相應(yīng)的偏移量。
[0097]所述步驟D包括以下步驟:
[0098]Dl、讀取新坐標(biāo)系下正射影像的坐標(biāo)數(shù)信息；
[0099]D2、將正射影像數(shù)據(jù)按照新的坐標(biāo)信息縮放、旋轉(zhuǎn)、移動(dòng)到新的坐標(biāo)文件指定的位置及范圍。
[0100]以廣西南寧市主城區(qū)2011年衛(wèi)星0.2m彩色影像為例，北京54坐標(biāo)系下兩點(diǎn)坐標(biāo)分別為乂() = 2507400、￥() = 515800，乂1 = 2540800、￥1 = 561080;西安80坐標(biāo)下兩點(diǎn)坐標(biāo)分別為父'0=2507340.277、Y'ο = 515726.960，X'I = 2540740.351、Y'I = 561007.487，則求得K =I.00000830743314,α = 00.927秒，Δ Χ = -78.234698，Δ Υ = -88.5935。假設(shè)某幅正射影像500000-2500000-1000-1000.jpg的原始北京54坐標(biāo)文件信息為東方向500000.1、北方向2500999.9、旋轉(zhuǎn)角為O秒、縮放因子為1、坐標(biāo)沿南北方向向南，東西方向向東；則新的西安80坐標(biāo)系下的東方向2499940.287、北方向499926.7956、旋轉(zhuǎn)角為00.927秒、縮放因子為1.00000830743314、坐標(biāo)沿南北方向向南，東西方向向東。其他原始北京54坐標(biāo)下的正射影像參也可以通過以上方法，求解四參數(shù)，實(shí)現(xiàn)全部正射影像文件的新舊坐標(biāo)轉(zhuǎn)換。
【主權(quán)項(xiàng)】
1.一種四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，其特征在于，所述的方法包括以下步驟: A、讀取舊坐標(biāo)文件的坐標(biāo)參數(shù)信息，包括北方向坐標(biāo)、東方向坐標(biāo)、北方向像元分辨率、東方向像元分辨率； B、計(jì)算不同坐標(biāo)系下正射影像的平面坐標(biāo)轉(zhuǎn)換四參數(shù)即豎軸北方向X坐標(biāo)的偏移量△X、橫軸東方向Y坐標(biāo)的偏移量A Y、縮放參數(shù)K、豎軸北方向旋轉(zhuǎn)角α; C、通過以上參數(shù)將舊的坐標(biāo)文件的坐標(biāo)參數(shù)進(jìn)行計(jì)算，得到新坐標(biāo)系下的對應(yīng)參數(shù)，將新信息寫入坐標(biāo)文件； D、讀取新坐標(biāo)系下正射影像的坐標(biāo)數(shù)信息，糾正新坐標(biāo)系下正射影像數(shù)據(jù)。2.根據(jù)權(quán)利要求1所述的四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，其特征在于，所述步驟A包括以下步驟: Al、檢查正射影像坐標(biāo)文件完整性，確保正射影像文件名稱與坐標(biāo)文件名稱一致； Α2、讀取坐標(biāo)文件信息，與文件名稱做一致性檢查，確保坐標(biāo)文件的北方向坐標(biāo)、東方向坐標(biāo)、北方向像元分辨率、東方向像元分辨率、坐標(biāo)走向等信息內(nèi)容準(zhǔn)確； A3、判定坐標(biāo)信息無誤，讀取坐標(biāo)文件信息。3.根據(jù)權(quán)利要求2所述的四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，其特征在于:所述的步驟Al中，坐標(biāo)信息錯(cuò)誤包括坐標(biāo)文件信息與文件名稱隱含的坐標(biāo)信息不一致的情形。4.根據(jù)權(quán)利要求1、2或3所述的四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，其特征在于，所述步驟B包括以下步驟: B1、通過坐標(biāo)轉(zhuǎn)換軟件，讀取不同坐標(biāo)系下同名點(diǎn)的兩套坐標(biāo)； Β2、計(jì)算不同坐標(biāo)系下正射影像的平面坐標(biāo)轉(zhuǎn)換四參數(shù)即豎軸北方向X坐標(biāo)的偏移量A X、橫軸東方向Y坐標(biāo)的偏移量Δ Y、縮放參數(shù)K、豎軸北方向旋轉(zhuǎn)角α，具體公式如下: 設(shè)(乂0、￥0)，(乂1、￥1)為舊坐標(biāo)系下兩點(diǎn)坐標(biāo)， (XO'、YO')，(XI'、Yl')為新坐標(biāo)系下的坐標(biāo)； DX'=Xl'-XO'， DY'=Yl'-YO'，DX = Xl-XO,DY = Yl-YO,Kl =K 氺 cosa，K2 = K>sina， 則， a = arcsin((DY'*DX-DX'*DY)/(DX'2+DY'2))，K = s qr t ((DY' 2+DX'2) / (DY2+DX2))， ΔΧ = Χ0'+Χ0*Κ1-Υ0*Κ2, ΔΥ = Υ0'+Υ0*Κ2+Χ0*Κ1； 其中，cos是余弦函數(shù)，sin是正弦函數(shù)，arcsin是反正弦函數(shù)，sqrt平方根計(jì)算函數(shù)。5.根據(jù)權(quán)利要求1、2或3所述的所述的四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，其特征在于，所述步驟C包括以下步驟: Cl、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)北方向偏移量△ X和東方向偏移量△ Y，將坐標(biāo)文件的起始點(diǎn)平移到新的坐標(biāo)系下； C2、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)中的縮放參數(shù)K，將坐標(biāo)文件的像元大小按照縮放參數(shù)進(jìn)行縮放； C3、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)豎軸北方向旋轉(zhuǎn)角α，將坐標(biāo)文件以起始點(diǎn)為圓心進(jìn)行旋轉(zhuǎn)，記錄轉(zhuǎn)換角度； C4、將新的坐標(biāo)文件的起始點(diǎn)、北方向像元大小、東方向像元大小、旋轉(zhuǎn)角α等信息寫入新的坐標(biāo)文件，具體計(jì)算公式如下: 設(shè)(X(KYO)為就坐標(biāo)，(XO'、Υ0')為新坐標(biāo)，則 XO'= Δ X+X0*K*cosa-Y0*K*sina， YO'= Δ Y+Y0*K*sina+X0*K*cosa， 其中，cos是余弦函數(shù)，sin是正弦函數(shù)。6.根據(jù)權(quán)利要求4所述的所述的四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，其特征在于，所述步驟C包括以下步驟: Cl、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)北方向偏移量△ X和東方向偏移量△ Y，將坐標(biāo)文件的起始點(diǎn)平移到新的坐標(biāo)系下； C2、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)中的縮放參數(shù)K，將坐標(biāo)文件的像元大小按照縮放參數(shù)進(jìn)行縮放； C3、讀取舊的坐標(biāo)，根據(jù)坐標(biāo)轉(zhuǎn)換參數(shù)豎軸北方向旋轉(zhuǎn)角a，將坐標(biāo)文件以起始點(diǎn)為圓心進(jìn)行旋轉(zhuǎn)，記錄轉(zhuǎn)換角度； C4、將新的坐標(biāo)文件的起始點(diǎn)、北方向像元大小、東方向像元大小、旋轉(zhuǎn)角a等信息寫入新的坐標(biāo)文件，具體計(jì)算公式如下: 設(shè)(Χ0、Υ0)為就坐標(biāo)，(XO'、Y0')為新坐標(biāo)，則 XO'= Δ X+X0*K*cosa-Y0*K*sina， YO'= Δ Y+Y0*K*sina+X0*K*cosa。7.根據(jù)權(quán)利要求1、2或3所述的四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，其特征在于，所述步驟D包括以下步驟: D1、讀取新坐標(biāo)系下正射影像的坐標(biāo)數(shù)信息； D2、將正射影像數(shù)據(jù)按照新的坐標(biāo)信息縮放、旋轉(zhuǎn)、移動(dòng)到新的坐標(biāo)文件指定的位置及范圍。8.根據(jù)權(quán)利要求4所述的四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，其特征在于，所述步驟D包括以下步驟: D1、讀取新坐標(biāo)系下正射影像的坐標(biāo)數(shù)信息； D2、將正射影像數(shù)據(jù)按照新的坐標(biāo)信息縮放、旋轉(zhuǎn)、移動(dòng)到新的坐標(biāo)文件指定的位置及范圍。9.根據(jù)權(quán)利要求5所述的四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，其特征在于，所述步驟D包括以下步驟: D1、讀取新坐標(biāo)系下正射影像的坐標(biāo)數(shù)信息； D2、將正射影像數(shù)據(jù)按照新的坐標(biāo)信息縮放、旋轉(zhuǎn)、移動(dòng)到新的坐標(biāo)文件指定的位置及范圍。10.根據(jù)權(quán)利要求6所述的四參數(shù)和坐標(biāo)文件的正射影像在不同坐標(biāo)系的轉(zhuǎn)換方法，其特征在于，所述步驟D包括以下步驟: D1、讀取新坐標(biāo)系下正射影像的坐標(biāo)數(shù)信息； D2、將正射影像數(shù)據(jù)按照新的坐標(biāo)信息縮放、旋轉(zhuǎn)、移動(dòng)到新的坐標(biāo)文件指定的位置及范圍。
【文檔編號】G06F17/30GK105893540SQ201610196807
【公開日】2016年8月24日
【申請日】2016年3月30日
【發(fā)明人】魏金占, 曹文海, 陳明輝, 凌建國, 岳國森, 張興國
【申請人】廣東中冶地理信息股份有限公司, 廣東中達(dá)規(guī)谷地信科技有限公司