タグ付けされた質問 「coordinate-system」

空間内のポイントの位置を2次元または3次元で定義するために使用される、ポイント、ライン、またはサーフェスのセット、およびルールのセットで構成される参照フレームワーク。このタグは、地図投影法と座標系変換も参照します。

3
QRコードの緯度/経度または経度/緯度のGEOタグですか
地理コンポーネントをQRコードに追加できることがわかりました。 また、これを使用してGoogle Chart APIで作成できることもわかりました。 https://chart.googleapis.com/chart?cht=qr&chs=400x400&chl=geo:30.3,-97.6 これにより、次のようなQRコードが作成されます。 ここまでは順調ですね。そして、QRコードリーダーアプリを使用してiPhoneでそれを読み返すと、正しい場所(オースティン、テキサス州)にあるGoogleマップに送られます。 問題は、これを作成するために使用される文字列がlat / longを使用することです。同じデータをKMLファイルに書き込んだ場合、long / latが予想されます。 問題は、QRコードにどのスペックを使用するかです。 QRスペックは実際に緯度/経度ですか? Google Chart APIは混乱していますか? QRリーダーアプリに問題がありますか? 私のお金は、QRリーダーがめちゃくちゃになっていますが、よくわかりません。誰もが決定的な答えを持っていますか?

3
pyprojとProj.4を使用してX、Y座標を緯度/経度に変換すると、誤った座標が返される
私は、x座標とy座標を含む複数のXMLファイルを読み取り、それらをすべて1つのcsvファイルに結合するpythonスクリプトを作成しています。緯度と経度はcsvの必須フィールドですが、オハイオ北部州平面usFtのX、Y座標をWGS84に変換するのが困難です。 >>> p = Proj(r'+proj=lcc +lat_1=41.7 +lat_2=40.43333333333333 +lat_0=39.66666666666666 +lon_0=-82.5 +x_0=600000 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=us-ft +no_defs') #Nad83 State Plane Ohio North US Feet Proj object using parameters >>> p(739400.91,2339327.3,inverse=True) (-80.138057868777224, 60.278230707978487) >>> p1 = Proj(init="epsg:3734") #Nad83 State Plane Ohio North US Feet Proj object using EPSG code >>> p1(739400.91,2339327.3,inverse=True) (-80.138057868777224, 60.278230707978487) …

1
非投影データと予測データを使用して行われた計算のエラー評価
この質問は、「投影データと非投影データからの流れ方向の計算と流域の描写」という題名の質問に基づいています。 投影対非投影DEMデータからの流れ方向の計算と流域の描写 ただし、これは完全に別の質問です。前述の質問では、球形/投影されていない地理座標系のデータでユークリッド距離を想定するアルゴリズム(ArcGISフロー方向など)の使用に問題があることが判明したためです。 地図の投影は、オレンジの皮を取り、それを机の上で平らにしようとするようなものであることがわかっています。地図の投影によって本質的にいくつかのエラーが発生します。ただし、特にデカルト/投影された平面サーフェスを想定する計算を実行している場合は、投影の利点により発生するエラーが相殺されるようです。この場合、私が興味を持っているアルゴリズムは、データが投影されていることを前提とするArcGIS Flow Directionアルゴリズムです(これは、私の研究に基づくほとんどのアプリケーションで採用されている前提です)。これは、距離の計算にユークリッドアプローチを使用しているためです。 私の質問は、非投影DEMデータ(地理座標系のDEMデータ)と投影データ(適切な投影などのDEMデータ)を使用して、特定の調査地域の流れの方向を計算することで発生する可能性のあるエラーをどのように定量化できるかです。 UTMまたは何かコンフォーマル)? 確かに、投影されていない同じDEMデータを使用して、流れ方向ラスターを導出できます。しかし、それではどうでしょうか。私たちの目標は、できる限り正確に地球の表面をモデル化することです(そして、元のDEMなどを作成するプロセスで発生する可能性のあるエラーには対処していません。これらは、私に関する限り一定です)。 ....投影されたDEMから導出された流れ方向データの方が優れていると想定し、2つのラスターの個々のセル値を比較して、(通常のD-8モデルのコンテキストで)方向の値が異なるセルを特定します)?これを行うには、投影されていないデータから派生したフロー方向ラスタを取得し、投影されたフロー方向ラスタで使用されているのと同じ投影法を適用する必要があると思います。 最も意味のあるものは何ですか。また、投影されていないDEMを精度のベンチマークとして何と比較する必要がありますか。 数学の方程式の細部に入ると、それを理解する人にとって、地上レベルでの証明が得られ、一部にとっては十分かもしれませんが、それだけでなく、数学を十分に理解しているが、危険なほど地理/ GISを知っているだけの場合は素晴らしいでしょう(理想的には、両方のレベルがハードコアな地理オタクと平均的なGISのやり手と共鳴するのが良いでしょう)。より高いレベルの人々にとって、証明は数学にあると言って、おそらく議論の余地があるかもしれません-私はもっと具体的なものを探しています(たとえば、政府のある種の非効率性にドルの数字を付けるようなものです)。 これをどのように定量化できるかについての考えやアイデアがあれば、大歓迎です。 トム


