地理情報システム

リーフレットマップに、最初のマップビューまたは特定の定義された位置にそれぞれズームするホームボタンのような可能性はありますか？

11 leaflet 

トラックのデータがあります（http://www.chorochronos.org/）。 このデータは、アテネのトラックの複数の軌跡のGPS座標です。 非常に類似しているものを削除するために、軌跡間の類似性を計算する必要があります！ 赤と緑は似ていますが、青、黒、（赤または緑）は異なる軌跡です。赤または緑の類似のものを削除したい。 データはポイント（geometry、latとlong、xとy）（座標gps）であり、画像は軌跡の例です

11 postgis  postgresql  similarity 

同じCRS（緯度と経度）を共有する2つのクラスがあります。 bolognaQuartieriMap：SpatialPolygonDataFrame都市のボロートのデータが含まれています。 crashPoints：SpatialPointsDataFrame事故のデータが含まれています。 彼らはよく使用してプロットされます： plot(bolognaQuartieriMap) title("Crash per quartiere") plot(crashPoints, col="red",add=TRUE) 必要なのはcrashPoints、構成する各ポリゴンのポイント数（）を取得することbolognaQuartieriMapです。使用するよう提案されましたover()が、成功しませんでした。

11 r  sp 

1バンドのGeoTIFFファイルがあります。バンドの色の解釈は灰色です。私は私の変換するために、カラーテーブルを必要Col.tiffとしますCol.jpegかCol.png。 私は.vrtを介してこのtiffにカラーテーブルを入れてみました： gdal_translate -of VRT Col.tiff Col.vrt 次に、.vrtを追加します。 <ColorInterp>Palette</ColorInterp> <ColorTable> <Entry c1="0" c2="0" c3="255" c4="240"/> <Entry c1="0" c2="0" c3="254" c4="240"/> <Entry c1="0" c2="0" c3="253" c4="240"/> <Entry c1="0" c2="0" c3="252" c4="240"/> </ColorTable> そしてもう一度翻訳してください： gdal_translate Col.vrt Col.tiff しかし、それは答えました： Warning 1: Unable to export color table to GeoTIFF file. Color tables can only …

11 gdal  gdal-translate  colortable 

曲線（小刻み）の線を同じ長さのセグメントに分割しようとしています。ある程度機能するが、私が求めている出力を提供しないGrassのv.splitを使用しました。 v.splitは、「最大セグメント長」または「セグメント内の頂点の最大数」に基づいてセグメント化された長さを抽出します。したがって、ラインを10メートルのセグメントに分割し、最大セグメント長を10にしたい場合、すべてのセグメントは、たとえば9.98または9.6のように短くなります。個々のライン、ウィグルまたはストレート。 代わりに、たとえば55mの線を指定します。一方の端から開始し、端に達するまで線に沿って10mごとにノード/ポイントを配置します。上記の場合、それは6つのセグメントに分割され、最後のセグメントは5mの長さで、他のすべてのセグメントは10mです。 より良いモジュールはありますか？

11 qgis  grass 

私は形を使ってポリラインの長さの非常に簡単な計算をしています： from shapely.geometry import LineString ... xy_list = [map(float,e) for e in xy_intm] line = LineString(xy_list) s = '%s,%s,%s' % (fr,to,line.length) 私の座標はWGS84です。shapelyの長さ属性に関する情報が見つからないようです。長さ属性の単位は何ですか？kmまたはメートルに変換する簡単な方法はありますか？

11 shapely  length  units 

私はGeoJSONをOpenLayersマップにレンダリングするアプリケーションに取り組んでいます。 アイデアは、外部サービス（OSMまたはBing）からダウンロードされたPNG画像ではなく、GeoJSONファイルに地理情報を持ち、それらをリアルタイムでレンダリングすることです。私の研究から、このテクノロジーは呼ばれているよう'vector tile'です。私はhttp://openstreetmap.us/~migurski/vector-datasource/ や http://mike.teczno.com/notes/postgreslessness-mapnik-vectiles.html（同じもの）のような彼らのウェブサイトをいくつか見ました 男）。 彼の仕事は私が達成したいもののようですが、私はそれがどのように機能するのか本当にわかりません。 私の理解では、OpenLayersを使用しているため、クライアントはGeoJSON形式のマップデータを取得し、OpenLayersのGeoJSONパーサーを使用してそれらをベクターに転送します。しかし、私が理解していないのは、このテクノロジーで使用されている用語だけです。 彼のブログで、彼は'MapNik Vector Tile'フォーマットと呼ばれるデータを提供するサーバーがあると述べました。ベクタータイルとGeoJSONの違いは何ですか？それは単に、境界ボックスに対応するタイル状のGeoJSONを意味するのでしょうか？それとも新しい画像タイプですか？ 彼の例を見ると、道路のラベルが道路のストロークに揃えられていることがわかります。OpenLayersのレンダラーコードを変更してOpenLayersでも同じことを試みましたが、ラベルを回転させることはできますが、道路に沿ってラベルを曲げることはできません（たとえば、道路に角がある場合、マップはラベルを「L」字形で表示します。* * + 非常によくレンダリングされます！ Mike（ブログの作成者）はOpenLayers以外のいくつかの異なるマップを使用しているようですが、彼がどのようにそれを行ったかはわかりません。質問がたくさんありますが、最初に上記を理解したいと思います。

