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

一部のシェープファイルには.prjファイルが関連付けられており、.prjファイルにはシェープファイルの投影情報がWKT形式で含まれています。WKTをproj4文字列に変換する必要がある場合もあれば、元に戻す必要がある場合もあります。 これを行うための既製のライブラリはありますか？

14 shapefile  coordinate-system  wkt  proj 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 2年前に閉店。 背景これは、異なる座標系で他のデータレイヤーと組み合わせて再表示するために、裸のラスターマップをジオリファレンスすることに関する2番目の質問です。前の質問は、メタデータを含まない任意のマップ画像をQGISプロジェクトに変換することです 問題私の目標は、このマップをジオリファレンスすることです。 これはPlate-Carréeではないようです。そこで、QGISで、いくつかの妥当なコントロールポイントを作成しました。完全を期すために、下部に追加しました[参照：1]。プロジェクトファイルEPSG：4326と同じターゲットSRSをQGIS Georeferencerに提供します。Helmertと多項式変換では非常に悪い結果が得られますが、薄いプレートスプラインでは妥当な画像が得られます（結果として得られる地理推定値がコントロールポイントを通過します）。ただし、この結果でさえ、たとえば高緯度では貧弱です（日本の北のロシア沿岸を参照）。これは、Natural Earthの背景を使用したQGIS画面のスクリーンショットです。 代替パスでの結果とコントロールポイントを参照してください。私はMapWarperではるかに簡単に使用できるツールと同様の行使しようとしたhttp://mapwarper.net/maps/758#Preview_Map_tab私はおそらく事実によるもの（貧しい結果を取得する場所を追加したコントロールポイントが少ないこと）。 一言で言えば質問 良いジオリファレンスを取得するのに欠けているトリックはありますか？ この投影はすぐに認識可能ですか？ で不明な古い図面上の座標系、gdaltransform実際に元のマップを生成するために使用される投影パラメータを暴くことを目標に、いくつかのターゲットSRSに、いくつかの座標点を変換することが示唆されます。私はこのような何かを試してみました：ポイントの私のQGISリストを保存した後、私は経由して、スペース区切りロング/ラッツのリストを取得するには、いくつかの文字列処理をしたcat eurasian-steppe-gcp.points | tail -n+2 | cut -d, -f1-2 | sed 's/,/ /'> tmp.txtとgdaltransformに結果のファイルを入力：gdaltransform -s_srs EPSG:3785 -t_srs EPSG:4326 < tmp.txtスイッチングs_srsとt_srsフラグ（プロジェクトはEPSG：4326を使用します）。幸運なことを望んで、暗闇で撮影していることを知っているので、出力を理解できなくても驚かなかった。誰かがこの方法を使用してソースマップの投影と投影パラメーターの最適な推定値を見つける方法を拡張できますか？この背後にある私の考えは、無数のコントロールポイントを配置して適切なジオリファレンスを作成するのではなく、すべての一般的な座標系をループするだけで、より少ないコントロールポイントでほぼ完璧なジオリファレンスを取得する方が簡単だと思いますか？テスト対象の各CRSについて、他のすべてのポイントに対する各ポイントの相互検証が含まれますか？ プロセスを自動化できるように、このアルゴリズムまたはジオリファレンスのいずれかを理解したいと思います---この問題に常に遭遇し、コンテンツ作成者がマップを統合されることのない一回限りの作成として扱うのをやめるまで他のコンテンツで、私は停止する予定はありません。 参照資料 [参照：1] QGIS GCPファイル： mapX,mapY,pixelX,pixelY,enable 142.632649100000009,54.453595900000003,505.941176470588232,-95.220588235293974,1 154.934252200000003,59.559921699999997,536.411764705882206,-52.779411764705742,1 80.080158100000006,9.657192300000000,291.558823529411711,-322.661764705882206,1 10.448442600000000,57.819128900000003,21.676470588235190,-103.926470588235134,1 34.007173000000002,27.761438299999998,101.117647058823422,-244.852941176470466,1 50.950890399999999,11.862196600000001,171.852941176470495,-313.955882352941046,1 29.713217199999999,60.024133200000001,90.779411764705799,-92.499999999999829,1 60.000000000000000,0.000000000000000,208.308823529411683,-362.382352941176350,1 69.867506500000005,66.639146199999999,224.088235294117567,-33.191176470588061,1 27.276107100000001,71.049154799999997,89.147058823529306,-21.764705882352814,1 140.000000000000000,0.000000000000000,536.955882352941217,-362.926470588235190,1 20.000000000000000,0.000000000000000,43.441176470588132,-362.926470588235190,1 …

