基于高清衛(wèi)星影像的地形圖一體化快速成圖方法
【專利摘要】本發(fā)明旨在提供一種基于高清衛(wèi)星影像的地形圖一體化快速成圖方法，該方法通過對衛(wèi)星影像的融合、分類識別，確定外業(yè)測繪的像控點，在野外同步完成數(shù)字化測圖、野外調(diào)繪及像控點采集；之后以配準后的彩色融合衛(wèi)星影像為底圖，補繪外業(yè)漏繪地物，完成地形圖的生成。本發(fā)明基于高清衛(wèi)星影像的地形圖一體化快速成圖方法將調(diào)繪、像控、補測等三個外業(yè)步驟融合為一，實現(xiàn)減少外業(yè)工作量的目標，與傳統(tǒng)的航測相比具有周期短的特點，與傳統(tǒng)人工測量相比，其外業(yè)工作量少、投入成本低且周期短。
【專利說明】基于高清衛(wèi)星影像的地形圖一體化快速成圖方法
【技術領域】
[0001]本發(fā)明涉及衛(wèi)星影像處理領域及地形圖成圖領域，具體涉及一種基于高清衛(wèi)星影像的地形圖一體化快速成圖方法。
【背景技術】
[0002]地形圖成圖多采用現(xiàn)場數(shù)字化測圖及攝影測量等常規(guī)手段，其中數(shù)字化測圖精度高但效率低，攝影測量受天氣條件限制等，當前前期規(guī)劃是對中大比例尺的地形圖需求周期非常短，一般勘測單位多采用人海戰(zhàn)術，以數(shù)字化測圖為主。
[0003]現(xiàn)有的GIS (地理信息系統(tǒng))平臺軟件，如ArcGIS，已經(jīng)具備比較完備的地理信息處理功能，包括不同格式影像圖的配準以及融合處理，對配準影像的裁剪和正確定位顯示，以及最后的地形圖生成功能。但工程應用中，地形圖的使用和采集多采用CAD平臺，而且ArcGIS的編輯采集功能較弱，很多復雜的大比例尺地形圖的符號、線性、字體及打印輸出比CAD弱且平臺應用范圍少，因此本發(fā)明必須結(jié)合生產(chǎn)實際，自動化部分步驟，節(jié)約大量的數(shù)據(jù)處理時間，比如影像圖在CAD的正確顯示、影像的自動裁剪等。
[0004]文獻“簡燦良，阮紅利，2007.利用高分辨率遙感衛(wèi)星影像更新1:1萬比例尺數(shù)字地形圖的研究[J].福建師范大學學報(自然科學版)，23(3):10-14”中公開了對IKONOS全色影像(Im)與IKONOS多光譜影像(4m)的配準、幾何糾正、輻射糾正的處理步驟，以及之后通過主成份變換法進行融合處理的步驟，實現(xiàn)對數(shù)字地形圖的更新。
[0005]這種模式僅僅是針對同一衛(wèi)星在一次采集中獲取的不同類型數(shù)據(jù)，是最為常見的融合模式。本發(fā)明主要解決的問題在于不同衛(wèi)星或數(shù)據(jù)來源的影像數(shù)據(jù)的融合。其主要解決的問題在于利用已有的彩色影像數(shù)據(jù)，解決當前僅有全色波段的影像色彩問題。以WorldViewl為例，其分辨率達到0.5米，采集能力強，但因其是全色影像(即黑白影像)，因此可嘗試將其他衛(wèi)星如同級別的GeoEyel或WorldView2的彩色影像色彩賦予WorldViewl影像，實現(xiàn)地物的高識別率。
[0006]文獻“徐穎，趙萍，黃亞萍，基于ArcGIS的標準分幅裁切新方法[J].現(xiàn)代測繪，29 (2): 43-45”中公開了一種基于ArcGIS中GRIDCLIP命令進行影像標準分幅的裁切方法，主要步驟包括:待裁切文件的預處理；確定標準分幅成果的圖幅號；確定每幅標準分幅圖幅的坐標范圍；運用GRIDCLIP命令完成裁切。
[0007]以上方法雖然可以實現(xiàn)影像的裁剪等問題，但存在如下問題:
1.海量影像的裁剪范圍自動計算問題
2.自動化分幅及自動裁剪
3.自動化的輸出與坐標文件的生成。
[0008]以上問題就如汽車的手動與自動的區(qū)別，本發(fā)明實現(xiàn)了海量影像的自動分幅裁剪及輸出功能。
[0009]發(fā)明人在實現(xiàn)本發(fā)明的過程中，發(fā)現(xiàn)現(xiàn)有技術存在以下缺點和不足:
1.數(shù)字化測圖成本高，效率低，工作周期長； 2.攝影測量技術在小范圍測圖是成本優(yōu)勢不明顯，工作周期較長；
3.中大比例尺地形圖需求范圍多為未開發(fā)區(qū)域，測區(qū)難以到達，基礎資料匱乏，上述兩種測圖模式無論成本還是作業(yè)周期都很難滿足客戶需求；
4.隨著衛(wèi)星影像分辨率的提高，其分辨率已經(jīng)滿足輔助測圖需求，但現(xiàn)有技術多停留在中小比例尺成圖應用方面，存在色彩信息缺失或應用范圍單一等缺陷，并且無法充分利用現(xiàn)有影像資源，成圖效率與質(zhì)量偏低。

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

