2024年6月18日 星期二

Some hint of run netcdf fortran at liunx

 不同版本linux distribution有所差異,自己記錄debug方向

1. yum or dnf 安裝

    dnf install netcdf netcdf-fortran netcdf-devel

2. 範例, 路徑自己改

fortran my_prog.f -o my_prog -I/usr/include -L/usr/lib64 -lnetcdff

2024年5月12日 星期日

Plot Shape file by GMT5 at CentOS 7.9

 因為工作上的需求,要把某些類型的shape file(如縣市邊界)資料用GMT繪製出來,

目前使用的Linux Distribution 是CentOS 7.9(我知道2024年中就EOL了,還沒空換新 Distribution)跟GMT 5(我知道這也舊了),先記錄一下,更新版本問題之後再說。

一開始查資料時,GMT好像無法直接畫Shape file這種向量(Vector)資料,需轉成光柵(Raster)格式,而不管是胡植慶老師的講義或是網路上的講法都是用ogr2org這個指令轉檔,但這不是GMT指令,需要安裝GDAL這個轉檔程式。而根據網路上的說法(多是5、6年前),要用source code compile proj跟GDAL,還要注意缺的套件,所幸到此時,yum install proj gdal就搞定!

接下來就是測試轉檔,這是用政府開放資料平台的縣市界線的shape file測試,為保證我抓下來的檔案是正確的,先用QGIS開檔。需注意要.shp、.shx跟.dbf一起存在,不然會開啟失敗。


嗯!在windows下用QGIS看起來正常,接下來將.shp、.shx跟.dbf一起傳到Linux上,目前因為已經裝好GDAL,所以可以用ogr2ogr。

用法範例 ogr2ogr -f  "GMT" test.gmt test.shp

接下來就可以用GMT畫出來了,成功!