點狀資訊KML檔之撰寫(csv2kml.py)

Table of contents

背景
- 前言
KML格式及內涵
座標轉換計算
多邊形的應用
- 基本
- 多邊形的特殊約定
程式說明
Reference

背景

點狀資訊在KML架構中是最單純的。KML提供了點狀性質(feature)的設定方式，因此使得單調的點，因著背景地圖的參考而變得豐富。
[csv2kml.py]的應用有很多，繪製軌跡圖、將點狀(與鏈結)資訊輸出到uMAP展現、貼在openTopo底圖(NCLonOTM)、ISC/AERMOD地理資訊的解讀展示等等。
如果要建立更複雜的KML，由熟悉點狀資訊的輸出，會是一個好的開始。

前言

csv是常見的資料表格式，KML則是google提供的語言介面，好讓使用者也可以在google map上貼圖，現已成為開放式地理空間協會地理資訊的編碼標準。因此csv2kml顧名思義就是嘗試把一個存有座標的csv檔案，寫成kml檔案，最終能貼在google map、openstreet map等網路地圖上，完成這項工作的python小程式。(成果如附圖)


csv2kml結果範例

輸入部分重點在於應用pandas模組來開啟csv檔案，呼叫csv的內容。其次則為處理命令列選項的模組argparse，這2者都是python資料處理常用的模組。
輸出部分是本程式的重點。主要參考KML官網的介紹，將要繪製的內容、點、線等元素，依據KML的順序寫出，以符合google map程式的規則。KML(Keyhole Markup Language)及其zip後的檔案格式KMZ，是google map/google earth平台上特有的檔案格式，詳細介紹可以參考官網及網友的說明。大致上KML語法和html相同，除了一般性的格式之外，在地圖上標示的元件以…做為起迄標示。而在本作業中乃是以LineString、Point 為繪圖形式，此外KML也有Path、Polygon(shaded色塊) 等形態。
處理計算部分google map用的是經緯度系統，台灣地區常用的是TWD97系統，如果使用後者，需要進行座標轉換。此處使用twd97模組。

KML格式及內涵

KML的內容可以參考其教學網站。如果只需要輸出單一分開的功能，可以考慮simplekml模組。範例圖檔(紅色部分)為如下共276行之kml檔案，kml並不需要跳行，純粹是為了閱讀解釋方便才加上跳行指令(‘\n’)，其內容說明如下：
起始
在此說明整個檔案會用到的宣告，包括
1. 開啟kml 及Document之宣告
2. 名稱name(會出現在google map的圖名標籤)(line 1)
3. 敘述說明descri(line 1)
4. 可能會用到的圖像定義Icon(line2~6)，是用style id 來定義並呼叫，範例呼叫了紅點reddotPlacemark。
icon選擇的考量：
- 大小：icon當點數增加的時候，icon會重疊在一起，點數太少或icon太小，icon會太分散不明顯。
- 顏色：具有鑑別能力意義：讓人一目瞭然其代表的意思
  點位段落
- 這一段(line7~184)繪圖範例圖檔中的「點」(Point)。
- 由於每一個點呼叫的style id可以不一樣，因此設計成每筆呼叫。 1. 每個點有獨立的座標coordinates，先給經度、再給緯度、高度(範例為0) 2. 各點的順序沒有特別的意義 3. 各點上可以有名稱(name)及敘述(desc)，範例為距離原點的小時數(hour=0 )、以及當時的年月日時標籤(ymd=2019062315 )