11 openlayers-2  openstreetmap  vector  vector-tiles 

SQL 2008に取得しようとしている数万のポイントで構成される数百の形状（polygonsとmultipolygons）があります。 残念ながら、インポートしようとしたシェイプは「右利き」です（各シェイプの周囲は、含まれているポイントの周りに時計回りに描かれています）。SQLサーバーは、少なくともgeography型については、「左回り」の形状（内部で反時計回り）を想定しています。これは、SQLが、自分の形状を除いて地球全体を選択しようとしていることを前提としていることを意味します。一部の人々は、これを「裏返し」の形状と表現しています。 MSDNから、どのリングの向きを使用するべきかについてイライラするほど言っていません： geographyデータタイプを使用して空間インスタンスを格納する場合、リングの方向を指定し、インスタンスの場所を正確に記述する必要があります。 SQL 2008で間違ったリング方向を使用すると、次のエラーでクラッシュします（強調は私のものです）。 ユーザー定義ルーチンまたは集計 "geography"の実行中に.NET Frameworkエラーが発生しました：Microsoft.SqlServer.Types.GLArgumentException：24205：指定された入力は、単一の半球を超えているため、有効なgeographyインスタンスを表していません。各geographyインスタンスは、単一の半球内に収まる必要があります。このエラーの一般的な理由は、ポリゴンのリングの向きが間違っていることです。 geometryではなくとしてシェイプをインポートすることはgeography問題ありませんが、geography可能であれば使用したいと思います。 SQL 2012では、この問題を修正するのはかなり簡単ですが、私は2008に縛られています。 形状を変換するにはどうすればよいですか？

11 sql  sql-server-spatial  geography-data-type 

私は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 

GeoJSON specは「WGS84データムを使用して、地理座標参照系であり、小数度の経度と緯度ユニットと、デフォルトのCRS」を言います。ポイントの標高/高度/ zの単位を指定していないようです。 GeoJSONのデフォルトの高度単位はありますか？そうでない場合、一般的な規則、または単位を指定する方法はありますか？

11 geojson 

私は主にQGISを使用しており、完了するには大きなプロジェクトがあります。完了したら、データをArcGISユーザーに送信する必要があります。データは主にポイントとポリゴンの形状ファイルで構成されます。 ArcGIS Desktopを使用する人にQGISプロジェクトを引き継ぐことがわかっている場合、最初からどのように設定しますか？

11 qgis  arcgis-desktop 

長い間、Windowsベースのソフトウェアを使用して、フォレストアプリケーションのLiDARデータセットを視覚化および分析してきました。最近、すべての作業をUbuntuプラットフォームに移行するようになりましたが、LiDARデータセットで動作する優れたソフトウェアをまだ探しています。 Windowsに関連して、私が体験した最高のフリーソフトウェアは次のとおりです。 視覚化のためだけに： a）クイックテレインリーダーは、事前に構築されたデジタル標高モデル（DEM）と点群を開くことができ、ユーザーが地形をすばやく直感的に自由に移動できるようにします。 b）FugroViewerは、ユーザーが地理空間データを最大限に活用できるように設計された、堅牢で使いやすいフリーウェアです。 c）PointVue LEは無料の3D LIDAR視覚化ツールで、LIDARデータをASPRS LASバージョン1.1形式で視覚化するために使用できます。 分析および視覚化の場合： a）FUSION / LDVは、LiDARおよび地形データセットへの高速で効率的かつ柔軟なアクセスを提供します。 Ubuntuに引っ越してから、慣れ親しんでいたソフトウェアで代用するのはかなり大変でした。今までのところ、いくつかの代替案を見つけましたが、Windowsのものほどよくありません。 視覚化のため： a）LiDAR視覚化は、コア外のマルチ解像度ポイントクラウドレンダラーとして実装されました。レンダラーは、固定サイズのメモリキャッシュを使用して、インタラクティブなフレームレートで、最大113億（11.3 * 109）個のサンプルポイントを含む、現在の最大のLiDARスキャンを視覚化できます。 b）LAGは、LiDAR点群の視覚化、検査、および分類のためのツールです。現在、LASおよびASCIIファイル形式をサポートしています。 残念ながら、FUSIONのような処理するものは何も見つかりませんでした。 LiDARファイル分析に関連するLinuxの優れた無料ソフトウェアを経験した人はいますか？ 他のソフトウェアを列挙するだけでなく、なぜそれを使用し、推奨するかを説明してください！