3
ArcGISは、非等間隔投影で2つのポイント間の距離をどのように計算しますか?
これは、私の前の質問に対するフォローアップの質問です。座標系の投影法について、よく書かれた紹介文を提案できますか? 私がCH1903マップ投影で作業していると仮定しましょう。これは、私が知っていることはすべて共形ですが、等距離ではありません。つまり、角度(形状)は保持されていますが、面積、距離、スケールは保持されていません。(少なくともこれらは正確に保存されていません)。ここまでは順調ですね。 2点間の距離を計算するときに、ArcGISがどのような計算を実行するのか疑問に思っています。ArcObjectsでは、次のようにIProximityOperatorインターフェイスを使用できます。 IPoint a = ..., b = ...; double distance = ((IProximityOperator)a).ReturnDistance(b); 質問:距離を正確に保持しないリファレンスシステムを使用している場合、(上記のように)2点間の距離をクエリするときにArcGISは何をしますか? それは単に距離を取得するためにいくつかのピタゴラスの数学(a 2 + b 2 = c 2)を実行しますか?つまり、返される距離は投影が許す限り正確になりますか? または、より正確な距離を取得するために、何らかの形の再投影など、より複雑なことを行いますか? (同じ質問ですが、より一般的には、そのジオメトリが投影された後、ArcGISはすべての計算をユークリッド空間で単純に実行しますか、それとも、使用されたマップ投影は、距離、角度、面積などの計算に影響しますか?)

6
オープンソースツールを使用してシェイプファイルにプロジェクションを追加する方法は?
プロジェクションなしでシェープファイルをいくつか与えられました。だから私はそれをグーグルマップに表示するために投影を追加する必要があります。必要な投影情報はhttp://spatialreference.org/ref/epsg/2066で見つかりました。 質問: オープンソースメソッドを使用して投影を追加するにはどうすればよいですか(ESRIのソフトウェアにアクセスできません)? ありがとう クリス


1
経度/緯度をマップ上のピクセルに変換します
ここから地図があります。任意のlon / latペアをマップ上のピクセルに簡単に変換できるようにしたい(逆変換も可能)。マップには.tfwファイルと投影情報が付属しています。 0.02222222222222 0.00000000000000 0.00000000000000 -0.02222222222222 -180.00000000000000 90.00000000000000 と投影情報: Projection: Plate Carree aka Geographic or "LatLong" Earth ellipsoid: Sphere, radius 6370997 m Datum: WGS84 Extent: 180 West to 180 East, 90 North to 90 South Size: 16,200 height samples wide x 8,100 high Resolution: 2.47 km/pixel 私は地図作成関連のことをまったく知らないので、私が理解している限り、まずWGS84(経度/緯度のペア)から地理的投影法に変換する必要があります(同じではありませんか?)。実際には同じように見えますが、上の投影情報の球の半径は6370997で、プレートカリー投影で見つけたspatialreference.comページとは異なります。とにかく、私はDotSpatial.Projectionsライブラリが次のコードでこれを行うことができることを発見しました: var kievCoordinates …