13 qgis  raster  coordinate-system  georeferencing 

代理店からプロジェクトのシェープファイルを受け取り、常にUTMを使用していると言われました。ただし、値の桁数が多すぎてUTMに投影できません。 min max x 32507319 32624343 y 5925974 6074881 ファイルをQGISに読み取ろうとしましたが、どの投影を試しても、世界地図を使用していても、形状は常にマップの境界の外側に表示されます。それらが表示される領域は、ドイツの最北端の州であり、UTMで次の最小値と最大値を持つシュレスヴィヒホルシュタインです。 min max x 426167.5 650147.9 y 5913377.7 6101334.8 簡単にわかるように、上位データに何か問題があります。これを修正するために使用できる投影法は誰にもわかりますか？

13 coordinate-system  germany 

NAD83 SC State Planeにフィーチャクラス（サウスカロライナ州の郡なので、かなり広い地理的領域）があります。2番目の投影（NAD83 UTM 17）に変換してから、元に戻す必要があります。これを実現するために、Esriのプロジェクトツールを使用します。 この二重変換は、ポリゴン座標の位置を、センチメートル、メートル、キロメートルのどれだけシフトさせることができますか？

13 arcgis-desktop  coordinate-system  united-states 

私はhttp://www.transpacificproject.comから太平洋を示すこの地図を見つけました： 地球の同じ地域について、同様の地図を作成したいと思います。また、大陸間の距離をよく示していると思うのと同じ投影法を使用したいと思います。 このプロジェクションの名前を検索しましたが、答えがまだわからないと思います。最初はロビソンだと思っていましたが、ウィキペディアによると： ロビンソン図法は、全世界を一度に表示する世界地図の地図図法です。 そうではないようです。（ここではヨーロッパとアフリカを見ることができません。） だから、私の質問は2つあります：この投影の名前は何ですか？以前に言ったことを説明するためにそれを使用する必要がありますか、それとも何か他のものを使用する必要がありますか？

13 coordinate-system  cartography 

これはグーグルアウトするのに本当に不便なトピックでした-私は明確化を求めています、私は物事を正しく理解しましたか？ メルカトル図法（3857）で短距離を測定することについて話しています。マップ単位で距離を計算できます（ピタゴラスの定理を使用）。赤道では、地球上の距離に等しい（スケールファクター= 1）。極に向かって移動すると、マップ単位と地球の距離は等しくなくなります-スケール係数は増加し始めます。 質問は：与えられた緯度のスケール係数を計算する方法ですか？ 私は正しく理解しましたfactor = 1 / cos(latitude)か、式は？

13 coordinate-system  distance  mercator 

このマップをロビンソン図法で投影したい： library(ggmap) world <- map_data("world") ggplot() + geom_path(data = world, aes(long, lat, group = group)) そして、投影を「ロビンソン」に変更したいと思います（前の質問への回答からのアドバイスに従ってください：ウィキペディアからの地球規模の気候地域マップはどの投影を使用していますか？ このプロジェクションのデフォルトの実装を見つけるのに苦労しました。proj4ライブラリを使用するために次のことを考えました。 library(proj4) robinson <- project(cbind(world$long, world$lat), proj = "+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs") 私は以下を含む多くのアプローチを試みました： # using ggmap::get.map() get_map("world", projection = mapprojection(robinson)) # using ggplot2::coord_map coord_map(projection = robinson) # and sp::coordinates: library(sp) …

