2012年10月20日土曜日

GPS トラックデータ

SimpleDEMViewer4.1でGPX形式のGPSデータをサポートしました。テストのために使わせていただいたデータはしまね自然歩道クラブのサイトのものです。目的は達したし、使わせていただいておいて文句を言う筋合ではないのですが、上の図のようにおかしなデータになっています。幾つかのログを接なぎ合せて一つのデータにしていると思うのですが、経路順になっていません。途中で突然ジャンプして逆走したりしています。データのテストはどうしているのでしょうね? テストしたプログラムでは線としての表示はせずに、点を並べるだけなのでしょうか。GPSへアップロードするルートデータを意識しているようですが、GPSって点で表示するものが普通なのでしょうか、そうとは思えないのですが? 私はGPSを使用したことがありません(^_^;
 地形は基盤地図情報の10mメッシュです。

5 件のコメント:

  1. 数値地図ビューアのユーザーです。SimpleDEMViewer4.1早速、使ってみました。GPSデータのサポートはとても便利だと思いました。それで、期待と質問なのですが、経路を描画した状態の画像をTIFFで保存する場合にワールドファイルを書き出す。にチェックしますと、TFWを書き出してくれますが、これをGEOTIFFとして一つの画像ファイルで書き出すようなことはできないでしょうか?

    返信削除
  2. コメントに気付かず遅くなりました。
     GeoTiffですか。実は「基盤地図標高変換」でもGeoTiffの方がファイル1個でいいよな、と思いながら実現できていません。どうやって書込んだらよいかまだ理解できていないのです。画像ファイルともなると、さらにハードルが上がります。一応念頭にはあるのですが、どうなるかは保証の限りでは有りません(^_^;

    返信削除
  3. 自分でも少しチャレンジしたのですが、やはり基本などが解らないただのユーザーでは無理でした。
    ・SDVに基盤地図標高と近所を歩き回って取ってきたGPSログを読み込んで
    ・SDVからTIFFとTFWを書き出して、
    ・Libgeotiffというユーティリティを使ってGeotiffに合成しようとしたのですが、
    http://freegeographytools.com/2007/a-new-gui-for-the-libgeotiff-worldfile-utilities
    (赤字の「at this link.」より)
    ・SDVから書き出されたTFWデータには座標値はあるのですが投影情報や座標系情報が無いみたいでGeoTiffとして使えない状態でした。
    ・何か、tfw(またはgtf)下記のようなデータ形式にならないといけないみたいですが、。

    GeoTiff_Information:
    Version: 1
    Key_Revision: 0.2
    Tagged_Information:
    ModelTiepointTag (2,3):
    0 0 0
    419610.38 3901689.53 0
    ModelPixelScaleTag (1,3):
    2.4384 2.4384 0
    End_Of_Tags.
    Keyed_Information:
    GTModelTypeGeoKey (Short,1): ModelTypeProjected
    GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
    ProjectedCSTypeGeoKey (Short,1): PCS_NAD27_UTM_zone_12N
    PCSCitationGeoKey (Ascii,25): “UTM Zone 12 N with NAD27″
    End_Of_Keys.
    End_Of_Geotiff.

    PCS = 26712 (name unknown)
    Projection = 16012 ()
    Projection Method: CT_TransverseMercator
    ProjNatOriginLatGeoKey: 0.000000 ( 0d 0′ 0.00″N)
    ProjNatOriginLongGeoKey: -111.000000 (111d 0′ 0.00″W)
    ProjScaleAtNatOriginGeoKey: 0.999600
    ProjFalseEastingGeoKey: 500000.000000 m
    ProjFalseNorthingGeoKey: 0.000000 m
    GCS: 4267/NAD27
    Datum: 6267/North American Datum 1927
    Ellipsoid: 7008/Clarke 1866 (6378206.40,6356583.80)
    Prime Meridian: 8901/Greenwich (0.000000/ 0d 0′ 0.00″E)

    Corner Coordinates:
    Upper Left ( 419610.380, 3901689.528)
    Lower Left ( 419610.380, 3885203.505)
    Upper Right ( 432770.425, 3901689.528)
    Lower Right ( 432770.425, 3885203.505)
    Center ( 426190.403, 3893446.517)

    返信削除
    返信
    1. gtfの知識はないですが、tfwには測地系を指定できないし、UTMの場合にゾーンも指定できません。そのためにgtfファイルというものを用意したのでしょう。ワールドファイルには左上コーナーアドレスと、縦横のデータピッチを指定する機能があるだけです。測地系はWGS84だと仮定して処理するしかないでしょうね。

      削除
  4. 試しに、サンプルダウンロードできる北海道の衛星データのGeoTiffからLibgeotiffを使ってtfwとgtfを抽出してみました。
    http://www.alos-restec.jp/staticpages/index.php/service-sampledata-10
    やはり、先に記載のような形式のデータが記述されているみたいです。


    「tfw」
    9.7171144950
    -2.3617125126
    -2.3617125190
    -9.7171144877
    521720.9613009879
    4823757.5360865006


    「gtf」
    Geotiff_Information:
    Version: 1
    Key_Revision: 1.0
    Tagged_Information:
    ModelTransformationTag (4,4):
    9.71711449 -2.36171252 0 521717.284
    -2.36171251 -9.71711449 0 4823763.58
    0 0 1 0
    0 0 0 1
    End_Of_Tags.
    Keyed_Information:
    GTModelTypeGeoKey (Short,1): ModelTypeProjected
    GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
    GTCitationGeoKey (Ascii,50): "Corrected AVNIR-2 Image Corresponding To RPC Data"
    ProjectedCSTypeGeoKey (Short,1): PCS_WGS84_UTM_zone_54N
    PCSCitationGeoKey (Ascii,12): "UTM / WGS84"
    End_Of_Keys.
    End_Of_Geotiff.

    PCS = 32654 (name unknown)
    Projection = 16054 ()
    Projection Method: CT_TransverseMercator
    ProjNatOriginLatGeoKey: 0.000000 ( 0d 0' 0.00"N)
    ProjNatOriginLongGeoKey: 141.000000 (141d 0' 0.00"E)
    ProjScaleAtNatOriginGeoKey: 0.999600
    ProjFalseEastingGeoKey: 500000.000000 m
    ProjFalseNorthingGeoKey: 0.000000 m
    GCS: 4326/WGS 84
    Datum: 6326/World Geodetic System 1984
    Ellipsoid: 7030/WGS 84 (6378137.00,6356752.31)
    Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)

    Corner Coordinates:
    Upper Left ( 521717.284, 4823763.576)
    Lower Left ( 502823.583, 4746026.660)
    Upper Right ( 592040.041, 4806671.862)
    Lower Right ( 573146.341, 4728934.946)
    Center ( 547431.812, 4776349.261)

    返信削除