[0010]本發(fā)明旨在提供一種基于高清衛(wèi)星影像的地形圖一體化快速成圖方法，將調(diào)繪、像控、補測等三個外業(yè)步驟融合為一，減少外業(yè)工作量、減少人力成本、加速成圖周期。
[0011]本發(fā)明的技術方案如下:一種基于高清衛(wèi)星影像的地形圖一體化快速成圖方法，包括以下步驟:
A、獲取最新的衛(wèi)星全色圖像和各類采集平臺、各種年份、各類分辨率的衛(wèi)星或者航拍彩色影像，利用各類分辨率的衛(wèi)星或者航拍彩色影像對衛(wèi)星或者航拍全色影像進行融合處理，得到彩色融合影像，對彩色融合影像進行分類識別，將識別得到的具有特定識別特征的地物作為像控點，進行外業(yè)數(shù)字化測繪，并補測彩色融合影像上遺失的地物；
B、對彩色融合影像進行配準處理，得到配準影像，以該影像為基準，讀取影像的四角坐標，根據(jù)四角坐標，計算出影像覆蓋的國家標準分幅圖幅，將落入國家標準分幅多邊形的的影像進行裁剪，對裁剪形成的單幅影像分幅命名，名稱中包含左下角東方向起始坐標和左下角北方向起始坐標，導出影像，生成坐標定位文件，坐標定位文件名稱與標準分幅影像同名；
C、根據(jù)步驟B中命名的單幅影像文件名稱，自動尋找對應坐標文件，讀取左下角東方向起始位置、左下角北方向起始位置、旋轉(zhuǎn)量、尺度變化量、東方向和北方向像元大小，根據(jù)影像的北方向和東方向像元個數(shù)，計算影像的正確定位并完成正確顯示；
D、循環(huán)進行步驟C，將各單幅影像依次正確定位顯示，生成初始地形圖；
E、以初始地形圖為底圖，根據(jù)外業(yè)測繪成果補繪重要地物和遺漏地物，對于次要地物以衛(wèi)星影像作為輔助成圖，生成成果地形圖。
[0012]其中，所述的配準是指:將測區(qū)的影像導入數(shù)據(jù)庫，以外業(yè)獲取的像控點為基準，進行精確配準；通常利用常見GIS平臺已有的配準功能，多采用均勻分布的NXN點模式，糾正的方法多采用四參數(shù)或七參數(shù)法，此配準方法可多次使用，其目的在于對于局部再次進行糾正提升整體影像配準精度。所述具有特定識別特征的地物指公路、河流、建筑等特征明顯的地物類別，上述地物同時也屬于地形圖的關鍵部分，需要作為像控點通過外業(yè)測繪進行確定與完善。在新建地形圖文件導入影像數(shù)據(jù)時，生成了國家規(guī)范的多個圖層，其中包括包含保密信息的控制點層、居民地層、交通層、水系層等。
[0013]所述步驟A中的融合處理包括以下步驟:
(1)將各類數(shù)據(jù)源的彩色影像與全色影像統(tǒng)一糾正到一個坐標下；
(2)將彩色影像重采樣，使得彩色影像與全色影像分辨率之比為1:1或1:4 ;
(3)將彩色影像的RGB波段分離出來；
(4)將分離出來的RGB波段影像與全色影像融合，得到彩色融合影像。[0014]所述步驟B中對配準影像的裁剪包括以下步驟:
(1)以配準衛(wèi)星影像為基準，讀取配準衛(wèi)星影像的四角坐標Xmin、Ymin、Xmax、Ymax，其中xmin、Yniin為左下角坐標，Xmax> Ymax為右上角坐標；
(2)根據(jù)四角坐標，計算出影像覆蓋的標準分幅圖幅，東方向和北方向的標準圖幅幅數(shù)計算公式如下:imax= int (XmaxX K/D) -1nt (XminX K/D)、jmax= int (YmaxX K/D) -1nt(YminX K/D)；
其中為東方向的標準圖幅幅數(shù)，Jmax為北方向的標準圖幅幅數(shù)，D為分幅比例尺，K為分幅因子，K=100/L，其中L為分幅板式大??；
(3)依次用標準分幅的圖框構(gòu)建多邊形，多邊形的四角坐標為(XJiXD/K ,Y0+jXD/K )、(X0+(i+l) XD/K，Y0+j XD/K)、(X0+(i+l) XD/K，Y0+(j+l) XD/K )、(X0+(i+l) XD/K，YfjXD/K)，其中i的取值范圍為零到東方向的標準圖幅幅數(shù)imax，j的取值范圍為零到北方向的標準圖幅幅數(shù) jmax，X。= int (XminX K/D) XD/K, Y0 = int (YminX K/D) XD/K；
(4)將標準分幅的多邊形與影像進行空間分析，將落入標準分幅多邊形的的影像進行裁剪；
(5)對裁剪形成的標準分幅影像按照標準分幅命名，后綴東北方向的范圍大小，具體為命名規(guī)則為TE-N-D/2-D/2，其中T為標準圖框分幅的含義，E為標準圖框左下角東方向坐標，N為標準圖框左下角北方向坐標，D為分幅比例尺；
(6)將標準分幅影像導出，導出格式為TIF或JPG，根據(jù)導出的影像，生成相應版本的坐標定位文件。
[0015]所述步驟C的中正確定位并完成正確顯示是通過以下方法實現(xiàn)的:
讀取需輸出的地形圖范圍內(nèi)各點的坐標(X，Y)，計算出對應的圖幅名稱的左下角坐標Xi, Yi,并進行各個標準分幅定位加載處理，其中Xi= int (XXK/D) XD/K, Yi= int (YXK/D) XD/K，對應的影像及坐標文件名稱為TX1-Y1-DA-D/^,其中T為標準圖框分幅的含義，D為分幅比例尺。
[0016]所述的基于高清衛(wèi)星影像的地形圖一體化快速成圖方法還包括以下步驟:
F.對成果地形圖進行保密信息處理。
[0017]所述對成果地形圖進行保密信息處理是這樣進行的:
調(diào)用成果地形圖中的各地物元素信息，判斷其所在圖層，若該地物元素信息處于控制點層，則讀取其編碼進行判斷，如果該編碼包含于保密數(shù)據(jù)庫中，則為保密地物元素，將該地物元素從成果地形圖中刪除；若該地物元素信息不處于控制點層或者該編碼沒有包含在保密數(shù)據(jù)庫中，則保留該地物元素不變。
[0018]以上所述的成圖方法用于生成1:2000地形圖或1:5000地形圖。
[0019]本發(fā)明基于高清衛(wèi)星影像的地形圖一體化快速成圖方法將調(diào)繪、像控、補測等三個外業(yè)步驟融合為一，實現(xiàn)減少外業(yè)工作量的目標，與傳統(tǒng)的航測相比具有周期短的特點，與傳統(tǒng)人工測量相比，其外業(yè)工作量少、投入成本低且周期短。
[0020]本發(fā)明融合處理的彩色數(shù)據(jù)源可采用各類采集平臺、各種年份、各類分辨率的彩色影像，降低數(shù)據(jù)采購門檻和成本，擴大了數(shù)據(jù)利用的范圍，有效利用各式各樣的衛(wèi)星影像資源；對融合后的融合影像進行分類識別，將無法識別的地物、作為像控點的具有重要定位信息的地物，進行外業(yè)數(shù)字化測繪，可增強地形圖對重要地物的重點呈現(xiàn)，也避免了對全圖范圍內(nèi)的外業(yè)測繪，減少了外業(yè)工作量，提高成圖效率；同時，融合后的彩色影像，既具有較高的分辨率，又具有色彩信息，使得地形圖的地物識別和采集效果更好。
[0021]對于標準分幅的裁剪分割，主要通過空間分析實現(xiàn)，適用于從1:500至1:5000比例尺等中小比例尺的標準分幅，其特點在于裁剪效率高，適用性強。
[0022]本發(fā)明還包括對保密元素進行相應處理的步驟，保密元素主要是存放在控制點層的控制點信息，控制點層在繪圖前已創(chuàng)建并初始化為空白層，在生成初始地形圖時將編碼寫入控制點元素并存放于控制點層；對保密元素的處理過程為:調(diào)用成果地形圖中的各地物元素信息，判斷其所在圖層，若該地物元素信息處于控制點層，則讀取其編碼進行判斷，如果該編碼包含于保密數(shù)據(jù)庫中，則為保密地物元素，將該地物元素從成果地形圖中刪除；若該地物元素信息不處于控制點層或者該編碼沒有包含在保密數(shù)據(jù)庫中，則保留該地物元素不變。
[0023]對保密元素進行刪除處理避免了機密信息的泄漏，保證了地形圖使用的穩(wěn)定性與安全性。還包括對特殊元素的處理，特殊元素是指特殊符號、字體和線型，特殊元素必須要有專有符號庫、字體庫和線型庫才能顯示，如丟失符號庫、字體庫和線型庫，則數(shù)據(jù)將出現(xiàn)亂碼，無法正常使用，為避免此問題，本發(fā)明對特殊元素進行了轉(zhuǎn)換，避免了因為客戶符號庫、字體庫和線型庫不同而影響地形圖的使用，具有更好的適用性與針對性。
[0024]本發(fā)明通過將數(shù)字化測圖與衛(wèi)星影像相結(jié)合，實現(xiàn)補測、調(diào)繪、像控一體化，將整個測繪周期大大縮短，完成了以往通過傳統(tǒng)數(shù)字化測圖或航測手段無法完成或難以完成的生產(chǎn)任務，既保證了成圖的質(zhì)量，又保證了成圖的效率提升及成本的降低。
【專利附圖】