11 geoprocessing  lidar  visualisation  linux  fusion 

あるフィーチャクラス属性を別の属性クラスにコピーすることで、同僚のワークフローをスピードアップするのに役立つPythonアドインボタンを作成しました。これは、arcpy.UpdateCursor関数を使用して、ターゲットフィーチャクラスの行を更新します。現在存在するため、このボタンスクリプトは編集モードに関係なく実行できます。明らかに編集セッションで実行する場合、ユーザーは編集を停止して変更を保存しないことを選択できますが、スクリプトが編集セッション外で実行される場合はそうではありません。 ArcMapが現在編集セッションにない場合にスクリプトの実行を停止するチェックをスクリプトに追加するにはどうすればよいですか？ これは、ArcMap 10および10.1に関係します。 また、他のArcMapユーザーに確認して、編集セッションに参加しない限り、テーブルの更新が通常許可されないことを確認したいと思います。 では、このスクリプトは編集セッションの外でどのように実行されるのでしょうか。 このスクリプトは、リストから2番目のフィーチャクラステーブルを更新したときに、たまたま機能するArcMapが実行する一見偶然のような選択順序に関する別の質問も表示しますが、それは別の日です。 これが現在動作するスクリプトです（10.1エディター実装なし）。 ユーザーが編集セッションにいることを確認するチェックを追加するにはどうすればよいですか？ def onClick(self): #Reference mxd mxd = arcpy.mapping.MapDocument("CURRENT") #Reference the main Data frame mm = arcpy.mapping.ListDataFrames(mxd, "MainMap")[0] #Reference the Water System Valve feature class waterValves = arcpy.mapping.ListLayers(mxd, "Water System Valve", mm)[0] #Reference the fire hydrant feature class fireHydrants = arcpy.mapping.ListLayers(mxd, "Water Hydrant", …

11 arcpy  arcmap  editing 

別のPostgresレイヤーに新しい要素（機能）を入力する場合、次の2つの方法で行うことができます。 まれに新しい要素を描画する（「機能の追加」を使用）または 別のPostgreレイヤー（ソースレイヤー）からいくつかの要素をコピー（またはカット）し、頻繁に行うターゲットレイヤーに貼り付けます 最初の例では、このレイヤーがpostgreデータベースシーケンスからgidを取得するため、編集の保存は正常に機能します* nextval（ 'layer_name_gid_seq' :: regclass）* 2番目の例では、ソースレイヤーからターゲットレイヤーに要素をコピーしているときにqgis がソースレイヤーから要素のgidをコピーしたため、編集の保存中にエラーが発生しました 。編集を保存しようとすると、次のエラーが返されます。 レイヤー「Cjevovodi」への変更をコミットできませんでした エラー：エラー：1つの機能が追加されていません。 プロバイダーエラー： 機能の追加中のPostGISエラー：エラー：重複するキー値が一意の制約 "cjevovodi_okill_pkey"に違反しています 詳細：キー（gid）=（5）は既に存在します。 * nextval（ 'layer_name_gid_seq' :: regclass）*をフィールドgidにコピーしようとしましたが、フィールドが数値として定義されているため、このシーケンスをフィールドgidに貼り付けることはできません。 （既存のgidで）ソースレイヤーから要素をコピーして新しいgidを割り当てる簡単な方法を知っている人はいますか？ ありがとう！

11 qgis  postgis 

熱帯雨林における道路網の拡大に関する研究のために、私はランドサット画像から道路を抽出しようとしています。道路が目ではっきりと見えるシャープで雲のないコンポジットがすでにありますが、ラインフィーチャにそれらを抽出するのは難しいことがわかっているので、Landsatが提供する大きな画像を処理できる優れたアルゴリズムまたは方法を誰かが知っているかどうか疑問に思いましたか？Grassのr.thinを試しましたが、うまくいかないようです。

11 raster  remote-sensing  landsat  road