$ cat -n trjguanshan2019062315.csv.kml
<?xml version="1.0" encoding="UTF-8"?><kml xmlns="http://www.opengis.net/kml/2.2"><Document><name>trjguanshan2019062315.csv</name><descri ption>trjguanshan2019062315.csv</description>
<Style id="normalPlacemark"><IconStyle><Icon><href>http://maps.google.com/mapfiles/kml/paddle/</href>wht-blank.png</Icon> </IconStyle></Style>
<Style id="highlightPlacemark"><IconStyle><Icon><href>http://maps.google.com/mapfiles/kml/paddle/red-stars.png</href> </Icon> </IconStyle></Style>
<Style id="dotPlacemark"><IconStyle><Icon><href>http://maps.google.com/mapfiles/kml/pal4/icon57.png</href> </Icon> </IconStyle></Style>
<Style id="blackdotPlacemark"><IconStyle><Icon><href>http://maps.google.com/mapfiles/kml/pal4/icon25.png</href> </Icon> </IconStyle></Style>
<Style id="reddotPlacemark"><IconStyle><Icon><href>http://maps.google.com/mapfiles/kml/pal4/icon49.png</href> </Icon> </IconStyle></Style>
<Placemark><name>hour=0</name><description>ymd=2019062315</description><styleUrl>#reddotPlacemark</styleUrl><Point><coordinates>
121.16166338764377,23.04499756226396,0</coordinates></Point></Placemark>
<Placemark><name>hour=1</name><description>ymd=2019062314</description><styleUrl>#reddotPlacemark</styleUrl><Point><coordinates>
121.16052210969468,23.06644148767801,0</coordinates></Point></Placemark>
<Placemark><name>hour=2</name><description>ymd=2019062313</description><styleUrl>#reddotPlacemark</styleUrl><Point><coordinates>
121.15857390654217,23.090818981273667,0</coordinates></Point></Placemark>
<Placemark><name>hour=3</name><description>ymd=2019062312</description><styleUrl>#reddotPlacemark</styleUrl><Point><coordinates>
121.17029256991134,23.110768605923834,0</coordinates></Point></Placemark>
...
<Placemark><name>hour=87</name><description>ymd=2019062000</description><styleUrl>#reddotPlacemark</styleUrl><Point><coordinates>
121.28943947883455,23.47345566850706,0</coordinates></Point></Placemark>
<Placemark><name>hour=88</name><description>ymd=2019061923</description><styleUrl>#reddotPlacemark</styleUrl><Point><coordinates>
121.28348218711278,23.475183559469606,0</coordinates></Point></Placemark>
<Placemark> <LineString> <coordinates>
121.16166338764377,23.04499756226396,0
121.16052210969468,23.06644148767801,0
121.15857390654217,23.090818981273667,0
121.17029256991134,23.110768605923834,0
121.19588315665499,23.134448977175694,0
121.23013239723322,23.161453958196724,0
...
121.28348218711278,23.475183559469606,0
</coordinates></LineString></Placemark>
</Document></kml>

直線段落

第3段(line185~275)為範例圖檔中的「線」(LineString)。「線」指令會將給定座標的第1點開始，按照順序連接所有的座標點，直到最後1點。

座標coordinates，先給經度、再給緯度、高度(範例為0)
各線段座標之間沒有分隔符號，空格或’\n’皆可，此處為’\n’
線段不能加註名稱或敘述
結尾
完成Document及kml之宣告。

df=read_csv(fname)#,encoding='big5')
col=df.columns
for i in range(len(df)):
  nam=df.loc[i,col[2]]
  desc=df.loc[i,col[3]]
  if TWD97:
    x,y=df.loc[i,col[0]],df.loc[i,col[1]]
    lat,lon=twd97.towgs84(x,y)
  if LL:
    lon,lat=df.loc[i,col[0]],df.loc[i,col[1]]

座標轉換計算

twd97
- 網友公開的twd97模組，其計算結果與學術單位網站服務所差小於1M
- 在119~121範圍內有效，其外結果則不保證。

多邊形的應用

基本

同樣2點以上的數據，可以用點(Point)或線(Line)形式來繪製KML的地標(Placemark)
4點以上，如果最後點與第一點重疊，則可以選擇以多邊形來展現地標，多邊形的特性就是除了線條之外，也可以針對填滿的顏色、透明度等進行設定。
等值圖KML檔之撰寫就是以具有順序的多邊形、來層疊出等值圖的特性。

多邊形的特殊約定

由於多邊形物件不限定點數，因此必須（至少）約定物件的開始點。
此處在df.desc內容中約定含有’p0’字串者，為該物件的第一筆，如此就可以正確切割csv檔案成為所要的多邊形群組，同時保有原來csv檔的單純性。