【附圖說明】
[0025]圖1為本發(fā)明基于高清衛(wèi)星影像的地形圖一體化快速成圖方法的流程示意圖 圖2為本發(fā)明基于高清衛(wèi)星影像的地形圖一體化快速成圖方法的優(yōu)選流程示意圖 圖3為本發(fā)明實施例1最后得到的Worldviewl衛(wèi)星的彩色融合影像的灰度圖。
【具體實施方式】
[0026]下面結(jié)合實施例詳細闡述本發(fā)明。
[0027]實施例1
本實施例基于高清衛(wèi)星影像的地形圖一體化快速成圖方法，包括以下步驟:
本實施例繪制1:2000地形圖。
[0028]Al.獲取2011年的Worldview2衛(wèi)星0.5m全色影像、2012年的Worldview2衛(wèi)星2m多光譜影像以及2012年的Worldviewl衛(wèi)星0.5m全色影像，對Worldview2衛(wèi)星0.5m全色影像、2012年的Worldview2衛(wèi)星2m多光譜影像進行融合處理,得到Worldview2衛(wèi)星的彩色融合圖像，再對Worldview2衛(wèi)星彩色融合圖像與Worldviewl衛(wèi)星0.5m全色影像進行融合處理，得到Worldviewl衛(wèi)星的彩色融合影像。
[0029]融合步驟具體如下:
(O首先將彩色影像、全色影像統(tǒng)一糾正到一個坐標下；
(2)將彩色影像重采樣，保證彩色影像與全色影像分辨率為1:1;
(3)將彩色影像的RGB波段分離出來；(4)將分離出來的RGB波段影像與全色波段融合，得到彩色融合影像；
A2.對彩色融合影像進行分類識別，將識別得到的具有特定識別特征的地物作為像控點，并將彩色融合影像打印后作為調(diào)繪底圖，對于具有重要地位、影像上遺失的地物及具有重要定位信息的地物進行實地數(shù)字化測繪；
B1.以步驟(A2)所得到的像控點為基準，利用Supermap平臺采用4X4點模式對彩色融合影像進行兩次配準處理，得到配準衛(wèi)星影像；
B2.將配準衛(wèi)星影像與國家標準分幅進行空間分析，并將落在標準分幅圖框內(nèi)的影像裁剪導出，同時以相應的圖幅編號為其命名，具體步驟如下:
(O以配準影像為基準，讀取配準影像的四角坐標Xmin、Ymin、Xmax、Ymax，其中Xmin、Ymin為左下角坐標，Xmax、Ymax為右上角坐標。
[0030](2)根據(jù)四角坐標，計算出配準影像覆蓋的國家標準分幅圖幅，東方向和北方向的標準圖幅幅數(shù)計算公式如下:imax= int (XmaxX K/D) -1nt (XminX K/D)、jmax= int (YmaxXK/D)-1nt (YminX K/D)；
其中為東方向的標準圖幅幅數(shù)，Jmax為北方向的標準圖幅幅數(shù)，D為分幅比例尺，K為分幅因子，K=100/L，其中L為分幅板式大小，當D為500、1000或2000時，L=50厘米;
(3)依次用國家標準分幅的圖框構(gòu)建多邊形，多邊形的四角坐標為(XJiXD/K , Y0+jXD/K )、(X0+(i + l) XD/K，Y0+jXD/K)、(X0+(i + l) XD/K，Y0+(j + 1) XD/K )、(X0+(i+l) XD/K ,YfjXD/K)，其中i的取值范圍為零到東方向的標準圖幅幅數(shù)imax，j的取值范圍為零到北方向的標準圖幅幅數(shù)」_，X0 = int (XminX K/D) XD/K，Y0 = int (YminXK/D) XD/K ；
(4)將標準分幅的多邊形與配準影像進行空間分析，將落入國家標準分幅多邊形的的影像進行裁到;
(5)對裁剪形成的標準分幅影像按照國家標準分幅命名，后綴東北方向的范圍大小，具體為命名規(guī)則為TE-N-D/2-D/2，其中T標準圖框分幅的含義，E為標準圖框左下角東方向坐標，N為標準圖框左下角北方向坐標，D為分幅比例尺；
(6)將標準分幅影像按照需求導出，導出格式為TIF，根據(jù)導出的影像，生成相應版本的坐標定位文件；
其中，導出的坐標定位文件名稱與標準分幅影像同名，后綴名與導出影像格式對應，TIF對應坐標文件后綴為TFW，文件內(nèi)容系統(tǒng)有六個要素:東方向像元大小、旋轉(zhuǎn)量、尺度變化量、北方向像元大小、左下角東方向起始位置、左下角北方向起始位置；
以南寧市六城區(qū)為例，Xmin=2456700、Ymin=475540、Xmax=2556500、Ymax=587730。當分幅比例尺D為2000時，則imax=100, Jmax=I 12,X0= 2456000,Y0 =475000，導出的標準圖框影像為 T475000-2456000-1000-1000.TIF 至 T587000-2556000-1000-1000.TIF 共計 11200 幅。
[0031](Cl)根據(jù)國家標準分幅的圖幅編號對各分幅影像進行定位，具體步驟為:
根據(jù)影像文件名稱，尋找對應坐標文件，讀取左下角東方向起始位置、左下角北方向起始位置、旋轉(zhuǎn)量、尺度變化量、東方向和北方向像元大小，根據(jù)影像的北方向和東方向像元個數(shù)，計算影像的正確定位并完成正確顯示；
C2.正確定位并完成正確顯示是通過以下方法實現(xiàn)的:
讀取需輸出的地形圖范圍內(nèi)各點的坐標(X，Y)，計算出對應的圖幅名稱的左下角坐標Xi;Yi，并進行各個標準分幅定位加載處理，其中Xi= int (XXK/D)XD/K, Yi= int (YXK/D) XD/K，對應的影像及坐標文件名稱為TX1-Y1-DA-D/^,其中T為標準圖框分幅的含義； D.循環(huán)進行步驟Cl與C2，將各單幅影像依次正確定位顯示，生成初始地形圖；
E.以初始地形圖為底圖，根據(jù)外業(yè)測繪成果，將采集的重要與遺漏地物補繪入影像圖中，對于次要地物如土路、植被邊界等以彩色融合影像作為輔助成圖，生成成果地形圖；
Fl.遍歷成果地形圖中的地物元素，判斷各地物元素是否為保密信息元素，并對保密信息元素進行處理；
對保密信息元素的處理方法具體步驟為；
調(diào)用成果地形圖中的各地物元素信息，判斷其所在圖層，若該地物元素信息處于控制點層，則讀取其編碼進行判斷，如果該編碼包含于保密數(shù)據(jù)庫中，則為保密地物元素，將該地物元素從成果地形圖中刪除；若該地物元素信息不處于控制點層或者該編碼沒有包含在保密數(shù)據(jù)庫中，則保留該地物元素不變。
[0032]F2.將需輸出的地形圖需求范圍線作為圍柵，對經(jīng)過Fl步驟處理的成果地形圖進行裁剪，將裁剪出的成果地形圖輸出，具體步驟為；
(O以客戶需求范圍為準，調(diào)用圖庫內(nèi)對應范圍內(nèi)的地形數(shù)據(jù)；
(2)裁剪數(shù)據(jù)前，將特殊數(shù)據(jù)轉(zhuǎn)換為客戶需求的格式，滿足裁剪需求；
(3)以客戶需求范圍為準，將需求范圍線作為圍柵，以剪切模式，對成果地形圖進行裁剪后輸出。
【權(quán)利要求】
1.一種基于高清衛(wèi)星影像的地形圖一體化快速成圖方法，其特征在于包括以下步驟: A、獲取最新的衛(wèi)星全色圖像和各類采集平臺、各種年份、各類分辨率衛(wèi)星影像或者航拍彩色影像，利用各類分辨率的衛(wèi)星或者航拍彩色影像對衛(wèi)星或者航拍全色影像進行融合處理，得到彩色融合影像，對彩色融合影像進行分類識別，將識別得到的具有特定識別特征的地物作為像控點，進行外業(yè)數(shù)字化測繪，并補測彩色融合影像上遺失的重要地物； B、對彩色融合影像進行配準處理，得到配準影像，以該影像為基準，讀取影像的四角坐標，根據(jù)四角坐標，計算出影像覆蓋的國家標準分幅圖幅，將落入國家標準分幅多邊形的的影像進行裁剪，對裁剪形成的單幅影像分幅命名，名稱中包含左下角東方向起始坐標和左下角北方向起始坐標，導出影像，生成坐標定位文件，坐標定位文件名稱與標準分幅影像同名； C、根據(jù)步驟B中命名的單幅影像文件名稱，自動尋找對應坐標文件，讀取左下角東方向起始位置、左下角北方向起始位置、旋轉(zhuǎn)量、尺度變化量、東方向和北方向像元大小，根據(jù)影像的北方向和東方向像元個數(shù)，計算影像的正確定位并完成正確顯示； D、循環(huán)進行步驟C，將各單幅影像依次正確定位顯示，形成初始地形圖； E、以初始地形圖為底圖，根據(jù)外業(yè)測繪成果補繪重要地物和遺漏地物，對于次要地物以衛(wèi)星影像作為輔助成圖，生成成果地形圖。
2.如權(quán)利要求1所述的基于高清衛(wèi)星影像的地形圖一體化快速成圖方法，其特征在于，所述步驟A中的融合處理包括以下步驟: (1)將各類數(shù)據(jù)源的彩色影像與全色影像統(tǒng)一糾正到一個坐標下； (2)將彩色影像重采樣，使得 彩色影像與全色影像分辨率之比為1:1或1:4 ; (3)將彩色影像的RGB波段分離出來； (4)將分離出來的RGB波段影像與全色影像融合，得到彩色融合影像。
3.如權(quán)利要求1所述的基于高清衛(wèi)星影像的地形圖一體化快速成圖方法，其特征在于，所述步驟B中對配準影像的裁剪包括以下步驟: (1)以配準衛(wèi)星影像為基準，讀取配準衛(wèi)星影像的四角坐標Xmin、Ymin、Xmax、Ymax，其中Xmin、Yniin為左下角坐標，Xmax> Ymax為右上角坐標； (2)根據(jù)四角坐標，計算出影像覆蓋的標準分幅圖幅，東方向和北方向的標準圖幅幅數(shù)計算公式如下:imax= int (XmaxX K/D)-1nt (XminX K/D)、jmax= int (YmaxX K/D) -1nt(YminX K/D)； 其中為東方向的標準圖幅幅數(shù)，Jmax為北方向的標準圖幅幅數(shù)，D為分幅比例尺，K為分幅因子，K=100/L，其中L為分幅板式大小； (3)依次用標準分幅的圖框構(gòu)建多邊形，多邊形的四角坐標為(XJiXD/K ,Y0+jXD/K )、(X0+(i+l) XD/K，Y0+jXD/K)、(X0+(i+l) XD/K，Y0+(j+l) XD/K )、(X0+(i+l) XD/K，YfjXD/K)，其中i的取值范圍為零到東方向的標準圖幅幅數(shù)imax，j的取值范圍為零到北方向的標準圖幅幅數(shù) jmax，X。= int (XminX K/D) XD/K, Y0 = int (YminX K/D) XD/K； (4)將標準分幅的多邊形與影像進行空間分析，將落入標準分幅多邊形的的影像進行裁剪； (5)對裁剪形成的標準分幅影像按照標準分幅命名，后綴東北方向的范圍大小，具體為命名規(guī)則為TE-N-D/2-D/2，其中T為標準圖框分幅的含義，E為標準圖框左下角東方向坐標，N為標準圖框左下角北方向坐標，D為分幅比例尺； (6)將標準分幅影像導出，導出格式為TIF或JPG，根據(jù)導出的影像，生成相應版本的坐標定位文件。
4.如權(quán)利要求1所述的基于高清衛(wèi)星影像的地形圖一體化快速成圖方法，其特征在于，所述步驟C的中正確定位并完成正確顯示是通過以下方法實現(xiàn)的: 讀取需輸出的地形圖范圍內(nèi)各點的坐標(X，Y)，計算出對應的圖幅名稱的左下角坐標Xi, Yi,并進行各個標準分幅定位加載處理，其中Xi= int (XXK/D) XD/K, Yi= int (YXK/D) XD/K，對應的影像及坐標文件名稱為TX1-Y1-DA-D/^,其中T為標準圖框分幅的含義，D為分幅比例尺。
5.如權(quán)利要求1所述的基于高清衛(wèi)星影像的地形圖一體化快速成圖方法，其特征在于，還包括以下步驟: F.對成果地形圖進行保密信息處理。
6.如權(quán)利要求5所述的基于高清衛(wèi)星影像的地形圖一體化快速成圖方法，其特征在于，所述對地形圖進行保密信息處理是這樣進行的: 調(diào)用成果地形圖中的各地物元素信息，判斷其所在圖層，若該地物元素信息處于控制點層，則讀取其編碼進行判斷，如果該編碼包含于保密數(shù)據(jù)庫中，則為保密地物元素，將該地物元素從成果地形圖中刪除；若該地物元素信息不處于控制點層或者該編碼沒有包含在保密數(shù)據(jù)庫中，則保留該地物元素不變。
7.如權(quán)利要求1飛任 意一項所述的基于高清衛(wèi)星影像的地形圖一體化快速成圖方法，其特征在于:所述的成圖方法用于生成1:2000地形圖或1:5000地形圖。
【文檔編號】G06T5/00GK103700063SQ201310386156
【公開日】2014年4月2日 申請日期:2013年8月30日 優(yōu)先權(quán)日:2013年8月30日
【發(fā)明者】魏金占, 陳大林, 黃炳強, 莫文通, 王文瑞 申請人:南寧市勘察測繪地理信息院