Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1009 1620 (N31E18 L=232)  1970    340        7 ( 897/1522-1100/1687) -15.5 E27-E09 N25-N34
      2  795 1573 (N22E29 L=221)  1095    202        3 ( 732/1496- 867/1654)  85.6 E30-E27 N16-N29
      3 1506 1321 (N26W22 L=272)  2351    546        1 (1372/1089-1572/1490) -48.0 W10-W29 N12-N34
      4  192 1223 (S11E65 L=186)   769     79        1 ( 162/1192- 222/1272)  71.6 E68-E60 S13-S07
      5  153 1102 (S19E67 L=183)  1987    250        2 (  82/1025- 234/1193)  24.0 E85-E58 S23-S11
      6  408  971 (S21E38 L=213)   723    211        1 ( 344/ 905- 507/1015)  42.3 E42-E30 S26-S16
      7 1780  958 (N12W51 L=301)   373     98        0 (1765/ 913-1792/1011) -80.5 W49-W51 N09-N15
      8  859  850 (S18E05 L=245)   349    153        1 ( 827/ 793- 909/ 895)  71.6 E06-E02 S20-S14
      9  431  853 (S28E34 L=216)   466     95        1 ( 386/ 838- 491/ 872)  12.5 E38-E29 S29-S26
     10 1324  781 (S10W23 L=274)   365    146        0 (1260/ 740-1379/ 807)   0.0 W19-W27 S11-S08
     11  750  387 (S51W00 L=251)  1615    249        3 ( 651/ 308- 846/ 478) -20.9 E10-W11 S54-S46
     12  948  244 (S55W27 L=277)   359     44        1 ( 925/ 224- 975/ 262)  45.0 W23-W29 S57-S52
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y