...
else: #Polygon case
  headP = '</styleUrl><Polygon><outerBoundaryIs><LinearRing><tessellate>1</tessellate><coordinates>'
  tailP = '</coordinates></LinearRing></outerBoundaryIs></Polygon></Placemark>'
  ps=[i[-2:] for i in list(df[col[3]])]
  npoly,i0,i1=1,[0],[len(df)]
  if 'p0' in ps:
    npoly=ps.count('p0')
    if npoly>1:
      i1=[]
      for i in range(npoly-1):
        i0.append(i0[-1]+1+ps[i0[-1]+1:].index('p0'))
        i1.append(i0[-1])
      i1.append(len(df))
  for ip in range(npoly):
#   level=int(np.random.rand()*10)
    level=5
    nam=df.loc[i0[ip],col[2]].replace('p0','')
    desc=df.loc[i0[ip],col[3]].replace('p0','')
    line.append('<Placemark><name>' + nam + '</name><description>'+desc+'</description><styleUrl>#level' + str(level) + headP)
    print (i0[ip],i1[ip])
    for i in range(i0[ip],i1[ip]):
      line.append(lonlat[i])
    line.append(tailP)
...

程式說明

csv2kml.py 程式碼

github.com/sinotec2

調用模組與引數之讀取

在命令列鍵入引數， python常用的讀取方式有2種，一者為簡易的sys.argv指令，另外最常使用的則是其內設的引數解析argparse模組。前者沒有help、不會解析引數的型態，一律為字串。相對而言，後者的功能就較為完整、足供應用了。使用argparser讓引數的種類可以多元化

本程式所需的引數有3個，分別為：

-f(csv檔案名稱)、
-n(點線型態的選擇、有N白色地點、H紅色地點、D細線圓點、R紅色圓點、B黑色暈點，再加L繪線、G則告訴程式csv的內容為多邊形 )
-g座標系統有TWD97及LL 等2個選項
- LL為經緯度
- TWD97為台灣地區2度分帶97年基準座標

程式語法以add_argument()來接收命令列訊息。

必要的選項有2，選項代碼(短名)、選項長名、
其他補充的選項包括是否必須(required)、引數的類型、以及說明內容。
在呼叫時以3個元素的tuple來承接(line 16)，其順序即為檔名、點線特性、以及座標系統種類。

$ cat -n ~/bin/csv2kml.py
#!/cluster/miniconda/envs/py37/bin/python
from pandas import *
import twd97,sys
     4
     5
def getarg():
    """ read the setting of plot from argument(std input)"""
    import argparse
    ap = argparse.ArgumentParser()
    ap.add_argument("-f", "--fname", required = True, type=str,help = "csv xy data")
    ap.add_argument("-n", "--NorH", required = True, type=str,help = "Normal/Highlight/Dot/Reddot/Blackdot w/o Line")
    ap.add_argument("-g", "--GEOG", required = True, type=str,help = "LL or TWD9")
    args = vars(ap.parse_args())
    return args['fname'],args['NorH'],args['GEOG']
    15

定義座標種類

標示方式可以是
- H：Hightlight Baloon
- N：Normal Baloon
- D: Dots
- R: Red Dots
- B: Black Corona

fname,NorH,geog=getarg()
NorH,geog=NorH.upper(),geog.upper()
LL,TWD97=False,False
if geog=='LL':
  LL=True
if geog=='TWD97':
  TWD97=True
if NorH[0] not in ['H', 'N', 'D','R','B'] :sys.exit('NorH not right')
NH=NorH.replace('L','')
    25
#B:target with black arona
#R:red target
#D:dot target
stl={'N':'normalPlacemark','H':"highlightPlacemark",'D':"dotPlacemark",'B':"blackdotPlacemark",'R':"reddotPlacemark"}
#csv with 4 columns,
#xp,yp,Hour,ymdh
#266567.0,2549275.0,hour=0,ymd=2019062315

pandas之應用

