地理情報システム

自己交差するポリゴンシェープファイルの修復について質問があります。自己交差機能を持つポリゴン（ZM）を含むシェープファイルがあります。「ジオメトリの修復」ツールを使用すると、いくつかの自己交差が検出され、修正されます。ほら！ しかし、このファイル（またはその選択）をエクスポートすると、このエクスポートには再び自己交差が含まれます。これは、次のような問題を処理するために行ったすべての試行で発生します。 失敗したツール「ジオメトリの修復」を使用してください。 すべての属性に基づいてポリゴンを溶解しますが、失敗しました。 ZMプロパティを破棄し、上記のすべてを実行すると、失敗します。 ArcCatalogを使用してこのシェープファイルをジオデータベースにインポートしても、正しくインポートされますが、データベースからエクスポートすると、自己交差が再び存在します。 したがって、この問題を修正する方法についてのアドバイスは大歓迎です。このタイプの動作が存在する理由を誰かが説明できれば、それは素晴らしいことです。 ArcGIS 9.3.1 build 4000（SP2）を使用しています。

13 arcgis-desktop  polygon  geometry  arcgis-9.3 

経度/緯度の座標セットがあります。ポリゴンを形成する少なくとも3つの順序付けられた座標が常に存在します（それらが平面に投影された場合）。有効な経度と緯度の範囲のセットとしてこれらの座標の最小境界を計算するにはどうすればよいですか？（「有効」とは、明示的にantimeridianを説明する範囲を意味します）。私が探しているものを説明するのは難しいので、写真があります。 写真のケースは簡単です。すべてのポイントの絶対的な最小値と最大値を見つけるだけです。ただし、これはすべての場合に機能するわけではありません。一般的な解決策はありますか？ 編集：「有効」の意味を明確にするために、データセットに-76、-135、および164の3つの経度値があったとしましょう。 -180および164から180まで。 さらにいくつかの説明。点は多角形を形成するため、特定の場合、必要な範囲は-180〜+180（つまり、360度全体）になります。 左側の画像は、地球の「半分」にある4つの座標の経度を示しています。北極（黒い点）を見下ろしているように想像してください。ピンクは、ポリゴンを含む最小の縦方向の範囲を示しています（ポリゴンは4つのポイントの間に紫色で表示されています）。左側のケースには、[-180〜-120]と[135〜180]（視覚的に推定）の2つの縦方向の範囲があります。 右側の画像は、ポイントが地球を一周する別のケースを示しています。この範囲は[-180〜180]です。

13 coordinates 

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 

ArcGIS Pythonコードとジオプロセシングツールを共有するための最適な組織構造は何ですか？または、コードの共有とツールの共有は別々の質問ですか？ Esriには、Arcgis 9.3および10.0向けに公開されているツール構造の配布方法があります。 しかし、他の場所Also do avoid distributing your code the way its done in Arc Scripts or Code Galleriesでは、ネイティブpython Distutilsを好むようなことを言っています。Esriには10.1（ref）に対応する配布ツールの記事がないようで、反論にある程度の重みを与えています。 GIS.seとは何ですか？ 更新：おそらく遅すぎるかもしれませんが、この質問の核心は、共有に使用されるツール（arcgis online、googleドライブ、dropbox、github、bitbucketなど）が登場する前のファイルおよびフォルダー構造のベストプラクティスに関するものです。演奏する。 Update2：そして、見かけ上孤児のdistutilsアプローチについて誰も発言しませんか？

13 arcgis-desktop  arcgis-10.1  arcpy  geoprocessing 

しばらくの間、TOCにレイヤーを追加することに取り組んでいます。追加したレイヤーを1つに管理しましたが、フォルダーをループしてすべてのレイヤーを追加する必要があります。シェープファイルをまったく追加できません。レイヤーファイルのみ。以下は私がいる場所のコードです： for layer in shp_List: mxd = arcpy.mapping.MapDocument("CURRENT") dataFrame = arcpy.mapping.ListDataFrames(mxd, "*")[0] outlayer = lyr layerfile = os.path.join(base_Folder, lyr + ".lyr") arcpy.MakeFeatureLayer_management(layer, outlayer) arcpy.SaveToLayerFile_management(outlayer, layerfile, "ABSOLUTE") addlayer = arcpy.mapping.Layer(layerfile) arcpy.mapping.AddLayer(dataFrame, addlayer, "BOTTOM") #addLayer = arcpy.mapping.Layer(layer) #arcpy.mapping.AddLayer(dataFrame, addLayer, "BOTTOM") #arcpy.RefreshTOC() #arcpy.RefreshActiveView() del addlayer, mxd 上記のコードは各レイヤーを目次に追加し、イライラさせてそれを削除し、次のレイヤーを追加します。スクリプト全体が終了すると、TOCには何も残りません。このスクリプトをモデルに追加し、派生した出力パラメーターを作成して、「表示に追加」をオンにしてモデルパラメーターとして追加しようとしました。また、ジオプロセシングオプションをチェックして、表示に追加するボックスがオンになっていることを確認しました。Arc 10を実行しています（サービスパックなし）。誰でも助けることができますか？

