Zoom d'une fonction à une autre en mode édition?

8

Je dois modifier manuellement un fichier de formes de quelques centaines de points.

Je voudrais un moyen rapide de passer d'un point à l'autre dans le sens de l'attribut et dans le sens visuel / spatial, c'est-à-dire que je voudrais en mode édition passer de l'ID d'objet 1 à l'ID d'objet 2 sans avoir à ouvrir la table des attributs , sélectionnez le point suivant, zoomez sur la sélection, etc.

Une sorte de bouton "Suivant" qui accélérerait le processus manuel.

user32882
la source
Avez-vous la capacité ArcObjects? J'ai remarqué cette lacune et j'ai écrit un outil pour cela. Je suis heureux de partager le code si vous pouvez l'utiliser.
Michael Stimson
Malheureusement, je ne le sais pas bien que je serais heureux d'apprendre ...
user32882
Vous devez installer le SDK à partir de votre support d'installation ArcGis, mais vous devez d'abord obtenir Visual Studio (Express). Jetez un œil à ma réponse gis.stackexchange.com/questions/62720/… pour les prérequis. Je dois mentionner que les fichiers de formes ont un problème lorsque le FID change lors de l'enregistrement, donc pour utiliser cet outil sur un fichier de formes, il est nécessaire de ne pas enregistrer jusqu'à la fin.
Michael Stimson
@ MichaelMiles-Stimson, avez-vous un complément de cet outil que vous pourriez partager pour que l'utilisateur n'ait pas à installer le SDK et VS?
artwork21
un complément serait très bien ....
user32882

Réponses:

6

La première partie est le complément, le vrai travail se fait sur un formulaire:

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

Lorsque vous créez un nouveau complément, la plupart de ces informations sont déjà là pour vous. Ajoutez ensuite un formulaire au projet (nom fFeatureInspector ou vous devrez le changer plusieurs fois dans le code).

entrez la description de l'image ici

Il est important d'obtenir les noms corrects ou vous devrez faire rechercher et remplacer dans le code du formulaire. La boîte à outils du formulaire possède tous les contrôles communs: bouton , case à cocher , zone de liste , zone de liste déroulante .

Comment cela fonctionne-t-il? L'outil obtient toutes les fonctionnalités sélectionnées et modifiables, copie leur nom et leur OID / FID dans la zone de liste, puis lorsque l'un est mis en surbrillance, il le sélectionnera (après avoir d'abord effacé la sélection) et zoomera dessus. Il y a un bouton d'enregistrement et de chargement pour enregistrer l'inspection, un en arrière et un en avant, un bouton d'enregistrement automatique et un bouton de réinitialisation. L'outil s'actualisera lorsqu'il sera chargé, mais vous pourrez ensuite l'actualiser à tout moment. L'enregistrement automatique n'est pas compatible avec l'édition de fichiers de formes car le FID n'est pas statique et est compressé lors de l'enregistrement.

Les points ont une étendue de largeur 0, il est donc important de définir une échelle min à quelque chose de réaliste; Le% de zoom est bien plus qu'un polygone / ligne que vous voulez voir autour de lui.

Voici le code du formulaire (désolé pour le manque de commentaires):

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

Autant que je n'aime pas partager le code compilé, voici le lien . Veuillez lire le document Esri sur «Partager et ajouter des compléments» .

Michael Stimson
la source
Merci beaucoup pour votre générosité. Ce sera très instructif et me fera gagner beaucoup de temps
user32882
7

Voici la version arcpy du zoom vers la fonctionnalité suivante. Vous pouvez exécuter ceci dans votre fenêtre python ArcMap:

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

Vous pouvez continuer à exécuter l' record = selectZoomNext(fc, 'FID', record)instruction pour continuer à sélectionner la fonction suivante dans le tableau et à zoomer dessus. Vous pouvez également inclure cet extrait dans un complément Python ou un outil de script Python. De plus, pour faciliter les choses lors de l'édition, vous pouvez désactiver les champs inutiles (dans les propriétés de la couche) et également ouvrir le panneau Attributs pour un accès rapide aux attributs.

oeuvre21
la source
Ça marche aussi. Je pense que c'est une bonne illustration de la différence entre ArcObjects et Python dans la longueur et la complexité du code. Notez que cela aura le même problème avec les fichiers de formes et ne fonctionnera pas sur les classes d'entités de géodatabase fichier / personnel / SDE car les valeurs OID / FID ne sont pas garanties d'être contiguës et basées sur 0 - c'est une bizarrerie du format de fichier de formes.
Michael Stimson
@ MichaelMiles-Stimson, oui, je recommanderais d'utiliser la méthode AddFieldDelimiters pour construire une expression de champ correcte entre les sources de fichier, .mdb et .gdb help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//…
artwork21
1
Ce n'est pas que, dans une classe d'entités de géodatabase, les valeurs OID n'ont pas besoin de commencer à 0 et même si elles le font, la suivante n'a pas besoin d'être 1. Vous devez obtenir les OBJECTID dans une liste ou un dictionnaire lorsque vous comptez avec le curseur, vous devez également utiliser arcpy.Describe (fc) .OIDFieldName comme champ car la géodatabase a OBJECTID et non FID.
Michael Stimson
3

Avez-vous l' extension Data Reviewer ? Le réviseur de données vous permet de «parcourir» toutes les fonctionnalités d'un simple clic sur le bouton «Suivant» (zoom à la fois sur l'emplacement spatial et sur l'enregistrement de la table d'attributs). En plus de cela, il y a beaucoup plus de fonctionnalités dans Data Reviewer (comme marquer les erreurs comme «corrigées», «marquées», etc. et exécuter des tâches par lots). Juste un outil standard, même si je suis sûr que votre outil @Michael est également très fantastique!

entrez la description de l'image ici

John
la source
On dirait exactement ce dont j'ai besoin, mais il semble que mon niveau de licence n'inclut pas le réviseur de données :(
user32882
J'ai écrit mon outil avant même d'entendre parler de PLTS (maintenant la cartographie de production), l'inspiration originale a été écrite dans un script AML parce que l'exigence est tellement basique mais nécessaire. Je ne suis pas surpris qu'Esri ait l'outil mais pourquoi n'est-il pas natif dans ArcMap? Le réviseur de données est très bon dans ce qu'il fait, mais j'ai constaté que j'avais déjà écrit la plupart des outils fournis avec quand je l'ai évalué.
Michael Stimson