python讀取ASCII資料一般也有2種方式，一者開啟檔案後直接讀取(with open(fname,’r’) as file:;data=[l for l in file])，讀出來的內容一律是字串(含跳行指令’\n’)，還必須另外分段、解析、給予型態定義。
另一方式則以pandas來讀取。在解析上方便太多。以一般具有表頭(header)、逗點分隔(delimitor)的csv檔而言，讀取指令最簡單如line 33所示。如果不是，也可以在一行內解決設定。
python 3必須指定coding，python 2 則不需要
pandas的資料結構(dataframe)具有陣列處理的功能，具有平行計算的能力，然本範例中仍然當成序列循續處理，也有淺顯易懂的好處。
df.columns：讀取表頭成為一序列
df.loc[i,j]：讀取資料表中i列(i為數字), j欄的內容(前述表頭欄位名稱之一)。

    33  df=read_csv(fname)#,encoding='big5')
    34  TITLE=fname

KML檔頭與標點

檔頭除了標示kml版本之外，也宣告了標點、線條的樣式

head0='<?xml version="1.0" encoding="UTF-8"?><kml xmlns="http://www.opengis.net/kml/2.2"><Document>'+ \
'<name>'+TITLE+'</name><description>'+TITLE+'</description>'
#store the lines
line=[head0]
    39
#icons may be used in KML file:
sty_pngs={'H':'paddle/red-stars.png','N':'paddle/white-stars.png','B':'pal4/icon25.png','D':'pal4/icon57.png','R':'pal4/icon49.png',}
for s in stl:
  stl_url='<Style id="'+stl[s]+'"><IconStyle><Icon><href>http://maps.google.com/mapfiles/kml/'+sty_pngs[s]+'</href> </Icon> </IconStyle></Style>'
  exec('styl'+s+"='"+stl_url+"'")
  exec('line.append(styl'+s+')')
    46

計算經緯度

使用twd97模組
輸出時經緯度及高度間以逗點隔開

col=df.columns
lonlat=[]
#plot the placemarks
for i in range(len(df)):
  nam=df.loc[i,col[2]]
  desc=df.loc[i,col[3]]
  sturl= '<styleUrl>#'+stl[NH]+'</styleUrl>'
  line.append('<Placemark><name>'+nam+'</name><description>'+desc+'</description>'+sturl+'<Point><coordinates>')
  if TWD97:
    x,y=df.loc[i,col[0]],df.loc[i,col[1]]
    lat,lon=twd97.towgs84(x,y)
  if LL:
    lon,lat=df.loc[i,col[0]],df.loc[i,col[1]]
  slonlat=str(lon)+','+str(lat)+',0'
  line.append(slonlat+'</coordinates></Point></Placemark>')
  lonlat.append(slonlat)
    63

附加軌跡線數據及檔尾

關鍵詞為
順序是軌跡線的經度、緯度、高度(0)，跳行指令不一定需要。

#if with L in NorH, connection the marks with line
if 'L' in NorH:
  line.append('<Placemark> <LineString> <coordinates>')
  for i in range(len(df)):
    line.append(lonlat[i]+'\n')
  line.append('</coordinates></LineString></Placemark>')
line.append('</Document></kml>')

輸出KML檔案

基本上KML是一個文字檔，因此使用open.write指令即可將其寫出

    71  with open(fname+'.kml','w') as f:
    72    [f.write(l) for l in line]

Reference

名詞解釋
- ASCII(American Standard Code for Information Interchange，美国信息交换标准代码) https://zh.wikipedia.org/wiki/ASCII
- KML(Keyhole Markup Language Keyhole公司標記語言),wiki
- TWD97, 大地座標系統漫談, sunriver.com
argparse
- Argparse 教學，官網
- 命令行选项、参数和子命令解析器,python官網
pandas
- Pandas 基礎教學
- 10 minutes to pandas
- 資料科學家的 pandas 實戰手冊：掌握 40 個實用數據技巧
KML and icons
- KML Samples - Google Developers
- KML Tutorial
- Google+Map讀取KML檔
- 在 Google 地球中匯入 KML 地圖資料 google support
- simplekml, pypi ;readthedoc
- google提供的icons
twd97
- pypi ; github