1
SRID 4326(WGS 84)から4269(NAD 83)SRIDへの再投影が必要ですか
MSSQL DBにSRID 4269として保存されている州の境界に対する米国国勢調査データがあります 。 そのデータセットに対する計算に使用される入力データは、SRID 4326として保存されます。 SpatialReference.orgから推測できる限り、SRID 4269はSRID 4326のサブセットにすぎず、座標の再投影は必要ありません。 入力データが常にSRID 4269の境界内にある場合、実際にSRID 4326からSRID 4269に再投影する必要がありますか?


3
このnetCDFラスターはどの座標系ですか?
netCDFラスターファイルを取得しましたが、ラスターが構築されている座標系の名前を取得するためのメタデータを取得できませんでした。ラスター自体には座標系が埋め込まれていません。WGS84だけだと思っていましたが、最初はそのように見えましたが、ArcMapをさらに調査したところ、それほど一般的ではないシステムであることがわかりました。表示方法は次のとおりです。 オレンジ色のラスターは、比較のためにここに挿入したWGS84の通常のラスターです。紫色のものは、座標系が不明なラスターです。これが何かの手掛かりはありますか? 一部の更新:ここでのnetCDFラスタがある:https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 私も画像提供者からいくつかのメタデータを取得することができました: netcdf sic_average_nclimate { dimensions: nlon = 361 ; nlat = 90 ; nseas = 4 ; variables: float SIC_Change(nlat, nlon) ; SIC_Change:Title = "Gridded Multi-Model Ensemble Mean Annual Mean Change in Ice Concentration 21C-20C" ; float SIC_Season_Change(nseas, nlat, nlon) ; SIC_Season_Change:Title = "Gridded Multi-Model Ensemble Mean …

3
PostGISのジオメトリSRIDの問題
PostGISで自分のテーブルを作成しました id_borne serial NOT NULL, num_borne character varying, shape_borne geometry を使用してテーブルのSRIDを正常に変更しました select UpdateGeometrySRID('public', 'borne', 'shape_borne', 26191) ; 以前のSRIDは0でした。 問題は、テーブルborneに新しい要素を追加しようとすると、次のメッセージエラーが発生することです。 geometry srid (0) does not match column srid (26191) それが私の問題です!SRIDは更新されていますが、新しいポイントを追加しようとすると、このメッセージが表示されます。

1
シェープファイルを再投影する方法は?
QGISプロジェクトを再投影したかったのですが、そのようなオプションはありません。回答を検索するときogr2ogr、各ベクターレイヤー(この例ではシェープファイル)を再投影することにより、このタスクの提案が表示されます。ただし、期待した結果が得られません。 C:\temp>ogrinfo -al -so misc.shp INFO: Open of `misc.shp' using driver `ESRI Shapefile' successful. Layer name: misc Geometry: Line String Feature Count: 10 Extent: (21.267388, 42.015857) - (21.270225, 42.017470) Layer SRS WKT: PROJCS["MGI_Balkans_zone_7_deprecated", GEOGCS["GCS_MGI", DATUM["Militar_Geographische_Institute", SPHEROID["Bessel_1841",6377397.155,299.1528128]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",21], PARAMETER["scale_factor",0.9999], PARAMETER["false_easting",7500000], PARAMETER["false_northing",0], UNIT["Meter",1]] id: Integer (10.0) C:\temp>ogr2ogr -t_srs …

2
Mapnikは「空白」のタイルのみをレンダリングします
私は魅力的なGISの世界に慣れていないため、この数日間、トピック全体を調査してきました。私の目標は、自分のサーバーでOSMデータを使用して地球規模のマップを提供することです。 私のスタックは: Postgres + Postgis -GISデータの保存、変換、処理用。 Mapnik-タイルをレンダリングするためのライブラリ Tilestache-タイルの提供とキャッシュ Nginx -Tilestacheへのリバースプロキシ リーフレット(Slippyマップクライアント) OSMデータをインポートしてスタイルを設定するために、私は PBFファイル(抽出/プラネット) osm2pgsql CartoCSS openstreetmap-carto Kosmtik-開発/デバッグ用のマップGUI 完全な惑星ファイルは大きいので、私は現在、サンパウロの抽出を使用しています。これは、テストするのに十分な大きさです。 私が抱えている問題に到達する前に、これまでに行った手順を説明します。 うまくいったこと PostgreSQLとPostGISが正しくインストールされている。前者には9.6、後者には2.3を使用しています。postgis_topologyとのhstore拡張機能も追加しました。これらは後で使用します。 osm2pgsql 0.92.0がインストールされています。サンパウロ抽出物をインポートするために、私は使用します osm2pgsql -G U <user> -d <db> -C 1000 -W --hstore --style openstreetmap-carto.style --tag-transform-script openstreetmap-carto.lua <pbf> hstore、styleおよびtag-transform-scriptとして引数は、OSMのスタイルの適切な利用のために必要とされるここで説明します。 openstreetmap-cartoのインストールで説明したように、カスタムインデックスを追加し、必要なシェープファイルとフォントをダウンロードしました(絵文字を除く)。 QGISを使用して、抽出データがデータベースに正しくロードされていることを確認しました。すべてのポイント、ポリゴン、ライン、道路を表示してクエリを実行できました。それだけです。 問題 次のステップはタイルのレンダリングです。そこで問題が発生します。mapnikとpython-mapnik 3.x、carto 0.18.2があります。 OSMのproject.mmlから、を使用して独自のproject.xmlを生成しましたcarto。結果に次の変更を加えました。 Postgresのユーザ名(パラメータの追加userにMAPNIKのXMLを)。残りの接続情報は〜/ .pgpassから取得されます。 相対ディレクトリdata/<shape>.shpを/full/path/data/<shape>.shp 相対ディレクトリsymbols/<symbol>.svgを/full/path/symbols/<symbol>.svg …

1
GDALを使用してPythonでスクリプトを再投影する
GDALで多くの問題を抱えています。ドキュメントが欠けていることもありますが、Pythonではほとんどサポートされていないようです。とにかく、これは次のように機能するバッチ再プロジェクトスクリプトです。 source_file = gdal.Open(filepath.encode("ascii")) source_wkt = source_file.GetProjectionRef() source_srs = osr.SpatialReference() source_srs.ImportFromWkt(source_wkt) reproj_file = gdal.AutoCreateWarpedVRT(source_file, source_wkt, dest_wkt) gdal.ReprojectImage(source_file, reproj_file, source_wkt, dest_wkt) reproj_attributes = reproj_file.GetGeoTransform() driver = gdal.GetDriverByName("GTiff") dest_file = driver.CreateCopy(outputpath.encode("ascii"), reproj_file, 0) dest_srsとdes_wktの部分はこのコードチャンクでは定義されていませんが、ループの外にあります(一度だけ定義する必要があるため)。それは一度動作するようです、私はそれから1つの見栄えの良いtifを得ることができ、それから私に「エラー6 WriteBlock()はサポートされていません」とpythonがクラッシュします。それらはすべて、同じ基本データを使用して同じ方法で作成された(時間だけが異なる)すべてのGeoTIFFです。 また、(GCSからPCSへの)再投影の性質により、AutoCreateWarpedVRTは多くの空白スペースを作成する傾向がありますが、値0を与えます。これは、実際のデータ値である可能性があるため問題です。代わりにnodata値を-99に設定する方法はありますか?

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.