地理情報システム

のレイヤーがありcoord. ref. : NAます。 私はどのように割り当てることができますcoord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0するにはcoord. ref. : NA？ class : RasterLayer dimensions : 127, 146, 18542 (nrow, ncol, ncell) resolution : 0.01, 0.01 (x, y) extent : -122.39, -120.93, 36.72, 37.99 (xmin, xmax, ymin, ymax) coord. ref. : NA data source : C:\2012Tif\2012ASC5min_ppt_spas1306_0001_20121222_0005_UTC.asc.tif …

16 raster  coordinate-system  r 

関数で構築した空間バッファー（30000バッファー）のリストがありますlapply。 buff.pts <- lapply(1:nrow(pts.prj), FUN=function(l){ buff <- gBuffer(pts.prj[l,], width=1000) ## 1km return(buff) })) > head(buff.pts) [[1]] class : SpatialPolygons features : 1 extent : 307941.8, 311941.8, 4994518, 4998518 (xmin, xmax, ymin, ymax) coord. ref. : +proj=tmerc +lat_0=0 +lon_0=-73.5 +k=0.9999 +x_0=304800 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs +towgs84=0,0,0 [[2]] class : SpatialPolygons …

16 polygon  r  sp  rgeos 

MacOSXにGDAL 2.0をインストールする最も簡単な方法を知っている人はいますか？ キングカオスのバイナリ配布とhomebrewの通常の容疑者はGDAL 1.11に固執しているようです 2.0に含まれていると思われるTopoJSON機能プロパティのサポートが必要です EDIT 2017/01/07： kingchaosバイナリディストリビューションがGDAL 2.1ディストリビューションを提供するようになりました

16 gdal  installation  ogr  mac  homebrew 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 先月閉鎖されました。 ナイジェリアのSentinel 2画像が必要です。CopernicusScihubにサインアップしましたが、テキストボックスにナイジェリアと入力しても結果は表示されません。私は何か間違ったことをしているに違いありません。彼らがナイジェリアのイメージを持っていない可能性はありますか？

16 data  remote-sensing  sentinel-2  nigeria 

別の列の値に基づいて、1つの列の値の合計を計算したいと思います。私の質問を説明するために、下の写真をアップロードしました。2番目の列は郵便番号で構成され、3番目の合計売上高です。郵便番号ごとの総売上はいくらですか？したがって、この場合： 列2の値の列3の合計は1023？出力はになるはずです4 (0+1+1+2+0)。その場合の合計は1024、同じロジックに従って11,5になります。 これは非常に大きなデータセットであるため、ZGISコードで合計売上を示す新しいフィールドをQGISフィールド計算機（またはその他のツール）で計算したいと思います。これまでのところ、私はそれを機能させることができないようで、ここの誰かが私を助けてくれることを望んでいました。 ご覧のとおり、いくつかの欠損値がありますが、これは計算上の問題ですか？おそらく言及しておくと良いのは、テーブルであることです。

16 qgis  field-calculator 

ArcGIS for Desktopで夜間にモデルが実行されている間、スリープ状態になりたいです。 ArcGIS for Desktopが終了したときにサウンドを再生して、私を目覚めさせたいと思います。 これは可能ですか？

16 arcgis-desktop  modelbuilder 

ここには、通りを表す2つのレイヤーがあり、アイデアは次のとおりです。最新のレイヤー（赤い線）を、他のレイヤー（明るい青）に保存されているデータの一部で更新しますが、これらのレイヤーは地図上の場所に応じて非常に近く、少し離れた地理的位置で類似しており、他の属性と同じ行がある場合に一致する属性は1つだけです（codlog x codlog5）。 たとえば、画像の中央には両方のレイヤーに正方形がありますが、青いレイヤーが緑のレイヤーと同じであることを認識して最近のレイヤーを更新するにはどうすればよいですか？ どうすればいいですか？それに役立つプラグインやツールはありますか？

16 qgis  layers  latitude-longitude 

私は座標系についてあまり知りません... 私のオフィスでは、考古学的なサイトからの空間データを処理するために使用しています。各サイトには独自のxyz座標系（GCS）があります。3つの単純な直交デカルト軸。ここ数年、特定の座標系を使用せずに、GISソフトウェア（ArcGIS）を介してこの空間データを管理しています（「未定義」のままにしておきます） 典型的なGCSのグリッドの歪みなしに、単純な直交直交軸を使用してそのようなデータセットを処理するように設計されたGCSが存在するかどうかを知りたいです。さらに、このシステムがオンラインマッピングアプリケーションでの使用に適しているかどうかを知りたいです。 ところで、2D（ArcMap）および3D（ArcScene）環境を管理し、長さの基本単位として「mm」を使用します。 そのようなものが存在しない場合、誰かがそれを作成する方法を知っている可能性があります。

16 arcgis-desktop  coordinate-system  epsg  srid 

米国の郡の境界を表す多数のポリゴンがあり、それらをできるだけ少ないポリゴンにマージする必要があります。 非商用ソフトウェア、またはリーズナブルな価格のソフトウェアでこれを行うことは可能ですか？ PostGISもJavaも利用できませんが、Cは利用できます。

16 ogr2ogr  aggregation 

さまざまな解像度と投影法を持つラスターデータレイヤーを扱うためのベストプラクティスに関する推奨事項を探しています。私が与えられたアドバイスは、分析を実行する前に常に最も低い解像度でレイヤーにリサンプリングすることですが、これは私にとって非常に正確な無駄のようです。 より高い解像度のグリッドに一致するようにリサンプリングするのが妥当なのはいつですか？また、より低い解像度へのリサンプリングと比較した場合の影響は何ですか？ これは状況に大きく依存する可能性が高いことを理解しています。私は主に一般的なガイドラインを探していますが、参照用の私の特定のシナリオは次のとおりです。 シナリオ：さまざまな環境層および社会経済層に基づいて、土地利用を予測する空間回帰モデルを構築しようとしています。私の土地利用マップはLandsatから派生したものであり、解像度は30mです。説明レイヤーの例には、SRTM DEM（3アーク秒、約90m）およびBioclim気候レイヤー（30アーク秒、約1km）が含まれます。

16 analysis  resolution  resampling 

現在、CANGRID気候データ（Surfer Grid ascii、「。grd」ファイルとして提供）をArcGISにインポートしています。グリッドのサイズは、95行125列です。メタデータは、緯度/経度の原点（左下隅）、セルサイズ（50km）に加えて、中央子午線（110度W）と原点の緯度（60度N）を備えた極立体図としてのノートの投影を提供します。 最初に.grdを.asciiおよび.fltとしてラスターに変換しようとして失敗した後、GDALを使用して範囲と投影を設定することができましたが、データセットは意図した領域の境界に正しく整列しません。下の画像をご覧ください。 このアライメントの欠如を説明することができる、ポーラーステレオグラフィックの受け入れられた地理変換はありますか？ たとえば、使用すべき特定の変換係数または回転はありますか？ データセットのサンプルファイルは次のとおりです。「t201113.grd」 ここに私が現在GDALで使用しているコードがあります ds = gdal.Open("t201113.grd") array = ds.ReadAsArray() x_rotation = 0 y_rotation = 0 xres = 1 yres = -1 llx = -129.8530 lly = 40.0451 ulx = -175.144 uly = 71.385 input_osr = osr.SpatialReference() input_osr.ImportFromWkt(ds.GetProjection()) wgs84_osr = osr.SpatialReference() wgs84_osr.ImportFromEPSG(4326) wgs_to_nps_trans = osr.CoordinateTransformation(wgs84_osr, input_osr) …

16 arcgis-desktop  python  gdal  coordinate-system 

以前のQ＆A（Esriの新しい「ArcGIS for Professionals」および「ArcGIS Pro」とは何ですか？）が、ArcPyコードがArcGIS Proで再利用価値を持つという安心感を提供しました。 ArcGIS Professionalの拡張性。ArcGIS Proの拡張性に関する現在の考えは次のとおりです。ArcGISでジオプロセシングツールを呼び出すPythonスクリプトを記述および実行し、ArcGIS Python API、ArcPyで使用可能な包括的なスクリプト機能スイートを使用して、GISタスクを自動化できます。 ただし、ArcGIS Proが、既存のModelBuilderと同等のGUIを使用してワークフローを自動化するモデル構築を可能にするか、または他のユーザーやIからの以前の拡張要求に対処するための大幅な改良については言及していません。 ArcGIS ProでModelBuilder（または少なくともモデルの構築）の道を聞いた人はいますか？

16 arcgis-desktop  modelbuilder  arcgis-pro  arcgis-10.3 

QGIS 2.2に読み込んだDEMに色を付ける必要があります。DEMの色が異なる標高を表すようにしたいのです。 https://dl.dropboxusercontent.com/u/26086662/color_ramp.jpg QGISに付属するカラーランプはどれも標高に関連するものではないため、QGISに読み込むことができる既存のランプを誰かに教えてもらえますか？ありがとう！

16 qgis  color-ramp 

ESRIジオデータベース形式の.gdbをPostGISに読み込むときに問題が発生しました。データ2.5GBのGDBファイルがあります。インターネットでいくつかのチュートリアルを行ったが、うまくいかないようだ。 データベース「SampleNY」を作成しました コンソールからこれを実行しました： ogr2ogr -f "PostgreSQL" PG： "dbname = SampleNY user = postgres" NYPluto / Pluto.gdb しかし、何も起こりません。エラーも操作も成功しませんでした。手順を逃しましたか？

16 postgis  ogr2ogr  file-geodatabase  translation 

一連の単一境界ポリゴンを使用して土壌ポリゴンのRで真のGISクリップを実行したいのですが、適切に実行するR関数が見つかりません。clipESRIのArcMapの関数と同じように機能するはずです。パッケージover内のメソッドを試しましたが、sppolys over polysでは機能しないようです。 1つの提案は、次のコードを使用して、gIntersectionin rgeosパッケージをクリップとして使用することです。 #------------------------------------ library(rgeos) library(maptools) #Read layers as SpatialPolygonsDataFrame (both the same Albers projection) Soils_poly = readShapePoly("Soils_polygons") #Note - Has 400 polygons clipper_poly = readShapePoly("clipper_polygon") #Note - Has 1 polygon #Try gintersection as clip Clipped_polys = gIntersection(Clipper_Tile_poly, Soils_poly) #----------------------------------- これには、実行に5分かかり（速度が遅すぎる）、次のエラーが発生します。 RGEOSBinTopoFunc（spgeom1、spgeom2、byid、id、drop_not_poly、 "rgeos_intersection"）のエラー：TopologyException：-721459.77681285271 2009506.5980877089で発信dirEdgeが見つかりません また、重複をチェックするためにこのコードを試しました。 gIntersects(Clipper_Tile_poly, Soils_poly) そして結果は真でした。clipESRI …

16 r  polygon  layers  vector  clip