Zoom dari satu fitur ke fitur lainnya dalam mode edit?

8

Saya harus mengedit secara manual shapefile beberapa ratus poin.

Saya ingin cara cepat untuk melompat dari satu titik ke titik berikutnya baik dalam arti atribut dan visual / spasial yaitu saya ingin sementara dalam mode edit melompat dari objek ID 1 ke Obyek ID 2 tanpa harus membuka tabel atribut , pilih titik berikutnya, zoom untuk pemilihan, dll.

Semacam tombol "Next" yang akan mempercepat proses manual.

pengguna32882
sumber
Apakah Anda memiliki kemampuan ArcObjects? Saya perhatikan kekurangan ini dan menulis alat untuk itu. Saya senang membagikan kode jika Anda dapat menggunakannya.
Michael Stimson
Sayangnya saya tidak meskipun saya akan senang untuk belajar ...
user32882
Anda perlu menginstal SDK dari media instal ArcGis Anda, tetapi pertama-tama Anda perlu mendapatkan Visual Studio (Express). Lihatlah jawaban saya gis.stackexchange.com/questions/62720/… untuk prasyaratnya. Saya perlu menyebutkan bahwa shapefile memiliki masalah di mana FID berubah pada save, jadi untuk menggunakan alat ini pada shapefile, tidak perlu menyimpan sampai akhir.
Michael Stimson
@ MichaelMiles-Stimson, apakah Anda memiliki alat tambahan yang bisa Anda bagikan sehingga pengguna tidak perlu menginstal SDK dan VS?
artwork21
sebuah tambahan akan sangat bagus ....
user32882

Jawaban:

6

Bagian pertama adalah AddIn, pekerjaan nyata dilakukan pada formulir:

Inherits ESRI.ArcGIS.Desktop.AddIns.Button
Private pForm As fFeatureInspector
Public Shared IsFormLoaded As Boolean = False

Public Sub New()

End Sub

Protected Overrides Sub OnClick()
    'My.ArcMap.Application.CurrentTool = Nothing
    If Not IsFormLoaded Then
        pForm = New fFeatureInspector
        pForm.pApp = CType(My.ArcMap.Application, ESRI.ArcGIS.ArcMapUI.IMxApplication)
        pForm.Show()
    Else
        pForm.sResetList()
    End If

End Sub

Protected Overrides Sub OnUpdate()
    Enabled = My.ArcMap.Application IsNot Nothing
End Sub

Ketika Anda membuat addin baru, sebagian besar sudah ada untuk Anda. Kemudian tambahkan formulir ke proyek (nama fFeatureInspector atau Anda perlu mengubahnya beberapa kali dalam kode).

masukkan deskripsi gambar di sini

Penting untuk mendapatkan nama yang benar atau Anda perlu menemukan & mengganti kode formulir. Toolbox untuk formulir memiliki semua kontrol yang umum: tombol , kotak centang , listbox , combobox .

Cara kerjanya adalah alat ini mendapatkan semua fitur yang dipilih dan dapat diedit, menyalin nama mereka dan OID / FID ke dalam kotak daftar dan kemudian ketika seseorang disorot akan memilihnya (setelah menghapus pilihan pertama) dan memperbesarnya. Ada tombol simpan dan muat untuk menyimpan inspeksi, kembali satu dan meneruskan satu, simpan otomatis centang dan tombol reset. Alat akan menyegarkan saat dimuat tetapi setelah itu Anda dapat menyegarkan kapan saja. Hemat otomatis tidak kompatibel dengan pengeditan shapefile karena FID tidak statis dan dikompresi saat disimpan.

Poin memiliki luas 0 lebar sehingga penting untuk mengatur skala min ke sesuatu yang realistis; Zoom% adalah berapa banyak lebih dari poligon / garis yang ingin Anda lihat di sekitarnya.

Ini adalah kode form (maaf karena tidak ada komentar):

Imports ESRI.ArcGIS.Framework
Imports ESRI.ArcGIS.ArcMapUI
Imports ESRI.ArcGIS.Carto
Imports ESRI.ArcGIS.Geometry
Imports ESRI.ArcGIS.Geodatabase
Imports ESRI.ArcGIS.esriSystem
Imports ESRI.ArcGIS.Editor
Imports ESRI.ArcGIS.Display