13 coordinate-system  r  proj  ggmap 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 4年前に閉鎖されました。 ArcGISのさまざまなシナリオで長さと面積を計算する方法を理解しようとしています。フィーチャクラスのフィールドで答えが見つからない理由はわかりませんが、何かを理解しておらず、いくつかの歴史があることを知っていない限り、正確な答えを見つけることができません。疑問符の記入を手伝ってもらえますか？または、なぜこれについてすべて間違っているのか教えてください。） GCS =地理座標系PCS =投影座標系 すべてのリンクは10.1ヘルプドキュメントにあります- フィーチャクラスのshape_lengthおよびshape_areaフィールド a。GCS-？ b。PCS-単純な平面の使用 c。シェープファイルを除き、常に自動更新されますか？はい ArcMap測定ツール a。GCS-デフォルトの測地線、代替案はLoxodromeとGreat Ellipticですが、平面ではありません。面積計算は利用できません！ b。PCS-デフォルトの平面、代替の測地線、ロクソドローム、グレートエリプティック http://resources.arcgis.com/en/help/main/10.1/index.html#//00s500000022000000 属性テーブル計算機 a。GCS-利用不可 b。PCS-プラナー http://resources.arcgis.com/en/help/main/10.1/index.html#//005s00000027000000 フィールド計算ツール（データ管理ツールボックス） a。GCS-測地線線形、利用可能な面積は疑わしい b。PCS-プラナー http://resources.arcgis.com/en/help/main/10.1/index.html#//00170000004m000000 バッファリングツール（およびその他のツール） a。GCS-測地線 b。PCS-平面またはGCS出力を指定 http://resources.arcgis.com/en/help/main/10.1/index.html#//000800000019000000 Javascript APIクライアントサイド a。GCS-測地線の面積と長さの関数 b。PCS-Webメルカトルから地理への変換（またはジオメトリサービスの使用） http://help.arcgis.com/en/webapi/javascript/arcgis/help/jsapi/namespace_geometry.htm Flex APIクライアントサイド a。GCS-測地線の面積と長さの関数、「長さ[または面積]は、カスタムの円柱等面積投影法を使用して計算されます」。これはjavascript apiには記載されていません!! b。PCS-Webメルカトルから地理http://resources.arcgis.com/en/help/flex-api/apiref/com/esri/ags/utils/GeometryUtil.htmlに変換できます ArcGIS Server REST API-ジオメトリサービス a。GCS-測地線 b。PCS-プラナー http://help.arcgis.com/en/webapi/javascript/arcgis/help/jsapi/geometryservice.htm 別の質問、測地線測定とは正確には何ですか？スフェロイドの3Dトリガー式を意味すると思いました（ハバーシン？）。また、面積の計算に使用するのは遅すぎるので、等しい面積の投影が使用されるのはなぜですか？ 長さと面積を決定するときの別の質問は、同じデータ、回転楕円体を使用した測地線計算よりも等面積投影法の方が正確ですか？そして、なぜ簡単に？

13 arcgis-10.1  arcgis-desktop  coordinate-system 

ArcGIS 10.1を使用していて、2つの既存のラスターに基づいて新しいラスターを作成したい。RasterToNumPyArrayは私が適応したい良い例があります。 import arcpy import numpy myArray = arcpy.RasterToNumPyArray('C:/data/inRaster') myArraySum = myArray.sum(1) myArraySum.shape = (myArray.shape[0],1) myArrayPerc = (myArray * 1.0)/ myArraySum newRaster = arcpy.NumPyArrayToRaster(myArrayPerc) newRaster.save("C:/output/fgdb.gdb/PercentRaster") 問題は、空間参照とセルサイズを削除することです。arcpy.envを実行する必要があると考えましたが、入力ラスターに基づいてそれらを設定するにはどうすればよいですか？私はそれを把握することはできません。 ルークの答えを取り上げると、これは私の暫定的な解決策です。 Lukeのソリューションはどちらも、空間参照、範囲、セルサイズを正しく設定しました。ただし、最初の方法では配列内のデータが正しく伝達されず、出力ラスタにはどこでもnodataが入力されます。彼の2番目の方法はほとんど動作しますが、nodataの大きな領域がある場合は、ブロック状のゼロと255で埋められます。これは、nodataセルをどのように処理したかに関係している可能性があります。私が話していることの画像を含めました。 #Setting the raster properties directly import arcpy import numpy inRaster0='C:/workspace/test0.tif' inRaster1='C:/workspace/test1.tif' outRaster='C:/workspace/test2.tif' dsc=arcpy.Describe(inRaster0) sr=dsc.SpatialReference ext=dsc.Extent ll=arcpy.Point(ext.XMin,ext.YMin) # sorry that i modify calculation …

13 arcpy  raster  arcgis-10.1  coordinate-system  numpy 

