Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  220 1310 (N35E65 L=150)   468     63        0 ( 183/1298- 254/1323)  -4.4 E72-E59 N34-N35
      2  684 1156 (N20E19 L=196)   550    222        0 ( 593/1127- 774/1187) -22.6 E25-E13 N17-N23
      3 1609  967 (S08W36 L=251)   378     97        0 (1575/ 927-1644/1018) -71.6 W33-W38 S10-S04
      4  801  702 (S09E19 L=196)   502    118        2 ( 757/ 661- 829/ 748)  68.2 E21-E16 S11-S06
      5  320  605 (S10E59 L=156)   838    145        2 ( 285/ 554- 357/ 673) -76.0 E61-E57 S13-S05
      6  372  629 (S09E53 L=163)   401     59        0 ( 346/ 604- 400/ 651)  18.4 E55-E49 S10-S06
      7 1691  574 (S36W43 L=259)  1374    198        1 (1601/ 539-1752/ 631)   0.0 W33-W50 S38-S33
      8 1455  531 (S33W21 L=236)   277     46        1 (1428/ 511-1478/ 547) -26.6 W18-W22 S34-S31
      9  521  504 (S18E45 L=170)   490    132        1 ( 473/ 479- 572/ 544) -38.7 E47-E40 S20-S14
     10 1046  451 (S30E10 L=206)   901    207        1 ( 984/ 377-1119/ 524)  31.6 E16-E03 S34-S25
     11  726  488 (S22E30 L=185)   356     71        0 ( 700/ 468- 763/ 513)  11.3 E31-E26 S23-S20
     12  575  335 (S31E51 L=165)  1879    236        3 ( 541/ 240- 620/ 428)  28.2 E69-E41 S39-S24
  

    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