專利名稱：一種復雜線狀要素的注記方法
技術領域：
本發明涉及地理信息系統技術領域，特別是指一種復雜線狀要素的注記方法。
背景技術：
線狀要素，例如境界線、海岸線、河流、鐵路、道路、等值線、等密度線等，是地理信息系統(Geographic Information System,GIS)中的一類重要地物。現有的線狀要素的注記方式根據地物類型和地圖用途的不同有較大的差別，具體可分為點定位注記配置方式和線定位注記配置方式(曲線注記方式)。其中，曲線注記方式可分為中心線注記配置方式和平行線注記配置方式。例如，道路和街道網一般采用中心線注記配置方式，在其中心線上選取合適的點位進行注記配置；線狀水系的注記一般采用平行線注記配置方式，先求出曲線的平行線，然后在平行線上選取適當的位置進行注記配置；等值線和等密度線則采用點定位注記配置方式，重點在于標注出等值線/等密度線間的拓撲關系，并在等值線/等密度線上選取合適的點位進行注記配置。相對復雜的平行線注記配置方式適用于諸如線狀水系等復雜線狀要素的注記，但是，在采用平行線注記配置方式注記復雜線狀要素時會出現平行線相對于線狀要素突出或回縮的問題，從而造成線狀要素曲線注記的自身沖突。
發明內容
有鑒于此，本發明的主要目的在于提供一種復雜線狀要素的注記方法，以解決平行線注記配置方式注記復雜線狀要素時出現的平行線相對于線狀要素突出或回縮的問題。本發明提供了一種復雜線狀要素的注記方法，包括以下步驟A、定義線狀要素的平行線，包括以下子步驟Al、構建用來定義平行線的參考線；A2、采用角平分線法求出所述參考線的平行線，其中，對于所述參考線的外側平行線作如下處理對于所述參考線上角度小于90°的夾角，利用該夾角平分線上在所述參考線外側的與該夾角的角點距離為一定值的新外側點代替外側平行線上原本與該角點對應的外側點；并在所述新外側點左右各加入一個與所述角點距離同樣為所述定值的輔助點，連接新外側點與左右輔助點構成一線段；最后將該線段與所述外側平行線連接；B、分段處理所述線狀要素的平行線；C、在每段線狀要素的平行線上進行注記定位。由上可以看出，通過上述方法可以解決在采用平行線注記配置方式注記復雜線狀要素時出現的平行線相對于線狀要素突出或回縮的問題，從而避免線狀要素曲線注記的自身沖突。在上述方法中，步驟A中所述參考線的構建包括掃描線狀要素的每一個點，依次以每個點作為角點，求其兩相鄰點與其共同形成的夾角α ；如果α <90°，則計算該角點與其兩相鄰點形成的兩條邊的長度如果這兩條邊的長度都小于d/tan(a ，則刪除該角點；如果在這兩條邊中只有一條邊的長度小于d/ tan( α/2)，則用另一條邊上的距所述角點的距離與較短邊的邊長相同的一點代替該角點， 其中d為參考線與其平行線的間距。由上可以看出，通過上述方法構建的參考線，要比中心線光滑，采用所述參考線適于構建相對于線狀要素不會突出或回縮的平行線。在上述方法中，在所述子步驟A2之后還包括檢查所述平行線中是否存在由相交線段形成的環路若存在，則刪除環路中的各個頂點，將環路起始點作為新頂點加入平行線坐標串中，并修改平行線的頂點數；輸出平行線新的頂點和線段。由上可以看出，采用上述方法可以解決影響注記的平行線自相交的問題。在上述方法中，在所述子步驟A2之后還包括從所述平行線的第一條線段開始，依次判斷每條線段與其后線段的距離是否大于 0且小于一設定值；如果是線段間彼此端點的距離大于0且小于一設定值，則連接所述端點，如果是線段間端點與線段的距離大于0且小于一設定值，則連接該端點和這個端點在其相近線段上的垂足；最后刪除端點間和/或端點與端點在其相近線段上的垂足間的其他點/線段。由上可以看出，采用上述方法可以解決影響注記的平行線自相近的問題。在上述方法中，在所述步驟A與B之間包括將所述線狀要素的平行線劃分為多條基礎線段，即所述線狀要素的平行線上每兩點間的線段；判斷所述線狀要素中是否存在長度大于Lg的基礎線段如果存在，則可以使用平均分割法對這部分基礎線段進行分段處理；其中，Lg為由兩條相鄰基礎線段合并成的合成基礎線段的長度閥值。在上述方法中，所述步驟B包括以下子步驟Bi、計算所述平行線中相鄰基礎線段的相對高度hk,k+1，并記錄下其中最小的相對聞度hmin ；B2、判斷最小的相對高度hmin是否大于相對高度的閥值H。。mp 如果大于H。。mp，則進入下一步；如果小于H。。mp，則合并相對高度等hmin的相鄰基礎線段以形成所述平行線的新曲線段，并求出所述新曲線段的長度Lk,k+1 = Lk+Lk+1 ；其中，如果Lk,k+1 < Ly則使用所述新曲線段替換原來的兩條相鄰基礎線段后，返回步驟Bl ；否則，從剩下的基礎線段中找出最小的相對高度來替換hmin，并重新判斷最小的相對高度hmin是否大于相對高度的閥值H。。mp ；B3、輸出所述平行線的曲線段；當最小的相對高度hmin大于相對高度閥值H。。mp時， 所述曲線段就是所述平行線的基礎線段；當最小的相對高度hmin大于相對高度閥值H。。mp時， 所述曲線段是由兩條相鄰基礎線段合并形成的、且合并長度小于Lg的新曲線段。由上可以看出，采用上述方法可以提取出相對光滑的曲線，相對于傳統的方位角法和單調分段法，上述方法對曲線彎曲度的適應程度較高，而且可以限制曲線段的合并長度，不對能夠進行注記的曲線段進行沒有必要的分段，從而提高了曲線分段的效率。
在上述方法中，所述步驟Bl包括假設所述平行線的相鄰基礎線段的高度為h，即當相鄰的兩條基礎線段合并為一條曲線時，該曲線上某所有點與該曲線兩端點所形成的直線距離中的最大值；設相鄰兩條基礎線段的長度分別為Ik和Ik+1，由該兩條基礎線段合并成的合成曲線兩端點之間的距離為d，則相鄰基礎線段的相對高度hk,k+1 = hd2/(lk2+lk+12)。在上述方法，所述步驟C包括以下子步驟Cl、根據注記字符的長度設定線狀要素注記的長度；C2、掃描線狀要素平行線的所有曲線段，對滿足注記長度的曲線段按照字符個數進行等分，其中，等分點就是注記定位點；C3、按照線注記的規則將每個注記字符放到注記定位點上。在上述方法中，在所述步驟A之前還包括壓縮線狀要素。由上可以看出，對復雜線狀要素的數據進行壓縮可以減少冗余數據對注記動態配置速度的影響。
圖1為本發明提供的一種復雜線狀要素的注記方法流程圖；圖2為角平分線法原理示意圖；圖3為改進的角平分線法原理示意圖；圖4為平行線自相交問題處理前后對照圖；圖5為平行線自相近問題處理前后對照圖；圖6為利用角平分線法求取線狀要素平行線的效果圖；圖7為利用改進的角平分線法求取線狀要素平行線的效果圖。
具體實施例方式下面結合附圖，以曲線線狀要素為例詳細介紹本發明提供的一種復雜線狀要素的注記方法。如圖1所示，所述復雜線狀要素的注記方法包括以下步驟步驟100 壓縮線狀要素。對于一個GIS系統來說，由于受軟硬件條件的限制，數據量的大小及質量將對整個系統的運行速度和效率產生直接的影響。如果GIS系統采用注記動態配置功能模塊對線狀要素進行注記，則需要不斷地動態響應GIS用戶的操作，此舉對系統速度要求較高。因此，精簡的數據對于GIS系統中注記操作來說顯得尤為重要。對于線狀要素的壓縮，一般主要集中在對于線狀要素的空間數據的壓縮上。而對于空間數據的壓縮主要包括無損壓縮和有損壓縮兩種方式。其中，無損壓縮可以在保證數據完整性的情況下一定程度上減少數據存儲量，但壓縮比較低，不能滿足注記配置對于數據壓縮的需要。本實施例在不改變原始數據的基礎上，采用了壓縮比更大的有損壓縮。具體為Douglas數據壓縮算法，包括基本步驟1、設定曲線數據壓縮的閥值。取一個基本上不會影響曲線視覺效果的值作為所述閥值。一般情況下，這個閥值不會太大。2、連接曲線的首尾兩點，并找出與兩點連線有最大偏移距離的點。
3、判斷該點的偏移距離是否小于上述閥值若是，則刪掉所述首尾兩點間所有的點；否則，以該點為分段點，將曲線分成一分為二，并返回執行第2步。由于所述Douglas數據壓縮算法為公知技術，在此不再贅述。在注記線狀要素前，執行線狀要素數據的壓縮不僅可以刪除冗余數據，減少數據存儲量，更重要的是通過壓縮線狀要素數據可以簡化線狀要素的輪廓線，這對于后期線狀要素的注記配置響應速度的提升有很大的幫助。步驟200 定義經過數據壓縮的線狀要素的平行線。由于復雜線狀要素的注記是在平行線上選取適當的位置進行配置的。因此，需要克服一切可能存在的誤差，準確地定義出線狀要素的平行線。以往一般采用如圖2示出的角平分線法來求取線狀要素的平行線。這種算法精度高，運算速度快。但通過這種方法生成的平行線具有如下缺點在中心線夾角α較小時，所求得的平行線夾角頂點(Α、Β)會離中心線非常遠(遠大于線狀要素與其平行線的間距d)。 如果平行線在中心線的外側，則平行線夾角頂點A會因遠離中心線而顯得過分突出；如果平行線在中心線的內側，則平行線夾角頂點B會因遠離中心線而顯得過分收縮。同時，如果中心線夾角α兩邊較短的話，平行線就會因平行線夾角頂點的過分外凸或收縮而改變其原有的形狀。在本實施例中，通過構建用來定義平行線的參考線，來對角平分線法進行改進，從而以解決中心線內側平行線因其平行線夾角頂點向內收縮而改變形狀的問題。如圖3所示，具體為掃描曲線上(相對于其兩側的平行線，也可將曲線稱作中心線)的每一個點，依次以每個點作為角點(是二維圖像亮度變化劇烈的點或曲線上曲率極大值的點)，求其兩相鄰點與其共同形成的夾角α 如果α <90°，則計算該角點與兩相鄰點形成的兩條邊的長度。如果這兩條邊的長度都小于d/tan(a/2)，則刪除該角點。如果在這兩條邊中只有一條邊的長度小于d/tan (α/2),則用另一條邊上的一點代替該角點，該點距所述角點的距離與較短邊(長度小于d/tan(a//2))的邊長相同。其中，d為參考線與其平行線的間距。圖 3示出的角點A、B分別與其相鄰點1、2和3、4構成的夾角α和β均小于90°，且角點Α、 B與其各自的相鄰點1、2和3、4之間形成的邊長均小于d/tan ( α/2),那么刪除角點Α、Β并將點1、2和3、4分別連接后，形成如圖3示出的參考線。最后再利用角平分線法求出所述參考線的內側平行線便可解決平行線夾角頂點向內收縮而改變形狀的問題。對于平行線在參考線的外側時因平行線夾角頂點會因遠離參考線而顯得過分突出的問題，可通過以下方法解決，即對于參考線上角度小于90°的夾角，利用該夾角平分線上在參考線外側的與該夾角的角點距離為d(可根據注記字符大小、間距來設置的定值) 的新外側點代替外側平行線上原本與該角點對應的外側點。然后在所述新外側點左右各加入一個與所述角點距離同樣為d的輔助點，連接新外側點與左右輔助點構成一線段。最后， 將該線段與外側平行線連接，構成新的外側平行線。從而解決了外側平行線因其外側點突出而改變形狀的問題。在完成所述平行線的定義后，如圖4、5所示，平行線可能存在自相交和/或自相近的的情況，由此導致注記重疊、顯示混亂等問題。在本實施例中，假設平行線N個坐標點(由該N個坐標點構成一坐標串)，則其包含的線段數量為N-1。其中，解決平行線自相交問題的具體方法如下如圖4所示，在平行線的N-I條線段中依次檢查所述平行線的是否存在相交線段，即是否存在由相交線段形成的環路若存在，則刪除環路中的各個頂點，將環路起始點a作為新頂點加入平行線坐標串中，并修改平行線的頂點數。最后輸出平行線新的頂點和線段。解決平行線自相近問題的具體方法如下如圖5所示，從平行線的第一條線段開始，依次判斷每條線段與其后線段的關系是否相近(以某一設定值為判斷依據，不大于該設定值且大于0則為相近)；如果是線段間的端點相近，則連接端點，如果是端點與線段相近，則連接該端點和這個端點在其相近線段上的垂足。最后刪除端點間和/或端點與端點在其相近線段上的垂足間的其他點/線段。由此，得到了消除自相近的平行線。下面參照圖6、圖7，對比分別利用角平分線法和本發明改進的角平分線法定義出的線狀要素平行線的效果，其中，Ll為經過數據壓縮的中心線，L2為通過角平分線法求出的Ll的平行線，L3為未經數據壓縮的中心線，“*”是平行線的分段點。如圖6所示，中心線Ll上的點a對應的外側平行線上的點遠離了中心線，且突出到了點b的位置。而中心線Ll上的點a對應的內側平行線上的點回縮到了點c的位置，使內側平行線的形狀發生了改變。而圖6中點d、e示出的位置則發生了平行線自相交的情況， 如果該平行線與中心線Ll的間距再小一點的話，那么點e示出的位置則會發生平行線自相近的情況。因為復雜線狀要素的注記是在平行線上選取適當的位置進行配置的，即注記字符的中心就在平行線上，在圖6示出的雜亂的平行線上顯然是很難進行注記。而圖7示出的由本發明改進的角平分線法定義的規整的平行線是很容易進行注記的。步驟300 分段處理所述線狀要素的平行線。對于長度較長的線狀要素，在注記之前需要對其平行線進行分段處理，然后再對每一段或其中的幾段進行注記，即在每一段或其中的幾段線狀要素的平行線上進行注記。 以曲線為例，現有的曲線分段方法主要有兩種一種是方位角法，具體為設定一個方位角變化閥值α ；計算曲線上每一條由曲線坐標點形成的線段的方位角；接著計算位于每個坐標點前后的線段的方位角之差，并將方位角之差的絕對值大于α的所有坐標點記錄下作為分段點；最后根據所述分段點對曲線進行分段。這種方法在一定程度上可以將曲線分成較為光滑的曲線段，但是對于一些復雜的曲線，特別是一些方位角變化比較大的曲線，容易將其劃分得太細，以至于無法進行注記。另一種是單調分段法，具體為提取出曲線中所有單調性發生變化的坐標點作為初始分段點；根據初始分段點對曲線進行分段，并提取出曲線單調段(由相鄰的且彼此間能夠同時投影在同一坐標軸上的曲線段合并而成)作為曲線的最終分段。這種方法可以保證曲線注記的通視性，但是不能將曲線劃分出相對平滑的曲線段。為了解決上述現有曲線分段方法存在的問題，本實施例基于逆向Douglas算法， 以曲線上所有的線段為基礎，逐步合并相對光滑的曲線段，并以最終輸出的曲線段作為注記曲線段。在利用逆向Douglas算法分段處理曲線之前，將所述曲線劃分為多條基礎線段 (曲線上每兩點間的線段)，設基礎線段的長度為lk。根據注記字符串的長度和字符的大小確定由兩條相鄰的基礎線段合并成的合成基礎線段(合成曲線)的相對高度的閥值為 H。。mp。設合成曲線的長度閥值為Ly判斷所述曲線中是否存在長度大于Lg的基礎線段如果存在，則可以使用平均分割法對這部分基礎線段進行分段處理；否則采用逆向Douglas算法分段處理所述曲線。所述逆向Douglas算法包括以下步驟步驟301 計算曲線中相鄰基礎線段的相對高度hk, k+1，并記錄下其中最小的相對高度hmin。其中，假設相鄰基礎線段的高度為h，即當相鄰的兩條基礎線段合并為一條曲線時，該曲線上某所有點與該曲線兩端點所形成的直線距離中的最大值；設相鄰兩條基礎線段的長度分別為Ik和lk+1，由該兩條基礎線段合并成的合成曲線兩端點之間的距離為d，則相鄰基礎線段的相對高度hk,k+1 = hd2/(lk2+lk+12)。步驟302 判斷最小的相對高度hmin是否大于相對高度的閥值H。。mp 如果大于H。。mp， 則進入下一步；如果小于H。。mp，則合并相對高度等于hmin的相鄰基礎線段以形成新曲線段， 并求出新曲線段的長度Lk, k+1 = Lk+Lk+1。其中，如果Lk,k+1 < Ly則使用新曲線段替換原來的兩條相鄰基礎線段后，返回步驟301 ；否則，從剩下的基礎線段中找出最小相對高度來替換hmin，并重新判斷最小的相對高度hmin是否大于相對高度的閥值H。。mp。步驟303 輸出曲線段。當最小的相對高度hmin大于相對高度閥值H。。mp時，所述曲線段就是所述曲線的基礎線段。當最小的相對高度hmin大于相對高度閥值H。。mp時，所述曲線段是由兩條相鄰基礎線段合并形成的、且合并長度小于Lg的新曲線段。步驟400 對每段線狀要素進行注記定位，即確定每段線狀要素數據上注記點的坐標以及該坐標上注記字符的字向。具體為首先根據字符的長度設定線狀要素注記的長度；掃描線狀要素的所有曲線段，對滿足注記長度的曲線段按照字符個數(注記點個數)進行等分，例如，注記有m個注記點，對滿足注記長度的曲線段按照其長度進行m+1等分，其中，等分點就是注記定位點；然后按照線注記的規則將每個注記字符放到注記定位點上；線注記應遵循以下規則1、沿線狀要素的平行線，一般采用離散的注記方式，屈曲字列，字向(字頭的朝向)與曲線的切線方向垂直；2、注記的字間距為一定范圍的等間隔；3、注記的排列方向沿著被注記線狀要素的方向，依據其傾角而定，采用從上往下， 從左往右注記的原則；4、注記不允許放在線狀要素的兩側；5、注記與線狀要素的距離要適當；6、若線狀要素較長，應分段重復注記；7、避免注記靠近線要素的端點，注記首末字到線狀要素兩端點的距離應相等且保持適當距離；8、避免壓蓋地圖要素，不壓蓋其他注記。最后校正已放到注記定位點上的注記字符的字向。例如，可通過縮短注記字符間距來校正注記首尾字符距離曲線段兩端點較近的問題。如果注記首尾字符距離曲線段兩端點過近，則可在調整注記字符間距的同時沿曲線段形成足夠長的折線，注記字符字向與所述折線兩端點連接線平行或垂直。以上所述僅為本發明的較佳實施例而已，并不用以限制本發明，凡在本發明的精神和原則之內，所作的任何修改、等同替換、改進等，均應包含在本發明的保護范圍之內。
權利要求
1.一種復雜線狀要素的注記方法，其特征在于，包括以下步驟A、定義線狀要素的平行線，包括以下子步驟 Al、構建用來定義平行線的參考線；A2、采用角平分線法求出所述參考線的平行線，其中，對于所述參考線的外側平行線作如下處理對于所述參考線上角度小于90°的夾角，利用該夾角平分線上在所述參考線外側的與該夾角的角點距離為一定值的新外側點代替外側平行線上原本與該角點對應的外側點；并在所述新外側點左右各加入一個與所述角點距離同樣為所述定值的輔助點，連接新外側點與左右輔助點構成一線段；最后將該線段與所述外側平行線連接；B、分段處理所述線狀要素的平行線；C、在每段線狀要素的平行線上進行注記定位。
2.根據權利要求1所述的方法，其特征在于，步驟A中所述參考線的構建包括掃描線狀要素的每一個點，依次以每個點作為角點，求其兩相鄰點與其共同形成的夾角ct ；如果α <90°，則計算該角點與其兩相鄰點形成的兩條邊的長度如果這兩條邊的長度都小于d/tan(a/2)，則刪除該角點；如果在這兩條邊中只有一條邊的長度小于d/ tan( α/2)，則用另一條邊上的距所述角點的距離與較短邊的邊長相同的一點代替該角點， 其中d為參考線與其平行線的間距。
3.根據權利要求2所述的方法，其特征在于，在所述子步驟A2之后還包括檢查所述平行線中是否存在由相交線段形成的環路若存在，則刪除環路中的各個頂點，將環路起始點作為新頂點加入平行線坐標串中，并修改平行線的頂點數；輸出平行線新的頂點和線段。
4.根據權利要求2所述的方法，其特征在于，在所述子步驟A2之后還包括從所述平行線的第一條線段開始，依次判斷每條線段與其后線段的距離是否大于0且小于一設定值；如果是線段間彼此端點的距離大于0且小于一設定值，則連接所述端點，如果是線段間端點與線段的距離大于0且小于一設定值，則連接該端點和這個端點在其相近線段上的垂足；最后刪除端點間和/或端點與端點在其相近線段上的垂足間的其他點/線段。
5.根據權利要求1所述的方法，其特征在于，在所述步驟A與B之間包括將所述線狀要素的平行線劃分為多條基礎線段，即所述線狀要素的平行線上每兩點間的線段；判斷所述線狀要素中是否存在長度大于Lg的基礎線段如果存在，則可以使用平均分割法對這部分基礎線段進行分段處理；其中，Lg為由兩條相鄰基礎線段合并成的合成基礎線段的長度閥值。
6.根據權利要求5所述的方法，其特征在于，所述步驟B包括以下子步驟Bi、計算所述平行線中相鄰基礎線段的相對高度hk, k+1，并記錄下其中最小的相對高度h . ·11Iiiin，B2、判斷最小的相對高度hmin是否大于相對高度的閥值H。。mp 如果大于H。。mp，則進入下一步；如果小于H。。mp，則合并相對高度等于hmin的相鄰基礎線段以形成所述平行線的新曲線段，并求出所述新曲線段的長度Lk,k+1 = Lk+Lk+1 ；其中，如果Lk,k+1 < Ly則使用所述新曲線段替換原來的兩條相鄰基礎線段后，返回步驟Bl ；否則，從剩下的基礎線段中找出最小的相對高度來替換hmin，并重新判斷最小的相對高度hmin是否大于相對高度的閥值H。。mp ；B3、輸出所述平行線的曲線段；當最小的相對高度hmin大于相對高度閥值H。。mp時，所述曲線段就是所述平行線的基礎線段；當最小的相對高度hmin大于相對高度閥值札_時，所述曲線段是由兩條相鄰基礎線段合并形成的、且合并長度小于Lg的新曲線段。
7.根據權利要求6所述的方法，其特征在于，所述步驟Bl包括假設所述平行線的相鄰基礎線段的高度為h，即當相鄰的兩條基礎線段合并為一條曲線時，該曲線上某所有點與該曲線兩端點所形成的直線距離中的最大值；設相鄰兩條基礎線段的長度分別為Ik和lk+1， 由該兩條基礎線段合并成的合成曲線兩端點之間的距離為d，則相鄰基礎線段的相對高度hk, k+i = hd2/(lk2+lk+12) ο
8.根據權利要求1所述的方法，其特征在于，所述步驟C包括以下子步驟 Cl、根據注記字符的長度設定線狀要素注記的長度；C2、掃描線狀要素平行線的所有曲線段，對滿足注記長度的曲線段按照字符個數進行等分，其中，等分點就是注記定位點；C3、按照線注記的規則將每個注記字符放到注記定位點上。
9.根據權利要求1所述的方法，其特征在于，在所述步驟A之前還包括壓縮線狀要
全文摘要
本發明提供了一種復雜線狀要素的注記方法，包括步驟A、定義線狀要素的平行線，包括構建用來定義平行線的參考線；采用角平分線法求出所述參考線的平行線，其中，對于所述參考線的外側平行線作如下處理對于所述參考線上角度小于90°的夾角，利用該夾角平分線上在所述參考線外側的與該夾角的角點距離為一定值的新外側點代替外側平行線上原本與該角點對應的外側點；并在所述新外側點左右各加入一個與所述角點距離同樣為所述定值的輔助點，連接新外側點與左右輔助點構成一線段；最后將該線段與所述外側平行線連接；B、分段處理所述線狀要素的平行線；C、在每段線狀要素的平行線上進行注記定位。以解決平行線相對于線狀要素突出或回縮的問題。
文檔編號G09B29/00GK102237010SQ20111021443
公開日2011年11月9日 申請日期2011年7月28日 優先權日2011年7月28日
發明者丁圣陶, 劉曉麗, 印潔, 葉關根, 孫偉, 張成成, 方馳宇, 李成名, 梁西騰, 殷勇, 沈濤, 焦孟凱, 趙占杰, 趙園春, 路文娟 申請人:中國測繪科學研究院