PostGISでの距離の計算について簡単な質問があります。 2つのジオメトリ間の距離を取得したいと思います。私はこのsidを使用しています：4269メートルで私が今していることはこれです：ST_Distance((a.geom,b.geom)) FROM ...しかし、私は度で結果を得ています。地理を扱うべきだと思いますが、どのようにジオムを地理にキャストできますか？（a.geom :: geography）で試しましたが、エラーが発生しています。 結果をメートル単位で取得するにはどうすればよいですか？ ありがとうございました

13 postgis  coordinate-system  postgresql  distance 

紙の地図を読んでいる場合、そのような情報がドキュメント上で明示的にラベル付けされていない場合、地図を描くためにどの投影法が使用されたかをどのように判断しますか？ オンラインでは、使用中のさまざまな予測を記述および分類するリソース（WikipediaやGIS Wikiなど）を簡単に見つけることができます。しかし、専門家ではない地図ユーザーとして、特定の地図を描くためにどの投影法が使用されたかを判断するために、建設の詳細を逆に使用することは難しいと感じています。地図ユーザーが使用中の投影法を決定するのに役立つフローチャートやその他の分類プロセスがあるかもしれないと思ったのでしょうか？ たとえば、私の地元の図書館とアーカイブには、1878年のDominion Lands Office（カナダ）からの地図のリトグラフのコピーがあります。この場合のポール）。これにより地図が円錐図法になると思いますが、よくわかりません。そして、それが正しければ、投影をより正確に分類できるかどうかを知ることに興味があります（たとえば、Lambert共形）。 誰かが私の質問に答えるのに役立つ場合に備えて、この1878マップの部分スキャンを作成しました（使用可能な機器で全体をスキャンするには大きすぎます）。

13 coordinate-system  paper 

無料のライブラリを使用して空間データを変換するにはどうすればよいですか？ たとえば、C＃Webアプリケーションのコード内でShapefileの投影を変更したい。それ、どうやったら出来るの？

13 shapefile  coordinate-system 

私の仕事では、昔から古いシステムから式がコピーされたステレオグラフィック投影法を実装すると主張する内部ソフトウェアのメンテナンスを行ってきました。 私が見たものから、それはEPSG（コード9809）からの「斜めと赤道のステレオグラフィック」を実装するようです。しかし、私が特定できなかった一歩があります。コードで使用される「R」値は、次のパラメーターを使用して計算されます。 (phi_n, lambda_n) - null distortion coordinates (PHI_nc, LAMBDA_nc) - conformal null distortion coordinates (PHI_o, LAMBDA_o) - conformal projection origin coordinates e - eccenttricity r_eq - earth radius at equator これを前提として、Rを取得する式は次のとおりです。 a = 1 + cos(PHI_nc) * cos(PHI_o) * cos(LAMBDA_nc - LAMBA_o) b = sin(PHI_nc) * sin(PHI_o) c = …

13 coordinate-system  datum 

私のQGISプロジェクトには2つのレイヤーがあり、1つのレイヤーはCRSが+ proj = merc + lon_0 = 0 + lat_ts = 0 + x_0 = 0 + y_0 = 0 + a = 6378137 +である静的なGoogleマップ画像を取得するGoogleレイヤープラグインによって取得されますb = 6378137 + units = m + no_defsで、別のベクターレイヤーはCRSが+ proj = longlat + ellps = WGS84 + datum = WGS84 + no_defsであるpostgisレイヤーです。 Googleレイヤーからいくつかのポリゴンフィーチャをデジタル化し、ベクターレイヤーに追加しました。ただし、QGISは自動CRS変換を行わないようです。投影を変更するにはどうすればよいですか？ ありがとう

13 qgis  coordinate-system 

閉まっている。この質問はトピック外です。現在、回答を受け付けていません。 この質問を改善したいですか？ 質問を更新して、地理情報システムスタック交換のトピックになるようにします。 2年前に閉店。 ArcPyからアクセスされているフィーチャクラスにポイントフィーチャがあります。ポイントは投影されますが、このポイントの未投影の緯度と経度を取得する効率的な手段を見つける必要があります。 再投影（投影解除）、新しいフィーチャクラスで検索カーソルを取得、フィーチャを検索、およびフィーチャの形状から緯度経度を取得する以外の方法はありますか？

13 arcpy  coordinate-system