I have access to ArcMap 10.0 (ArcInfo/Workstation) on Windows 7 SP1.

My usual scenario is to produce a map book having 3 shapes:

  • 1 polyline (a road track about 200-300 km without any useful attribute data)
  • 1 polygon (boundaries with all the county data like county names, numbers etc. in the attribute table)
  • 1 polygon (road track based Strip Map Index Feature (SMIF) to get the DDP running with attributes generated by ArcMap)

All works fine and ArcMap produces around 130-150 ISO A4-PDF pages per map book at a given scale.

But how to get the county names fieldbased on each single page without some manual action?

I know that ArcMap only uses the SMIF-Layer with fixed attribute data as a valid source for DDP, so no Joins/Relates to the attribute table are allowed. And adding new fields manually to SMIF table and filling them actually works, I did that last time, but, you know, it stinks somehow.

So, how to get the attributes from boundaries.shp into SMIF.shp to provide dynamic data for each page in DDP?

  • 1
    To my knowledge, you have to put the county information into the layer that defines your DDP. – Fezter Nov 2 '14 at 23:28
  • Well, Fezter...yes. – mapcat Nov 3 '14 at 11:36

I finally figuerd it out. It may seem a little awkward since it is no arcpy solution, but it worked for me.

  1. Create the regular trackbased SIMPLE_SMIF.shp
  2. Run INTERSECT from toolbox with SIMPLE_SMIF.shp county.shp as input layers. Note that the rank (which of these layers is put above the other) in the dialogue window will influence the result. I placed SIMPLE_SMIF.shp above county.shp. The result layer is INTER_SMIF.shp with all the required county data 'stamped' in it, but it is cracked in far too many pieces, so...
  3. Run DISSOLVE from toolbox with INTER_SMIF.shp as input having the "wanted" fields 'PageNumber', 'CountyInfo1', 'CountyInfo2' as dissolve fields. The result layer DISS_SMIF.shp looks much more handy now with less records, but still there are multiple values for each PageNumber.
  4. To get unique PageNumbers, Add Field 'Area', Type 'Double', Edit, run Calculate Geometry, set Area/sq km. Now you've got a 'decider' on which PageNumber to take. DISS_SMIF.shp
  5. Now run SUMMARIZE STATISTICS from toolbox with DISS_SMIF.shp as input. It will produce an Info Table 'maxgroup' as a result with a new field 'MAX_Area' by calculating the maximum value from a PageNumber group having set these parameters in the dialogue before: Statistics Field = Area, Method = MAX, Case Field = PageNumber.
  6. Now JOIN DISS_SMIF.shp with 'maxgroup' Info Table by 'MAX_Area' and 'Area' keeping only matching records. Export this layer into DISS_SUM_SMIF.shp to have the attributes fixed to the new feature. It should have unique PageNumber values and required county infos as well. enter image description here
  7. Finally JOIN SIMPLE_SMIF.shp with DISS_SUM_SMIF.shp by 'PageNumber' keeping all records. This helps providing the previous settings for overlapping and prevents from creating empty pages in the PDF, f.e. because of missing county.shp of neighborhood. Export this result into VALUED_SMIF.shp and yippieh,
  8. Choose VALUED_SMIF.shp for main dialogue of DDP setup.
| improve this answer | |

To get the attributes from boundaries.shp into SMIF.shp I would use the Union (Analysis) tool:

Computes a geometric union of the input features. All features and their attributes will be written to the output feature class.

Once you have that you will probably see that some of your SMIF polygons fall within more than one of your county polygons so you will need to decide how you wish them to be allocated. I would suggest using either the county that it is mostly in by area, or to list all counties that it falls within comma separated.

All of the above should be easily automated using ArcPy.

| improve this answer | |
  • I'll have a try on the Intersect method as well as mentioned in the ArcToolbox Window here: "Use Intersect when you want to overlay a layer with the polygons in another layer so that the resulting output layer a) has the combined attribute data of the features in the two inputs, and b) only contains features that fall within the spatial extent of the overlay polygons. In this way, you can find just those features that overlap and 'stamp' the attributes of the overlay polygons in the second layer onto the features in the first layer." It's the stamp-attribute-thing I like best in both meths. – mapcat Nov 3 '14 at 11:35

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.