Public Class fFeatureInspector
    Const FormCaption As String = "Feature Inspector (22 Feb 10)"
    Const FormName As String = "fFeatureClass"
    Public pApp As IApplication
    Private pDoc As IMxDocument
    Private pMap As IMap

    Dim pEd As IEditor2
    Dim pID As UID = New UID
    Dim pFeatFrom() As String
    Dim pFeatWS As IFeatureWorkspace
    Dim pWS As IWorkspace
    Dim pFeatOID() As Long
    Dim pFeatCnt As Long
    Dim pInRefresh As Boolean
    Dim pPointExtent As IEnvelope
    Dim pSaveEdits As ICommandItem
    Dim pLoadTime As Long
    Dim pNow As Date
    Dim pStartIndex As Long

    Dim vStartTime As DateTime
    Dim vCurrentTime As DateTime

    Private Sub fFeatureInspector_Disposed(ByVal sender As Object, ByVal e As System.EventArgs) Handles Me.Disposed
        StartFeatureInspector.IsFormLoaded = False
    End Sub
    Private Sub fFeatureInspector_Load(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles MyBase.Load
        pEd = pApp.FindExtensionByName("Esri Object Editor")
    End Sub
    Private Sub form1_Move(ByVal sender As Object, ByVal e As System.EventArgs) Handles Me.Move
        Dim pOutFile As Integer
        Dim pTempDir As String = Environ("Temp")

        If Me.Visible Then
            pOutFile = FreeFile()
            FileOpen(pOutFile, pTempDir & "\" & FormName & ".xy", OpenMode.Output)
            WriteLine(pOutFile, Me.Left & "," & Me.Top)
            FileClose(pOutFile)
        End If
    End Sub
    Private Sub Form1_Shown(ByVal sender As Object, ByVal e As System.EventArgs) Handles Me.Shown
        StartFeatureInspector.IsFormLoaded = True
        Me.Text = FormCaption
        pDoc = CType(pApp.Document, IMxDocument)

        Dim pTempDir As String = Environ("temp")
        Dim pInFile As Integer
        Dim pReadString As String = ""
        Dim pResyk As String = ""
        Dim pXpos As Integer = 0
        Dim pYpos As Integer = 0

        fZoomPercent.Items.Add(110)
        fZoomPercent.Items.Add(150)
        fZoomPercent.Items.Add(200)
        fZoomPercent.Text = "110"

        fPointScale.Items.Add(200)
        fPointScale.Items.Add(500)
        fPointScale.Items.Add(1000)
        fPointScale.Items.Add(2500)
        fPointScale.Items.Add(10000)
        fPointScale.Text = "1000"

        If My.Computer.FileSystem.FileExists(Environ("Temp" & "\" & FormName & ".xy")) Then
            pApp.StatusBar.Message(0) = "Loading position"
            pInFile = FreeFile()
            FileOpen(pInFile, pTempDir & "\" & FormName & ".xy", OpenMode.Input)
            pReadString = LineInput(pInFile)
            pReadString = Mid(pReadString, 2, Len(pReadString) - 2)
            pApp.StatusBar.Message(0) = pReadString

            pResyk = Microsoft.VisualBasic.Left(pReadString, InStr(pReadString, ",") - 1)
            pApp.StatusBar.Message(0) = pResyk
            pXpos = CInt(pResyk)

            pApp.StatusBar.Message(0) = "Xposition " & pXpos
            pResyk = Microsoft.VisualBasic.Right(pReadString, Len(pReadString) - InStr(pReadString, ","))
            pApp.StatusBar.Message(0) = pResyk
            pYpos = CInt(pResyk)
            pApp.StatusBar.Message(0) = "Yposition " & pYpos
            FileClose(pInFile)
            Me.Left = pXpos
            Me.Top = pYpos
        End If
        sResetList()
        pID.Value = "{59D2AFD2-9EA2-11D1-9165-0080C718DF97}"
        Dim pComBars As ICommandBars = pApp.Document.CommandBars
        pSaveEdits = pComBars.Find(pID, False, False)

    End Sub
    Private Sub fSaveButton_Click()
        Dim pOutfile As Integer
        Dim cnt As Long

        pOutfile = FreeFile()
        FileOpen(pOutfile, (Environ("Temp") & "\" & "FeatInspect"), OpenMode.Output, OpenAccess.Write)

        Print(pOutfile, pFeatCnt & vbNewLine)
        For cnt = 0 To pFeatCnt - 1
            Print(pOutfile, pFeatFrom(cnt) & "|" & pFeatOID(cnt) & vbNewLine)
        Next cnt
        Print(pOutfile, fFeatureList.SelectedIndex & vbNewLine)
        Print(pOutfile, fZoomPercent.Text & vbNewLine)
        Print(pOutfile, fPointScale.Text & vbNewLine)
        FileClose(pOutfile)
    End Sub
    Private Sub fSaveButton_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles fSaveButton.Click
        fSaveButton_Click()
    End Sub
    Private Sub fLoadButton_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles fLoadButton.Click
        Dim cnt As Long
        Dim pInFile As Integer
        Dim pReadString As String
        Dim pSplitString() As String

        pInRefresh = True
        fFeatureList.Items.Clear()
        pInFile = FreeFile()
        FileOpen(pInFile, (Environ("Temp") & "\" & "FeatInspect"), OpenMode.Input, OpenAccess.Read)
        pReadString = LineInput(pInFile)
        pFeatCnt = pReadString
        ReDim pFeatFrom(pFeatCnt)
        ReDim pFeatOID(pFeatCnt)

        fFeatureList.Items.Clear()

        For cnt = 0 To pFeatCnt - 1
            pReadString = LineInput(pInFile)
            pSplitString = Split(pReadString, "|")

            pFeatFrom(cnt) = pSplitString(0)
            pFeatOID(cnt) = pSplitString(1)
            fFeatureList.Items.Add(pFeatFrom(cnt) & " - " & pFeatOID(cnt))
        Next cnt
        pInRefresh = False
        pReadString = LineInput(pInFile)
        fFeatureList.SelectedIndex = pReadString
        pReadString = LineInput(pInFile)
        fZoomPercent.Text = pReadString
        pReadString = LineInput(pInFile)
        fPointScale.Text = pReadString
        FileClose()
        pStartIndex = fFeatureList.SelectedIndex
        pNow = Now()
        pLoadTime = (Hour(pNow) * 3600) + (Minute(pNow) * 60) + Second(pNow)

    End Sub

    Private Sub fBackButton_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles fBackButton.Click
        If fFeatureList.SelectedIndex = 0 Then
            MsgBox("But you're already at the start!")
            Exit Sub
        End If
        fFeatureList.SelectedIndex = fFeatureList.SelectedIndex - 1
    End Sub
    Private Sub fGoDown_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles fGoDown.Click
        If fFeatureList.SelectedIndex = fFeatureList.Items.Count - 1 Then
            MsgBox("That's all there is, there isn't anymore.")
            Exit Sub
        End If
        fFeatureList.SelectedIndex = fFeatureList.SelectedIndex + 1
    End Sub
    Private Sub bReset_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles bReset.Click
        pInRefresh = True
        sResetList()
        pInRefresh = False
        fFeatureList_Change()
        pLoadTime = (Hour(pNow) * 3600) + (Minute(pNow) * 60) + Second(pNow)
    End Sub
    Public Sub sResetList()
        Dim pEnumFeat As IEnumFeature
        Dim pFeature As IFeature
        Dim pFeatClass As IFeatureClass

        If pEd.EditState = esriEditState.esriStateNotEditing Then Exit Sub
        pEnumFeat = pEd.EditSelection

        pFeature = pEnumFeat.Next
        If pFeature Is Nothing Then
            MsgBox("Nothing selected", vbCritical)
            Exit Sub
        End If
        pFeatCnt = pEd.SelectionCount - 1

        ReDim pFeatFrom(pFeatCnt)
        ReDim pFeatOID(pFeatCnt)

        pFeatCnt = 0

        fFeatureList.Items.Clear()

        Do Until pFeature Is Nothing
            pFeatClass = pFeature.Class
            pFeatFrom(pFeatCnt) = pFeatClass.AliasName
            pFeatOID(pFeatCnt) = pFeature.OID
            fFeatureList.Items.Add(pFeatFrom(pFeatCnt) & " - " & pFeatOID(pFeatCnt))
            pFeatCnt = pFeatCnt + 1
            pFeature = pEnumFeat.Next
        Loop

        pEd.Map.ClearSelection()
        fFeatureList.SelectedIndex = 0
        pNow = Now()
        pLoadTime = (Hour(pNow) * 3600) + (Minute(pNow) * 60) + Second(pNow)

    End Sub

    Private Sub fFeatureList_SelectedIndexChanged(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles fFeatureList.SelectedIndexChanged
        fFeatureList_Change()
    End Sub
    Private Sub fFeatureList_Change()
        Dim pFeatClass As IFeatureClass
        Dim pFeature As IFeature
        Dim pSelection As ISelection
        Dim pLayer As ILayer
        Dim pEnumLayer As IEnumLayer
        Dim pEditLayers As IEditLayers
        Dim pFeatLayer As IFeatureLayer
        Dim pEnv As IEnvelope2
        Dim pDispTran As IDisplayTransformation
        Dim pThisTime As Long
        Dim pAvVis As Single
        Dim pTotTime As Long
        Dim pDeltaVis As Long
        Dim pLeft As Long
        Dim pTPF As Long
        Dim pAvVisStr As String
        Dim pETAstr As String
        Dim pPoint As ESRI.ArcGIS.Geometry.IPoint
        Dim cnt As Long

        If pInRefresh Then Exit Sub
        If fFeatureList.SelectedIndex < 0 Then Exit Sub

        If pEd.EditState = esriEditState.esriStateNotEditing Then
            MsgBox("This tool only works on EDIT features" & vbNewLine & "Please start editing", vbCritical)
            Exit Sub
        End If
        pFeatWS = pEd.EditWorkspace
        pFeatClass = pFeatWS.OpenFeatureClass(pFeatFrom(fFeatureList.SelectedIndex))
        On Error Resume Next
        pFeature = pFeatClass.GetFeature(pFeatOID(fFeatureList.SelectedIndex))
        If pFeature Is Nothing Then
            MsgBox("Feature not found", vbCritical) ' comment this out if you don't want to see errors
        End If
        pEd.Map.ClearSelection()
        pID.Value = "{6CA416B1-E160-11D2-9F4E-00C04F6BC78E} "
        pEnumLayer = pEd.Map.Layers(pID, True)
        pEditLayers = pEd

        pLayer = pEnumLayer.Next
        Do Until pLayer Is Nothing
            If TypeOf pLayer Is IFeatureLayer Then
                If pEditLayers.IsEditable(pLayer) And pLayer.Visible = True Then
                    pFeatLayer = pLayer
                    If pFeatLayer.Selectable Then
                        If pFeatLayer.FeatureClass.AliasName = pFeatFrom(fFeatureList.SelectedIndex) Then
                            pEd.Map.SelectFeature(pLayer, pFeature)
                            If pFeatLayer.FeatureClass.ShapeType = esriGeometryType.esriGeometryPoint Then
                                pEnv = New Envelope
                                pEnv.PutCoords(pDoc.ActiveView.Extent.XMin, pDoc.ActiveView.Extent.YMin, pDoc.ActiveView.Extent.XMax, pDoc.ActiveView.Extent.YMax)
                                pEnv.SpatialReference = pFeature.Shape.SpatialReference
                                If pEnv.SpatialReference.FactoryCode <> pDoc.FocusMap.SpatialReference.FactoryCode Then pEnv.Project(pDoc.FocusMap.SpatialReference)
                                pPoint = New ESRI.ArcGIS.Geometry.PointClass
                                pPoint = pFeature.ShapeCopy
                                pEnv.CenterAt(pPoint)
                                pDoc.ActiveView.Extent = pEnv
                                pDispTran = pDoc.ActiveView.ScreenDisplay.DisplayTransformation
                                If Len(fPointScale.Text) > 0 Then
                                    pDispTran.ScaleRatio = Int(fPointScale.Text)
                                Else
                                    pDispTran.ScaleRatio = 1000
                                End If
                                If fSaveOnNext.Checked Then
                                    pSaveEdits.Execute()
                                    fSaveButton_Click()
                                End If 'If fSaveOnNext.Checked Then
                                pDoc.ActiveView.Refresh()
                            Else

                                If Not pFeature.Shape.Envelope.IsEmpty Then
                                    pEnv = New Envelope
                                    pEnv.PutCoords(pFeature.Extent.XMin, pFeature.Extent.YMin, pFeature.Extent.XMax, pFeature.Extent.YMax)
                                    pEnv.SpatialReference = pFeature.Shape.SpatialReference
                                    If pEnv.SpatialReference.FactoryCode <> pDoc.FocusMap.SpatialReference.FactoryCode Then pEnv.Project(pDoc.FocusMap.SpatialReference)
                                    If Len(fZoomPercent.Text) > 0 Then
                                        pEnv.Expand(Int(fZoomPercent.Text) / 100, Int(fZoomPercent.Text) / 100, True)
                                    End If
                                    pDoc.ActiveView.Extent = pEnv
                                    If fPointScale.Text.Length > 0 Then
                                        If pDoc.ActiveView.ScreenDisplay.DisplayTransformation.ScaleRatio < Int(fPointScale.Text) Then pDoc.ActiveView.ScreenDisplay.DisplayTransformation.ScaleRatio = Int(fPointScale.Text)
                                    End If
                                    If fSaveOnNext.Checked Then
                                        pSaveEdits.Execute()
                                        fSaveButton_Click()
                                    End If 'If fSaveOnNext.Checked Then
                                    pDoc.ActiveView.Refresh()
                                End If 'Not pFeature.Shape.Envelope.IsEmpty
                            End If 'pFeatLayer.FeatureClass.ShapeType = esriGeometryType.esriGeometryPoint Then
                        End If 'pFeatLayer.FeatureClass.AliasName = pFeatFrom(fFeatureList.SelectedIndex) Then
                    End If 'pFeatLayer.Selectable Then
                End If 'pEditLayers.IsEditable(pLayer) And pLayer.Visible = True Then
            End If
            pLayer = pEnumLayer.Next
        Loop 'Until pLayer Is Nothing

        pNow = Now()
        pThisTime = (Hour(pNow) * 3600) + (Minute(pNow) * 60) + Second(pNow)
        pTotTime = pThisTime - pLoadTime
        pDeltaVis = fFeatureList.SelectedIndex - pStartIndex
        If pDeltaVis <= 0 Then
            fProgressLabel.Text = "Unable to Calculate"
            Exit Sub
        Else
            pAvVis = pTotTime / pDeltaVis
            pLeft = fFeatureList.Items.Count - fFeatureList.SelectedIndex + 1
            pETAstr = fLongTime_to_TimeString(pLeft * pAvVis)

            fProgressLabel.Text = pDeltaVis & " Inspected of " & fFeatureList.Items.Count & ". ETA " & pETAstr
            Me.Update()

        End If
    End Sub

    Private Function fLongTime_to_TimeString(ByVal pLongTime As Long) As String
        Dim pRemainder As Long
        Dim pHour As Integer
        Dim pMin As Integer
        Dim pSec As Integer

        pRemainder = pLongTime Mod 3600
        pHour = pLongTime - pRemainder
        If pHour > 0 Then pLongTime = pLongTime - pHour
        pRemainder = pLongTime Mod 60
        pMin = pLongTime - pRemainder
        If pMin > 0 Then pLongTime = pLongTime - pMin
        pSec = pLongTime

        pHour = pHour / 3600
        pMin = pMin / 60

        fLongTime_to_TimeString = pHour & ":" & pMin & ":" & pSec
    End Function


End Class

Sebanyak yang saya tidak suka berbagi kode kompilasi, Inilah tautannya . Silakan baca dokumen Esri pada 'Berbagi dan menambahkan Addins' .

Michael Stimson
sumber
Terima kasih banyak atas kemurahan hati Anda. Ini akan sangat instruktif dan akan menyelamatkan saya banyak waktu
user32882
7

Berikut ini adalah versi perbesaran untuk fitur selanjutnya. Anda dapat menjalankan ini di jendela python ArcMap Anda:

mxd = arcpy.mapping.MapDocument("CURRENT") # currently opened map doc
df = arcpy.mapping.ListDataFrames(mxd, "Layers") [0]

# define layer you want to iterate and zoom on
for lyr in arcpy.mapping.ListLayers(mxd):
    if lyr.name == 'myTOCLayerNameHere':
        fc = lyr

# get total record count of fc
with arcpy.da.SearchCursor(fc, ["FID"]) as cursor:
    for row in cursor:
        totalCount+=1

def selectZoomNext(fc, field, record):
    if record > totalCount:
        record = 0 # reset to first feature
    expression = '{} = {}'.format(field, record)
    arcpy.SelectLayerByAttribute_management (fc, "NEW_SELECTION", expression)
    df.zoomToSelectedFeatures()
    nextRecord = record + 1
    return nextRecord

record = 0
record = selectZoomNext(fc, 'FID', record) # second argument is the field name, this could be OBJECTID too

Anda dapat terus menjalankan record = selectZoomNext(fc, 'FID', record)pernyataan untuk terus memilih fitur berikutnya dalam tabel dan memperbesarnya. Anda juga bisa memasukkan potongan ini ke alat python add-in atau skrip python. Selain itu, untuk mempermudah saat mengedit, Anda dapat mematikan bidang yang tidak perlu (di properti layer) dan juga membuka panel Atribut untuk akses atribut cepat.

artwork21
sumber
Itu juga berhasil. Saya pikir ini adalah ilustrasi yang bagus tentang perbedaan antara ArcObjects dan Python dalam panjang dan kompleksitas kode. Perhatikan bahwa ini akan memiliki masalah yang sama dengan shapefile dan tidak akan berfungsi pada kelas fitur file / personal / SDE geodatabase karena nilai OID / FID tidak dijamin bersebelahan dan berbasis 0 - itu adalah kekhasan dari format shapefile.
Michael Stimson
@ MichaelMiles-Stimson, ya saya akan merekomendasikan menggunakan metode AddFieldDelimiters untuk membangun ekspresi bidang yang benar antara file, .mdb, dan .gdb sumber help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//…
artwork21
1
Bukan itu, dalam kelas fitur geodatabase nilai-nilai OID tidak perlu mulai dari 0 dan bahkan jika mereka melakukan yang berikutnya tidak perlu 1. Anda perlu memasukkan OBJECTIDs ke dalam daftar atau kamus ketika Anda menghitung dengan kursor, Anda juga perlu menggunakan arcpy.Description (fc) .OIDFieldName sebagai bidang karena geodatabase memiliki OBJECTID dan bukan FID.
Michael Stimson
3

Apakah Anda memiliki ekstensi Peninjau Data ? Data Reviewer memungkinkan Anda untuk "Menelusuri" melalui semua fitur dengan mengklik sederhana tombol "Next" (memperbesar lokasi spasial dan catatan tabel atribut). Ada lebih banyak fungsi untuk Data Reviewer selain dari ini (seperti menandai kesalahan sebagai "diperbaiki", "ditandai" dll. Dan menjalankan pekerjaan batch). Hanya alat off-the-shelf, meskipun saya yakin alat Anda @Michael juga sangat fantastis!

masukkan deskripsi gambar di sini

John
sumber
Sepertinya hanya apa yang saya butuhkan, Namun sepertinya doo tingkat lisensi saya tidak termasuk peninjau data :(
user32882
Saya menulis alat saya bahkan sebelum saya mendengar tentang PLTS (sekarang pemetaan produksi), inspirasi asli ditulis ke dalam skrip AML karena persyaratannya sangat mendasar tetapi perlu. Saya tidak terkejut Esri memiliki alat tetapi mengapa bukan yang asli di ArcMap? Peninjau data sangat baik dalam hal apa yang dilakukannya tetapi saya menemukan bahwa saya sudah menulis sebagian besar alat yang datang ketika saya mengevaluasinya.
Michael Stimson