タグ付けされた質問 「proj」

0〜360の経度の座標系を持つジオティフラスターイメージがあります。画像の水平方向の中心は180経度です。下の画像をご覧ください： -180 180経度範囲のEPSG：4326 SRSに変換したい。そして、画像の中心がグリニッジ子午線（0）になるようにします。このsrsは非常に広く使用されていると思います。結果は次のようになります。 そこで、gdalwarpコマンドを使用して再投影します。 gdalwarp -s_srs '+proj=latlong +datum=WGS84 +pm=180dW' -t_srs EPSG:4326 test_col.tif test_4326.tif しかし、より大きなサイズ（より多くのピクセル）とEPSG：4326メタデータを持つTIFFのみを取得します。画像自体は、最初のものと同じに見えます。しかし、私はそれが半球を交換することを期待しています。 質問は-経度0の中心で厳密に-180 180 EPSG：4326になるように画像をgdalwarpする方法です。 これは私の初期ファイルのgdalinfoです： Origin = (-0.102272598067084,89.946211604095552) Pixel Size = (0.204545196134167,-0.204423208191126) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( -0.1022726, 89.9462116) ( 0d 6' 8.18"W, 89d56'46.36"N) Lower Left ( -0.1022726, -89.9462116) ( …

31 coordinate-system  gdal  geotiff-tiff  proj  gdalwarp 

QvGISにNTv2変換を統合する正しい方法は何だろうと思っています。 NTv2変換をQGISに正常に統合しましたが、意図したとおりに実行されたかどうかはわかりません。 NTGISをQGISに統合するために、グリッドの名前をsrs.dbのテーブルtbl_datum_transformのp1列に保存する必要があるといういくつかの情報を見つけました。 https://github.com/qgis/QGIS/pull/1502#issuecomment-50580199 https://dragons8mycat.com/2014/03/03/gis-tips-getting-ostn02-to-work-with-qgis-2-2/ https://web.archive.org/web/20140224182157/http://sourcepole.ch/ntv2-transformations-with-qgis http://isticktoit.net/?p=371 （http://www.gkg-kassel.de/pdf/Praezise_Transformation_im_QGIS_mit_Batchmodus.pdf（ドイツ語）のスクリーンショット） それが正しいかどうか疑問に思うのは、列に数値フィールドタイプがあるように見えるという事実です。 7パラメーター変換の最初のパラメーターはQGISのtbl_datum_transformテーブルのp1列に格納されていませんか？ sqlite FAQによると、sqliteの数値列に文字列を保存することは可能ですが（https://sqlite.org/faq.html#q3）、文字列を数値列に保存するハックや回避策のように思えます。これは本当にやるべきことですか？ tl; dr：QvGISにNTv2変換を統合する正しい方法は何ですか？

15 qgis  coordinate-system  proj  epsg 

ポリゴンのシェープファイルと、（Lat、Lng）ペアとしてのポイントのリストを含む別のCSVファイルがあります。 CSVファイルの各（lat、lng）ペアをチェックして、どのポリゴンに含まれるかを確認します。 シェープファイルが投影され、projファイルは次のようになります。 PROJCS["Transverse_Mercator",GEOGCS["GCS_OSGB 1936", DATUM["D_OSGB_1936",SPHEROID["Airy_1830",6377563.396,299.3249646]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",49],PARAMETER["central_meridian",-2],PARAMETER["scale_factor",0.9996012717],PARAMETER["false_easting",400000],PARAMETER["false_northing",-100000],UNIT["Meter",1]] 私の計画は次のとおりです。 readShapePolyR MapToolsパッケージの関数を使用してシェープファイルを読み取ります。 CSVファイルからポイント座標をデータフレームに読み取り、SpatialPointsDataFrameに変換します over関数を使用して、内部にあるポリゴンを判別します。 そのためにはproj4string、ステップ1でシェープファイルをロードするときに指定し、ステップ3で関数をspTransform適用する前にover、ポイントとポリゴンが必要であるため、関数を使用してCSVファイルから座標を同じ投影システムに変換する必要があります同じ投影システムの下にある。 上に示したprojファイルの内容の正しい値は何でしょうか？

14 coordinate-system  shapefile  r  proj 

proj 4.8.0がインストールされたQGISがあります。次のように定義されていると思われるWinkel Tripelプロジェクションを追加したい +proj=wintri しかし、その定義文字列は機能しません。これは私のprojバージョンのバグですか、それとも正しい定義文字列は何ですか？

14 qgis  coordinate-system  proj 

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

14 shapefile  coordinate-system  wkt  proj 

ライブラリの問題を解決しましたが、PROJ.4ライブラリがロードされていないライブラリの問題と漠然と関連したものだけが残ります。私は尋ねたフォローアップの質問をこれらのために。 GDALパッケージとPROJ.4ライブラリが破損しているため、ファイルを再投影しようとするとエラーが発生します（以下の古い質問を参照）。GDALの別のインスタンスを誤ってインストールしたときに以前にこの問題が発生したかどうかはわかりません。ただし、アンインストールしてQGISを再インストールしましたが、問題は残ります。 再インストールのたびに、QGISは以前と同じプラグインを使用するため、エラーの原因となっている部分はインストール/アンインストールプロセスによって変更されません。 私が試したもの： sudo apt-get remove build-essential python-all-devそしてsudo apt-get remove gdal-bin手動でインストールGDALをパージします Ubuntu Software CenterでのQGISの（アン）インストール sudo apt-get --purge autoremove dpkg -l | grep '^rc' | awk '{print $2}' | sudo xargs dpkg --purgeおよびdpkg --list | grep "^rc" | cut -d " " -f 3 | xargs sudo dpkg --purgeコンフィギュレーションファイルを削除するには まだ運がありません。関連しているように見えるものを見つけましたが、それはWindows用であり、Ubuntu …

13 installation  ogr2ogr  ubuntu  proj 

このマップをロビンソン図法で投影したい： 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 

WGS84緯度/経度ポイントを変換しようとしています 11d42'32.10629" E 5d12'56.75371" S 3つの異なるソフトウェアパッケージ（Proj4、GeoTrans、Leica GeoOffice）で提供されていますが、結果には非常に大きな違いがあります-XとYで数百メートル前後！これらのバリエーションは、3ではなく7つのパラメーターでのみ表示されることに注意してください。 Proj4 cs2cs + proj = longlat + ellps = WGS84 + datum = WGS84 + no_defs + to + proj = utm + zone = 32 + ellps = clrk80 + towgs84 = 178.3,316.7,131.5、-5.278、-6.077、-10.9792、-19.166 +南+ units = m + no_defs GeoTrans Delta X …

12 coordinate-system  proj  datum 

私はgdalwarpを使用して、静止からランバート正角に再マッピングしようとしています。入力データはnetcdfにあり、地理座標（度）にあります。再マップしたデータをnetcdfに出力したいと思います。入力netcdfデータに対応するvrtファイルを作成しました。Gdalwarpはnetcdfファイルを出力しますが、出力データはすべてゼロであり、次のエラーを受け取ります。 Creating output file that is 5120P x 5120L. Processing input file netcdf.vrt. ERROR 1: Too many points (441 out of 441) failed to transform, unable to compute output bounds. Warning 1: Unable to compute source region for output window 0,0,5120,5120, skipping. 0...10...20...30...40...50...60...70...80...90...100 - done. 私は次のコマンドを試みました： /usr/bin/gdalwarp -s_srs "+proj=geos +h=35785831 …

11 gdal  coordinate-system  proj  gdalwarp  netcdf 

ドキュメントに記載されているように、proj4 foramtを使用して、QGISでオルテリウス楕円投影とニコロージ球状投影を定義しようとしていますが、proj4ライブラリでサポートされている場合でも、QGISはこれらの投影に投影できません...古いものをさらに分析するには、これらの投影が必要ですマップ。 誰が何が問題になるのか知っていますか？

10 qgis  proj 

私は、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) …

10 python  coordinate-system  wgs84  proj 

ローカル座標系を基準とした座標を持つ点群があります。GPS値のある地上管理ポイントもあります。PROJ.4または他のライブラリを使用して、これらのローカル座標をグローバル座標系に変換できますか？ 上記の問題に関するPythonのコードは、非常に役立ちます。

9 python  coordinate-system  gps  proj  ground-control 

QGISコース内で、proj.4パラメーターについて少しお話ししたいと思います。ただし、「+ no_defs」パラメーターについてはわかりません。proj.4のWebサイト（http://proj4.org/parameters.html#parameter-list）には次のように書かれています。 + no_defs /usr/share/proj/proj_def.datデフォルトファイルを使用しない これが何を意味するのか誰かが私に説明できますか？（/usr/share/proj/proj_def.datデフォルトファイルを使用するとはどういう意味ですか？）

9 qgis  proj 

世界地図の新しい正確な投影►Equal Earth。 そのproj4は知っていますが、QGISに追加できません（以下を参照）。QGISの更新が必要ですか？

8 qgis  coordinate-system  proj  qgis-3.2  qgis-3.4 

私は、オープンソースソフトウェアがオーストラリアの新しいデータを適切に処理できるように支援しています。GDA2020プロジェクトの詳細については、ICSMのWebサイトを参照してください。 現在、QGISにはGDA2020の定義が含まれていますが、GDALを介して理解しています。 GDA2020座標参照系の例は次のとおりです。 +proj=utm +zone=55 +south +ellps=GRS80 +units=m +no_defs GDA94 CRSを見ると、次のように定義されています。 +proj=utm +zone=55 +south +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs ご覧のとおり、これらは非常によく似ています。 現在、2つのCRSはまったく同じように定義されていますが、GDA94からGDA2020では、北東に約1.5 mの座標のシフトがあります。（NTv2形式のグリッドシフトファイルがあり、間もなく準備が整い、正確な変換が可能になりますが、それはこの質問の目的ではありません。） あなたはGDA94とGDA2020の間で変換する場合でも、今 QGISを使用して、座標系に変化はありません。それは本質的にそれを異なって分類するだけです。 GDA94とGDA2020の間のデフォルトの変換（ただし不完全）である、Proj.4または他のオープンソースツールに実装された単純な7パラメーター変換が必要ですか？ それとも、ツールが常に変化しないということですか？ これはどのように処理する必要がありますか？ （そして、グリッドを使用した変換が理想的であり、これはこのQGISプラグインを含むいくつかの方法で処理されることを再度述べておきます。）

8 qgis  coordinate-system  gdal  proj  gda2020