13 arcgis-desktop  arcgis-10.0  arcpy  shapefile  layers 

私は主にベクターの男ですが、現在のプロジェクトでは、スキャンした古い紙の地図（ロンドンの場合はww2爆弾被害マップ、興味があるなら！） マップをスキャンしてジオリファレンスし、Webサイトで提供するためのタイル合成レイヤーを作成したいと考えています。明らかに境界線を切り取りますが、ここでは問題ではありません。 問題は、マップシート間で見苦しい視覚的な色と明るさの違いがあることです。一貫性のある視覚的な外観を与えるためにそれらを均等化する方法について私は少し迷っています。ヒストグラムイコライズについて調べましたが、現在のツールボックス（Manifold GIS、GDAL、GeoServer）に必要な機能がないようです。 すでにジオリファレンスされた4つのスキャンの例：

13 gdal  imagery  optimization  paper 

PostGISでは、ST_length()関数はメートル単位ではなく、いくつかの単位でラインストリングの長さを返します。 ST_Length()メートル単位のような機能はありますか？

13 postgis  length 

多くのGISチームまたは部門が持っていると確信している問題があります。それは、複数のソースからの複数の異なるデータセットを見つけ、追跡し、管理する方法です。2年前に尋ねられた同様の質問の回答を読みました 大量の地理空間データを管理していますか？ 他のソフトウェアソリューションがあるかどうか疑問に思っていました。VoyagerGISについて知っています。これは見栄えが良いですが、ESRIとリンクしており、同様のオープンソースバージョンがあるかどうか疑問に思いました。 現在の状況の概要を説明するには： 6週間ごとに変わるテラバイトの画像、歴史を残したい 1テラバイトのMapInfo TAB ESRI Shapefilesは6週間ごとに再び変更されます 100 GbpsのCSV、GeoJSON、KML、SQLite、Excelスプレッドシート、JPEG、TIFファイル（一部はジオリファレンス、一部はワールドファイル）、これらは複数の異なるソースから提供されます。 私の完璧なソリューションは、メタデータの更新、ファイル名や属性に基づく検索、視覚画像（VoyagerGISで見られるような）、メタデータ、属性のリストの表示を可能にするGUIインターフェイスだと思います。ただし、データをすべてPostGISに読み込むのではなく、データの元の形式を保持したいと思います。

13 qgis  postgis  open-source-gis  data-management 

プロジェクトのベクターレイヤーから新しいラスターファイルを生成しています。プロジェクトまたはレイヤーファイルと同じディレクトリに保存したいと思います。 PyQGISを使用してこのパスを見つけるにはどうすればよいですか？

13 pyqgis  layers 

ヘルプが必要です... QGISに問題があります。ブールラスターを取得するために再分類したにもかかわらず、null値を持つように見えるラスターがいくつかあります:( NULL値を0に変更するにはどうすればよいですか？

13 qgis  raster  reclassify  nodata 

ラインレイヤーの矢印を作成していました。そのため、三角形を単一のマーカーとして線に追加していました。ここまでは順調ですね。しかし、独自のシングルマーカーシンボルまたはsvgマーカーを作成したい場合、これらの独自のシンボルをQGISに追加するにはどうすればよいですか？これらの独自のシンボルはどこに保存されますか？それらが単一のマーカーまたはsvgマーカー選択リストに追加されるように？ .svgファイルを検索してもフォルダーが見つかりませんでしたか？;（独自の単一マーカーまたはsvgマーカーを選択リストに配置する方法は何ですか？

13 qgis 

GDALを使用してGeoTiffsのラスター属性テーブルを作成するPythonコードがいくつかあります。 ArcGISでは、これらの属性テーブルをうまく見ることができ、（分類として）好きなように動作しますが、QGISでは属性テーブルの概念全体が欠落しているようです。 QGISでRATにアクセスする方法はありますか？

13 qgis  raster  attribute-table 

第一次世界大戦（1914-1918年）の境界線ポリゴンを探しています。この種の情報を保持するリソースはありますか？

13 data  borders 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 2年前に閉店。 簡単なマッピングアプリケーションを開発したいと思います。ユーザーがライン、ポリゴン、ポイントをデジタル化し、フィーチャを属性付けできる場所。属性は、後で評価するためにデータベースに保存されます。 リーフレットライブラリを使用してこのようなアプリを見たことがありますか？

13 leaflet 

QGISでは、特定のレイヤーのポリゴンとラインの頂点の数をカウントし、これらの値を含む別のフィールドを作成する方法はありますか？ 私はそれがフィールド計算機のどこかの機能だと思うが、私はそれを見つけることができません。

13 qgis  vertices  count