Ruby~Python^gdal GeoTIFF

はじめに

コンピュータを使いはじめて、30年になる。
初めは、ash-awk(2年)、次にksh-awk(3年)、次にperl(10年)、そして、ruby(15年)を使ってきた。
昨年からbs4を使うため、pythonも始めた。
でも、rubyの書きやすさは、pythonよりもずっと好きだ。

しかし、gdalライブラリを使いたくなった。(とりあえず、geotiffを解釈したいだけだが。)
すると、rubyより、pythonのドキュメントが圧倒的に豊富なのだ～(残念)

PYTHON　anaconda-gdalのインストール

で、pythonをまた使うことに。しかし、python の使い方もいろいろだ。

前は、pyenvとかpipとか使ったが、
今度は、anacondaを使う。なんか、複雑だ。(後で思えば、検索するとそんなサイトが出てきただけで、pipでもやれたのかもしれない。)
①Python3×地理空間データ 地理空間データプログラミングの流れ | Mapping Life
②OSGeo4W 64ビット版 (OSGeo4W64) のインストール（Windows 上）
③【超初歩】AnacondaにGDALをインストールしてみた - Qiita
④Anaconda

いろんなサイトを見て、anaconda-pronpt(DOSプロンプトに環境変数入れただけ)で、
gdalinfoを使うために、libopenssl-3-x64.dll、libcrypto-3-x64.dllを別途入れて、
なんとかたどり着いた。半日以上かかってしまった。
Aw3d-30mのデータの例を下記する。

???\Aw3d-30m\id\N000E110\ALPSMLC30_N000E110_DSM.tif
Driver: GTiff/GeoTIFF
Files:???\Aw3d-30m\id\N000E110\ALPSMLC30_N000E110_DSM.tif
Size is 3600, 3600
Coordinate System is:
GEOGCRS["WGS 84",
    DATUM["World Geodetic System 1984",
        ELLIPSOID["WGS 84",6378137,298.257223563,
            LENGTHUNIT["metre",1]]],
    PRIMEM["Greenwich",0,
        ANGLEUNIT["degree",0.0174532925199433]],
    CS[ellipsoidal,2],
        AXIS["geodetic latitude (Lat)",north,
            ORDER[1],
            ANGLEUNIT["degree",0.0174532925199433]],
        AXIS["geodetic longitude (Lon)",east,
            ORDER[2],
            ANGLEUNIT["degree",0.0174532925199433]],
    USAGE[
        SCOPE["unknown"],
        AREA["World"],
        BBOX[-90,-180,90,180]],
    ID["EPSG",4326]]
Data axis to CRS axis mapping: 2,1
Origin = (110.000000000000000,1.000000000000000)
Pixel Size = (0.000277777777778,-0.000277777777778)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_XRESOLUTION=1
  TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( 110.0000000,   1.0000000) (110d 0' 0.00"E,  1d 0' 0.00"N)
Lower Left  ( 110.0000000,   0.0000000) (110d 0' 0.00"E,  0d 0' 0.01"N)
Upper Right ( 111.0000000,   1.0000000) (111d 0' 0.00"E,  1d 0' 0.00"N)
Lower Right ( 111.0000000,   0.0000000) (111d 0' 0.00"E,  0d 0' 0.01"N)
Center      ( 110.5000000,   0.5000000) (110d30' 0.00"E,  0d30' 0.00"N)
Band 1 Block=3600x1 Type=Int16, ColorInterp=Gray
<||

(open-sslサイトから、lightの64bit版をインストールすれば、簡単だったのだろうが、
誤って、32bit版を入れたのが大失敗か？)

※OSGeoサイトからのインストールも並行してやってしまったが、無意味だったようだ。
※anaconda,conda install -c conda-forge gdal 

cygwin にaucondaの設定を適用しようとするとまた工夫が必要なようだ。
私はcygwin 上で python のgdalライブラリ使いたいのだ。
上記URLの④が、使えそうだ。